Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load factor”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Urban wash-off of tire wear particles

Tire wear particles (TWPs) are an important class of microplastics due to their toxicity and abundance. Because most TWPs are generated on impervious road surfaces, urban wash-off is the critical first phase of waterborne transport from their zone of production to stormwater drainage. However, little is known about the driving factors behind their mobilization. In this study, we use a rainfall simulator to investigate how surface roughness, rainfall intensity, and surface slope affect wash-off behaviors of TWPs. We also analyze how the size and shape of mobilized TWPs change over the course of simulated storm events. We found that low surface roughness, high rainfall intensity (most significant factor), and low slope result in the most rapid conveyance of TWP load. On average, large particles (>1000 µm) travelled faster than small particles (<125 µm). Particle shape explained a very small amount of variance in TWP wash-off velocity but was found to be more important under higher surface roughness conditions. In addition to wash-off velocity, we found similar conditions controlled the percent mobilization of TWPs. Low surface roughness and high rainfall intensity resulting in higher TWP wash-off rates is consistent with mineral sediment wash-off behavior. Conversely, low surface slope and large particle size leading to faster conveyance is directly opposed to mineral sediment wash-off. Our findings suggest drag-dominated flow and that sufficient runoff depth is the most important parameter governing TWP wash-off. These findings are important first steps to understanding wash-off behaviors of TWPs and informing future modeling efforts and mitigation strategies.

13 HYDRO ENERGY↗
DOE OSTI2021

A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.

divertor physics↗
DOE OSTI2024

Understanding early HIV-1 rebound dynamics following antiretroviral therapy interruption: The importance of effector cell expansion

Most people living with HIV-1 experience rapid viral rebound once antiretroviral therapy is interrupted; however, a small fraction remain in viral remission for an extended duration. Understanding the factors that determine whether viral rebound is likely after treatment interruption can enable the development of optimal treatment regimens and therapeutic interventions to potentially achieve a functional cure for HIV-1. We built upon the theoretical framework proposed by Conway and Perelson to construct dynamic models of virus-immune interactions to study factors that influence viral rebound dynamics. We evaluated these models using viral load data from 24 individuals following antiretroviral therapy interruption. The best-performing model accurately captures the heterogeneity of viral dynamics and highlights the importance of the effector cell expansion rate. Our results show that post-treatment controllers and non-controllers can be distinguished based on the effector cell expansion rate in our models. Furthermore, these results demonstrate the potential of using dynamic models incorporating an effector cell response to understand early viral rebound dynamics post-antiretroviral therapy interruption.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Studies of optical and chemical properties of aged and fresh biomass burning absorbing aerosols for climate models

Particulate emissions from burning fuels from sub-Saharan Africa is currently poorly characterized. With increasing population and industrialization, emissions are expected to grow. Optical methods for measuring light absorption in the ultraviolet (UV) and visible (vis) regions of the spectrum are crucial for understanding the optical properties of BB aerosols (Bond and Bergstrom, 2006; Lack and Cappa, 2010; Moschos et al., 2021). One such instrument is the dual-spot AE33 model aethalometer is a common filter-based instrument that can be used to measure aerosol absorption properties (Hansen et al., 1984). This instrument has many figures of merit; however, the method's limitations include potential biases from scattering within the filter material, multiple scattering effects of the particles, filter loading (i.e. particle shadowing), and the need for wavelength-specific correction factors (Cλ), which depend on accurate intercomparison with other in situ techniques (Bond et al., 1999; Collaud Coen et al., 2010). This work aims to: 1. Optically characterize BB aerosol generated from African fuel 2. Correct the aethalometer at multiple wavelengths 3. Parametrize these correction factors. 4. Combusting the same fuels with two difference tube furnace systems assess the general applicability of this flexible means of producing BB aerosol. This was accomplished by conducting laboratory studies at Los Alamos National Laboratory (LANL) and North Carolina A&T State University (NCA&T), where aethalometer measurements were compared to in situ optical measurements. Mass cross sections of BB aerosol from these under-studied fuels will be determined along with their wavelength dependence.

09 BIOMASS FUELS↗
DOE OSTI2026

Ultra-High Vacuum Outgassing Characterization of Thermally Processed Low-Carbon Steel for Advanced Particle Accelerator and Gravitational Wave Detector Applications

