North American Plate Research Articles

New insights into the petrogenesis of the Puerto Vallarta Batholith, Mexico: Evidence from petrology, zircon petrochronology, and phase equilibrium modeling

At subduction zones, continental magmatic arcs are produced by differentiation processes. The Puerto Vallarta Batholith, at the Mexican Pacific margin, is the final product of the subduction of oceanic Farallon plate beneath the North American plate during the Cretaceous. Although this Cordilleran Batholith has been subject to different studies, its petrogenetic evolution remains unclear. To better constrain the magmatic processes that gave origin to the granitic rocks of the Puerto Vallarta Batholith, we combined detailed fieldwork, zircon petrochronology, and conventional thermobarometry coupled with phase equilibrium modeling. We present the crystallization and emplacement evolution, and its relationship to magmatic processes of the Puerto Vallarta Batholith. In the granitic rocks, we observed calc-alkaline and metaluminous signatures, as well as enrichment in light-rare earth elements (LREE) relative to the heavy-rare earth elements (HREE). The early crystallization stages of the Puerto Vallarta Batholith started at ~729–776 °C and ~3.5–5.3 kbar, with the last melt portion crystallized at ~675–685 °C. Fieldwork observations and Zr saturation thermometry reveal mixing of melts during the last crystallization stages. Geochemical signatures of zircons from granitic rocks suggest crystallization after minor plagioclase fractionation with likely assimilation of crustal components. Our results from geochemical whole-rock data further indicate partial melting of a heterogenous source in the mantle wedge and the differentiation of those primary melts as the principal processes that gave origin to the Puerto Vallarta Batholith. We show that the magmatic processes that originated the granitic rocks of this Cordilleran Batholith occurred in the middle crust during the Late Cretaceous. Our results provide new evidence to better constrain the origin of this magmatic system.

Cenozoic Exhumation History of the Eastern Margin of the Northern Canadian Cordillera

AbstractNew low‐temperature thermochronology data from clastic sedimentary rocks in the northern Richardson Mountains, Canada, indicate significant exhumational cooling during late Eocene–early Oligocene time. Apatite (U‐Th‐Sm)/He (AHe) data were collected from 19 Proterozoic–Paleocene rocks across a 115 km transect. Eighty‐eight single‐grain AHe dates range from 16–300 Ma and are generally younger than stratigraphic ages, indicative of thermal resetting by burial. Additionally, zircon (U‐Th)/He (ZHe) dates from two Proterozoic–Cambrian rocks range from 49–123 Ma and suggest burial to >160°C. In contrast, ZHe dates from Jurassic sandstones are older than the stratigraphic age, which limits maximum burial to <160°C. Thermal history modeling reveals three phases of cooling, during the Paleocene–early Eocene (>65–50 Ma), late Eocene–early Oligocene (40–30 Ma), and late Oligocene–early Miocene (30–15 Ma). Most samples cooled during the first and second phases, whereas the third phase is less well constrained. In general, most rocks were below the sensitivity of AHe analysis since the early–middle Miocene. The results suggest a previously unrecognized phase of inferred deformation in the northern Richardson Mountains between 40–30 Ma. Our findings contribute to previous work that recognizes Late Cenozoic deformation along the eastern margin of the Northern Cordillera. We further investigated the potential mechanisms of this widespread deformation and suggest exhumation may relate to kinematic changes of the North American plate relative to structural trends along the margin of the Northern Cordillera.

New Paleomagnetic Data on Late Cretaceous Chukotka Volcanics: the Chukotka Block Probably Underwent Displacements Relative to the North American and Eurasian Plates after the Formation of the Okhotsk-Chukotka Volcanic Belt?

New Paleomagnetic Data on Late Cretaceous Chukotka Volcanics: the Chukotka Block Probably Underwent Displacements Relative to the North American and Eurasian Plates after the Formation of the Okhotsk-Chukotka Volcanic Belt?

Influence of Mantle Dynamic Topographical Variations on US Mid‐Atlantic Continental Margin Estimates of Sea‐Level Change

