Tectonics and Structure Commons^™

Crustal And Upper-Mantle Structure Beneath Ice-Covered Regions In Antarctica From S-Wave Receiver Functions And Implications For Heat Flow, C. Ramirez, Andrew A. Nyblade, S. E. Hansen, Douglas A. Wiens, Sridhar Anandakrishnan, Richard C. Aster, Audrey D. Huerta, Partick Shore, Terry Wilson

Geological Sciences Faculty Scholarship

S-wave receiver functions (SRFs) are used to investigate crustal and upper-mantle structure beneath several ice-covered areas of Antarctica. Moho S-to-P (Sp) arrivals are observed at ~6–8 s in SRF stacks for stations in the Gamburtsev Mountains (GAM) and Vostok Highlands (VHIG), ~5–6 s for stations in the Transantarctic Mountains (TAM) and the Wilkes Basin (WILK), and ~3–4 s for stations in the West Antarctic Rift System (WARS) and the Marie Byrd Land Dome (MBLD). A grid search is used to model the Moho Sp conversion time with Rayleigh wave phase velocities from 18 to 30 s period to estimate crustal …

Upper Mantle Structure Of Central And West Antarctica From Array Analysis Of Rayleigh Wave Phase Velocities, David S. Heeszel, Douglas A. Wiens, Sridhar Anandakrishnan, Richard C. Aster, Ian W.D. Dalziel, Audrey D. Huerta, Andrew A. Nyblade, Terry J. Wilson, J. Paul Winberry

Geological Sciences Faculty Scholarship

The seismic velocity structure of Antarctica is important, both as a constraint on the tectonic history of the continent and for understanding solid Earth interactions with the ice sheet. We use Rayleigh wave array analysis methods applied to teleseismic data from recent temporary broadband seismograph deployments to image the upper mantle structure of central and West Antarctica. Phase velocity maps are determined using a two–plane wave tomography method and are inverted for shear velocity using a Monte Carlo approach to estimate three-dimensional velocity structure. Results illuminate the structural dichotomy between the East Antarctic Craton and West Antarctica, with West Antarctica …

Hazard Identification And Coastal Stratigraphy In Crescent Harbor, Northeast Whidbey Island, Washington, Brian Ostrom

All Master's Theses

Crescent Harbor marsh, on northeastern Whidbey Island, records evidence of co-seismic land-level change 1825 to 1925 cal. yrs. BP. The lithostratigraphy and diatom microfossil assemblages reveal a marsh peat abruptly overlain by intertidal mud, indicating rapid subsidence. Analysis of the modern-day position of depositional facies indicates subsidence from a high marsh to a tidal-flat environment representing an estimated 1.7 m elevation change. The timing of subsidence fits within the dates of a rupture found on the nearby Utsalady Point fault between 1,100 and 2,200 years BP (Johnson et al. 2004). Likely, the stratigraphy at Crescent Harbor records the same event …

Investigating Taphonomic Changes Of Deposits And Modeling Of The 2010 Earthquake And Tsunami In South-Central Chile, Alexandra Carranco Ruiz

All Master's Theses

South-central Chile has an extensive written catalog of historic earthquakes and tsunamis, but such records can be subject to inconsistencies. Dated tsunami deposits are more objective data that provide hard evidence of past tsunamis. The inland extent of deposits from past tsunamis (paleodeposits) can be used in tsunami modeling to reveal characteristics of the source earthquake, but these deposits may have undergone taphonomic processes since initial deposition. Therefore, to determine how tsunami deposits change during burial and preservation and the potential limitations of using paleodeposits in modeling, I investigated the modern 2010 M_w 8.8 Maule earthquake and tsunami as …

Spatiotemporal Slip Rate Variations Along Surprise Valley Fault In Relation To Pleistocene Pluvial Lakes, Brian N. Marion

All Master's Theses

Using mapped paleoshoreline features with high-resolution topographic data and obtained radiocarbon dates on paleoshoreline tufas, I documented precise fault offsets of dated features over the last 25 ka along the Surprise Valley Fault (SVF). Fault offset measured in three lake sections within Surprise Valley ranged from 3.6 m in the southern section to 14.4 m in the central section. The offset paleoshorelines are dated to the late Pleistocene (<22 >ka) and were formed during the latest impoundment of pluvial Lake Surprise since the last glacial maximum. Slip rates vary along strike, assuming a fault dip of 68° with 0.25 ± …