This dissertation investigated AISI 1020 low-carbon steel as an alternative vacuum chamber material to conventional stainless steel for ultra-high vacuum (UHV) and extreme-high vacuum (XHV) applications. After a 400 °C/48 h bake, AISI 1020 tube chambers achieved a hydrogen outgassing rate of 2.4 × 10¿¹6 Torr·L·s¿¹·cm¿², approximately 2,300 times lower than the prebaked 316L stainless-steel comparator, among the lowest hydrogen outgassing rates ever reported for an uncoated metallic vacuum chamber. Bare and magnetite-coated AISI 1020 chambers were then compared using throughput and rate-of-rise methods. The magnetite coating yielded 5× lower water outgassing at room temperature, but this advantage disappeared after 80 °C baking. After full thermal conditioning (400 °C/48 h prebake followed by 150 °C/96 h and 200 °C/110 h), bare steel achieved 25× lower hydrogen outgassing than the magnetite-coated chamber (9.6 × 10¿¹6 Torr·L·s¿¹·cm¿²) and >99% H2 purity with carbon species below RGA detection. Monte Carlo molecular flow simulations of a CEBAF photogun beamline (96 scenarios) showed that replacing 304L stainless steel with AISI 1020 reduces equilibrium H2 pressure by a factor of 833; a single 304L electrode contributes 98.8% of the gas load despite occupying only 9.1% of the internal surface area. A 500-m Einstein Telescope beampipe screening showed that corrugated bellows contribute 18% of the gas load from only 0.7% of the surface area. A five-model adsorption isotherm framework applied to 22 pumpdown datasets (164 fits with AR(1)-GLS correction) established that the experimental protocol, not the material, controls isotherm identifiability: Dubinin–Radushkevich wins isothermal pipe pumpdowns; Langmuir wins thermally dominated chamber bakes. Cross-dataset joint fitting of the AISI 1020 pipe pumpdowns yielded an H2 diffusion activation energy Ed = 7.24 ± 1.28 kcal·mol¿¹, consistent with trap dominated diffusion in commercial low-carbon steels. Two companion innovations were developed: a Variable Conductance Device (VCD, patent pending IDF-00723) for XHV outgassing measurement, and VacuumDesignerPro (VDP), a MATLAB-based design tool validated against LIGO benchmarks.

Al-Allaq, Aiman H [Old Dominion University]↗
DOE OSTI2026

Assessing Geospatial and Seasonal Influences on Energy and Cost-Efficiency of Drayage Trucks

The electrification of heavy-duty vehicles is a critical pathway toward improved energy efficiency in the freight sector. The current battery electric truck technology poses several challenges to commercial vehicle operations, such as limited driving range, sensitivity to climate conditions, and long recharging times. Estimating the energy consumption of heavy-duty electric trucks is crucial to assessing the feasibility of fleet electrification and its impact on the electric grid. This article focuses on developing a model-based simulation approach to predict and analyze the energy consumption of electric trucks by considering the impact of weather and geographical conditions on vehicle road load and auxiliary components power consumption, as well as the impact these factors have on driving range. Specifically, drayage trucks employed in logistics around maritime ports are used as a case study, with consideration of seasonal climate variations and geographical characteristics at different locations. The article includes results for three major container ports within the United States, providing region-specific insights into the energy requirements and driving range of the electric drayage trucks in these regions, which will inform decision-makers in integrating electric trucks into the existing drayage operations and plan investments for electric grid development.

Sujan, Vivek [ORNL] (ORCID:0000000269882342)↗
DOE OSTI2021

The Building Science Advisor: A Web-Based Tool to Assess the Durability of Building Envelope Components