AbstractSpatial analysis of discrepancies in sea‐level estimates derived from “backstripping” Mid‐Atlantic margin cores reveals a coherent signal that can be fit with a 103 km wavelength, 45 m amplitude sinusoid that moved across the margin at a rate and direction of motion generally opposite to that of the North American Plate. This signal we observe suggests topographical uplift occurred over much of the Mid‐Atlantic region since 35 Ma and may be superimposed upon a longer‐wavelength signal of Cenozoic subsidence associated with the subducted Farallon plate passing beneath the Mid‐Atlantic margin. Our statistical modeling of Mid‐Atlantic margin strata suggests that: (a) Cenozoic subsidence is likely to have occurred, but is very unlikely to have exceeded 100 m in magnitude; and (b) variations in ocean basin volume are likely to have contributed to 39 ± 24 m (1σ) of global‐mean geocentric sea‐level fall over the past 55 million years.

The great tectonic discontinuity of Bahía de Banderas, Mexico

Detachment of the Baja California peninsula from the North American plate altered the crustal region between them, including an extension of the crust between the peninsula and Mexico's mainland. In the SE portion of this region Bahía de Banderas, a 60-km indentation in the littoral of Mexico, is a boundary region across which significant changes in topography, seismicity, seismic attenuation, trench subduction angle, continental platform width, geomagnetic field, and crustal thickness have been documented. Here we identify two additional, major discontinuities: one in gravity and the second in the regional electrical resistivity to 50-km depth, both associated with the region in which the crust changes from extended to un-extended. Under the former, we ratify and extend the presence of the oceanic slab obtained elsewhere with seismological determinations. Plate rollback is also reported in the boundary region. Gravity and electrical resistivity models locate the NW and SE limits of the rolled-back portion of the oceanic slab; within this region, we identify upwelling asthenospheric material. Bahía de Banderas region represents a major tectonic discontinuity, implying the existence of a continental sliver prolonging ~160 km into the San Blas basin north of the bay, the Tres Marías sliver, independent from the Jalisco block, in which crustal fragmentation processes appear to be still in progress. We argue that this sliver is tectonically different from the Jalisco block since at least the late Miocene, when extensional tectonics stretched and thinned its corresponding crust and lithosphere.

Present‐Day Deformation Mechanism of the Northeastern Mina Deflection Revealed by the 2020 Mw 6.5 Monte Cristo Range Earthquake

AbstractThe 15 May 2020 Monte Cristo Range Mw 6.5 earthquake occurred in the northeast of the Mina deflection, which accommodates approximately a quarter of the relatively dextral motion between the Pacific and North American plates. The Monte Cristo Range event provides an opportunity to study the present‐day regional deformation mechanism and active tectonics. In this study, we investigate the source rupture process of the event using joint inversion of interferometric synthetic aperture radar and broadband seismic data. We find that the rupture propagates almost simultaneously on two main segments. The fault motion changes from the predominantly sinistral slip near the epicenter on the eastern segment to the oblique slip on the western segment, with a maximum coseismic slip of 0.8 m. Our results suggest that the accommodation of slip transfer localized in the northeastern Mina deflection tends to transform from the wrench‐ to extension‐dominated transtension.

Two Cases of COVID-19 in Adults who Vape

In 2019, a severe acute respiratory disease with coupled multiple organ pathologies, namely, Coronavirus Disease 2019 (COVID-19), emerged. The lungs and other organs in this disease exhibit glass-like (hyaline) tissue pathologies with iron oxide deposits, which mimics rock-related mineralization. The impact of significant geological and geophysical perturbations on aberrant mineralization in living organisms has not been explored. Many Indigenous Knowledge systems posit that unseen forces emitted by orogenic structures and modulated by geophysical dynamics play a significant role in health. Here, it is proposed that carbon dioxide-rich water-lithospheric ferromagnesian silicates interactions generate aberrant lithospheric long-wavelength magnetic anomalies (LWMA) via serpentinization, during geologic periods with polar ice melt-induced hydrosphere redistribution, increased atmospheric carbon dioxide, and weakened geomagnetic field intensity. Here, it is also proposed that COVID-19 in humans (and animals) is a pathogenic manifestation of aberrant LWMA-induced magnetic catalysis of silicate/carbonate-like-iron oxides minerals from tissues (carbon-based molecules), gases (dioxygen and carbon dioxide) and metals (iron, magnesium, silicon), and coronavirus from the genome (endogenous viruses), with the virus capable of replication and limited transmission. These COVID-19-LWMAs are generated in tectonic plates with degraded Proterozoic bedrocks and are linked to the refertilization of iron oxides-silicate/carbonate rocks. Thus, severe COVID-19 outbreaks are/will predominately occur on the Eurasian and North American tectonic plates and governed by the diurnal, semiannual (equinoctial maxima, solstitial minima) and quasi-biannual oscillations of the geomagnetic magnetic field intensity (analogous to a vibrating drum head). Furthermore, abnormally high ferromagnetic iron content in humans is the unifying determinant for COVID-19-induced morbidity and mortality.