Hidden Intrabasin Extension: Evidence For Dike-Fault Interaction From Magnetic, Gravity, And Seismic Reflection Data In Surprise Valley, Northeastern California, Noah D. Athens, Jonathan M.G. Glen, Simon L. Klemperer, Anne E. Egger, Valentina C. Fontiveros

All Faculty Scholarship for the College of the Sciences

The relative contributions of tectonic and magmatic processes to continental rifting are highly variable. Magnetic, gravity, and seismic reflection data from Surprise Valley, California, in the northwest Basin and Range, reveal an intrabasin, fault-controlled, ~10-m-thick dike at a depth of ~150 m, providing an excellent example of the interplay between faulting and dike intrusion. The dike, likely a composite structure representing multiple successive intrusions, is inferred from modeling a positive magnetic anomaly that extends ~35 km and parallels the basin-bounding Surprise Valley normal fault on the west side of the valley. A two-dimensional high-resolution seismic reflection profile acquired across the …

Slip Pulse And Resonance Of The Kathmandu Basin During The 2015 Gorkha Earthquake, Nepal, John Galetzka, Walter Szeliga

All Faculty Scholarship for the College of the Sciences

Detailed geodetic imaging of earthquake ruptures enhances our understanding of earthquake physics and associated ground shaking. The 25 April 2015 moment magnitude 7.8 earthquake in Gorkha, Nepal was the first large continental megathrust rupture to have occurred beneath a high-rate (5-hertz) Global Positioning System (GPS) network. We used GPS and interferometric synthetic aperture radar data to model the earthquake rupture as a slip pulse ~20 kilometers in width, ~6 seconds in duration, and with a peak sliding velocity of 1.1 meters per second, which propagated toward the Kathmandu basin at ~3.3 kilometers per second over ~140 kilometers. The smooth slip …

The Mantle Transition Zone Beneath West Antarctica: Seismic Evidence For Hydration And Thermal Upwellings, E. L. Emry, A. A. Nyblade, J. Juliá, S. Anandakrishnan, R. C. Aster, D. A. Wiens, Audrey D. Huerta, T. J. Wilson

All Faculty Scholarship for the College of the Sciences

Although prior work suggests that a mantle plume is associated with Cenozoic rifting and volcanism in West Antarctica, the existence of a plume remains conjectural. Here we use P wave receiver functions (PRFs) from the Antarctic POLENET array to estimate mantle transition zone thickness, which is sensitive to temperature perturbations, throughout previously unstudied parts of West Antarctica. We obtain over 8000 high‐quality PRFs using an iterative, time domain deconvolution method filtered with a Gaussian width of 0.5 and 1.0, corresponding to frequencies less than ∼0.24 and ∼0.48 Hz, respectively. Single‐station and common conversion point stacks, migrated to depth using the …

Estimation Of Offsets In Gps Time-Series And Application To The Detection Of Earthquake Deformation In The Far-Field, Jean-Philippe Montillet, S. D. P. Williams, A. Koulali, S. C. Mcclusky

All Faculty Scholarship for the College of the Sciences

Extracting geophysical signals from Global Positioning System (GPS) coordinate time-series is a well-established practice that has led to great insights into how the Earth deforms. Often small discontinuities are found in such time-series and are traceable to either broad-scale deformation (i.e. earthquakes) or discontinuities due to equipment changes and/or failures. Estimating these offsets accurately enables the identification of coseismic deformation estimates in the former case, and the removal of unwanted signals in the latter case which then allows tectonic rates to be estimated more accurately. We develop a method to estimate accurately discontinuities in time series of GPS positions at …

Rapid Middle To Late Miocene Slip Along The Zanskar Normal Fault, Greater Himalayan Range, Nw, India: Constraints From Low-Temperature Thermochronometry, Brett L. Shurtleff

All Master's Theses

The Zanskar normal fault (ZF) is a NW-striking, moderately NE-dipping, normal fault that bounds the northern flank of the Greater Himalaya Range, NW India. The ZF is the far west continuation of the South Tibetan Detachment System (STDS), a major arc-parallel normal sense shear zone that spans the length of the Himalayan orogen. Detailed new zircon and apatite (U-Th)/He (ZHe and AHe) and apatite fission-track (AFT) thermochronometric data from high-grade (amphibolite-migmatite) Greater Himalayan Sequence (GHS) metamorphic rocks, exposed in the footwall immediately adjacent to the ZF, provide constraints on the middle Miocene to present exhumation history of the footwall. The …