The Department of Energy (DOE) and Oak Ridge National Laboratory (ORNL) have developed an innovative new design tool that will put building science expert advice in the hands of every building designer who cares to manage moisture risk in high-R envelope components. As modern buildings become increasingly more airtight, constructed with modern labor-saving materials, and equipped with air-conditioning, they have become less forgiving to moisture intrusion. The Building Science Advisor (BSA) provides building science knowledge and advice based on expert experience, field measurements, laboratory tests, and computer simulations.BSA users are prompted to enter relevant information about the building location, design, and material selection options like cladding, structural system, and insulation. Influential factors such as climate, building air tightness, material properties, and internal moisture loads are also considered to estimate and compare the moisture durability performance of several design options. If the BSA deems the assembly’s performance unsatisfactory, it will provide the reasons why, and suggest necessary changes to wall design to ensure more robust performance.Using the BSA tool enables building designers to confidently select assembly design characteristics that achieve their design goals with the least moisture durability risk. Links to design-specific guidance are also be provided to help users manage any remaining risk. This tool will further enable DOE’s Building Technologies Office (BTO) to meet its long-term energy goal of a 50% reduction in building energy consumption by reducing builder concern about using highly energy efficient wall systems.

Desjarlais, Andre Omer↗
DOE OSTI2023

Experimental and computational studies of load ratio effect on fatigue crack growth rates in Alloy 709 at elevated temperatures

Abstract Experimental testing and finite element simulations of fatigue crack growth were performed in austenitic 20Cr‐25Ni (Alloy 709) steel at different load ratios and elevated temperatures. The experimental tests were performed using compact tension specimens, and crack growth rates were measured at stress intensity factors ranging between 5 and 35 MPa√m. Fractographic analysis using scanning electron microscopy indicated crack surface roughness and secondary cracking depending on testing temperature and load ratio. Finite element simulations of fatigue crack growth were performed to compute plasticity‐induced crack opening loads and predict crack growth rates. Predictions of fatigue crack growth rates using finite element simulations were performed using the computed crack‐tip opening loads, and they are shown to match well the experimental measurements.

Engineering↗
DOE OSTI2021

Steady-State Thermal-Hydraulic Analysis of the LEU-Fueled Dalat Nuclear Research Reactor

This paper presents results of steady-state thermal-hydraulic analysis for the designed working core of the Dalat Nuclear Research Reactor (DNRR) using the PLTEMP/ANL code. The core was designed to be loaded with 92 low-enriched uranium (LEU) VVR-M2 fuel bundles (FBs) and 12 beryllium rods surrounding a neutron trap at the core center, for replacement of the previous core with 104 high-enriched uranium (HEU) VVR-M2 FBs. Before using this code for thermohydraulic analysis of the designed LEU working core, it was validated by comparing calculation results with experimental data collected from the HEU working core of the DNRR. The discrepancy between calculated results and measured data was at the maximum about 0.8°C and 1.5°C of fuel cladding and outlet coolant temperatures, respectively. In the design calculation, thermohydraulic safety was confirmed through evaluation of the fuel cladding and coolant temperatures, as well as of other safety parameters such as Departure from Nucleate Boiling Ratio (DNBR) and Onset of Nucleate Boiling Ratio (ONBR). The calculation results showed that, in normal operation conditions at full nominal thermal power of 500 kW without uncertainty parameters, the maximum fuel cladding temperature of the hottest FB was about 90.4°C, which is lower than its limit value of 103°C, the minimum DNBR was 32.0, which is much higher than the recommended value of 1.5, and the minimum ONBR was 1.43, which is higher than the recommended value of 1.4 for VVR-M2 LEU fuel type. When the global and local hot channel factors were taken into account, the maximum temperature of fuel cladding at the hottest FB was about 98.4 °C, for global only, and 114.3°C, for global together with local hot channel factors. The calculation results confirm the safety operation of the designed LEU core loaded with 92 fresh VVR-M2 FBs.

Nguyen, Kien-Cuong↗
DOE OSTI2020

Sensitivity Effects of High Temperature Overhead Conductors to Line Rating Variables