Focused fluid flow and methane venting along the Queen Charlotte fault, offshore Alaska (USA) and British Columbia (Canada)

AbstractFluid seepage along obliquely deforming plate boundaries can be an important indicator of crustal permeability and influence on fault-zone mechanics and hydrocarbon migration. The ∼850-km-long Queen Charlotte fault (QCF) is the dominant structure along the right-lateral transform boundary that separates the Pacific and North American tectonic plates offshore southeastern Alaska (USA) and western British Columbia (Canada). Indications for fluid seepage along the QCF margin include gas bubbles originating from the seafloor and imaged in the water column, chemosynthetic communities, precipitates of authigenic carbonates, mud volcanoes, and changes in the acoustic character of seismic reflection data. Cold seeps sampled in this study preferentially occur along the crests of ridgelines associated with uplift and folding and between submarine canyons that incise the continental slope strata. With carbonate stable carbon isotope (δ13C) values ranging from −46‰ to −3‰, there is evidence of both microbial and thermal degradation of organic matter of continental-margin sediments along the QCF. Both active and dormant venting on ridge crests indicate that the development of anticlines is a key feature along the QCF that facilitates both trapping and focused fluid flow. Geochemical analyses of methane-derived authigenic carbonates are evidence of fluid seepage along the QCF since the Last Glacial Maximum. These cold seeps sustain vibrant chemosynthetic communities such as clams and bacterial mats, providing further evidence of venting of reduced chemical fluids such as methane and sulfide along the QCF.

Eocene intra-plate shortening responsible for the rise of a faunal pathway in the northeastern Caribbean realm.

Intriguing latest Eocene land-faunal dispersals between South America and the Greater Antilles (northern Caribbean) has inspired the hypothesis of the GAARlandia (Greater Antilles Aves Ridge) land bridge. This landbridge, however, should have crossed the Caribbean oceanic plate, and the geological evolution of its rise and demise, or its geodynamic forcing, remain unknown. Here we present the results of a land-sea survey from the northeast Caribbean plate, combined with chronostratigraphic data, revealing a regional episode of mid to late Eocene, trench-normal, E-W shortening and crustal thickening by ∼25%. This shortening led to a regional late Eocene-early Oligocene hiatus in the sedimentary record revealing the location of an emerged land (the Greater Antilles-Northern Lesser Antilles, or GrANoLA, landmass), consistent with the GAARlandia hypothesis. Subsequent submergence is explained by combined trench-parallel extension and thermal relaxation following a shift of arc magmatism, expressed by a regional early Miocene transgression. We tentatively link the NE Caribbean intra-plate shortening to a well-known absolute and relative North American and Caribbean plate motion change, which may provide focus for the search of the remaining connection between 'GrANoLA' land and South America, through the Aves Ridge or Lesser Antilles island arc. Our study highlights the how regional geodynamic evolution may have driven paleogeographic change that is still reflected in current biology.

Geochronology and Geochemistry of Gabbros from Moa‐Baracoa Ophiolitic Massif, Eastern Cuba: Implication for Early Cretaceous SSZ Magmatism

