Search for: All records

We undertook Zr isotope measurements on zircon, titanite, biotite, amphibole, and whole rocks from the La Posta pluton (Peninsular Ranges, southern California) together with trace element analyses and U-Pb age measurements to understand the controls on Zr isotope fractionation in igneous rocks, including temperature, crystallization sequence, and kinetic effects. We find large (>0.6‰) Zr isotope fractionations (expressed as δ94/90Zr) between titanite and zircon forming at approximately the same temperature. Using equilibrium fractionation factors calculated from ionic and ab initio models, we infer the controls on Zr isotope evolution to include the relative order in which phases appear on the liquidus, with titanite fractionation resulting in isotopically lighter melt and zircon fractionation resulting in isotopically heavier melt. While these models of Zr fractionation can explain δ94/90Zr variations in zircon of up to ∼1.5‰, crystallization order, temperature and presence of co-crystallizing phases do not explain all aspects of the intracrystalline Zr isotopic distribution in zircons in the La Posta pluton or the large range of Zr isotopic values among zircons (>2‰). Without additional constraints, such as knowledge of co-crystallizing phases and a better understand of the true causes of Zr isotope fractionation, Zr isotopes in zircon remains an ambiguous proxy of magmatic evolution. 

Semi-competing risks refer to the time-to-event analysis setting, where the occurrence of a non-terminal event is subject to whether a terminal event has occurred, but not vice versa. Semi-competing risks arise in a broad range of clinical contexts, including studies of preeclampsia, a condition that may arise during pregnancy and for which delivery is a terminal event. Models that acknowledge semi-competing risks enable investigation of relationships between covariates and the joint timing of the outcomes, but methods for model selection and prediction of semi-competing risks in high dimensions are lacking. Moreover, in such settings researchers commonly analyze only a single or composite outcome, losing valuable information and limiting clinical utility—in the obstetric setting, this means ignoring valuable insight into timing of delivery after preeclampsia has onset. To address this gap, we propose a novel penalized estimation framework for frailty-based illness–death multi-state modeling of semi-competing risks. Our approach combines non-convex and structured fusion penalization, inducing global sparsity as well as parsimony across submodels. We perform estimation and model selection via a pathwise routine for non-convex optimization, and prove statistical error rate results in this setting. We present a simulation study investigating estimation error and model selection performance, and a comprehensive application of the method to joint risk modeling of preeclampsia and timing of delivery using pregnancy data from an electronic health record.

 

Tectonic deformation can influence spatiotemporal patterns of erosion by changing both base level and the mechanical state of bedrock. Although base-level change and the resulting erosion are well understood, the impact of tectonic damage on bedrock erodibility has rarely been quantified. Eastern Tibet, a tectonically active region with diverse lithologies and multiple active fault zones, provides a suitable field site to understand how tectonic deformation controls erosion and topography. In this study, we quantified erosion coefficients using the relationship between millennial erosion rates and the corresponding channel steepness. Our work shows a twofold increase in erosion coefficients between basins within 15 km of major faults compared to those beyond 15 km, suggesting that tectonic deformation through seismic shaking and rock damage significantly affects eastern Tibet erosion and topography. This work demonstrates a field-based, quantitative relationship between rock erodibility and fault damage, which has important implications for improving landscape evolution models. 

Stability of a solitary wave disturbed by a submerged flat sill is investigated experimentally. For sills narrow compared with the solitary wave, the transmitted waves are found to be unaffected in waveform and amplitude. A wider sill disturbs the solitary wave resulting in the formation of a dispersive wavetrain following the transmitted wave. In some cases, the wave amplitude recovers, despite being perturbed, to the state of an unobstructed solitary-wave state at a certain distance beyond the sill. Wider sills cause wave breaking that occurs over the sill or, in some cases, after the wave passes through the sill. Details of waveform transformation leading toward the breaking and subsequent energy dissipation are discussed. 

Metamorphic rocks from the Connecticut Valley Trough (CVT), Vermont, and Massachusetts, have been examined using quartz‐in‐garnet (QuiG) and conventional thermobarometry, thermodynamic reaction modelling, diffusion modelling, and⁴⁰Ar/³⁹Ar thermochronology to constrain theirP–T–tpaths during Acadian metamorphism and subsequent exhumation. Numerous samples, collected in the vicinity of the Acadian domes, contain garnet porphyroblasts that display cloudy zones characterized by numerous fluid inclusions and modified garnet compositions associated with the replacement of the original garnet by biotite±muscovite±plagioclase±quartz±lowX_grs/enrichedX_sps. QuiG and conventional thermobarometry constrain both the conditions of garnet nucleation and peakP–Tconditions to have occurred at ~0.85–1.05 GPa, ~550–600°C. Most notably, QuiG barometry was performed on inclusions adjacent to these reaction zones in conjunction with Gibbs method reaction modelling to reveal that these dissolution–reprecipitation reactions occurred during nearly isothermal decompression from the peakP–Tconditions to around ~0.3 GPa, 550°C. Diffusion modelling reveals that the Mn zoning profiles created during garnet resorption that accompanied decompression formed in less thanc. 3 Ma, which constrains the tectonic exhumation to have occurred at 8–10 mm/year. Subsequent cooling to 500°C occurred rapidly at a rate of 100°C/Ma, followed by slower cooling reaching 1.7°C /Ma by the mid Carboniferous. This is the first reported example of QuiG barometry revealing a multi‐stage metamorphic history and highlights the utility of this method for unravelling complex metamorphic terranes.

 

The time of origin of the geodynamo has important implications for the thermal evolution of the planetary interior and the habitability of early Earth. It has been proposed that detrital zircon grains from Jack Hills, Western Australia, provide evidence for an active geodynamo as early as 4.2 billion years (Ga) ago. However, our combined paleomagnetic, geochemical, and mineralogical studies on Jack Hills zircons indicate that most have poor magnetic recording properties and secondary magnetization carriers that postdate the formation of the zircons. Therefore, the existence of the geodynamo before 3.5 Ga ago remains unknown.