Geologic Mapping In The Black Mountain Area, Northern Eastern California Shear Zone: Testing A Kinematic And Geometric Fault Slip Transfer Model, Kevin M. Delano

All Master's Theses

New geologic mapping, structural, kinematic, and ⁴⁰Ar/³⁹Ar geochronology studies in the Black Mountain area, northern eastern California shear zone, are used to test a kinematic fault slip transfer model for the Owens Valley fault-Mina deflection transition. In the Black Mountain area, range bounding ~NNW- to ~NS-striking and lesser NW- to NE-striking normal faults cut Mesozoic, Miocene (22.42 ± 0.05 Ma), Pliocene (3.53 ± 0.06 to 3.29 ± 0.02 Ma), and early-middle Pleistocene (1936 ± 12.7 to 766 ± 3.1 ka) rocks. Palinspastically restored cross-sections show that offset Pliocene markers record 1.5 +0.7/-0.6 km of ~ENE-WSW extension since …

Thermodynamic Model For Energy-Constrained Open-System Evolution Of Crustal Magma Bodies Undergoing Simultaneous Recharge, Assimilation And Crystallization: The Magma Chamber Simulator, Wendy A. Bohrson, Frank J. Spera, Mark S. Ghiorso, Guy Brown, Jeffrey Creamer, Aaron Mayfield

All Faculty Scholarship for the College of the Sciences

The Magma Chamber Simulator quantifies the impact of simultaneous recharge, assimilation and crystallization through mass and enthalpy balance in a multicomponent–multiphase (melt + solids ± fluid) composite system. As a rigorous thermodynamic model, the Magma Chamber Simulator computes phase equilibria and geochemical evolution self-consistently in resident magma, recharge magma and wallrock, all of which are connected by specified thermodynamic boundaries, to model an evolving open-system magma body. In a simulation, magma cools from its liquidus temperature, and crystals ± fluid are incrementally fractionated to a separate cumulate reservoir. Enthalpy from cooling, crystallization, and possible magma recharge heats wallrock from its …

Clockwise Rotation Of The Brahmaputra Valley Relative To India: Tectonic Convergence In The Eastern Himalaya, Naga Hills, And Shillong Plateau, ‪Philippe Vernant, R. Bilham, Walter Szeliga, D. Drupka, S. Kalita, A. K. Bhattacharyya, V. K. Gaur

All Faculty Scholarship for the College of the Sciences

GPS data reveal that the Brahmaputra Valley has broken from the Indian Plate and rotates clockwise relative to India about a point a few hundred kilometers west of the Shillong Plateau. The GPS velocity vectors define two distinct blocks separated by the Kopili fault upon which 2–3 mm/yr of dextral slip is observed: the Shillong block between longitudes 89 and 93°E rotating clockwise at 1.15°/Myr and the Assam block from 93.5°E to 97°E rotating at ≈1.13°/Myr. These two blocks are more than 120 km wide in a north‐south sense, but they extend locally a similar distance beneath the Himalaya and …

Seismic And Geodetic Evidence For Grounding-Line Control Of Whillans Ice Stream Stick-Slip Events, Martin J. Pratt, J. Paul Winberry, Douglas A. Wiens, Sridhar Anandakrishnan, Richard B. Alley

All Faculty Scholarship for the College of the Sciences

The tidally modulated, stick‐slip events of Whillans Ice Stream in West Antarctica produce seismic energy from three locations near the grounding line. Using ice velocity records obtained by combining time series from colocated broadband seismometers and GPS receivers installed on the ice stream during the 2010–2011 and 2011–2012 austral summers, along with far‐field seismic recordings of elastic waves, we locate regions of high rupture velocity and stress drop. These regions, which are analogous to “asperities” in traditional seismic fault studies, are areas of elevated friction at the base of the ice stream. Slip events consistently initiate at one of two …

Tidal Pacing, Skipped Slips And The Slowdown Of Whillans Ice Stream, Antarctica, J. Paul Winberry, Sridhar Anandakrishnan, Richard B. Alley, Douglas A. Wiens, Martin J. Pratt

All Faculty Scholarship for the College of the Sciences

