Construction of late Pleistocene ice history model that is consistent with Composite Rheology and 3D material property variations
In our research, the response of the solid earth and its oceans to the glaciation and deglaciation events since the last Ice Age is modeled numerically. This response of the Earth is called the Glacial Isostatic Adjustment (or GIA) process and numerical modeling of such process is very important both practically and academically (see ‘significance’ below). Previous models of the GIA process assume that the flow law or relation between stress and strain in the mantle is linear and that material properties of the earth only vary in the radial direction. However, creep experiments on mantle rocks shows that the flow law in the mantle is likely to be nonlinear. Also, seismic tomography shows that earth material properties vary strongly in the lateral direction. What is special in our model is that material properties are allowed to vary both in the radial and lateral directions as described by seismic tomography. Also a combination of both linear and nonlinear creep is implemented in our model (i.e. Composite Rheology). The ABAQUS finite element software package is used in this study. Our finite element model has spatial resolution at the surface of 1x1 degree grid covering the whole mantle with more than 20 layers. This result in a finite element model of at least 1.3 million finite elements. Thus such computation requires high performance computer with at least 6.4 GB of RAM to run and 200 GB of storage for the results of each job.