Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geoscience”

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.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development of Cr, Se, U, Sb, and Te isotopes as indicators of redox reactions, contaminant fate, and contaminant transport in aqueous systems: A review

Cr, Se, U, Sb, and Te are toxic, redox-active elements that are more mobile and environmentally problematic in their oxidized forms, and less mobile and bioavailable in their reduced forms. This chapter reviews the development of Cr, Se, U, Sb, and Te isotope ratio measurements as new indicators of redox reactions and contaminant migration. Reliable analytical methods exist, but are still evolving. Understanding of isotopic fractionation induced by various (bio)geochemical processes has been explored in dozens of publications, yet is far from complete: Reduction reactions, the major driver of isotopic variation, have been relatively well studied. However, the magnitude of fractionation is variable and the systematics of that variation are still being explored. Isotopic fractionation induced by oxidation reactions is not well understood. Non-redox reactions, which involve smaller changes in bonding of these elements, tend to induce less isotopic fractionation, but can nonetheless cause significant isotopic shifts. Field applications of Cr, Se, U isotope ratios have demonstrated that they are useful as indicators of reduction in natural systems. A few studies suggest they are also useful as indicators of oxidation and contaminant sources. The physical and chemical complexity of groundwater systems hinders accurate quantitative interpretation of Cr, Se, U isotope data using simple models. Numerical models have been developed that capture the behavior of complex, coupled systems and enable the most effective extraction of information from field data sets.

58 GEOSCIENCES↗
DOE OSTI2021

Appropriate “marginal” farmlands for second‐generation biofuel crops in North Carolina

Abstract Current research on bioenergy crops shows that perennial grasses can yield substantial amounts of dry biomass with relatively low inputs of water and fertilizer. In order to minimize competition with food production, it has been suggested that bioenergy crops could be directed to land areas less suitable for commodity crops, commonly referred to as “marginal” lands. These are land units with inherent limitations to vegetative growth and production, which may be due to several factors (soil physical and chemical properties, climatic conditions, etc.). However the term “marginal” is an adjective with imprecise meaning, and objective criteria for determining “marginal” lands for siting bioenergy crops are necessary. Here we propose that such criteria may be based on soil survey classifications and realistic yield estimates, and we show an example of its use to justify site selection for bioenergy crops in different regions of North Carolina.

09 BIOMASS FUELS↗
DOE OSTI2026

Co‐Electrolysis of CO 2 and H 2 O to Syngas on Bimetallic Pd x Cu 1‐ x Catalysts for Tandem Thermochemical Conversion to Carbon Nanofibers

Electrification of chemical production using renewable energy and abundant feedstocks offers a promising pathway for decarbonizing the chemical industry. Current efforts on CO 2 valorization largely focus on making chemicals and fuels. Here, to help achieve net-negative emissions through long-term carbon storage, this study aims to develop efficient electrocatalysts for a tandem electrochemical-thermochemical process to convert CO 2 into carbon nanofibers (CNFs). CO 2 and water are first electrochemically reduced in a membrane electrode assembly (MEA) electrolyzer to produce syngas (CO + H 2 ), which is subsequently fed into a thermochemical packed bed reactor to facilitate CNF growth. This work systematically evaluated Pd x Cu 1-x bimetallic electrocatalysts to assess the effect of Pd–Cu alloying on enhancing syngas production while reducing Pd loading. Transmission electron microscopy and Raman spectroscopy confirmed the formation of high-purity, crystalline CNFs, regardless of the syngas composition from the MEA. In situ X-ray absorption spectroscopy and X-ray diffraction measurements revealed that increasing Cu content in the Pd x Cu 1-x alloy progressively inhibited palladium hydride formation, consistent with DFT calculations on the stability of Pd x Cu 1-x under reducing electrochemical potentials.

58 GEOSCIENCES↗
DOE OSTI2025

Enhancing Acidic Oxygen Evolution Activity by Controlling Oxidation State of Iridium

