Crustal Viscosity Structure Estimated from Multi-Phase Mixing Theory

Estimates of lower crustal viscosity are typically constrained by analyses of isostatic rebound, post seismic creep, and laboratory-derived flow laws for crustal rocks and minerals. Here we follow a new approach for calculating the viscosity structure of the lower continental crust. We use Perple_X to calculate mineral assemblages for different crustal compositions. Effective viscosity is then calculated using the rheologic mixing model of Huet et al. (2014) incorporating flow laws for each mineral phase. Calculations are performed along geotherms appropriate for the Basin and Range, Tibetan Plateau, Colorado Plateau, and the San Andreas Fault. To assess the role of crustal composition on viscosity, we examined two compositional gradients extending from an upper crust with ~67 wt% SiO2 to a lower crust that is either: (i) basaltic with ~53 wt% SiO2 (Rudnick and Gao, 2003), or (ii) andesitic with ~64% SiO2 (Hacker et al., 2011). In all cases, the middle continental crust has a viscosity that is 2-3 orders of magnitude greater than that inferred for wet quartz, a common proxy for mid-crustal viscosities. An andesitic lower crust results in viscosities of 1020-1021 Pa-s and 1021-1022 Pa-s for hotter and colder crustal geotherms, respectively. A mafic lower crust predicts viscosities that are an order of magnitude higher for the same geotherm. In all cases, the viscosity calculated from the mixing model decreases less with depth compared to single-phase estimates. Lastly, for anhydrous conditions in which alpha quartz is stable, we find that there is a strong correlation between Vp/Vs and bulk viscosity; in contrast, little to no correlation exists for hydrous conditions.