Geophysics and Seismology Commons^™

Subduction Initiation Recorded In The Dadeville Complex Of Alabama And Georgia, Southeastern United States, Naomi A. Becker, Freya R. George, George L. Guice, James L. Crowley, Wendy R. Nelson, Joseph F. Browning-Hanson, Supratik Roy, Daniel R. Viete

Geosciences Faculty Publications and Presentations

The Dadeville Complex of Alabama and Georgia (southeastern United States) represents the largest suite of exposed mafic-ultramafic rocks in the southern Appalachians. Due to poor preservation, chemical alteration, and tectonic reworking, a specific tectonic origin for the Dadeville Complex has been difficult to deduce. We obtained new whole-rock and mineral geochemistry coupled with zircon U-Pb geochronology to investigate the magmatic and metamorphic processes recorded by the Dadeville Complex, as well as the timing of these processes. Our data reveal an up-stratigraphic evolution in the geochemistry of the volcanic rocks, from forearc basalts to boninites. Our new U-Pb zircon crystallization data—obtained …

A Multidisciplinary Approach To Resolving The End-Guadalupian Extinction, Christopher R. Fielding, Scott E. Bryan, James L. Crowley, Tracy D. Frank, Michael T. Hren, Chris Mays, Stephen Mcloughlin, Jun Shen, Peter J. Wagner, Arne Winguth, Cornelia Winguth

Geosciences Faculty Publications and Presentations

The transition from the middle to late Permian (Guadalupian–Lopingian) is claimed to record one or more extinction events that rival the ‘Big Five’ in terms of depletion of biological diversity and reorganization of ecosystem structure. Yet many questions remain as to whether the events recorded in separate regions were synchronous, causally related, or were of a magnitude rivaling other major crises in Earth's history. In this paper, we survey some major unresolved issues related to the Guadalupian–Lopingian transition and offer a multidisciplinary approach to advance understanding of this under-appreciated biotic crisis by utilizing records in Southern Hemisphere high-palaeolatitude settings. We …

How Low Can You Go?: Widespread Challenges In Measuring Low Stream Discharge And A Path Forward, Erin C. Seybold, Anna Bergstrom, C. Nathan Jones, Amy J. Burgin, Sam Zipper, Sarah E. Godsey, Walter K. Dodds, Margaret A. Zimmer, Margaret Shanafield, Thibault Datry, Raphael D. Mazor, Mathis L. Messager, Julian D. Olden, Adam Ward, Songyan Yu, Kendra E. Kaiser, Arial Shogren, Richard H. Walker

Geosciences Faculty Publications and Presentations

Low flows pose unique challenges for accurately quantifying streamflow. Current field methods are not optimized to measure these conditions, which in turn, limits research and management. In this essay, we argue that the lack of methods for measuring low streamflow is a fundamental challenge that must be addressed to ensure sustainable water management now and into the future, particularly as climate change shifts more streams to increasingly frequent low flows. We demonstrate the pervasive challenge of measuring low flows, present a decision support tool (DST) for navigating best practices in measuring low flows, and highlight important method developmental needs.

Geophysics-Informed Hydrologic Modeling Of A Mountain Headwater Catchment For Studying Hydrological Partitioning In The Critical Zone, Hang Chen, Qifei Niu, Aida Mendieta, John Bradford, James Mcnamara

Geosciences Faculty Publications and Presentations

Hydrologic modeling has been a useful approach for analyzing water partitioning in catchment systems. It will play an essential role in studying the responses of watersheds under projected climate changes. Numerous studies have shown it is critical to include subsurface heterogeneity in the hydrologic modeling to correctly simulate various water fluxes and processes in the hydrologic system. In this study, we test the idea of incorporating geophysics-obtained subsurface critical zone (CZ) structures in the hydrologic modeling of a mountainous headwater catchment. The CZ structure is extracted from a three-dimensional seismic velocity model developed from a series of two-dimensional velocity sections …

Electrofacies Analysis Using A Geostatistical Approach, Northern Iraq Case Study, Hussein S. Hussein

Polytechnic Journal

The distribution of petrophysical parameters is governed by lithology, hence understanding the spatial variation in lithology is essential for reservoir characterisation. This study points out the use of well logs to estimate the lithology of carbonate rocks (limestone, marly limestone, dolomite, and dolomitic limestone) found in Upper Cretaceous to Miocene formations at the Khabbaz Oil Field in Northern Iraq. Applying the multivariate regression technique to neutron, and density data enabled accurate lithology prediction. By using independent values (well log data) to predict a dependent value (lithology), this technique is a progression of regression analysis. The formations under study comprised a …

