Search for: All records

Ganderia and the Southeastern New England Avalon terrane are both terranes that rifted from Gondwana and accreted to North America in the early to mid-Paleozoic. Accretion of the Avalon terrane was accompanied by plutonism, deformation, and metamorphism including partial melting within the Nashoba terrane, the trailing edge of Ganderia, and may be interpreted as indicators for mid- to lower-crustal channel flow. Channel flow describes the flow of weak, partially molten material between more competent crust as a result of pressure gradients in the mid- to lower crustal levels. Such flow should typically result in seismic anisotropy due to the crystallographic preferred orientations of minerals and shape preferred orientations at various scales. Here, we present first results for the crustal anisotropic structure beneath the Nashoba terrane that were produced with a newly developed approach from currently collected data in the region. To investigate the hypothesis of crustal flow during the orogenic history of Southeastern New England, we deployed a dense profile of 6 broadband seismic stations crossing the Nashoba terrane. We analyze the harmonic variation of amplitudes in teleseismic P-Receiver Functions (RFs) to identify interfaces of isotropic and anisotropic contrasts within the crust. In the case of particularly prominent anisotropic features that have significantly larger amplitudes than other signals, it is feasible to derive quantitative constraints on the strength and orientation of the anisotropy. However, with growing complexity, a classical forward modelling or grid search approach becomes unfeasible. These difficulties can be mitigated by applying Bayesian inversion, which infers values of model parameters from a probabilistic perspective. Applying a Bayesian inversion to the harmonically decomposed RFs has the potential to infer complex anisotropic seismic structures. We find evidence for two crustal anisotropic layers with confined properties to the geologic units of the Ganderia, Nashoba and Avalon terranes that might be related to episodes of lower crustal flow. In addition, we identify anisotropy at shallow mantle levels beneath the Nashoba and Avalon Terranes possibly indicating the upper interface of a shallow asthenosphere in the region. 

Seismic waves with different propagation and oscillation directions can exhibit different velocities when going through a medium with some directional properties; this phenomenon is called seismic anisotropy. Seismic anisotropy observed beneath eastern North America is often attributed to present-day flow in the upper mantle. The mantle flow causes shear waves oscillating in the direction of flow (e.g., in the direction of North America plate motion) to travel faster than those that travel in other directions. However, this pattern does not hold true for some regions along the Appalachian orogen, suggesting that past tectonic events can result in long-lived, ‘frozen-in’ anisotropy in the lithosphere, which modifies the predicted anisotropic behavior beneath these regions. In this study, we investigate sources of seismic anisotropy beneath southern New England using a method based on directionally dependent variations of P-wave to S-wave conversions at interfaces with contrasts in anisotropy. This method can separate signals caused by different anisotropic features and constrain the depth distribution of anisotropy. Within the crust there are multiple features that may be correlated with stratification in the Hartford Basin, faults in the Taconic thrust belt, shear zones formed during Salinic/Acadian terrane accretion events, and orogen-parallel crustal flow in the Acadian orogenic plateau. We apply a Bayesian inversion method to obtain quantitative constraints on the direction and strength of intra-crustal anisotropy beneath the Hartford Basin. In the upper mantle, we identify a fossil shear zone possibly formed during oblique subduction of Rheic Ocean lithosphere. We also find evidence for a plate motion-parallel flow zone in the asthenosphere that is likely disturbed by mantle upwelling near the southern margin of the Northern Appalachian Anomaly in the eastern part of the study area. 

Ganderia and the Southeastern New England Avalon terrane, are both terranes that rifted from Gondwana and accreted to North America in the early to mid-Paleozoic. Accretion of the Avalon terrane was accompanied by plutonism, deformation, and metamorphism including partial melting within the Nashoba terrane, the trailing edge of Ganderia and may be interpreted as indicators for mid- to lower-crustal channel flow. Channel flow describes the flow of weak, partially molten material between more competent crust as a result of pressure gradients in the mid- to lower crustal levels. Such flow should typically result in seismic anisotropy due to the crystallographic preferred orientations of minerals and shape preferred orientations at various scales. Here, we present a new method designed to analyze the crustal anisotropic structure beneath the Nashoba terrane and provide insight into its capabilities in a first application to permanent stations in the area and currently collected data. To investigate the hypothesis of crustal flow during the orogenic history of Southeastern New England, we deployed a dense profile of 6 broadband seismic stations crossing the Nashoba terrane. We analyze the harmonic variation of amplitudes in teleseismic P-Receiver Functions (RFs) to identify interfaces of isotropic and anisotropic contrasts within the crust. In the case of particularly prominent anisotropic features that have significantly larger amplitudes than other signals, it is feasible to derive quantitative constraints on the strength and orientation of the anisotropy. However, with growing complexity, a classical forward modelling or grid search approach becomes unfeasible. These difficulties can be mitigated by applying Bayesian inversion, which infers values of model parameters from a probabilistic perspective. Here we use a Bayesian framework to invert for the anisotropic model that best fits the observed constant and harmonic terms. Applying a Bayesian inversion to the harmonically decomposed RFs instead of full RF waveforms has the potential to infer the anisotropic seismic model more faithfully, without attempting to fit unrelated signals and artifacts. 