Iridium oxides with high oxidation states have been reported to be effective in enhancing the acidic oxygen evolution reaction (OER) performance. In this work, we develop ultrasmall IrO x nanoparticles (NPs) over titanium nitride (TiN), which undergoes surface oxidation under oxidative conditions to form oxygen‐modified TiN (oxi‐TiN), enabling the formation of highly oxidized Ir δ+ (δ > 4). This IrO x /oxi‐TiN catalyst delivers higher Ir mass activity than commercial IrO 2 , while comparable stability is maintained. The superior OER activity of IrO x /oxi‐TiN is further demonstrated in a proton exchange membrane water electrolyzer (PEMWE), requiring only 1.88 V to reach 3 A cm −2 , achieving the U.S. Department of Energy 2025 target (1.90 V at 3 A cm −2 ). In situ X‐ray absorption spectroscopy (XAS) confirms that the superior OER activity of IrO x /oxi‐TiN originates from highly oxidized Ir δ+ under OER conditions. Density functional theory (DFT) calculations reveal a general correlation between the oxidation state of Ir and OER overpotential. Specifically, the introduction of interfacial oxygen at the Ir/TiN interface increases the oxidation state of deposited Ir δ+ from δ < 4 to δ > 4, decreasing the OER overpotential. This study highlights the critical role of high oxidation states of Ir δ+ in enhancing OER activity, providing guidance for the development of advanced acidic OER catalysts.

58 GEOSCIENCES↗
DOE OSTI2024

Rejection of low-molecular weight neutral organics is highly sensitive to reverse osmosis system design and operation

A computational model was developed to investigate the significance of system design and operating conditions on the rejection of neutral, low-MW organics by reverse osmosis for potable reuse. Here, the model demonstrated that the decrease in local rejection as net driving pressure decreases is substantially greater for moderately rejected compounds than for highly rejected compounds. At recovery values less than 70%, the local permeate concentration can exceed the pressure vessel feed concentration for moderately rejected compounds. System-level rejection of moderately rejected compounds is likewise substantially more sensitive to operating conditions than highly rejected compounds. The findings highlight a drawback of relying on rejection results from bench-scale testing that operates at low recovery, which invariably has higher rejection than full-scale systems operating at similar pressure. The analysis demonstrates a trade-off in which the low-pressure, high-recovery operation desired for potable reuse systems can be detrimental to the removal of low-MW neutral organics. The removal of low-MW neutral organics can be improved if organics rejection is explicitly evaluated during the design process.

42 ENGINEERING↗
DOE OSTI2024

Comparative clumped isotope temperature relationships in freshwater carbonates

Lacustrine, riverine and spring carbonates represent archives of terrestrial climates and their geochemistry has been used to study palaeoenvironments. Clumped isotope thermometry is an emerging tool that has been applied to freshwater carbonates. Limited work has been done to evaluate comparative relationships between clumped isotopes and temperature in different types of modern freshwater carbonates. This study assembles an extensive calibration data set with 135 samples of modern freshwater carbonates from 96 sites and constrains the relationship between independent observations of water temperature and the clumped isotopic composition of carbonates (denoted by Δ 47 ), including new measurements, and recalculates published data in accordance with current community-defined standard values. For temperature reconstruction, the study reports a composite freshwater calibration and material-specific calibrations for biogenic carbonates (freshwater gastropods and bivalves), fine-grained carbonate (e.g. micrites), biologically mediated carbonates (microbialites and tufas) and travertines. Material-specific calibration trends show a convergence of slopes that are in agreement with recently published syntheses, but statistically significant differences in intercepts occur between some materials (e.g. some biogenics, fine-grained carbonates). These differences may arise due to unresolved seasonal biases, kinetic isotope effects and/or varying degrees of biological influence. The impact of different calibrations is shown through application to new data for glacial and deglacial age travertines from Austria and published data sets. While material-specific calibrations may yield more accurate results for biogenic and fine-grained carbonate samples, the use of material-specific and the composite freshwater calibrations generally produces values within 1.0–1.5°C of each other, and typically fall within calibration uncertainty given limitations of precision.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Plant responses to elevated CO 2 under competing hypotheses of nitrogen and phosphorus limitations