Revealing The Core Dynamics Of Vesta: Insights From Experimental Investigations Of Electrical Resistivity And Thermal Conductivity, Oluwasanmi A. Orole Mr

Electronic Thesis and Dissertation Repository

Insights from high pressure and temperature experiments involving in-situ measurements of the electrical resistivity of Fe-5wt%Ni at temperatures of up to 2000 K, under pressures of 2-5 GPa in a 1000-ton cubic-anvil press have been used to reveal Vesta’s core dynamics. The Wiedemann–Franz law was used to calculate the thermal conductivity from the measured electrical resistivity data. Comparing the findings of this study with prior investigations on both pure Fe and Fe-10wt%Ni indicates that an increase in Ni ranging from 0-10wt% has negligible effect on the electrical resistivity of Fe alloys. By comparing the range of estimated heat flux through …

Evaluating 3 Decades Of Precipitation In The Upper Colorado River Basin From A High-Resolution Regional Climate Model, William Rudisill, Alejandro Flores, Rosemary Carroll

Geosciences Faculty Publications and Presentations

Convection-permitting regional climate models (RCMs) have recently become tractable for applications at multi-decadal timescales. These types of models have tremendous utility for water resource studies, but better characterization of precipitation biases is needed, particularly for water-resource-critical mountain regions, where precipitation is highly variable in space, observations are sparse, and the societal water need is great. This study examines 34 years (1987–2020) of RCM precipitation from the Weather Research and Forecasting model (WRF; v3.8.1), using the Climate Forecast System Reanalysis (CFS; CFSv2) initial and lateral boundary conditions and a 1 km × 1 km innermost grid spacing. The RCM is centered …

Astronomical And Tectonic Influences On Climate And Deposition Revealed Through Radioisotopic Geochronology And Bayesian Age-Depth Modeling Of The Early Eocene Green River Formation, Wyoming, Usa, Benjamin T. Bruck, Brad S. Singer, Mark D. Schmitz, Alan R. Carroll, Stephen Meyers, Andrew P. Walters, Brian R. Jicha

Geosciences Faculty Publications and Presentations

The Wilkins Peak Member (WPM) of the Green River Formation in Wyoming, USA, comprises alternating lacustrine and alluvial strata that preserve a record of terrestrial climate during the early Eocene climatic optimum. We use a Bayesian framework to develop age-depth models for three sites, based on new ⁴⁰Ar/³⁹Ar sanidine and ²⁰⁶Pb/²³⁸U zircon ages from seven tuffs. The new models provide two- to ten-fold increases in temporal resolution compared to previous radioisotopic age models, confirming eccentricity-scale pacing of WPM facies, and permitting their direct comparison to astronomical solutions. Starting at ca. 51 Ma, the median …

U–Pb Zircon Ages, Mapping, And Biostratigraphy Of The Payette Formation And Idaho Group North Of The Western Snake River Plain, Idaho: Implications For Hydrocarbon System Correlation, Renee L. Love, Reed S. Lewis, Spencer H. Wood, Dennis M. Feeney, Mark D. Schmitz

Geosciences Faculty Publications and Presentations

Sedimentary deposits north of the western Snake River Plain host Idaho’s first and only producing oil and gas field. They consist of the lower to middle Miocene Payette Formation, the middle to upper Miocene Poison Creek and Chalk Hills Formations, and the Pliocene to lower Pleistocene Glenns Ferry Formation. Using new geochronology, palynomorph biostratigraphy, and geologic mapping, we connect updip surface features to subsurface petroleum play elements. The Payette Formation is a likely main source of the hydrocarbons, and acts as one of the reservoirs in the unnamed basin. Here, we redefine the Payette Formation as 0 to ~3,500 ft …

Empirical Nonlinear Site Response Applicable To Greater Vancouver, British Columbia, Natalia Gomez Jaramillo

Electronic Thesis and Dissertation Repository