AbstractA series ophiolitic massifs in the northeastern part of Cuba are generally considered to be relics of proto‐Caribbean or Caribbean lithosphere, emplaced during Paleogene collision of the NW Caribbean and North American plates (Iturralde‐Vinent et al., 2016). However, neither their age nor their geodynamic settings are well constrained. They are generally considered to be Creataceous in age but they have been interpreted as having formed in a back‐arc basin, forearc basin or plume setting (e.g. Giunta et al., 2009; Marchesi et al., 2016; Kerr et al., 1999). The Moa‐Baracoa ophiolitic massif consists of mantle harzburgite, a Moho transition zone, layered and isotropic gabbro, and pillow basalts. The isotropic gabbros, collected west of Moa (20°35′ 17.99″N, 75°02′56.58″W), are composed mainly of subhedral plagioclase and subhedral to euhedral clinopyroxene with minor primary brown amphibole. The plagioclase grains are fresh, whereas the clinopyroxene very weakly altered and partly replaced by secondary amphibole. Cathodoluminescence (CL) images of zircon grains separated from isotropic gabbro display oscillatory zoning and broad zoning, analogous to typical igneous zircon in gabbro. Ten zircon grains with oscillatory zoning in CL images yielded a mean U‐Pb SHRIMP age of 136.7±1.8 Ma. Isotropic gabbro samples have low SiO2 (47.22 wt%–49.25 wt%), TiO2 (0.13wt%–0.26 wt%) and K2O (0.04 wt%–0.08 wt%) and high MgO (8.22 wt%–16.72 wt%), Cr (1337–2065 μg/g) and Ni (192–507 μg/g) contents. Their CaO, Al2O3 and Na2O contents vary from 13.56 wt% to 15.78 wt%, 14.98 wt% to 21.98 wt%, and 0.71 wt% to 1.87 wt%, respectively. In the SiO2 versus Zr/TiO2∗0.0001 diagram, the samples plot in the sub‐alkaline basalt field but show calc‐alkaline affinities in the FeOT/MgO vs. SiO2 diagram. The gabbros contain low REE contents (ΣREE = 4.71–7.71), and their chondrite normalized REE patterns show very weak depletion in light REE (LREE, (La/Yb)N = 0.34–0.50) and flat heavy REE (HREE, (Gd/Yb)N = 1.02–1.15) profiles, with slightly positive Eu anomalies (Eu/Eu∗ = 1.55–1.76) suggesting magmatic accumulation of plagioclase. Furthermore, all samples are depleted in Nb and Th, and enriched in Ba, Pb and Sr in the primitive mantle‐normalized trace element diagram. Geochemical features, such as high values of Ba/La (30–38) and low Th/Yb (0.02–0.03), medium values of (Ta/La)PM (0.30–0.35) and (Hf/Sm)PM (0.58–0.80), indicate enrichment by slab‐derived fluids and subduction‐related metasomatism. These isotropic gabbros are geochemically comparable with those of the Kermanshah and Neyriz ophiolites of Iran (Monsef et al., 2018) that are thought to have originated in a SSZ environment. Therefore, we suggest that the isotropic gabbros in the Moa‐Baracoa massif crystallized from early Cretaceous SSZ magmas that were most likely derived from a mantle source metasomatized by slab fluids in a suprasubduction zone environment.

Petrology and U–Pb geochronology of high-grade metavolcano-sedimentary rocks from central Xolapa Complex, southern Mexico

The Xolapa Complex (XC) is a high-grade metamorphic belt exposed along the southern margin of the North American plate in Mexico. Its evolution includes an episode of widespread anatexis related to crustal thickening and subsequent uplift to upper crustal levels, which culminated in the Paleocene. On the basis of field and petrographic work, this study integrates a petrological modelling approach with LA-ICP-MS U–Pb zircon and monazite geochronology to elucidate the tectonothermal evolution of this metamorphic belt. The study area in central XC includes a suite of alternating migmatitic paragneiss and garnet-bearing mafic schist. Petrological evidence and phase equilibria calculations demonstrate that a staurolite–kyanite grade metamorphism occurred at ~640–670 °C and 8–9 kbar before extensive migmatization took place. Differential anatexis of fertile and refractory layers occurred progressively via biotite and amphibole dehydration-melting reactions, and continued during peak metamorphism at granulite-facies conditions of ~800–820 °C and ~5–7 kbar. According to Ti-in-zircon thermometry, melt generation/mobilization progressed during cooling until crystallize at ~700 °C. The age of prograde metamorphism was likely recorded by monazite cores that give mid- to Late Cretaceous ages. Cooling and melt crystallization could be constrained by leucosome zircon dated at 61.8±0.6 Ma. Both zircon and low-Y monazite from whole rock and melanosome portions, respectively, are consistent with an episode of growth around 60.9±0.5 Ma. Monazite, however, would have continued its crystallization process for a period of ca. 10 Ma. Thin zircon and high-Y monazite rims show ages between ~34–25 Ma and are interpreted to reflect a stage of reheating by magmatic advection related to large Oligocene plutonism in the region. Age and nature of central XC basement show good correlation with other Early-Cretaceous volcano-sedimentary basins of central and southern Mexico. A period of crustal thickening would have buried the sequence at depths of ~30 km, causing Barrovian-type metamorphism during mid- to Late Cretaceous time and eventually, anatexis that evolved at granulite-facies conditions during orogenic collapse in the Paleocene (≥62 Ma). The time span from leucosome crystallization to late-stage reheating in central XC would imply a protracted high-temperature evolution indicating a middle- to upper-crustal residence time of at least ca. 30 Ma.