We summarize new observations of the deceleration and stick–slip motion of Whillans Ice Stream (WIS), Antarctica. We refine the location of the large sticky spots that resist motion between slip events, the locations of which are controlled by the patterns of subglacial water flow. Our examination of the long-term velocity time series for the ice stream reveals that the decadal-scale deceleration is not occurring at a steady rate, but varies at the sub-decadal timescale. This unsteady deceleration modulates the temporal evolution of a broad (~50 km across) surface-elevation bulge forming at the junction between the relatively narrow upstream portion of …

Thermobarometry And Geochronology In The Dulan Region, North Qaidam Ultrahigh-Pressure Terrane: Resolving Spatial Variation Of Ages, Temperatures, And Pressures, Megan Eileen Regel

All Master's Theses

The Dulan ultrahigh-pressure (UHP) terrane exposes 449-422 Ma eclogites in the east and 449-418 Ma high-pressure granulites to the west, but the relationship between the units is poorly understood. Zirconium-in-rutile temperatures show an east-west 3- 5°C/km gradient in four UHP eclogites: 662-716±10°C; garnet-omphacite-phengite thermobarometry of peak conditions yields 659-730±82°C, 25.3-27.1±3.2 kbars. Zircon U-Pb ages of four eclogites are the same within uncertainty, averaging 434±2.9 Ma. Tiin- zircon thermometry and rare earth element patterns indicate these ages represent eclogite facies metamorphism. One high-pressure/high-temperature (HP/HT) granulite yields a minimum rutile temperature of 773±4°C. Previous thermobarometry of highpressure granulites yields ≥800°C and 14-17 …

Spatial Relationship Between Gps Slip And Seismic Tremor During Cascadia Slow Slip Events, Hillary Lynn Goodner

All Master's Theses

We model GPS deformation and timing of seismic tremor associated with transient deformation in Cascadia to test the hypothesis that tremor and slip occur synchronously but are spatially offset. For the period 2010–2013, we use seismic tremor data with a duration-moment relationship to predict GPS time series and compare them to observations. We find that observed GPS displacements are best predicted when tremor locations on the plate interface are shifted 15 km up-dip of their published epicenter. To test whether the spatial offset of tremor and slip is due to systematic mislocation of published epicenters, we attempt to identify individual …

Constraints On Interplate Locking Within The Makran Subduction Zone Using Gps, Elyse Erin Frohling

All Master's Theses

We use GPS measurements from 20 stations located in Iran, Pakistan, and Oman along with International Seismological Centre hypocentral locations to create a fault coupling model for the Makran subduction zone. The 3-D geometry of the megathrust fault is estimated from five profiles landward of the Makran trench. Continuous GPS on the Arabian plate show convergence of ~27 mm/year, and interseismic deformation is calculated using the back-slip method of Savage [ 1983]. We find the entire subduction zone appears to be accumulating stress uniformly from ~15-60 km depth and find no evidence of along-strike segmentation in strain accumulation, suggesting the …

Proxy-Based Reconstructions Of Earthquakes And Tsunamis At Quidico, South-Central Chile, Isabel J. Hong

All Master's Theses

We utilized geomorphic, microfossil, sedimentological, and stratigraphic methods to investigate the paleotsunami history at Quidico, Chile (38.1° S, 73.2° W). A combination of pits, cores, and riverbank exposures provide stratigraphic evidence of 8 sand layers at Quidico, including tsunami deposits from 2010 and 1960. The sands are laterally continuous but display landward thinning, disappearing completely from the stratigraphy 1.2 km inland from the coast. Intervening sediment layers consist of organicrich silts and peat. Grain size analysis shows the sand units are bimodal in distribution whereas the organic-rich silts and peat have a broader, bimodal distribution. Lithological analysis showed similarities among …

Controls On Fault Geometry During Early Stages Of Extension In The Larkspur Hills, Northwest Basin And Range, Diana Jean Strickley

All Master's Theses

Detailed analyses of normal faults in the Larkspur Hills, CA-NV, northwest Basin and Range, offer insight into factors controlling normal fault initiation, growth, and distribution. N-trending faults in the southern portion of the study area share trends of major range-bounding structures and Pliocene linear volcanic vents; in contrast, NNW- and NNE- trending faults dominate further north and into south-central Oregon. Stress analyses and comparison with experimental and field data suggest that preexisting structures control faults in the northern Larkspur Hills, while faults form perpendicular to σ3 in the southern hills. The change in fault orientations is abrupt, occurring across a …