In Greater Vancouver, strong earthquake shaking has not been recorded to quantify expected nonlinear soil response. Earthquake horizontal to vertical spectral ratios (eHVSR) are evaluated at seismic stations located at the edges or deeper center of sedimentary basins in Mexico City, the Kanto basin and Anchorage with available weak-to-strong motion recordings that have similar linear site response to Greater Vancouver to constrain the expected nonlinear soil response for future large earthquakes affecting Greater Vancouver. The empirical results confirm that the very low fundamental frequency (f₀ < 0.6 Hz) at deep Fraser River delta sites will be less affected by strong shakings and suggest that for Fraser River delta edge sites the amplification at the predominant frequency (f₁) could increase up to 1.77 times with a strong event. This thesis’ …

Reducing Uncertainty In Sea-Level Rise Prediction: A Spatial-Variability-Aware Approach, Subhankar Ghosh, Shuai An, Arun Sharma, Jayant Gupta, Shashi Shekhar, Aneesh Subramanian

I-GUIDE Forum

Given multi-model ensemble climate projections, the goal is to accurately and reliably predict future sea-level rise while lowering the uncertainty. This problem is important because sea-level rise affects millions of people in coastal communities and beyond due to climate change's impacts on polar ice sheets and the ocean. This problem is challenging due to spatial variability and unknowns such as possible tipping points (e.g., collapse of Greenland or West Antarctic ice-shelf), climate feedback loops (e.g., clouds, permafrost thawing), future policy decisions, and human actions. Most existing climate modeling approaches use the same set of weights globally, during either regression or …

Seabed Reflection Coefficient Variability In The Mono Lake Using Seismic Reflection Sparker Data, Francis Nnamdi Okeifufe

Earth Sciences Theses and Dissertations

Acoustic reflections from a lakebed provide valuable information about the dynamic boundary between water and the lakebed, where sediments and biota in the water column eventually settle. In this study, we present a comprehensive analysis of reflected sound wave amplitudes to gain insights into the physical properties of a water-sediment boundary. We decided to carry out this analysis within the unique ecosystem of the Mono Lake basin due to the fact that the Mono Lake, a meromictic and hypersaline waterbody is situated within one of the most active yet understudied volcanic regions in the United States. The fact that the …

Geophysical Survey In Central North Dakota, University Of North Dakota. Energy And Environmental Research Center

EERC Brochures and Fact Sheets

Fact sheet about a 2023–2024 Energy & Environmental Research Center (EERC) geophysical survey in central North Dakota to evaluate the feasibility of developing safe and permanent CO₂ storage in the project area. Describes the purpose of the survey and potential impact on local landowners.

The Influence Of Volcano Topographic Changes On Infrasound Amplitude: Lava Fountains At Mt. Etna In 2021, Adriana Iozzia, Leighton M. Watson, Massimo Cantarero, Emanuela De Beni, Giuseppe Di Grazia, Gaetana Ganci, Jeffrey B. Johnson, Eugenio Privitera, Cristina Proietti, Mariangela Sciotto, Andrea Cannata

Geosciences Faculty Publications and Presentations

Infrasound signals are used to investigate and monitor active volcanoes during eruptive and degassing activity. Infrasound amplitude information has been used to estimate eruptive parameters such as plume height, magma discharge rate, and lava fountain height. Active volcanoes are characterized by pronounced topography and, during eruptive activity, the topography can change rapidly, affecting the observed infrasound amplitudes. While the interaction of infrasonic signals with topography has been widely investigated over the past decade, there has been limited work on the impact of changing topography on the infrasonic amplitudes. In this work, the infrasonic signals accompanying 57 lava fountain paroxysms at …

Assessing Controls On Ice Dynamics At Crane Glacier, Antarctic Peninsula, Using A Numerical Ice Flow Model, Rainey Aberle, Ellyn M. Enderlin, Hans-Peter Marshall, Michal Kopera, Tate G. Meehan

Geosciences Faculty Publications and Presentations

The Antarctic Peninsula's widespread glacier retreat and ice shelf collapse have been attributed to atmospheric and oceanic warming. Following the initial post-collapse period of retreat, several former tributary glaciers of the Larsen A and B ice shelves have been slowly re-advancing for more than a decade. Here, we use a flowline model of Crane Glacier to gauge the sensitivity of former tributary glaciers to future climate change following this period of long-term dynamic adjustment. The glacier's long-term geometry and speed changes are similar to those of other former Larsen A and B tributaries, suggesting that Crane Glacier is a reasonable …

Early Ordovician Seamounts Preserved In The Canadian Cordillera: Implications For The Rift History Of Western Laurentia, Rose N. Cobbett, Luke P. Beranek, Stephen J. Piercey, James L. Crowley, Maurice Colpron