Multicriteria seismic hazard assessment in Puerto Vallarta metropolitan area, Mexico

Puerto Vallarta, a medium-size tourist city, located in the Pacific Coast of Mexico, in a similar way as many other coastal cities, combines human activity with the potential occurrence of natural hazard events. In this way, the use of new tools to evaluate the impact of such events seems imperative. Puerto Vallarta is located within a tectonic setting where the Rivera microplate subducts beneath the North American plate and is affected by seismic activity. We performed a seismic hazard assessment by implementing a GIS-based multicriteria evaluation model. The seismic microzonation map of Puerto Vallarta was performed using a criteria set of six thematic layers, i.e., peak ground acceleration values, soil, bedrock, slope gradient, curvature, and flow accumulation. We performed the integration of the criteria set by implementing the Analytical Hierarchy Process to assign a weight to each criterion according to its contribution to the seismic hazard, i.e., PGA (0.38), soil (0.25), rock (0.14), curvature (0.10), slope (0.08), and flow accumulation (0.07). The thematic maps were integrated using GIS according to the normalized weights. We classified the seismic hazard microzonation of Puerto Vallarta into five hazard levels, i.e., low (18%), low-medium (28%), medium (22%), medium–high (20%), and high (12%). The map shows heterogeneous distribution over the territory. However, the study area can be divided into three zones, i.e., the northern mountainous area, the Ameca River Valley, and the southern mountainous area. There is an overall increment of seismic hazard from south to north. However, the highest seismic hazard levels dominate the Rio Ameca Valley showing that it is more susceptible to deposits of soft sediment and thus can be affected in the occurrence of a major earthquake. The main objective of this paper was to implement a technique to quickly estimating seismic hazards levels using available data when there is no sophisticated geophysical and engineering analysis. Using the GIS-based multicriteria techniques in seismic hazard assessment allows to elucidated areas where factors influencing surface response to earthquakes interact and raise the soil amplification susceptibility.

Mantle flow distribution beneath the California margin

Although the surface deformation of tectonic plate boundaries is well determined by geological and geodetic measurements, the pattern of flow below the lithosphere remains poorly constrained. We use the crustal velocity field of the Plate Boundary Observatory to illuminate the distribution of horizontal flow beneath the California margin. At lower-crustal and upper-mantle depths, the boundary between the Pacific and North American plates is off-centered from the San Andreas fault, concentrated in a region that encompasses the trace of nearby active faults. A major step is associated with return flow below the Eastern California Shear Zone, leading to the extrusion of the Mojave block and a re-distribution of fault activity since the Pleistocene. Major earthquakes in California have occurred above the regions of current plastic strain accumulation. Deformation is mechanically coupled from the crust to the asthenosphere, with mantle flow overlaid by a kinematically consistent network of faults in the brittle crust.

Imaging the rupture process of the 10 January 2018 MW7.5 Swan island, Honduras earthquake*