For traditional overhead transmission lines, the maximum allowable conductor operating temperature on ASCR, AAC, or AAAC lines for determining the static line ratings can be quite low, typically well under 100 C. At these low temperatures, the primary driver for determining static line ratings are assumptions made about the wind speed and direction and the ambient temperature. Assumptions for solar loading and emitted thermal radiation, driven by line emissivity and absorptivity, are secondary factors, and only provide small changes to the static ratings. However, when considering newer high temperature conductors such as ACSS, ACCR, ZTACIR, and ACCC® the maximum allowable conductor temperatures can be significantly higher, up to 250° C. This higher temperature shifts the importance of the emissivity assumptions of the overhead transmission line, due to the nature of the T4 dependence of the radiative heat loss compared to the T1 dependence of the convective heat loss, as well as shifting the dependency on local wind conditions. Typical assumptions in the United States for the emissivity/absorptivity of overhead transmission lines in determining the static ratings are to use a value of 0.5 for both parameters or set a value in the range of 0.7 to 0.9. Recent experimental tests at EPRI have shown that higher emissivity aging assumptions may not be valid in some regions of the US, and that actual values may be much lower (0.25-0.45) than older studies have predicted (0.8-0.9). Here, the assumptions of different US Regional Transmission Operators (RTOs) and different utilities within RTOs are examined for the static ratings for traditional lower operating temperatures up to 100° C as well as with the higher maximum temperature conductors associated with newer conductor designs at 250° C. The sensitivity of the static ratings is examined at both lower and higher temperature conductors with respect to overhead line emissivity and absorptivity, wind speed, wind direction, and solar loading, and ambient temperature. It is shown that the effects of the conductor absorptivity, solar flux and ambient temperature on the static rating are increased when the transmission line has a lower maximum conductor temperature. The effects of the conductor emissivity, wind direction and wind speed have increased effects on the static rating when the transmission line has a higher maximum conductor temperature. In addition, example weather data is used to calculate the conductor temperature with different line emissivity and absorptivity assumptions. This shows the effect of assumptions made for emissivity only has a minor impact for low temperature lines, but for high temperature transmission lines can cause a 150-200° C temperature swing.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Benchmarking Demand Flexibility in Commercial Buildings and Flattening the Duck – Addressing Baseline and Commissioning Challenges

With the transition from our traditional electric grid to a cleaner grid with renewable power generation, there is a need to enable building loads to be flexible. Load shedding and shifting will be essential for flattening the “Duck” for decarbonization. This paper explored the trend in the timing of DR events as a reflection of the grid’s needs using recent four years of event data from 203 retail stores in 11 states. The events are becoming significantly shorter with 2-hour duration being the most popular; shifting to late afternoon and early evening is another trend beyond California. Benchmarking will be essential for accounting DF as a reliable grid resource. This paper addresses a challenging aspect of benchmarking – inaccuracies in counterfactual baseline methods can introduce significant DF metrics variations in addition to weather and building characteristics related factors. The conventional “10/10” with adjustment baseline method has inherent limitation by design for load shifting applications. Therefore, it is imperative to identify alternative methods. This study compared three hourly regression baseline methods with “10/10” methods using two groups of commercial buildings that participated in DR programs: (1) 121 big-box retail stores, and (2) 11 office buildings in CA. The 14-day hourly outdoor temperature regression method was found to produce least error in the tested datasets and is promising for load shifting. The paper also pointed out that commissioning issues can also be a significant barrier for achieving consistent DF performance, which building managers and utilities should be aware of.

Liu, Jingjing↗
DOE OSTI2025

Heat load measurements for the PIP-II pHB650 cryomodule