Geosciences Faculty Publications and Presentations

The breakup of the supercontinent Rodinia and development of the western Laurentian rifted margin are in part recorded by Neoproterozoic to mid-Paleozoic igneous and sedimentary rock successions in the Canadian Cordillera. New bedrock mapping and volcanic facies analysis of Early Ordovician mafic rocks assigned to the Menzie Creek Formation in central Yukon allow reconstruction of the depositional environment during the volcanic eruptions, whole-rock geochemical data constrain the melting depth and crust-mantle source regions of the igneous rocks within the study area, and zircon U-Pb age studies provide determination of the precise timing of submarine eruptions. Menzie Creek Formation volcanic rocks …

Snowpack Relative Permittivity And Density Derived From Near-Coincident Lidar And Ground-Penetrating Radar, Randall Bonnell, Daniel Mcgrath, Andrew R. Hedrick, Ernesto Trujillo, Tate G. Meehan, Keith Williams, Hans-Peter Marshall, Graham Sexstone, John Fulton, Michael J. Ronayne, Steven R. Fassnacht, Ryan Webb, Katherine E. Hale

Geosciences Faculty Publications and Presentations

Depth-based and radar-based remote sensing methods (e.g., lidar, synthetic aperture radar) are promising approaches for remotely measuring snow water equivalent (SWE) at high spatial resolution. These approaches require snow density estimates, obtained from in-situ measurements or density models, to calculate SWE. However, in-situ measurements are operationally limited, and few density models have seen extensive evaluation. Here, we combine near-coincident, lidar-measured snow depths with ground-penetrating radar (GPR) two-way travel times (twt) of snowpack thickness to derive >20 km of relative permittivity estimates from nine dry and two wet snow surveys at Grand Mesa, Cameron Pass, and Ranch Creek, Colorado. …

Application Of Model-Based Time Series Prediction Of Infrared Long-Wave Radiation Data For Exploring The Precursory Patterns Associated With The 2021 Madoi Earthquake, Jingye Zhang, Ke Sun, Junqing Zhu, Ning Mao, Dimitar Ouzounov

Mathematics, Physics, and Computer Science Faculty Articles and Research

Taking the Madoi M_S 7.4 earthquake of 21 May 2021 as an example, this paper proposes using time series prediction models to predict the outgoing long-wave radiation (OLR) anomalies and study short-term pre-earthquake signals. Five time series prediction models, including autoregressive integrated moving average (ARIMA) and long short-term memory (LSTM), were trained with the OLR time series data of the aseismic moments in the 5° × 5° spatial range around the epicenter. The model with the highest prediction accuracy was selected to retrospectively predict the OLR values during the aseismic period and before the earthquake in the area. It …

Archean (3.3 Ga) Paleosols And Paleoenvironments Of Western Australia, Gregory J. Retallack, Mark D. Schmitz

Geosciences Faculty Publications and Presentations

The Pilbara craton of northwestern Australia is known for what were, when reported, the oldest known microfossils and paleosols on Earth. Both interpretations are mired in controversy, and neither remain the oldest known. Both the microfossils and the paleosols have been considered hydrothermal artefacts: carbon films of vents and a large hydrothermal cupola, respectively. This study resampled and analyzed putative paleosols within and below the Strelley Pool Formation (3.3 Ga), at four classic locations: Strelley Pool, Steer Ridge, Trendall Ridge, and Streckfuss, and also at newly discovered outcrops near Marble Bar. The same sequence of sedimentary facies and paleosols was …

Two New Absolute Gravity Base Stations Established In Lincoln, Ne, Kris Guthrie, Irina Filina

Transactions of the Nebraska Academy of Sciences and Affiliated Societies

Two new absolute gravity base stations were established on the University of Nebraska-Lincoln City Campus in April 2022. The values were transferred with a relative gravity meter from four existing stations within 238 miles and were later validated by the National Oceanic and Atmospheric Administration. The final values at 22.5 cm above marker height are 980177.781 ± 0.010 mGal at the Lincoln Bell Tower Station and 980178.079 ± 0.010 mGal at the Lincoln Museum station with a gravity gradient of -3.105 µGal/cm. This paper describes the process of transferring gravity values and outlines recommendations for similar projects in the future.

Constraints On Near-Ridge Magmatism Using 40Ar/39Ar Geochronology Of Enriched Morb From The 8°20' N Seamount Chain, Molly K. Anderson, Michael R. Perfit, Leah E. Morgan, Daniel J. Fornari, Michael Cosca, V. Dorsey Wanless

