SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This paper will discuss the objectives, approach, and anticipated results of a study of earthquakes in Costa Rica. GPS measurements will be taken and field surveys will be made. Assessments of seismic strain accumulation and post-seismic deformation will be made in an effort to understand the effect these processes have on regional tectonic models.
Most of the Zagros deformation resulting from the convergence of Arabia and Eurasia takes place in the Southeast Zagros. To apply the SAR interferometry geodetic technique, a few ERS 1 & 2 satellite images were used to map this continuing deformation proven by GPS. Interferograms over 7 years show surprisingly high coherence. The unwrapped phases display a high correlation with topography reflecting atmospheric noise in addition to the desired tectonic signal. We estimate two simple linear trends and remove them from interferograms. The preliminary results show local uplift rates with a likely minimum of 1-2 mm/yr. These early crude results will be tested by more data in project No. 3174.
Europa's surface shows signs of extension, which is revealed as lithospheric dilation expressed along ridges, dilational bands and ridged bands. Ridges, the most common tectonic feature on Europa, comprise a central crack flanked by two raised banks a few hundred meters high on each side. Together these three classes may represent a continuum of formation. In Tufts' Dilational Model ridge formation is dominated by daily tidal cycling of a crack, which can be superimposed with regional secular dilation. The two sources of dilation can combine to form the various band morphologies observed. New GPS data along a rift on the Ross Ice Shelf, Antarctica is a suitable Earth analog to test the framework of Tufts' Dilational Model. As predicted by Tufts' Dilational Model, tensile failures in the Ross Ice Shelf exhibit secular dilation, upon which a tidal signal can be seen. From this analog we conclude that Tufts' Dilational Model for Europan ridges and bands may be credible and that the secular dilation is most likely from a regional source and not tidally driven.
Pluto and its largest satellite, Charon, currently orbit in a mutually synchronous state; both bodies continuously show the same face to one another. This orbital configuration is a natural end-state for bodies that have undergone tidal dissipation. In order to achieve this state, both bodies would have experienced tidal heating and stress, with the extent of tidal activity controlled by the orbital evolution of Pluto and Charon and by the interior structure and rheology of each body. As the secondary, Charon would have experienced a larger tidal response than Pluto, which may have manifested as observable tectonism. Unfortunately, there are few constraints on the interiors of Pluto and Charon. In addition, the pathway by which Charon came to occupy its present orbital state is uncertain. If Charon's orbit experienced a high-eccentricity phase, as suggested by some orbital evolution models, tidal effects would have likely been more significant. Therefore, we determine the conditions under which Charon could have experienced tidally-driven geologic activity and the extent to which upcoming New Horizons spacecraft observations could be used to constrain Charon's internal structure and orbital evolution. Using plausible interior structure models that include an ocean layer, we find that tidally-driven tensile fractures would likely have formed on Charon if its eccentricity were on the order of 0.01, especially if Charon were orbiting closer to Pluto than at present. Such fractures could display a variety of azimuths near the equator and near the poles, with the range of azimuths in a given region dependent on longitude; east-west-trending fractures should dominate at mid-latitudes. The fracture patterns we predict indicate that Charon's surface geology could provide constraints on the thickness and viscosity of Charon's ice shell at the time of fracture formation.
Dynamic interactions between atmosphere, hydrosphere, lithosphere, and biosphere, driven by solar energy and modulated by atmospheric greenhouse gasses, volcanism, and plate tectonics produce substantial variations in Earth's climate over geologic time. Major Quaternary drivers of climate change include variations in Earth's orbital cycles, global temperature, natural atmospheric greenhouse gasses, and continental ice mass. Multiple abrupt climate perturbations with onsets of a century to decades have also occurred. Since the 1970s, geologically unprecedented rising levels of anthropogenic greenhouse gasses from fossil fuel combustion, industry, and agriculture have become the dominant contributors to rising global temperature, glacier and ice sheet losses, and sea level rise.
Granitic inselbergs are enduring landforms that punctuate planation surfaces in tectonically stable settings, yet the mechanisms and timescales of their persistence remain debated. Here, in this study, we present a multinuclide cosmogenic dataset ( 10 Be, 26 Al and 21 Ne) from bedrock, fluvial sediments and planation surfaces in the semi-arid Central Ceará Domain, north-eastern Brazil, to quantify denudation rates and reconstruct the long-term dynamics of residual relief formation. This triple-nuclide approach enables independent validation of denudation contrasts and long-term exposure histories, providing quantitative constraints on inselberg evolution that a single-nuclide approach could not resolve. Bedrock denudation rates (1–4 m Myr −1 ) are systematically lower than basin-wide rates (8–14 m Myr −1 ), revealing a regional pattern of vertical differential denudation. When combined with summit elevations and surface exposure ages of up to 5 Myr, these contrasts suggest that cumulative relief growth has been operating over timescales of tens of millions of years (1–57 Myr). Morphological and structural observations further indicate localized mass wasting and heterogeneous erosion styles, modulated by rock fabric and joint density. Cosmogenic nuclide inventories in detrital sediments reveal partial reworking from older sediment stocks, including surface pebbles with exposure ages up to 5 Myr, which record multistage of burial and re-exposure histories and reflect episodic remobilization in low-connectivity, semi-arid catchments. Rather than static relicts of ancient surfaces, inselbergs in this cratonic setting emerge as dynamic landforms actively sustained by slow but persistent differential denudation. These results provide new constraints on long-term cratonic landscape evolution and underscore the interplay between lithological resistance, structural inheritance and sediment transfer processes in sustaining topographic relief. They also call for a critical reassessment of classical inselberg models, suggesting that hybrid mechanisms, rather than single-process paradigms, more accurately capture the complexity of residual landform development in cratonic interiors.
Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).
Serpentine minerals have received a lot of attention because of their unique crystal structures, their wide occurrence in orogenic belts and their potential role in contributing seismic anisotropy in subducting slabs. Several studies have investigated crystal preferred orientation (CPO) in high temperature antigorite serpentinites from Japan, the Alps, Spain, Cuba and Tibet, documenting significant crystal alignment. However, only a limited number of lower grade serpentines have been explored to date. Mainly because of submicroscopic microstructural heterogeneities CPO cannot be measured with conventional methods such as optical microscopy and EBSD. In this study 15 serpentinites from different tectonic settings in California, the Central Alps and Northern Spain have been investigated, mainly with high energy synchrotron X-ray diffraction, to quantify bulk crystal alignment. We find that CPO is strong on sheared surfaces of fractured blocks and secondary veins but the bulk of most serpentinite samples, except high-grade recrystallized antigorite serpentinite, show only weak crystal alignment. Correspondingly calculated seismic anisotropy based on CPO is not very significant. This is supported by very heterogeneous microstructures as documented with SEM and TEM analyses.
Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.
Early isotopic studies of sulfate in carbonate minerals (carbonate associated sulfate; CAS) suggested that carbonates can provide a reliable, well-dated archive of the marine sulfur cycle through time. However, subsequent research has shown that diagenetic alteration can impose highly heterogeneous CAS sulfur isotopic compositions (δ 34 S CAS ) among different carbonate phases within sediments. Such alteration necessitates targeted sampling of well-preserved, primary carbonate phases. Here, we present a new record of Carboniferous and Early Permian brachiopod δ 34 S CAS generated from over 130 measurements of microsampled brachiopod shells. Our record refines existing brachiopod δ 34 S CAS records and confirms a large, ~6.5‰ δ 34 S CAS decrease in the Early Carboniferous. Importantly, the record also features a novel 3–5‰ increase in δ 34 S CAS near the Serpukhovian-Bashkirian boundary (323.4 Ma) that coincides with carbonate δ 13 C and δ 18 O increases. Variability in δ 34 S CAS is minor both within (≤0.3‰) and among (≤2‰) individual co-depositional brachiopod specimens. A taxon-specific δ 34 S CAS offset is present one species (Composita subtilita) that also exhibits a δ 13 C offset, supporting the existence of biological “vital effects” on δ 34 S CAS . Geologic evidence and mathematical modeling of the Permo-Carboniferous carbon and sulfur cycles suggest that changes in the burial ratio of organic carbon to pyrite sulfur (R C:S ) are insufficient to explain the observed mid-Carboniferous δ 34 S CAS record. We find that changes in the 34 S depletion of pyrite relative to seawater sulfate ( 34 ε) or in the δ 34 S of the input to the ocean (δ 34 S in ) are also needed. Large additions of O 2 from organic carbon burial during the Permo-Carboniferous cannot be entirely compensated for with sulfur cycle changes; lower than modern late Visean pO 2 and/or additional O 2 sinks are needed to keep pO 2 at plausible levels. Based on the geologic context surrounding our record's mid-Carboniferous δ 34 S CAS increase, we advocate for simultaneous changes in pyrite burial, 34 ε, and δ 34 S in , driven by sea level or tectonically induced changes in environments of sulfur burial, as a viable mechanism to produce rapid seawater δ 34 S changes.
The Pennsylvania Anthracite Fields are in a complex tectonic and metamorphic terrain, historically hosting one of the largest concentrations of coal mining in the USA. Anthracite mining now largely consists of the surface mining of the pillars remaining from the prime years of underground mining. The geochemical study of the sized coal products and the refuse (largely rock) from three preparation plants (breakers) demonstrates that Principal components analysis (PCA) of select major oxide, minor element, and rare earth elements illustrates some differentiation among the products from the individual plants. The rock samples, with abundant quartz and metamorphic Al-Si minerals and with a lower ash-basis REE concentration than the coals, were distinctly separated from the coal samples on the PCA plots. Plots of Gd N /Gd N * vs. Eu N /Eu N * and Ce N /Ce N * vs. Eu N / Eu N * showed differentiation between the individual suites of coals showed that the refuse samples had distinct REE distributions compared to the associated coals. Further, several minor and trace elements show enrichments in the coal samples. Lithium, with concentrations of up to 314 ppm on an ash basis, is among the most promising of the critical elements, exceeding the enrichment of the REY and Sc.
The ability to precisely constrain the physical–chemical conditions of serpentinization, such as temperature and fluid source, is limited by the accuracy of calibrations for oxygen isotope fractionation between serpentine and water – i.e., 1000 lnα( serpentine-water ) – which disagree by up to 20‰ when extrapolated to T<200°C. In this study, we present a new empirical calibration of 1000 lnα( serpentine-water ) aiming to improve applications of oxygen isotope geochemistry to very low-T serpentinization (T<100°C). We used the high-spatial resolution capabilities of Secondary Ion Mass Spectrometry (SIMS) to analyze oxygen isotope ratios in mineral pairs of calcite + serpentine, quartz + serpentine and talc + serpentine co-crystallized at scales≤50μm in six serpentinite samples from the Samail ophiolite (Oman). SIMS analysis shows that the mineral pairs are relatively homogenous in oxygen isotope ratios with variability in δ 18 O values ≤2‰ (2s). Clumped isotope thermometry and petrological constraints indicate crystallization temperatures from 25 to 100°C for the investigated samples. These independent constraints on temperature allowed us to derive 1000 lnα( serpentine-water ) by combining Δ 18 O( quartz-serpentine ), Δ 18 O( talc-serpentine ) and Δ 18 O( calcite-serpentine ) measured by SIMS in the investigated samples with Δ 18 O( quartz-water ), Δ 18 O( talc-water ) and Δ 18 O( calcite-water ) calibrations available in the literature. Our empirical calibration of 1000 lnα( serpentine-water ) = 1.04±0.20 (2SE) × 10 6 /T 2 (T in K), from T = 25 to 100°C, is within uncertainty of former high-temperature empirical calibrations extrapolated to T<100°C and matches experimental calibrations when extrapolated to T=250°C. In conclusion, the new serpentine-water calibration enables more accurate reconstructions of fluid-rock interactions occurring during low-temperature serpentinization processes in various tectonic settings.
The islands of the Galápagos archipelago are of different ages due to the combination of a geostationary volcanic hotspot and the eastward movement of the Nazca tectonic plate. Taking advantage of this, a soil chronosequence (1.5–1070 ka) has recently been established covering four of the islands. This was used to investigate the temporal development of soil microstructure and micromineralogy using micromorphological methods in combination with synchrotron-based mineralogical and geochemical microanalyses. Along the chronosequence, the microstructure changed from intergrain microaggregates in the young soils to angular blocky microaggregates in the older soils. In the young soils, a high amount of vitric scoria and weatherable minerals occur, but they quickly disappear with age, and micromass increased fast. The micromorphological studies as well as micro-XRD and micro-XRF revealed increasing heterogeneity of soil microstructure and micromineralogy due to weathering and pedogenic processes until intermediate soil age. In the oldest soils, the microstructure became more uniform again, de-mixing pedofeatures like clay coatings have completely disappeared, and kaolinite, hematite and gibbsite dominate the mineralogical composition. Our study shows that the development of soil microstructure mirrors macroscopic soil profile development showing increasing differentiation until intermediate developmental stages and receding heterogeneity in highly weathered soils.
Understanding the cause of the November 2017 Pohang main shock is of considerable importance to the future of the geothermal industry because of its large magnitude compared to prior expectations based on case histories of other projects involving underground fluid injection. Of the three possibilities – induced, triggered or natural, “induced” can be ruled out based on the disproportionately large seismic moment of the main shock. Furthermore, whether natural or triggered, the source of seismic hazard at Pohang was tectonic strain accumulation, not fluid injection. Arguably, the most timely indicator of seismic hazard and risk in the environs of Pohang was the September 2016 Mw 5.4 Gyeongju earthquake, which was natural and located about 40 km south of Pohang along the same active fault system.
Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.
There is a growing sociopolitical desire to develop cleaner energy sources in the United States and maintain energy security. Regardless of politics, many coal-fired electric plants have already been shut down and many utilities are vowing to retire their current coal-fired assets within the next two decades. Replacement power assets require consideration of appropriate siting. A geographic information system (GIS)-based multicriteria decision analysis approach is useful to assist utility and energy companies, as well as policymakers, to evaluate potential areas for siting new plants in the contiguous United States. A GIS-based framework is simply a database of location information that allows for mapping, querying, modeling, and analyzing data based on location. The spatial output can be structured to be visual, allowing for easier analysis of location data. The need to site additional power assets, including renewable resources and clean power sources, such as nuclear, led to the development of the Oak Ridge Siting Analysis for power Generation Expansion (OR-SAGE) tool discussed in this paper. The tool takes inputs such as population growth, water availability, environmental indicators, and tectonic and geological hazards to provide an in-depth visual analysis for siting options. Energy companies and other stakeholders can use OR-SAGE to procure feedback quickly and effectively on land suitability based on technology specific inputs. Policymakers can use OR-SAGE to analyze the impacts of future energy technology decisions, while balancing competing resource use. Overall, this paper discusses the recent use of OR-SAGE for these purposes and plans for future development.
The Gulf of Elat-Aqaba (GEA), located in the southern part of the Dead Sea Transform (DST), is one of the most active segments of the Dead Sea fault system. Yet, some fundamental details concerning its evolution in space and time are still not fully understood. To better constrain the tectonic history of this region, we study a succession of alluvial, lacustrine, and fluvial deposits in the north-western edge of the GEA, which belong to the Quaternary Eilot and Garof Fms. The Eilot Fm. represents low-energy shallow freshwater bodies and sabkhas, which predate the development of the deep part of the basin. Later, in response to the subsidence of the basin, the Garof Fm. was deposited in alluvial fans that built up following the increase in the depositional energy. Cosmogenic burial ages combined with magnetostratigraphy show that these Fms. were deposited as early as 3.03 Ma and as late as 1.65 Ma, depending on the applied age model. Considering the present depth of the northern GEA, this age indicates a maximum subsidence rate of the head of the GEA of ~1.98mm/yr. We suggest that faults along the onshore margin of the deep basin were active during the deposition of the Garof Fm. and subsequently became inactive, while the faulting activity migrated toward the center of the basin. In conclusion, closed Holocene basins north of the present seashore, such as the Elat and Evrona playas, which may represent a modern analogue to the Plio-Pleistocene morphology, raise the possibility that the GEA head is still propagating northward.
This paper presents modeling of mechanical anisotropy in argillaceous rocks using an irregular lattice modeling approach, namely the rigid-body-spring network. To represent the mechanical anisotropy, new schemes are implemented in the modeling framework. The directionality of elastic deformation is resolved by modifying the element formulation with anisotropic elastic properties. The anisotropy of strength and failure characteristics is facilitated by adopting orientation-dependent failure criteria into the failure model. The verification of the improved modeling procedures is performed against theoretical model predictions for unconfined compression tests with various bedding orientations. Furthermore, excavation damage and fracturing processes in rock formations are simulated for different geomechanical configurations, such as rock anisotropy and tectonic heterogeneity. The simulated excavation damage characteristics are realistic and comparable with the actual field observation at a tunnel located in an argillaceous clay formation. The simulation results provide insights into the excavation damage zone phenomena with an explicit representation of fracturing processes.