This study presents a brief overview of the 1st and 2nd phases and an in-depth analysis of the 3rd phase heat load testing performed on the pHB650 (prototype High Beta 650 MHz) cryomodule at PIP2IT (PIP-II Injector Test Facility), with a focus on both the results and the methodological advancements that have improved testing efficiency and accuracy. A key challenge identified in the testing campaign is the higher-than-expected heat loads observed in the first PIP-II (Proton Improvement Plan II) prototype cryomodules (pSSR1 and pHB650) tested at PIP2IT. Elevated heat loads are concerning given the fixed capacity of the PIP-II cryoplant that is currently being installed at Fermilab. However, understanding the sources of these elevated heat loads offers a critical opportunity to implement effective heat load mitigations on upcoming PIP-II cryomodules to stay within the available capacity of the PIP-II cryoplant. The study includes a summary of test results, descriptions of measurement procedures, and key observations on parameters directly and indirectly related to heat load measurements. Direct observations include measured heat loads and the effectiveness of JT heat exchanger under varying conditions, while indirect observation analyze factors such as the temperature distribution on the two-phase pipe and relief piping under varying conditions. Thermal acoustic oscillations (TAO) were identified during testing, which was mitigated by replacing the original G10 stem with a stainless steel stem equipped with wipers for the cryomodule cooldown valve. A major innovation during pHB650 Phase 3 testing was the development of an automated Python script to streamline data acquisition, analysis, and reporting of heat load results. This script automatically retrieved data from ACNET (Accelerator Control Network), performed heat load calculations, and generated detailed reports featuring plots and tables. This advancement significantly reduced manual labor and enhanced the thoroughness of data analysis compared to earlier campaigns. The heat load test reports were promptly uploaded to the electronic logbook shortly after each test, enabling rapid feedback and collaboration between the SRF and cryogenic teams. The heat load measurements included various components: HTTS (high-temperature thermal shield), LTTS (low-temperature thermal shield), 2K isothermal and non-isothermal heat loads. Results were recorded both within the cryomodule and between the bayonet can supply and return. Measurements were conducted under different operating conditions such as "standard", "linac", and "simulated dynamic". Additionally, HTTS and LTTS heat loads were calculated in real time, allowing for the tracking of thermal stability and identification of changes during testing, both in steady-state and transient conditions. The results of this testing campaign not only provide valuable insights into the performance of the pHB650 cryomodule but also highlight best practices and lessons learned that will inform future cryomodule testing at PIP2IT. These include adopting automated tools for data analysis, refining real-time measurement capabilities, and emphasizing detailed pre-test planning. The framework established in this campaign aims to set an improved standard for cryomodule testing and heat load reporting in future cryomodule test campaigns.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.]↗
DOE OSTI2020

Power Factor revisited

Power factor is a term whose value can be found by dividing the power in a load to the product of the voltage and the current. It has been used as a rule of thumb to estimate the cost of the losses incurred in delivering reactive power. But it has been extremely troublesome. The same load might be penalized by some definitions, and might not by others. That is still the situation today. There have been two efforts at solution. In 1919 a committee of experts was formed to solve the problem of its definition. The effort did not succeed at producing a definition; two alternatives were recommended in the committee’s final report, but both were adopted. In 2010, IEEE issued Standard 1459 to give definitions. In this standard, definitions of power factor appear under no less than seven different headings. The problem demands solution. The matter is of considerable economic importance to power producers, consumers, and instrument makers. This paper examines the history of the topic, and explores reasons for the difficulty of definition in terms of measurement theory. It is shown that power factor is not suited to purpose. In fact, no definition can exist that will define power factor usefully for its intended economic purpose. Two alternatives are proposed.

power factror, definitions, measurement theory, IE↗
DOE OSTI2024

Self‐Standing Chiral Covalent Organic Framework Thin Films with Full‐Color Tunable Guest‐Induced Circularly Polarized Luminescence

Abstract Exploring self‐standing chiral covalent organic framework (COF) thin films with controllable circularly polarized luminescence (CPL) is of paramount significance but remains a challenging task. Herein, we demonstrate the first example of self‐standing chiral COF films employing a polymerization‐dispersion‐filtration strategy. Pristine, low‐quality chiral COF films were produced by interfacial polymerization and then re‐dispersed into COF colloidal solutions. Via vacuum assisted assembly, these COF colloids were densely stacked and assembled into self‐standing, pure chiral COF films ( L‐ / D‐ CCOF‐F) that were transparent, smooth, crack‐free and highly crystalline. These films were tunable in thicknesses, areas, and roughness, along with strong diffuse reflectance circular dichroism (DRCD) and cyan CPL signals, showing an intrinsic luminescence asymmetric factor ( g lum ) of ~4.3×10 −3 . Furthermore, these COF films served as host adsorbents to load various achiral organic dye guests through adsorption. The effective chiral transfer and energy transfer between CCOF‐F and achiral fluorescent dyes endowed the dyes with strong chirality and tunable DRCD, resulting in intense, full‐color‐tunable solid‐state CPL. Notably, the ordered arrangement of dye guest molecules within the preferentially oriented chiral pores of CCOF‐F contributed to an amplified | g lum | factor of up to 7.2×10 −2 , which is state‐of‐the‐art for COF‐based CPL materials. This work provides new insights into the design and fabrication of self‐standing chiral COF films, demonstrating their great potential for chiroptical applications.

Tang, Xihao↗
DOE OSTI2021

Isotopic analysis of sub-nanogram neodymium loads using new ATONA™ amplifiers