The Crustal Thickness Of West Antarctica, Julien A. Chaput, Richard C. Aster, Audrey D. Huerta, X. Sun, Andrew Lloyd, Douglas A. Wiens, Andrew Nyblade, Sridhar Anandakrishnan, J. Paul Winberry, Terry Wilson

Geological Sciences Faculty Scholarship

P-to-S receiver functions (PRFs) from the Polar Earth Observing Network (POLENET) GPS and seismic leg of POLENET spanning West Antarctica and the Transantarctic Mountains deployment of seismographic stations provide new estimates of crustal thickness across West Antarctica, including the West Antarctic Rift System (WARS), Marie Byrd Land (MBL) dome, and the Transantarctic Mountains (TAM) margin. We show that complications arising from ice sheet multiples can be effectively managed and further information concerning low-velocity subglacial sediment thickness may be determined, via top-down utilization of synthetic receiver function models. We combine shallow structure constraints with the response of deeper layers using a …

Glacier Slip And Seismicity Induced By Surface Melt, Peter L. Moore, J. Paul Winberry, Neal R. Iverson, Knut A. Christianson, Sridhar Anandakrishnan, Miriam Jackson, Mark E. Mathison, Denis Cohen

All Faculty Scholarship for the College of the Sciences

Many of the key processes governing fast glacier flow involve interaction between a glacier and its basal hydrological system, which is hidden from direct observation. Passive seismic monitoring has shown promise as a tool for remotely monitoring basal processes, but lack of glacier-bed access prevents clear understanding of the relationships between subglacial processes and corresponding seismic emissions. Here we describe direct measurements of basal hydrology, sliding, and broadband seismicity made in a unique subglacial facility in Norway during the onset of two summer melt seasons. In the most pronounced of these episodes, rapid delivery of surface meltwater to the bed …

Alluvial Sediment Or Playas: What Is The Dominant Source Of Sand And Silt In Desert Soil Vesicular A Horizons, Southwest Usa, Mark R. Sweeney, Eric V. Mcdonald, Christohper E. Markley

All Faculty Scholarship for the College of the Sciences

Vesicular A (Av) soil horizons form beneath desert pavements from the accretion of aeolian sediment (dust) commonly thought to be derived primarily from desiccating pluvial lakes and playas, with contributions from ephemeral washes and alluvial fans. Particle size distributions of Av horizons are typically bimodal with primary modes of very fine silt and fine sand, suggesting that the horizon matrix is derived from multiple sources. Here we conduct detailed chemical and physical analysis of both Av horizon soil samples and potential sources of aeolian sediment to better constrain the relative contributions of dust associated with the development of Av horizons. …

2012 Haida Gwaii Quake: Insight Into Cascadia's Subduction Extent, Walter Szeliga

Geological Sciences Faculty Scholarship

The limits of Cascadia were first defined to contain nearly the entire margin of the Pacific Northwest, from Cape Mendocino through the Alaska Panhandle [Schuchert, 1910; Schuchert and Barrell, 1914]. Since that time, the boundary of Cascadia has shrunk to become essentially synonymous with the region where the Juan de Fuca plate subducts beneath the North American plate. As a consequence, seismic hazard assessments in the Pacific Northwest have conventionally focused on the potential for large megathrust earthquakes along the interface of the Juan de Fuca and North American plates.

Nucleation And Seismic Tremor Associated With The Glacial Earthquakes Of Whillans Ice Stream, Antarctica, J. Paul Winberry, Sridhar Anandakrishnan, Douglas A. Wiens, Richard B. Alley

All Faculty Scholarship for the College of the Sciences

The ability to monitor transient motion along faults is critical to improving our ability to understand many natural phenomena such as landslides and earthquakes. Here, we usedata from a GPS and seismometer network that were deployed to monitor the regularly repeating glacial earthquakes of Whillans Ice Stream, West Antarctica to show that a unique pattern of precursory slip precedes complete rupture along the bed of the ice stream. Additionally, we show that rupture can be independently tracked by increased levels of microseismic activity, including harmonic tremor, that are coincident with the onset of slip at any location, thus providing a …

Pliocene Sinistral Slip Across The Adobe Hills, Eastern California-Western Nevada: Kinematics Of Fault Slip Transfer Across The Mina Deflection, Sarah Nagorsen-Rinke, Jeffrey Lee, Andrew Calvert

All Faculty Scholarship for the College of the Sciences