Abstract The future ecosystem carbon cycle has important implications for biosphere‐climate feedback. The magnitude of future plant growth and carbon accumulation depends on plant strategies for nutrient uptake under the stresses of nitrogen (N) versus phosphorus (P) limitations. Two archetypal theories have been widely acknowledged in the literature to represent N and P limitations on ecosystem processes: Liebig's Law of the Minimum (LLM) and the Multiple Element Limitation (MEL) approach. LLM states that the more limiting nutrient controls plant growth, and commonly leads to predictions of dramatically dampened ecosystem carbon accumulation over the 21st century. Conversely, the MEL approach recognizes that plants possess multiple pathways to coordinate N and P availability and invest resources to alleviate N or P limitation. We implemented these two contrasting approaches in the E3SM model, and compiled 98 in situ forest N or P fertilization experiments to evaluate how terrestrial ecosystems will respond to N and P limitations. We find that MEL better captured the observed plant responses to nutrient perturbations globally, compared with LLM. Furthermore, LLM and MEL diverged dramatically in responses to elevated CO 2 concentrations, leading to a two‐fold difference in CO 2 fertilization effects on Net Primary Productivity by the end of the 21st century. The larger CO 2 fertilization effects indicated by MEL mainly resulted from plant mediation on N and P resource supplies through N 2 fixation and phosphatase activities. This analysis provides quantitative evidence of how different N and P limitation strategies can diversely affect future carbon and nutrient dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

California Trees Seasonally Use Augmented Water Sources: Water Isotope Tracking in a Groundwater‐Dependent Ecosystem

Sustainable groundwater management must account for the needs of groundwater dependent ecosystems. To understand the relationship of ecosystems and seasonal water use, we studied the stable isotope composition (δ 18 O and δ 2 H) of water in streamside trees in a semi-arid streamside environment (Livermore, California, USA). We sampled seven trees at two sites every other month from April 2024 through April 2025 for tree xylem stable water isotope signatures. These data were compared to potential source waters: precipitation, imported surface water, soil water and regional groundwaters. Large daily precipitation events were found to be isotopically similar to regional groundwater and were thus treated as one water source. A Bayesian mixing model using stable water isotopes was used to determine the ratios of these three potential source waters (small daily precipitation events, groundwater/large daily precipitation events and imported water) present in tree xylem water. On average, tree water sources include 32% imported water (SD = 9%), 35% small daily precipitation events (SD = 10%) and 33% groundwater (SD = 3%), with significant seasonal variation (t summer-winter = 30.8, p < 0.01), particularly drawing more imported water (more than 55%) in the summer. While small daily precipitation events contribute only 10% of the total precipitation in our dataset, it represents a third of water used by trees. In addition, while these ecosystems are designated as groundwater dependent ecosystems, the trees use approximately one third imported water and even more during dry summer months. This approach provides water managers with a practical tool for quantifying ecosystem water needs, supporting data-driven decisions and regulatory compliance. While California's Sustainable Groundwater Management Act emphasizes supporting groundwater dependent ecosystems, it allows flexibility in demonstrating benefits. In conclusion, our methodology provides a way to document that management actions (such as managed aquifer recharge with imported water) deliver measurable co-benefits to GDEs.

Geosciences↗
DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗
DOE OSTI2021

Variability of near-fault seismic risk to reinforced concrete buildings based on high-resolution physics-based ground motion simulations

Broadband physics-based simulated earthquake ground motions are utilized to characterize the regional-scale seismic risk to modern reinforced concrete (RC) structures. A highly dense dataset of ground motions covering a 100-km urn:x-wiley:00988847:media:eqe3413:eqe3413-math-0001 40-km domain was generated using kinematic fault rupture models with varying rupture characteristics to represent shallow crustal earthquakes and resolved up to frequencies of 5 Hz. Over 40,000 nonlinear response history simulations of short- and mid-rise RC special moment frame buildings were conducted using simulation models that are capable of representing nonlinear behavior and component deterioration effects. The spatial variability of structural risk within a single earthquake scenario and between different rupture scenarios is examined, and the regions of strongest directivity effects and highest structural demands are identified. The structural demands may vary by factor of up to 8.0 at very short distances from the fault, and the large dispersion in the demands decreases significantly beyond a distance of 15 km. The interstory drift and member rotation demands are substantially impacted by important features of the geological structure and the characteristics of the rupture scenarios, particularly the presence of localized high-slip regions. The frequency characteristics of the structures are found to play an important role in determining the effects of near-fault ground motions on the structural response and expected damage. The results of this study suggest that the simulated ground motions, particularly those generated using the hybrid rupture approach, may offer reasonable structural risk estimates for low-frequency structures and conservative estimates for high-frequency structures.