The long-term gains of the ATONA™ amplifiers are shown to have less than 1 ppm variability. Repeat measurements of the JNdi-1 reference material demonstrate the ability of the ATONA™ amplification technology to make measurements of 143 Nd/ 144 Nd ratios with 23 ppm external reproducibility on 100-pg loads. The effect of increasing integration time on analytical reproducibility is also displayed as increasing integration time from 10 to 30s reduced the external measurement uncertainty from 37 to 23 ppm. These measurements represent an improvement of more than a factor of 3 in external measurement reproducibility relative to previously published 143 Nd/ 144 Nd measurements of 100-pg loads. Furthermore, this new technology will allow for the measurement of smaller samples for precise isotope ratios and open new avenues of research in the geological and nuclear forensic communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Direct ink writing of shear exfoliated two-dimensional nanomaterial- elastomeric multifunctional nanocomposite

Direct ink writing (DIW) of polymer nanocomposites with high loadings of two-dimensional (2D) nanofillers (graphene and hexagonal boron nitride (hBN)) is challenging because of potential clogging, use of hazardous solvents, and agglomeration. Here, in this work, a shear exfoliation and sieving method to prepare DIW ink with high loading of nanofillers produced from low-cost bulk layered materials such as graphite and bulk hBN powder for successful DIW printing without the use of any solvents, binders, or plasticizers. The single-step exfoliation technique resulted in a composite with substantial layer reduction along the c-axis, as confirmed by SEM, TEM, XRD, and Raman analysis. Incorporating exfoliated graphene (40 wt%) increased viscosity by ∼6 orders of magnitude due to enhanced particle–matrix interactions, leading to pronounced yield stress behavior and a yield stress of approximately 1598 Pa, which enabled excellent shape retention during extrusion. Using the DIW technique, porous structures such as desalination membranes, self-sensing bone scaffolds, thermal management coating, and serpentine strain sensors were fabricated. When tested in a direct contact membrane distillation setup, the fabricated membrane demonstrated a promising permeate flux of 21.85 Lm −2 h −1 and a salt rejection of 74.3 %. The fabricated serpentine sensor exhibited stable signal variations under cyclic tensile loading, with a working range of 0–200 % strain and a maximum gauge factor of 43,735. A cell culture test using the printed bone scaffold demonstrated promising cell attachment and proliferation. The DIW printed hBN nanocomposite exhibited reversible shape change under heat, demonstrating potential 4D printing capability and efficient thermal management when exposed to high heat or flame.

Desalination↗
DOE OSTI2020

Seasonal and Long‐Term Groundwater Unloading in the Central Valley Modifies Crustal Stress

Abstract Changes in terrestrial water content cause elastic deformation of the Earth's crust. This deformation is thought to play a role in modulating crustal stress and seismicity in regions where large water storage fluctuations occur. Groundwater is an important component of total water storage change in California, helping to drive annual water storage fluctuations and loss during periods of drought. Here we use direct estimates of groundwater volume loss during the 2007–2010 drought in California's Central Valley obtained from high resolution Interferometric Synthetic Aperture Radar‐based vertical land motion data to investigate the effect of groundwater volume change on the evolution of the stress field. We show that GPS‐derived elastic load models may not capture the contribution of groundwater to terrestrial water loading, resulting in an underestimation of nontectonic crustal stress change. We find that groundwater unloading during the drought causes Coulomb stress change of up to 5.5 kPa and seasonal fluctuations of up to 2.6 kPa at seismogenic depth. We find that faults near the Valley show the largest stress change and the San Andreas fault experiences only ~40 Pa of Coulomb stress change over the course of a year from groundwater storage change. Annual Coulomb stress change peaks dominantly in the fall, when the groundwater level is low; however, some faults experience peak stress in the spring when groundwater levels are higher. Additionally, we find that periods of increased stress correlate with higher than average seismic moment release but are not correlated with an increase in the number of earthquakes. This indicates groundwater loading likely contributes to nontectonic loading of faults, especially near the Valley edge, but is not a dominant factor in modulation of seismicity in California because the amplitude of stress change declines rapidly with distance from the Valley. By carefully quantifying and spatially locating groundwater fluctuations, we will improve our understanding of what drives nontectonic stress and forces that modulate seismicity in California.

58 GEOSCIENCES↗