Glacial Loading and Tectonics

Modeling the Influence of the Last Glacial Maximum Ice Load on the Tectonics of southeast Alaska

Lauren Wheeler (a), Peter O. Koons (a), Phaedra Upton (b), Sean Birkel (a)

a School of Earth and Climate Sciences, University of Maine b GNS Science, Lower Hutt, New Zealand

Abstract

Using three-dimensional mechanical modeling of plate boundary deformation and dynamic ice sheet modeling of the LGM (Last Glacial Maximum) ice extent, we tested the sensitivity of local kinematics to rheological structure of southeast Alaska and glacial loading during the LGM. The coast of southeast Alaska is defined by a subduction boundary to the north, a strike-slip fault to the east (Fairweather fault), and the Transition fault to the southwest. The Transition fault separates the thin (10-12 km) Pacific oceanic crust from the thickened continental/oceanic crust of the Yakutat block (~30 km). Seismic profiles and kinematic models imply that the Transition fault is currently active, however, a relative lack of data makes it difficult to constrain the style of deformation taken up along the Transition fault. We incorporate modern topography, a modeled LGM glacial load, and strength profiles based on the velocity and porosity interpretations from seismic profiles within the Gulf of Alaska (Christeson et al. 2010; Worthington et al. 2012; Van Avendonk et al. 2013). Two sets of simulations were run, one without weak zones along the major fault boundaries, referred to as the ‘unfaulted’ model. The second with weak zones along the faults, referred to as the ‘faulted’ model. For both the unfaulted and models we find that LGM glacial loading and unloading changes the pattern and magnitude of strain. The faulted model matches observations best with higher strain along the faulted boundary as less stress is required to overcome the strength along the fault. After glacial unloading there is high strain in the YZ shear strain increment in the Yakutat bay region. The modeled high strain corresponds well with the epicenters of the 1899 earthquakes, a combined tectonic and glacial isostatic response to Little Ice Age unloading (Elliott et al. 2010). These models show that the spatial and temporal distribution of glaciers influences the distribution of strain and that observations of the tectonic response to these loads can be used to constrain the rheology of a region.

Funding sources: The thesis was primarily funded by STEEP (EAR-0409162), salary during the Summer and Fall of 2013 by EAPSI (OISE-1310847) and GeoPRISMS (OCE-1249909)

• 

  Figure 1: STEEP01 and STEEP02 refraction profiles (Worthington et al. 2012; Christeson et al. 2010). A. Location map of local tectonics, Kayak Island Zone (KIZ), Pamplona Fault Zone (PFZ), Dangerous River Zone (DRZ), the Fairweather fault...read more

• 

  Figure 2: Maximum ice extent during the Last Glacial Maximum (at 18,900 years BP) modeled using the University of Maine Ice Sheet Model (UMISM) for the last glacial (~116 ka). Base precipitation and temperature grids are hybrid, derived from...read more

• 

  Figure 3: Topographic maps of A) New Zealand and the regions of B) Wairoa, C) Central Otago, D) Marlborough, E) Canterbury...read more

• 

  Figure 4: A) Dendritic surface created from a stream power-based landscape evolution model created on a uniform substrate...read more

• 

  Figure 5: A) Sine wave surface used as a simplified replication of synform-antiform pairs found in Central Otago. B) Anisotropy...read more

• 

  Figure 6: A) Low relief trough bisecting high relief plateau; a simplified replication of fault damage influenced topography...read more