The 10 January 2018 MW7.5 Swan island, Honduras earthquake occurred on the Swan island fault, which is a transform plate boundary between the North American and Caribbean plates. Here we back-project the rupture process of the earthquake using dense seismic stations in Alaska, and find that the earthquake ruptured at least three faults (three stages) for a duration of ~40 s. The rupture speed for the longest fault (stage 3) is as fast as 5 km/s, which is much faster than the local shear wave velocity of ~4 km/s. Supershear rupture was incidentally observed on long and straight strike-slip faults. This study shows a supershear rupture that occured on a strike-slip fault with moderate length, implying that supershear rupture might commonly occur on large strike-slip earthquakes. The common occurrence of supershear rupture on strike-slip earthquakes will challenge present understanding of crack physics, as well as strong ground motion evaluation in earthquake engineering.

Plate boundary trench retreat and dextral shear drive intracontinental fault-slip histories: Neogene dextral faulting across the Gabbs Valley and Gillis Ranges, Central Walker Lane, Nevada

Abstract The spatial-temporal evolution of intracontinental faults and the forces that drive their style, orientation, and timing are central to understanding tectonic processes. Intracontinental NW-striking dextral faults in the Gabbs Valley–Gillis Ranges (hereafter referred to as the GVGR), Nevada, define a structural domain known as the eastern Central Walker Lane located east of the western margin of the North American plate. To consider how changes in boundary type along the western margin of the North American plate influenced both the initiation and continued dextral fault slip to the present day in the GVGR, we combine our new detailed geologic mapping, structural studies, and 40Ar/39Ar geochronology with published geologic maps to calculate early to middle Miocene dextral fault-slip rates. In the GVGR, Mesozoic basement is nonconformably overlain by a late Oligocene to Miocene sequence dominated by tuffs, lavas, and sedimentary rocks. These rocks are cut and offset by four primary NW-striking dextral faults, from east to west the Petrified Spring, Benton Spring, Gumdrop Hills, and Agai Pah Hills–Indian Head faults. A range of geologic markers, including tuff- and lava-filled paleovalleys, the southern extent of lava flows, and a normal fault, show average dextral offset magnitudes of 9.6 ± 1.1 km, 7.0 ± 1.7 km, 9.7 ± 1.0 km, and 4.9 ± 1.1 km across the four faults, respectively. Cumulative dextral offset across the GVGR is 31.2 ± 2.3 km. Initiation of slip along the Petrified Spring fault is tightly bracketed between 15.99 ± 0.05 Ma and 15.71 ± 0.03 Ma, whereas slip along the other faults initiated after 24.30 ± 0.05 Ma to 20.14 ± 0.26 Ma. Assuming that slip along all four faults initiated at the same time as the Petrified Spring fault yields calculated dextral fault-slip rates of 0.4 ± 0.1–0.6 ± 0.1 mm/yr, 0.4 ± 0.1–0.5 ± 0.1 mm/yr, 0.6 ± 0.1 mm/yr, and 0.3 ± 0.1 mm/yr on the four faults, respectively. Middle Miocene initiation of dextral fault slip across the GVGR overlaps with the onset of normal slip along range-bounding faults in the western Basin and Range to the north and the northern Eastern California shear zone to the south. Based on this spatial-temporal relationship, we propose that dextral fault slip across the GVGR defines a kinematic link or accommodation zone between the two regions of extension. At the time of initiation of dextral slip across the GVGR, the plate-boundary setting to the west was characterized by subduction of the Farallon plate beneath the North American plate. To account for the middle Miocene onset of extension across the Basin and Range and dextral slip in the GVGR, we hypothesize that middle Miocene trench retreat drove westward motion of the Sierra Nevada and behind it, crustal extension across the Basin and Range and NW-dextral shear within the GVGR. During the Pliocene, the plate boundary to the west changed to NW-dextral shear between the Pacific and North American plates, which drove continued dextral slip along the same faults within the GVGR because they were fortuitously aligned subparallel to plate boundary motion.

A 36-Year Record of Rock Avalanches in the Saint Elias Mountains of Alaska, With Implications for Future Hazards

