North American Crustal Motion and Glacial Isostatic Adjustment Model Predictions

A suite of forward GIA model predictions, spanning a wide range of layered mantle viscosity and lithospheric thickness values, is compared to observed horizontal crustal motions in North America to discern optimal model parameters in order to minimize a root-mean-square (RMS) measure of the velocity residuals. To obtain the Earth model response, a combination of the full normal mode analysis and the collocation method is implemented. It provides a means to determine the surface loading response automatically and robustly to 1-dimensional (radially varying) Earth models, while retaining as much of the physics of the normal mode method as numerically feasible, given documented issues with singularities along the negative inverse-time axis in the Laplace transform domain. This method enables the exploration across a wide parameter range (for the lower mantle, transition zone, asthenosphere, and thickness of the elastic lithosphere) to find optimal combinations to explain horizontal crustal motion in North America. The analysis utilizes crustal motion rates from approximately 300 GNSS sites in central North America (Canada and United States) provided by the Nevada Geodetic Laboratory.  Preliminary results indicate that as the lithospheric thickness increases, from 60 km to 240 km, the horizontal motion residuals gradually decrease with no minimum apparent for the thicknesses thus far considered. The residual velocities for the best-fitting models appear to carry a remaining signal, confirming previous inferences of limitations to spherically symmetric Earth models in modeling horizontal crustal motions in North America.