The New England Appalachians provide a fascinating window into a host of fundamental geological problems. These include the modification of crustal and mantle lithospheric structure via orogenesis, terrane accretion, and continental rifting, the evolution of individual terranes through processes such as channel flow and ductile extrusion, and the causes and consequences of the Northern Appalachian Anomaly (NAA), a prominent geophysical anomaly in the upper mantle. Recent and ongoing deployments of dense seismic arrays in New England are providing images of the crust and upper mantle in unprecedented detail, allowing us to address both new and longstanding science questions. These deployments include the Seismic Experiment for Imaging Structure beneath Connecticut (SEISConn, 2015-2019), the New England Seismic Transects (NEST, 2018-present), and the GEology of New England via Seismic Imaging Studies (GENESIS, 2022-present) arrays. Here we present results from these experiments that are shedding new light on the tectonic evolution of New England and the ways in which structures and processes in the upper mantle can affect the structure of the overlying lithosphere. These include detailed new images of crustal architecture beneath central and southern New England, including a sharp transition from thick (~48 km) crust Laurentia terranes to thin (~32 km) crust beneath Appalachian terranes. The character of this offset beneath the SEISConn and NEST arrays suggests an overlap of two Moho boundaries, forming an overthrust-type structure that may have resulted from reactivation of faults during the compression and shortening associated with the formation of the hypothesized Acadian Altiplano. Beneath SEISConn, there is evidence for multiple relict structures preserved in the lithosphere from past episodes of terrane accretion and suturing, as well as anisotropic layering that constrains the kinematics of past lithospheric deformation events. Beneath the NEST line in central New England, we infer a relatively shallow (~80 km) lithosphere-asthenosphere boundary above the NAA upper mantle geophysical anomaly, providing evidence for lithospheric thinning above a presumed asthenospheric upwelling. Finally, preliminary results suggest layered crustal anisotropy beneath the GENESIS array, perhaps corresponding to a past episode of channel flow in the mid-crust. 

Neural-network quantum molecular dynamics (NNQMD) simulations based on machine learning are revolutionizing atomistic simulations of materials by providing quantum-mechanical accuracy but orders-of-magnitude faster, illustrated by ACM Gordon Bell prize (2020) and finalist (2021). State-of-the-art (SOTA) NNQMD model founded on group theory featuring rotational equivari- ance and local descriptors has provided much higher accuracy and speed than those models, thus named Allegro (meaning fast). On massively parallel super- computers, however, it suffers a fidelity-scaling problem, where growing number of unphysical predictions of interatomic forces prohibits simulations involving larger numbers of atoms for longer times. Here, we solve this problem by com- bining the Allegro model with sharpness aware minimization (SAM) for enhanc- ing the robustness of model through improved smoothness of the loss landscape. The resulting Allegro-Legato (meaning fast and “smooth”) model was shown to elongate the time-to-failure tfailure, without sacrificing computational speed or accuracy. Specifically, Allegro-Legato exhibits much weaker dependence of time- to-failure on the problem size, t_failure = N^−0.14 (N is the number of atoms) compared to the SOTA Allegro model (t_failure ∝ N^−0.29), i.e., systematically delayed time-to-failure, thus allowing much larger and longer NNQMD simulations without failure. The model also exhibits excellent computational scalabil- ity and GPU acceleration on the Polaris supercomputer at Argonne Leadership Computing Facility. Such scalable, accurate, fast and robust NNQMD models will likely find broad applications in NNQMD simulations on emerging exaflop/s computers, with a specific example of accounting for nuclear quantum effects in the dynamics of ammonia to lay a foundation of the green ammonia technology for sustainability.