Glacial retreat and mountain-permafrost degradation resulting from rising global temperatures have the potential to impact the frequency and magnitude of landslides in glaciated environments. Several recent events, including the 2015 Taan Fiord rock avalanche, which triggered a tsunami with one of the highest wave runups ever recorded, have called attention to the hazards posed by landslides in regions like southern Alaska. In the Saint Elias Mountains, the presence of weak sedimentary and metamorphic rocks and active uplift resulting from the collision of the Yakutat and North American tectonic plates create landslide-prone conditions. To differentiate between the typical frequency of landsliding resulting from the geologic and tectonic setting of this region, and landslide processes that may be accelerated due to changes in climate, we used Landsat imagery to create an inventory of rock avalanches in a 3700 km2 area of the Saint Elias Mountains. During the period from 1984-2019, we identified 220 rock avalanches with a mean recurrence interval of 60 days. We compared our landslide inventory with a catalog of M ≥ 4 earthquakes to identify potential coseismic events, but only found three possible earthquake-triggered rock avalanches. We observed a distinct temporal cluster of 41 rock avalanches from 2013 through 2016 that correlated with above average air temperatures (including the three warmest years on record in Alaska, 2014-2016); this cluster was similar to a temporal cluster of recent rock avalanches in nearby Glacier Bay National Park and Preserve. The majority of rock avalanches initiated from bedrock ridges in probable permafrost zones, suggesting that ice loss due to permafrost degradation, as opposed to glacial thinning, could be a dominant factor contributing to rock-slope failures in the high elevation areas of the Saint Elias Mountains. Although earthquake-triggered landslides have episodically occurred in southern Alaska, evidence from our study suggests that area-normalized rates of non-coseismic rock avalanches were greater during the period from 1964 to 2019, and that the frequency of these events will continue to increase as the climate continues to warm. These findings highlight the need for hazard assessments in Alaska that address changes in landslide patterns related to climate change.

Geochemical study of Cenozoic mafic volcanism in the west-central Great Basin, western Nevada, and the Ancestral Cascades Arc, California

AbstractProcesses linked to shallow subduction, slab rollback, and extension are recorded in the whole-rock major-, trace-element, and Sr, Nd, and Pb isotopic compositions of mafic magmatic rocks in both time and space over southwestern United States. Eocene to Mio-Pliocene volcanic rocks were sampled along a transect across the west-central Great Basin (GB) in Nevada to the Ancestral Cascade Arc (ACA) in the northern Sierra Nevada, California (∼39°–40° latitude), which are interpreted to represent a critical segment of a magmatic sweep that occurred as a result of subduction from east-northeast convergence between the Farallon and North American plates and extension related to the change from a convergent to a transform margin along the western edge of North America.Mafic volcanic rocks from the study area can be spatially divided into three broad regions: GB (5–35 Ma), eastern ACA, and western ACA (2.5–16 Ma). The volcanic products are dominantly calc-alkalic but transition to alkalic toward the east. Great Basin lavas erupted far inland from the continental margin and have higher K, P, Ti, and La/Sm as well as lower (Sr/P)pmn, Th/Rb, and Ba/Nb compared to ACA lavas. Higher Pb isotopic values, combined with lower Ce/Ce* and high Th/Nb ratios in some ACA lavas, are interpreted to come from slab sediment. Mafic lavas from the GB and ACA have overlapping 87Sr/86Sr and 143Nd/144Nd values that are consistent with mantle wedge melts mixing with a subduction-modified lithospheric mantle source. Eastern and western ACA lavas largely overlap in age and elemental and isotopic composition, with the exception of a small subset of lavas from the westernmost ACA region; these lavas show lower 87Sr/86Sr at a given 143Nd/144Nd. Results show that although extension contributes to melting in some regions (e.g., selected lavas in the GB and Pyramid Lake), chemical signatures for most mafic melts are dominated by subduction-related mantle wedge and a lithospheric mantle component.

Integration of tectonic geomorphology and crustal structure across the active obliquely collisional zone on the island of Hispaniola, northeastern Caribbean

