Along your requirements we can generate a realistic planet starting from:
- Height field
- Water Table
- Chemical Ratios and Other Coefficients
GIMP/Photoshop can actually pull off every necessary transformation of that data set using built-in filters if your willing to give up accuracy a little bit and think out the connections you want to model. The one's it can't model you could script in or use other software (one example being measuring insolation. I don't know of a way in standard GIMP to make a shadow map from a height map. QGIS can automatically do this for you though.).
TimB's great series of questions on world-building will probably end up covering more than you care to model but:
- Height Map: This give you pretty much everything you need. Especially if it's filled out with erosion channels for rivers, fault-lines, and a water table.
- Heat Map: Necessary but difficult if you wan't it coupled with everything tightly.
- Shadow Map: Your going to need it for the heat.
- Water Table/Moisture Map: This is going to drive fertility, erosion, rain, a lot.
- Wind Map: This is actually an output of the other maps mixed together.
- Soil Map: type, thickness, etc. These are important for pollution and moisture and other things.
In general gaussian blur is from dispersion and bleed effects (its the solution to heat.) Gradients are used to set up constants (insolation, etc.) Averaging and blurring can mix multiple season and years into one overall climate despite periodic flux. Smearing in a direction (NOT directional blur because that goes in front of the mountain as well) can give you rain shadow effects. A Dijkstra Map/Distance Transform/Grassfire/Eikonal Solution/etc. can get odd thing like distance around a mountain or distance to the nearest set of _____. I'm assuming we already have a height map. If we don't then Voronoi Regions/Stained Glass/Mural/Cobblestone Effects give us segmented regions. Actual basins are simply Voronoi regions mapped to mountain ranges and fractally nested inside each other. To get the crinkly look from actual plate movement you'd need Distortion/Displacement/Glass Maps Primary rain erosion can be modeled by throwing random noise scaled by moisture intensity, etc. River carving and a whole mess of other things can be modeled using the first derivatives of your maps (which is the gradient). etc.
QGIS was mentioned and goes with GRASS GIS and a lot of other open-source GIS. It doesn't do a whole lot along these lines but it can be interfaced with R. So there's probably a program out there for everything you can think of and some you can't.
If you're not in the mood to code something from scratch but know how to run a compiler you can try doing searches for Fortran, Astrophysics, and Geo-etc. as most of the physics have remained the same, so all the good programs are old University code.
Blender and Rhino (and Houdini?) can be used in a similar way to GIMP/Photoshop but to greater effect. (We drop some error as we move from a flat representation to a globe. We can also actually just outright do some things like place the sun and model an elliptic orbit).
EdGCM was apparently developed as an educational but full functional climate-modelling tool by NASA (it supposedly also doesn't kill your computer like some other newer climate-modelling software).
Gocad looks amazing although I haven't tried it.