skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An optimization tool for identifying multiple-diffusion domain model parameters
Abstract. The multiple-diffusion domain (MDD) model empirically describes the diffusive behavior of noble gases in some terrestrial materials and has been commonly used to interpret 40Ar/39Ar stepwise degassing observations in K-feldspar. When applied in this manner, the MDD model can be used to test crustal exhumation scenarios by identifying the permissible thermal paths a rock sample could have undergone over geologic time, assuming the diffusive properties of Ar within the mineral are accurately understood. More generally, the MDD model provides a framework for quantifying the temperature-dependent diffusivity of noble gases in minerals. However, constraining MDD parameters that successfully predict the results of step-heating diffusion experiments is a complex task, and the assumptions made by existing numerical methods used to quantify model parameters can bias the absolute temperatures permitted by thermal modeling. For example, the most commonly used method assumes that no domains lose more than 60 % of their gas during early heating steps (Lovera et al., 1997). This assumption is unverifiable, and we show that the Lovera et al. (1997) procedure may bias predicted temperatures towards lower values when it is violated. To address this potential bias and to provide greater accessibility to the MDD model, we present a new open-source method for constraining MDD parameters from stepwise degassing experimental results, called the “MDD Tool Kit” (https://github.com/dgorin1/mddtoolkit, last access: 11 October 2024). This software optimizes all MDD parameters simultaneously and removes any need for user-defined Ea or regression fitting choices used by other tools. In doing so, this new method eliminates assumptions about the domain size distribution. To test the validity of our thermal predictions, we then use the MDD Tool Kit (https://github.com/dgorin1/mddtoolkit) to interpret 40Ar/39Ar results from the Grayback Fault, AZ, USA. Although the resulting thermal histories are consistently ∼ 60–75 °C higher than those found in previous studies, they agree with independent observations from apatite fission track, zircon fission track, and (U-Th)/He.  more » « less
Award ID(s):
1935907
PAR ID:
10560781
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
European Geophysical Union / Copernicus Publications
Date Published:
Journal Name:
Geochronology
Volume:
6
Issue:
4
ISSN:
2628-3719
Page Range / eLocation ID:
521 to 540
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The 40Ar/39Ar dating method is among the most versatile of geochronometers, having the potential to date a broad variety of K-bearing materials spanning from the time of Earth’s formation into the historical realm. Measurements using modern noble-gas mass spectrometers are now producing 40Ar/39Ar dates with analytical uncertainties of ∼0.1%, thereby providing precise time constraints for a wide range of geologic and extraterrestrial processes. Analyses of increasingly smaller subsamples have revealed age dispersion in many materials, including some minerals used as neutron fluence monitors. Accordingly, interpretive strategies are evolving to address observed dispersion in dates from a single sample. Moreover, inferring a geologically meaningful “age” from a measured “date” or set of dates is dependent on the geological problem being addressed and the salient assumptions associated with each set of data. We highlight requirements for collateral information that will better constrain the interpretation of 40Ar/39Ar data sets, including those associated with single-crystal fusion analyses, incremental heating experiments, and in situ analyses of microsampled domains. To ensure the utility and viability of published results, we emphasize previous recommendations for reporting 40Ar/39Ar data and the related essential metadata, with the amendment that data conform to evolving standards of being findable, accessible, interoperable, and reusable (FAIR) by both humans and computers. Our examples provide guidance for the presentation and interpretation of 40Ar/39Ar dates to maximize their interdisciplinary usage, reproducibility, and longevity. 
    more » « less
  2. The Mesozoic–Cenozoic convergent margin history of southern Alaska has been dominated by arc magmatism, terrane accretion, strike-slip fault systems, and possible spreading-ridge subduction. We apply 40Ar/39Ar, apatite fission-track (AFT), and apatite (U-Th)/He (AHe) geochronology and thermochronology to plutonic and volcanic rocks in the southern Talkeetna Mountains of Alaska to document regional magmatism, rock cooling, and inferred exhumation patterns as proxies for the region’s deformation history and to better delineate the overall tectonic history of southern Alaska. High-temperature 40Ar/39Ar thermochronology on muscovite, biotite, and K-feldspar from Jurassic granitoids indicates postemplacement (ca. 158–125 Ma) cooling and Paleocene (ca. 61 Ma) thermal resetting. 40Ar/39Ar whole-rock volcanic ages and 45 AFT cooling ages in the southern Talkeetna Mountains are predominantly Paleocene–Eocene, suggesting that the mountain range has a component of paleotopography that formed during an earlier tectonic setting. Miocene AHe cooling ages within ~10 km of the Castle Mountain fault suggest ~2–3 km of vertical displacement and that the Castle Mountain fault also contributed to topographic development in the Talkeetna Mountains, likely in response to the flat-slab subduction of the Yakutat microplate. Paleocene–Eocene volcanic and exhumation-related cooling ages across southern Alaska north of the Border Ranges fault system are similar and show no S-N or W-E progressions, suggesting a broadly synchronous and widespread volcanic and exhumation event that conflicts with the proposed diachronous subduction of an active west-east–sweeping spreading ridge beneath south-central Alaska. To reconcile this, we propose a new model for the Cenozoic tectonic evolution of southern Alaska. We infer that subparallel to the trench slab breakoff initiated at ca. 60 Ma and led to exhumation, and rock cooling synchronously across south-central Alaska, played a primary role in the development of the southern Talkeetna Mountains, and was potentially followed by a period of southern Alaska transform margin tectonics. 
    more » « less
  3. Abstract Our study used zircon (U-Th)/He (ZHe) thermochronology to resolve cooling events of Precambrian basement below the Great Unconformity surface in the eastern Grand Canyon, United States. We combined new ZHe data with previous thermochronometric results to model the <250 °C thermal history of Precambrian basement over the past >1 Ga. Inverse models of ZHe date-effective uranium (eU) concentration, a relative measure of radiation damage that influences closure temperature, utilize He diffusion and damage annealing and suggest that the main phase of Precambrian cooling to <200 °C was between 1300 and 1250 Ma. This result agrees with mica and potassium feldspar 40Ar/39Ar thermochronology showing rapid post–1400 Ma cooling, and both are consistent with the 1255 Ma depositional age for the Unkar Group. At the young end of the timescale, our data and models are also highly sensitive to late-stage reheating due to burial beneath ∼3–4 km of Phanerozoic strata prior to ca. 60 Ma; models that best match observed date-eU trends show maximum temperatures of 140–160 °C, in agreement with apatite (U-Th)/He and fission-track data. Inverse models also support multi-stage Cenozoic cooling, with post–20 Ma cooling from ∼80 to 20 °C reflecting partial carving of the eastern Grand Canyon, and late rapid cooling indicated by 3–7 Ma ZHe dates over a wide range of high eU. Our ZHe data resolve major basement exhumation below the Great Unconformity during the Mesoproterozoic (1300–1250 Ma), and “young” (20–0 Ma) carving of Grand Canyon, but show little sensitivity to Neoproterozoic and Cambrian basement unroofing components of the composite Great Unconformity. 
    more » « less
  4. Deep-time thermochronology by the zircon (U-Th)/He (ZHe) method is an emerging field of study with promise for constraining Precambrian rock thermal and exhumation histories. The Grand Canyon provides an opportunity to further explore this method because excellent geologic constraints can be integrated with multiple thermochronometers to address important questions about the spatial variability of basement erosion below the sub-Cambrian Great Unconformity composite erosional surface. In this study, we synthesize new ZHe results (n = 26) and published (n = 77) ZHe data with new K-feldspar 40Ar/39Ar data and models (n = 4) from Precambrian basement rocks of the Grand Canyon, USA. We use HeFTy and QTQt thermal history modeling to evaluate the ability of the individual ZHe and K-feldspar 40Ar/39Ar thermochronometric data sets to resolve Precambrian thermal histories and compare those results with jointly modeled data using the QTQt software. We also compare Precambrian basement thermal histories of the eastern and western Grand Canyon, where the eastern Grand Canyon has ∼4 km of Grand Canyon Supergroup strata deposited and preserved, and the western Grand Canyon, where the Supergroup was either never deposited or not preserved. In all locations, models constrained only by ZHe data have limited resolving power for the past ∼600 m.y., compared to models that combine K-feldspar 40Ar/39Ar and ZHe data, which extends the recorded history into the Mesoproterozoic. Our model results suggest that two regional basement unroofing events occurred. A ca. 1350−1250 Ma cooling event is interpreted to record basement exhumation from depths of ∼10 km, and a second cooling episode (∼200−100 °C total) records exhumation from a depth of ∼3 km to 7 km to near-surface conditions between ca. 600 Ma and 500 Ma. Easternmost Grand Canyon models suggest that the preserved maximum ∼4 km thickness of the Grand Canyon Supergroup (with burial heating at ∼100 °C) approximates the total original Mesoproterozoic and Neoproterozoic stratal thickness. Whether these Supergroup rocks were present and then eroded in the western Grand Canyon, as suggested by regional geologic studies, or were never deposited is not constrained by thermochronological data. 
    more » « less
  5. Measurements of cosmic-ray-produced beryllium-10, neon-21, and helium-3 in quartz in a soil profile from a forested landscape in the Oregon Coast Range show that the cosmogenic noble gases 21Ne and 3He are depleted relative to 10Be in the shallow subsurface. The noble gases are mobile in mineral grains via thermally activated diffusion and 10Be is not, implying that noble gas depletion is the result of surface heating by wildfires and subsequent mixing of partially degassed quartz downward into the soil. Cosmogenic noble gas depletion by wildfire heating of soils is a potential means of estimating wildfire intensity and/or frequency over pre-observational timescales. 
    more » « less