Abstract Active tectonic deformation and seismicity of Hispaniola define a 250 km-wide, oblique collisional zone between the Bahamas, the island arc of Hispaniola and the Caribbean Large Igneous Province (CLIP). To reflect how collision is accommodated within Hispaniola, we calculate river normalized steepness and terrain surface roughness to reveal the areas of the most active uplift within central and western Hispaniola compared to eastern Hispaniola. We use gravity modelling to show thickness variations in the main crustal types in the obliquely convergent zone: (1) 33–45 km-thick arc crust in central and western Hispaniola; (2) 15–25 km-thick oceanic crust beneath the Bahamas north of Hispaniola; (3) 5–8 km-thick Atlantic oceanic crust NE of Hispaniola; and (4) 6–16 km-thick CLIP south of Hispaniola. Intermediate to deep earthquakes beneath eastern Hispaniola indicate active southwestward subduction of normal oceanic crust and northward subduction of the CLIP. We interpret that the west-to-east geomorphological and crustal variations within Hispaniola to be the result of an along-strike transition from crustal shortening without subduction between the Bahamas and arc crust in central and western Hispaniola to subduction of the North American and Caribbean plates beneath eastern Hispaniola. Crustal shortening in central and western Hispaniola produces thrust-fault-bounded basins with sufficient clastic sedimentary infill to produce hydrocarbon maturity.

Regionally continuous Miocene rhyolites beneath the eastern Snake River Plain reveal localized flexure at its western margin: Idaho National Laboratory and vicinity

The eastern Snake River Plain (ESRP) is a northeast-trending topographic basin interpreted to be the result of the time-transgressive track of the North American plate above the Yellowstone hotspot. The track is defined by the age progression of silicic volcanic rocks exposed along the margins of the ESRP. However, the bulk of these silicic rocks are buried under 1 to 3 kilometers of younger basalts. Here, silicic volcanic rocks recovered from boreholes that penetrate below the basalts, including INEL-1, WO-2 and new deep borehole USGS-142, are correlated with one another and to surface exposures to assess various models for ESRP subsidence. These correlations are established on U/Pb zircon and 40Ar/39Ar sanidine age determinations, phenocryst assemblages, major and trace element geochemistry, δ18O isotopic data from selected phenocrysts, and initial εHf values of zircon. These data suggest a correlation of: (1) the newly documented 8.1 ± 0.2 Ma rhyolite of Butte Quarry (sample 17KS03), exposed near Arco, Idaho to the upper-most Picabo volcanic field rhyolites found in borehole INEL-1; (2) the 6.73 ± 0.02 Ma East Arco Hills rhyolite (sample 16KS02) to the Blacktail Creek Tuff, which was also encountered at the bottom of borehole WO-2; and (3) the 6.42 ± 0.07 Ma rhyolite of borehole USGS-142 to the Walcott Tuff B encountered in deep borehole WO-2. These results show that rhyolites found along the western margin of the ESRP dip ~20º south-southeast toward the basin axis, and then gradually tilt less steeply in the subsurface as the axis is approached. This subsurface pattern of tilting is consistent with a previously proposed crustal flexural model of subsidence based only on surface exposures, but is inconsistent with subsidence models that require accommodation of ESRP subsidence on either a major normal fault or strike-slip fault.

Natural Time Analysis of Seismicity within the Mexican Flat Slab before the M7.1 Earthquake on 19 September 2017.

One of the most important subduction zones in the world is located in the Mexican Pacific Coast, where the Cocos plate inserts beneath the North American plate. One part of it is located in the Mexican Pacific Coast, where the Cocos plate inserts beneath the North American plate with different dip angles, showing important seismicity. Under the central Mexican area, such a dip angle becomes practically horizontal and such an area is known as flat slab. An earthquake of magnitude M7.1 occurred on 19 September 2017, the epicenter of which was located in this flat slab. It caused important human and material losses of urban communities including a large area of Mexico City. The seismicity recorded in the flat slab region is analyzed here in natural time from 1995 until the occurrence of this M7.1 earthquake in 2017 by studying the entropy change under time reversal and the variability β of the order parameter of seismicity as well as characterize the risk of an impending earthquake by applying the nowcasting method. The entropy change ΔS under time reversal minimizes on 21 June 2017 that is almost one week after the observation of such a minimum in the Chiapas region where a magnitude M8.2 earthquake took place on 7 September 2017 being Mexico’s largest quake in more than a century. A minimum of β was also observed during the period February–March 2017. Moreover, we show that, after the minimum of ΔS, the order parameter of seismicity starts diminishing, thus approaching gradually the critical value 0.070 around the end of August and the beginning of September 2017, which signals that a strong earthquake is anticipated shortly in the flat slab.