58 GEOSCIENCES↗
DOE OSTI2021

Engineering evaluation of the EQSIM simulated ground-motion database: The San Francisco Bay Area region

Ground-motion simulations for infrequent earthquake scenarios are gaining increasing interest in the engineering community for their potential to inform performance-based structural design and assessment, particularly in regions where major earthquakes are expected, but recordings from consistent historical events are not available. However, the absence of empirical data makes the evaluation of such scenarios and the interpretation of the simulation results very challenging. In this context, this paper presents the evaluation of the first EQSIM (v.1.1.0) ground-motion simulations generated for the San Francisco Bay Area (SFBA) region. The current database, which is at the first stages of development, includes eight realizations of a Hayward Fault M w 7 event. The evaluation of the simulated motions is first conducted on the average ground-motion amplitudes through comparison against the NGA-W2 GMPEs and a selected population of real records. A methodology for selecting a population of real records consistent with the simulated event is proposed. The objective of the proposed evaluation is twofold: (1) build confidence in the realistic character of the simulated motions for their use in engineering applications and (2) offer a critical physics-based interpretation of the simulation results that can help improve key features of the simulation models. To further investigate the implications of using the simulated ground motions for site-specific structural assessments, the ASCE7-16 approach is employed for the analysis of two hazard-consistent near-field sites in the SFBA. Results are discussed. Furthermore, this study offers a critical review of the EQSIM (v.1.1.0) SFBA ground-motion simulations and suggestions for improving the earthquake rupture models.

58 GEOSCIENCES↗
DOE OSTI2022

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2024

Quantification of modeling uncertainty in the Rayleigh damping model

Understanding and accurately characterizing energy dissipation mechanisms in civil structures during earthquakes is an important element of seismic assessment and design. The most commonly used model is attributed to Rayleigh. This paper proposes a systematic approach to quantify the uncertainty associated with Rayleigh's damping model. Bayesian calibration with embedded model error is employed to treat the coefficients of the Rayleigh model as random variables using modal damping ratios. Through a numerical example, we illustrate how this approach works and how the calibrated model can address modeling uncertainty associated with the Rayleigh damping model.

42 ENGINEERING↗
DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗
DOE OSTI2021

Hurricanes alter 10 Be concentrations in tropical river sediment but do not change regional erosion rate estimates

Tropical islands, including many in island arcs, are subjected to recurring disturbances from extreme storms such as tropical cyclones. To test whether such storms influence cosmogenic nuclide concentrations such that they do not reflect long-term rates of erosion, we measured meteoric and in situ 10Be in river sediment samples from Dominica, an andesitic island in the Caribbean, before and after category five Hurricane Maria (2017). Populations of before- and after-storm concentrations are statistically indistinguishable (n=7 pairs for in situ 10 Be, n=11 pairs for meteoric 10 Be). 10 Be concentrations vary from -138% to +73% within before-after sample pairs relative to the mean of the pair. These new data suggest that the effects of extreme storms on the depth and amount of near-surface erosion on Dominica vary spatially. Our data support calculations of Niemi, et al. (2005) and Yanites, et al. (2009) suggesting that basin-by-basin comparisons of erosion rates based on cosmogenic nuclides should be approached with caution in small (<~100 km 2 ) watersheds affected by mass movements and extreme storms. Erosion rates determined from in situ 10Be in Dominica (geometric mean=0.102 mm y-1, n=12) are low compared to similarly steep and wet areas globally and correlate positively with the spatial density of mass movements.

58 GEOSCIENCES↗