There is one problem with ice age modelling in general (and I think we should focus on that, earthquake and volcanic events are already there and as Koshling said most other tectonic changes have a far larger time scale): The map scripts have different approaches and currently yield a post ice age map with only plot and terrain information, no height map or climate map (although several use that internally).
So either those are changed to provide height and climate maps for a pre ice age map or there has to be an inverse approach to ice age modelling. Instead of computing the ice age effect as it happens, you start with the post ice age map from the map scripts, then apply the inverse ice age effects, then when the actual turns happen, you slowly undo them until you end up with the post ice age map from the map scripts (which you have stored separately on the plots).
On the other hand going for the full change to the map scripts so you always have the height, temperature and climate information for all plots would allow for quite some interesting dynamic simulations with the property system.
So either those are changed to provide height and climate maps for a pre ice age map or there has to be an inverse approach to ice age modelling. Instead of computing the ice age effect as it happens, you start with the post ice age map from the map scripts, then apply the inverse ice age effects, then when the actual turns happen, you slowly undo them until you end up with the post ice age map from the map scripts (which you have stored separately on the plots).
On the other hand going for the full change to the map scripts so you always have the height, temperature and climate information for all plots would allow for quite some interesting dynamic simulations with the property system.