Geosciences Faculty Publications and Presentations

Our understanding of the spatial-temporal-compositional relationships between off-axis magmatism and mid-ocean ridge spreading centers is limited. Determining the ⁴⁰Ar/³⁹Ar ages of mid-ocean ridge basalt (MORB) lavas erupting near mid-ocean ridges (MOR) has been a challenge due to the characteristically low K₂O contents in incompatible element-depleted normal MORB (NMORB). High-precision ⁴⁰Ar/³⁹Ar geochronology is used here to determine ages of young, basaltic lavas erupted along the 8°20' N seamount chain west of the East Pacific Rise (EPR) axis that have a range of incompatible element enrichments (EMORB) suitable for ⁴⁰Ar/³⁹Ar geochronology …

U-Pb Zircon Dates From North American And British Avalonia Bracket The Lower–Middle Cambrian Boundary Interval, With Evaluation Of The Miaolingian Series As A Global Unit, Ed Landing, Mark D. Schmitz, Stephen R. Westrop, Gerd Geyer

Geosciences Faculty Publications and Presentations

High-precision U-Pb zircon ages on SE Newfoundland tuffs now bracket the Avalonian Lower–Middle Cambrian boundary. Upper Lower Cambrian Brigus Formation tuffs yield depositional ages of 507.91 ± 0.07 Ma (Callavia broeggeri Zone) and 507.67 ± 0.08 Ma and 507.21 ± 0.13 Ma (Morocconus-Condylopyge eli Assemblage interval). Lower Middle Cambrian Chamberlain’s Brook Formation tuffs have depositional ages of 506.34 ± 0.21 Ma (Kiskinella cristata Zone) and 506.25 ± 0.07 Ma (Eccaparadoxides bennetti Zone). The composite unconformity separating the Brigus and Chamberlain’s Brook formations is constrained between these ages. An Avalonian Lower–Middle Cambrian boundary between 507.2 ± 0.1 …

Is Macroporosity Controlled By Complexed Clay And Soil Organic Carbon?, Aaron N. Koop, Daniel R. Hirmas, Sharon A. Billings, Li Li, Alejandro Cueva, Xi Zhang, Hang Wen, Attila Nemes, Lígia F. T. Souza, Hoori Ajami, Alejandro N. Flores, Aoesta K. Rudick, Annalise Guthrie, Lola M. Klamm, Micah Unruh, Pamela L. Sullivan

Geosciences Faculty Publications and Presentations

Multi-scale evidence of rapid, climate-induced soil structural changes occurring at yearly to decadal timescales is mounting. As a result, it has become increasingly important to identify the properties and mechanisms controlling the development and maintenance of soil structure and associated macroporosity. This is especially relevant since macroporosity has disproportionate effects on saturated hydraulic conductivity (K_sat) which strongly influences water storage and flux, thus, affecting the water cycle. In this study, we use decision trees and piecewise linear regression to assess the influence of soil and climate properties on effective porosity (EP; a proxy of macroporosity) in both …

Applications Of Observational Seismology: Insights Into Volcanic And Near Surface Processes, Justin T. Wilgus

Earth and Planetary Sciences ETDs

The field of observational seismology has made tremendous progress in the past two decades. This progress has been multi-faceted in form, but significant contributions emanated from 1) increases in both the quality and the quantity of seismic data 2) advances in computational power 3) advances in algorithmic capability, including machine learning. In this dissertation I report on three distinctly different seismic applications made possible by the aforementioned progress and discuss the insights these applications have provided in understanding volcanic and near surface processes of the Earth.

In the first chapter titled, “Shear Velocity Evidence of Upper Crustal Magma Storage Beneath …

Petrophysical Properties Of The Hed Meteorites And The Bushveld Igneous Complex: Implications For Future Asteroid Prospecting Missions, Lukas Louwerse

Electronic Thesis and Dissertation Repository

In this work, I make a quantitative and mineralogical comparison of Howardite-EuctriteDiogenite (HED) meteorites with terrestrial samples from the Bushveld Igneous Complex (BIC) to assess how geophysical surveys of 4 Vesta or of V-type asteroids could detect BIClike ore deposits. Petrophysical methods in this work include bulk and grain density, derived porosity (using helium pycnometry, the Archimedean method, and X-ray micro-computed tomography), magnetic susceptibility, induced polarization and resistivity, and P-wave velocity. Mineralogical methods in this work include optical microscopy, micro X-ray diffraction, micro X-ray fluorescence element mapping, and electron beam methods. Howardite NWA 15199 was classified as part of this …