The Adobe Hills region (California and Nevada, USA) is a faulted volcanic field located within the western Mina deflection, a right-stepping zone of faults that connects the northern Eastern California shear zone (ECSZ) to the south with the Walker Lane belt (WLB) to the north. New detailed geologic mapping, structural studies, and ⁴⁰Ar/³⁹Ar geochronology in the Adobe Hills allow us to calculate fault slip rates and test predictions for the kinematics of fault slip transfer into the Mina deflection. The Adobe Hills are dominated by Pliocene tuffaceous sandstone, basaltic lavas that yield ⁴⁰Ar/³⁹Ar ages …

Paleoseismologic Evidence For Holocene Activity On The Pinto Mountain Fault, Twentynine Palms, California, Ana Maria Cadena

All Master's Theses

Excavations across the Pinto Mountain fault in Twentynine Palms, California exposed faulted strata across a 32-m wide zone. Trench wall exposures revealed clear evidence for five ground-rupturing events during the Holocene, and two additional events in the late Pleistocene. Optically stimulated luminescence ages from alluvial sediments suggest that the most recent event occurred between 1.7-2.9 ka B.P. and the penultimate event between 2.7-4.2 ka B.P.. Prior to the penultimate event, there were five ground-rupturing earthquakes on the eastern Pinto Mountain fault between 3.5-13.6 ka B.P.. The average recurrence interval since 13.3-13.6 ka B.P. is 1510-1680 years, and 1200-1500 years in …

Interseismic Strain Accumulation Along The Western Boundary Of The Indian Subcontinent, Walter Szeliga, Roger Bilham, Din Mohammad Kakar, Sarosh H. Lodi

All Faculty Scholarship for the College of the Sciences

Despite an overall sinistral slip rate of ≈3 cm/yr, few major earthquakes have occurred in the past 200 years along the Chaman fault system, the western boundary of the India Plate with the Eurasia Plate. GPS and InSAR data reported here indicate sinistral shear velocities of 8–17 mm/yr across the westernmost branches of the fault system, suggesting that a significant fraction of the plate boundary slip is distributed in the fold and fault belt to the east. At its southernmost on‐land segment (≈26°N), near the triple junction between the Arabia, Eurasia, and India Plates, we find the velocity across the …

Wilson Cycles, Tectonic Inheritance, And Rifting Of The North American Gulf Of Mexico Continental Margin, Audrey D. Huerta, Dennis L. Harry

Geological Sciences Faculty Scholarship

The tectonic evolution of the North American Gulf of Mexico continental margin is characterized by two Wilson cycles, i.e., repeated episodes of opening and closing of ocean basins along the same structural trend. This evolution includes (1) the Precambrian Grenville orogeny; (2) formation of a rift-transform margin during late Precambrian opening of the Iapetus Ocean; (3) the late Paleozoic Ouachita orogeny during assembly of Pangea; and (4) Mesozoic rifting during opening of the Gulf of Mexico. Unlike the Atlantic margins, where Wilson cycles were first recognized, breakup in the Gulf of Mexico did not initially focus within the orogen, but …

Size And Exhumation Rate Of Ultrahigh-Pressure Terranes Linked To Orogenic Stage, Andrew R.C. Kylander-Clark, Bradley R. Hacker, Chris G. Mattinson

All Faculty Scholarship for the College of the Sciences

A growing set of data indicates a stark contrast between the evolution of two types of ultrahigh-pressure (UHP) terranes: large terranes that evolved slowly (over 10–30 Myr), and small terranes that formed and were exhumed on timescales of < 10 Myr. Here we compare the characteristics – area, thickness, formation rate, exhumation rate, age, and tectonic setting – of these two endmember types of UHP terrane worldwide. We suggest that the two UHP terrane types may form during different orogenic stages because of variations in the buoyancy and traction forces due to different proportions of subducting crust and mantle lithosphere or to different rates of subduction. The initial stages of continent collision involve the subduction of thin continental crust or microcontinents, and thus tectonic forces are dominated by the density of the oceanic slab; subduction rates are rapid and subduction angles are initially steep. However, as collision matures, thicker and larger pieces of continental material are subducted, and the positive buoyancy of the down-going slab becomes more prominent; subduction angles become gentle and convergence slows. Assessing the validity of this hypothesis is critical to understanding the physical and chemical evolution of Earth's crust and mantle.

Included here is the post-print copy of this article. The final publication is available via ScienceDirect at http://www.sciencedirect.com/science/article/pii/S0012821X11007564