Partial Auto-Correlation Of Low Magnitude Earthquakes From The 2016 Iris Aray In Grant County, Oklahoma, Alex Christopher Eddy

Open Access Theses & Dissertations

The Anadarko Basin in is a seismically active basin that spans central and western Oklahoma into Kansas. In June and July of 2016, the IRIS Community Wavefield Experiment array was deployed for approximately a month, recording over 300 earthquakes. Using 4 of these earthquakes and an east â?? west array of 129 Fairfield Nodal Z-land 3C nodes, I present a novel application of correlation. These four earthquakes with a ML of 2.8-3.0 represent a varied azimuthal distribution and are between 9 and 40km from the array and have a frequency peak at approximately 10 Hz.

As a proof of concept, …

Seismic Analysis Of The Upper & Lower Falls Of The Yellowstone River, Loring Schaible

Earth and Planetary Sciences ETDs

Twelve years of concurrent hydrologic and continuous seismic data along with temporary seismic data demonstrate that the Upper and Lower Falls of the Grand Canyon of the Yellowstone River comprise a highly localized source of 0.5-5 Hz seismic energy that overwhelms anthropogenic contributions. In aggregate, seismic amplitude from 2008-2019 is linearly related to discharge with a correlation coefficient of 0.96. Repeated deviations from this linear relationship persist for 1-2 weeks prior to the date that Yellowstone Lake becomes clear of winter ice coverage. Seismic efficiency increases by ~50-250% during this period of ice-breakup, during which lake ice flows into the …

Seismic Shadow Zone Investigation In The Upper Magma Reservoir Of The Yellowstone Caldera, Sarah Nolt-Caraway

Earth and Planetary Sciences ETDs

This study investigates whether the Yellowstone Caldera has enough melt to mute S-waves, creating a seismic shadow zone. Using a dense nodal deployment of ~650 stations, 7-9 earthquakes during the nodal deployment, and 21 broadband stations with 3,000-4,000 events per station; amplitude and noise maps, seismograms, and automatic phase picks probabilities from a deep learning model were analyzed to assess the potential role of melt, crustal attenuation, and noise in affecting body waves, particularly S-phases. The results are inconclusive, with unclear evidence whether observed amplitude decay is normal signal decay due to distance, noise-related, melt, or from scattering and intrinsic …

Evaluating The Feasibility Of Using Strain Measured During Sinusoidal Rate Pumping Tests To Characterize An Aquifer, Riley Blais

All Theses

Pumping tests with sinusoidal variation in pumping rate have been proposed as a method for improving aquifer characterization. These tests can interrogate a larger aquifer volume than slug tests and they can be more sensitive to small variations in drawdown. Current methods of using sinusoidal variations of rate are based on measuring pressure signals from the reservoir or aquifer, which requires access to monitoring wells. An alternative approach has been developed that measures the strain in the vadose zone instead of pressure in the reservoir. An instrument has been developed at Clemson University that can measure small strains using optical …

The Temporal Evolution Of Subduction Initiation In The Samail Ophiolite: High-Precision U–Pb Zircon Petrochronology Of The Metamorphic Sole, Matthew Rioux, Joshua M. Garber, Michael Searle, James L. Crowley, Sally Stevens, Mark Schmitz, Andrew Kylander-Clark, Kayla Leal, Tyler Ambrose, Andrew J. Smye

Geosciences Faculty Publications and Presentations

High-precision dating of the metamorphic sole of ophiolites can provide insight into the tectonic evolution of ophiolites and subduction zone processes. To understand subduction initiation beneath a young, well-preserved and well-characterized ophiolite, we performed coupled zircon laser-ablation inductively coupled mass spectrometry trace element analyses and high-precision isotope dilution-thermal ionization mass spectrometry U–Pb dating on 25 samples from the metamorphic sole of the Samail ophiolite (Oman-United Arab Emirates). Zircon grains from amphibolite- to granulite-facies (0.8–1.3 GPa, ~700–900°C), garnet- and clinopyroxene-bearing amphibolite samples (n = 18) show systematic trends of decreasing heavy rare earth element slope (HREE; Yb/Dy) with decreasing Yb …