Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity profiles”

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 199 records · Page 11

MicroFisher: Fungal taxonomic classification for metatranscriptomic and metagenomic data using multiple short hypervariable markers

AbstractProfiling the taxonomic and functional composition of microbes using metagenomic (MG) and metatranscriptomic (MT) sequencing is advancing our understanding of microbial functions. However, the sensitivity and accuracy of microbial classification using genome– or core protein-based approaches, especially the classification of eukaryotic organisms, is limited by the availability of genomes and the resolution of sequence databases. To address this, we propose the MicroFisher, a novel approach that applies multiple hypervariable marker genes to profile fungal communities from MGs and MTs. This approach utilizes the hypervariable regions of ITS and large subunit (LSU) rRNA genes for fungal identification with high sensitivity and resolution. Simultaneously, we propose a computational pipeline (MicroFisher) to optimize and integrate the results from classifications using multiple hypervariable markers. To test the performance of our method, we applied MicroFisher to the synthetic community profiling and found high performance in fungal prediction and abundance estimation. In addition, we also used MGs from forest soil and MTs of root eukaryotic microbes to test our method and the results showed that MicroFisher provided more accurate profiling of environmental microbiomes compared to other classification tools. Overall, MicroFisher serves as a novel pipeline for classification of fungal communities from MGs and MTs.

Wang, Haihua↗

Effects of 63 MeV proton irradiation on the performance of MWIR InGaAs/InAsSb nBn photodetectors

An investigation into the effects of 63 MeV proton irradiation on high-sensitivity mid-wave infrared InGaAs/InAsSb nBn devices is performed. Three different structures with various absorber region doping profiles are irradiated and characterized to assess their impact on performance degradation. Minority carrier lifetime is measured using time-resolved photoluminescence and lifetime damage factors are assessed. The majority carrier concentration is determined via capacitance–voltage measurements and dopant introduction rates are calculated. An analysis of dark current density is performed using these material parameters, revealing a reduction in mobility with proton fluence and the emergence of a proton-induced trap energy level. Quantum efficiency is calculated at each proton fluence, and quantum efficiency damage factors show that the graded doping structure exhibits the least reduction of quantum efficiency with dose, attributed to its effective mobility enhancement. Conclusively, detector sensitivity, assessed via shot-noise limited noise-equivalent irradiance, shows that the graded doping structure is the least susceptible to high energy proton irradiation-induced performance degradation.

Physics↗

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗

Behavior and mechanisms of Doppler wind lidar error in complex terrain: stable flow case study at Perdigão

A numerical experiment is carried out investigating the magnitude of biases in ground-based lidar measurements in complex flow conditions. Biases assessed include those arising from flow curvature and from the interaction of turbulence with the wind field reconstruction (WFR) algorithms used by a WindCube lidars and anemometers. RANS-CFD and WRF-LES simulations were performed for the Perdig˜ao Field Experiment site for a range of atmospheric conditions. Virtual anemometer and lidar data were generated for four locations: two near exposed ridge tops and two in low-speed regions in the valley. The LES data at these four locations show that the scalar inflation terms (the relation between scalar and vector averaged wind speed) for virtual lidar and virtual cups agree very well with predictions using perturbation theory. While the lidar errors vary greatly with location and height, the contribution from the flow curvature tends to be larger than the differences arising from scalar inflation. For one lidar/mast pair near the ridge top, comparisons between simulations and measurements are carried out for a resonant mountain wave event on June 14th, 2017, and for the whole duration of the Perdigão campaign for winds perpendicular to the ridges. The lidar error during the mountain wave, a period of strong stability and low inversion height, is significantly larger than the campaign average. The sensitivity of the lidar error to atmospheric stability is confirmed by the RANS simulations, which suggests strong sensitivity of flow curvature error to stability conditions and to the shape of the wind speed profile near the top of the boundary layer.

17 WIND ENERGY↗

Impact Analysis of Data Integrity Attacks on FACTS-based Wide-Area Voltage Control System

Energy management system (EMS) consists of several wide-area control applications that serve as a backbone for security, stability, and reliability of the power system. Wide-area voltage control system (WAVCS), one of the critical wide-area applications, operates in coordination with local Flexible AC Transmission System (FACTS) devices to provide voltage security and optimal management of active and reactive power resources. Since the WAVCS relies on wide-area communication and data sharing devices, possible cybersecurity vulnerabilities have to be addressed to ensure the closed-loop operation of WAVCS. In this paper, we present a methodology for performing an impact analysis of cyber-attacks in WAVCS cybersecurity. In particular, different types of data integrity attacks, such as malicious tripping, fault replay, and signal altering attacks, are considered, and detailed impact analysis is conducted in a testbed environment using the Kundur's four machine two-area system. For performing an impact analysis, the transient voltage stability of the sensitive bus voltage is studied, followed by the quantitative assessment and severity ranking using the voltage profile index. Our experimental evaluation reveals that the data integrity attacks on control signals exhibit a higher attack severity than on the measurement signals. Further, the severity of these attacks varies with nature (static or dynamic), location, and types of attacks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Hybrid chemical characterization of latent images in EUV resist with 12 nm half-pitch features

With the advancement of high numerical aperture extreme ultraviolet (EUV) lithography, the new platform will enable chipmakers to achieve critical dimensions of 8 nm. However, resist materials face significant challenges in delivering increased sensitivity while managing rising stochastic variations. We aim to develop comprehensive techniques to characterize the chemical profile of latent images, stored in EUV resists after exposure and postexposure baking, which is essential for understanding the origin of stochastic effects. Infrared photo-induced force microscopy (IR PiFM) is a bimodal atomic force microscopy technique combined with an infrared light source, allowing for simultaneous sub-5 nm topographic and chemical characterization within a localized environment. Critical-dimension resonant soft X-ray scatterometry (CD-RSoXS) provides statistical data that reveal structural and chemical information for comparative analysis. For the first time, IR PiFM has been used to chemically map the latent images of EUV resists (after exposure and postexposure baking) at a record high resolution of 12 nm half-pitch, enabling nondestructive analysis of patterns prior to development. Furthermore, CD-RSoXS offers direct experimental observation and comparison of exposed, postexposure baked, and developed patterns, which align with the IR PiFM results. We demonstrate that the IR PiFM technique offers valuable insights into both high spatial resolution and local chemical information simultaneously. In addition, CD-RSoXS provides statistical results that support our main findings. This hybrid metrology approach leverages a multifaceted dataset by combining the most reliable information from each source, which is essential for a comprehensive understanding of the stochastic effects in EUV lithography processes.

O’Reilly, Padraic↗

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗

Compatibility of Fe-Cr-Al and Fe-Cr-Al-Mo oxide dispersion strengthened steels with static liquid sodium at 700 °C

We report the corrosion behavior and tensile properties of oxide dispersion strengthened (ODS) FeCr, FeCrAl, FeCrAlMo steels and 316H stainless steel were evaluated after 1000 h at 700 °C in static sodium with and without an oxygen getter (Zr foil). Without Zr foil to getter the dissolved oxygen in the liquid sodium, combined oxidation (formation of Cr, Fe-rich or Al-rich oxides) and dissolution of the studied materials was observed. For the ODS FeCr, Cr- and Fe-rich oxide formation and dissolution induced Cr depletion in the underlying alloy and led to the formation of laves phases that were observed to be detrimental to the room temperature tensile properties. For ODS FeCrAlMo and FeCrAl Al-rich oxide formation and dissolution resulted in Al depletion in the underlying alloys but did not impact tensile properties at room temperature. For 316H temperature induced sensitization combined with oxidation decreased room temperature tensile properties. With Zr foil added as oxygen getter, roughened surfaces and Al depletion profiles suggested the dissolution of ODS FeCr and FeCrAl into the liquid sodium while limited dissolution was observed for ODS FeCrAlMo and 316H. For the 316H and ODS FeCrAlMo steels, precipitates were formed at the grain boundaries after exposure. Tensile properties of ODS FeCr, FeCrAl and FeCrAlMo were not affected by the 1000 h corrosion exposure with Zr foil added to the sodium but for 316H temperature-induced grain boundary sensitization negatively affected the tensile properties.

36 MATERIALS SCIENCE↗

TATB thermal decomposition: Expanding the molecular profile with cryo-focused pyrolysis GC-MS

Understanding the molecular composition of high explosives during thermal decomposition is vital for predicting the sensitivity, safety, and performance of explosive materials. The thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has been linked to the formation of furazans through a series of dehydration reactions of the NO2 and NH2 groups on the phenyl ring, along with breakdown into small molecules (≤120 amu). Molecular identification of compounds formed in this transformation of the furazans to light gases has been lacking. To address this, we have applied a pseudo-confined sampling system in a cryo-focused pyrolysis gas chromatography-mass spectrometry (pyGC-MS) system to molecularly identify these intermediates. By design, sublimation of TATB, which has complicated MS analyses of thermal degradation, was significantly reduced and additional compounds were identified with potential structural information. In addition to the known furazan compounds, one of these compounds forms from the loss of oxygen from benzo-trifurazan (F3) and produces an open ring structure that may be the first step in the formation of lower molecular weight furazan breakdown products. The loss of a nitro group from benzo-monofurazan (F1) was also discovered and implicates the formation of oxidizing NO2 gas in the thermal decomposition mechanism. So these findings are vital for understanding the proper heat flow from energetic materials on a molecular level, necessary when measuring enthalpy and developing decomposition models based on kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Production of hydrogen peroxide in an intra-meander hyporheic zone at East River, Colorado.

The traditionally held assumption that photo-dependent processes are the predominant source of H 2 O 2 in natural waters has been recently questioned by an increasing body of evidence showing the ubiquitiousness of H 2 O 2 in dark water bodies and in groundwater. In this study, we conducted field measurement of H 2 O 2 in an intra-meander hyporheic zone and in surface water at East River, CO. On-site detection using a sensitive chemiluminescence method suggests H 2 O 2 concentrations in groundwater ranging from 6 nM (at the most reduced region) to ~ 80 nM (in a locally oxygen-rich area) along the intra-meander transect with a maxima of 186 nM detected in the surface water in an early afternoon, lagging the maximum solar irradiance by ~1.5 h. Our results suggest that the dark profile of H 2 O 2 in the hyporheic zone is closely correlated to local redox gradients, indicating that interactions between various redox sensitive elements could play an essential role. Due to its transient nature, the widespread presence of H 2 O 2 in the hyporheic zone indicates the existence of a sustained balance between H 2 O 2 production and consumption, which potentially involves a relatively rapid succession of various biogeochemically important processes (such as organic matter turnover, metal cycling and contaminant mobilization). More importantly, this study confirmed the occurrence of reactive oxygen species at a subsurface redox transition zone and further support our understanding of redox boundaries on reactive oxygen species generation and as key locations of biogeochemical activity.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Predicting Frigid Mixed-Phase Clouds for Pristine Coastal Antarctica

Supercooled water is common in the clouds near coastal Antarctica and occasionally occurs at temperatures at or below -30°C. Yet the ice physics in most regional and global numerical models will glaciate out these clouds. This presents a challenge for the simulation of highly supercooled clouds that were observed at McMurdo, Antarctica during the Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) project during 2015–2017. The polar optimized version of the Weather Research and Forecasting model (Polar WRF) with the recently developed two-moment P3 microphysics scheme was used to simulate observed supercooled liquid water cases during March and November 2016. Nudging of the simulations to observed rawinsonde profiles and Antarctic automatic weather station observations provided increased realism and much greater cloud water amounts. Furthermore, sensitivity tests that adjust the ice physics for extremely low ice nucleating particle (INP) concentrations decrease cloud ice and increases the cloud liquid water closer to observed amounts. In these tests, a liquid layer near cloud top is simulated, in agreement with observations. Accurate representation of INP concentrations appears to be critical for the simulation of coastal Antarctic clouds.

54 ENVIRONMENTAL SCIENCES↗

Trajectory of a flying plasma mirror traversing a target with density gradient

It has been proposed that laser-induced relativistic plasma mirror can accelerate if the plasma has a properly tailored density profile. Such accelerating plasma mirrors can serve as analog black holes to investigate Hawking evaporation and the associated information loss paradox. Here we reexamine the underlying dynamics of mirror motion in a graded-density plasma to provide an explicit trajectory as a function of the plasma density and its gradient. Specifically, a decreasing plasma density profile (down-ramp) along the direction of laser propagation would in general accelerate the mirror. In particular, a constant-plus-exponential density profile would generate the Davies–Fulling trajectory with a well-defined analog Hawking temperature, which is sensitive to the plasma density gradient but not to the density itself. Furthermore, we show that without invoking nano-fabricated thin-films, a much lower density gas target at, for example, ~1 × 10 17 cm -3 , would be able to induce an analog Hawking temperature, k B T H ~3.1 × 10 -2 eV, in the far-infrared region. We hope that this would help to better realize the experiment proposed by Chen and Mourou.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Timing characterization of fast hCMOS sensors

We describe a method of analyzing gate profile data for ultrafast x-ray imagers that allows pixel-by-pixel determination of temporal sensitivity in the presence of substantial background oscillations. With this method, systematic timing errors in gate width and gate arrival time of up to 1 ns (in a 2 ns wide gate) can be removed. In-sensor variations in gate arrival and gate width are observed, with variations in each up to 0.5 ns. This method can be used to estimate the coarse timing of the sensor, even if errors up to several ns are present.

Benedetti, L. R. (ORCID:0000000160502983)↗

Realization of a gas puff imaging system on the Wendelstein 7-X stellarator

A system for studying the spatiotemporal dynamics of fluctuations in the boundary of the W7-X plasma using the “Gas-Puff Imaging” (GPI) technique has been designed, constructed, installed, and operated. This GPI system addresses a number of challenges specific to long-pulse superconducting devices, such as W7-X, including the long distance between the plasma and the vacuum vessel wall, the long distance between the plasma and diagnostic ports, the range of last closed flux surface (LCFS) locations for different magnetic configurations in W7-X, and management of heat loads on the system’s plasma-facing components. The system features a pair of “converging–diverging” nozzles for partially collimating the gas puffed locally ≈135 mm radially outboard of the plasma boundary, a pop-up turning mirror for viewing the gas puff emission from the side (which also acts as a shutter for the re-entrant vacuum window), and a high-throughput optical system that collects visible emission resulting from the interaction between the puffed gas and the plasma and directs it along a water-cooled re-entrant tube directly onto the 8 × 16 pixel detector array of the fast camera. The DEGAS 2 neutral code was used to simulate the Hα (656 nm) and HeI (587 nm) line emission expected from well-characterized gas-puffs of H2 and He and excited within typical edge plasma profiles in W7-X, thereby predicting line brightnesses used to reduce the risks associated with system sensitivity and placement of the field of view. Operation of GPI on W7-X shows excellent signal-to-noise ratios (>100 at 2 Mframes/s) over the field of view for minimally perturbing gas puffs. The GPI system provides detailed measurements of the two-dimensional (radial and poloidal) dynamics of plasma fluctuations in the W7-X edge and scrape-off layer and in and around the magnetic islands outside the LCFS that make up the island divertor configuration employed on W7-X.

Instruments & Instrumentation↗

Observation of a cross-section enhancement near the production threshold in $\sqrt{s} = 13$ TeV pp collisions with the ATLAS detector

A measurement of $t\bar{t}$ production is presented in the invariant-mass region near the pair production threshold, $m_{t\bar{t}} \sim 345$ GeV, in final states with two charged leptons and multiple jets. The measurement is based on $140\,\mathrm{fb}^{-1}$ of proton–proton collision data collected at $\sqrt{s} = 13$ TeV with the ATLAS detector at the Large Hadron Collider. The data are compared to two models of $t\bar{t}$ production: a baseline model including only perturbative quantum chromodynamics (pQCD) predictions for the hard process at approximate next-to-next-to leading order accuracy in the strong coupling, and an extended model that, in addition, incorporates non-relativistic QCD simulations that also include the formation of colour-singlet quasi-bound-states near the $t\bar{t}$ threshold. The agreement between the data and the models is quantified via a profile-likelihood fit to the reconstructed $m_{t\bar{t}}$ distributions, in bins of two angular observables sensitive to spin-correlations in the $t\bar{t}$ system. An excess of events is observed over the baseline pQCD prediction, with an observed significance over 8 standard deviations. This excess is consistent with the formation of colour-singlet and spin-singlet S-wave quasi-bound $t\bar{t}$ states, as predicted by non-relativistic QCD, and corresponds to an observed cross-section of $9.3^{+1.4}_{-1.3}$ pb.

Large Hadron Collider↗

Frequentist cosmological constraints from full-shape clustering measurements in DESI DR1

We present a frequentist analysis of clustering measurements from Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI) using the standard profile likelihood method. While Bayesian inferences for effective field theory models of galaxy clustering can be highly sensitive to prior choices for extended cosmological models, frequentist inferences are not susceptible to such effects. We compare frequentist and Bayesian constraints for the parameter set {σ 8 , H 0 , Ω m , w 0 , w a } using the full-shape power spectrum multipoles, post-reconstruction baryon acoustic oscillation (BAO) measurements, and external datasets from the CMB and type Ia supernovae measurements. The frequentist confidence intervals are significantly shifted relative to the Bayesian credible intervals for the w 0 w a CDM model, unless supernovae data are included. When DESI full-shape and BAO data are fit jointly, we obtain the following 1σ frequentist confidence intervals for ΛCDM (w 0 w a CDM): σ 8 = 0.863 +0.048 -0.040 , H 0 = 68.96 +0.81 -0.80 km s -1 Mpc -1 , Ω m = 0.3034 ± 0.0110 (σ 8 = 0.782 +0.060 -0.036 , H 0 = 63.7 +4.2 -2.0 km s -1 Mpc -1 , Ω m = 0.378 +0.024 -0.047 , w 0 = -0.16 +0.10 -0.50 , w a = -3.0 +1.7 ), corresponding to 0.8σ, 0.3σ, 0.7σ (2.1σ, 4.1σ, 6.5σ, 6.3σ, 6.6σ) shifts between the maximum likelihood estimate and the Bayesian posterior mean for ΛCDM (w 0 w a CDM) respectively.

Bayesian reasoning↗

Cherenkov and scintillation separation in water-based liquid scintillator using an LAPPDTM

Abstract This manuscript describes measurements of water-based liquid scintillators (WbLS), demonstrating separation of the Cherenkov and scintillation components using a low energy $$\beta $$ β source and the fast timing response of a Large Area Picosecond Photodetector (LAPPD). Additionally, the time profiles of three WbLS mixtures, defined by the relative fractions of scintillating compound, are characterized, with improved sensitivity to the scintillator rise-time. The measurements were made using both an LAPPD and a conventional photomultiplier tube (PMT). All samples were measured with an effective resolution $$O\left( 100~\text {ps}\right) $$ O 100 ps , which allows for the separation of Cherenkov and scintillation light (henceforth C/S separation) by selecting on the arrival time of the photons alone. The Cherenkov purity of the selected photons is greater than 60% in all cases, with greater than 80% achieved for a sample containing 1% scintillator. This is the first demonstration of the power of synthesizing low light yield scintillators, of which WbLS is the canonical example, with fast photodetectors, of which LAPPDs are an emerging leader, and has direct implication for future mid- and large-scale detectors, such as Theia , ANNIE, and AIT-NEO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-frequency and Vertical Variations of Soil Moisture and their Impact on Ecosystem Dynamics

The vertical distribution of subsurface moisture is a major control on ecosystem dynamics and climate. The Mediterranean climate of California, with wet winters and dry summers, provides a testbed for understanding how transpiration is sustained through extended dry periods. The project developed and tested a new subsurface hydrology model that includes a parameterization of stochastic hydraulic conductivity to capture fracture flow, a diagnostic model of rooting profile that depends on the subsurface hydrology, and a representation of hydraulic redistribution by plant roots. A sensitivity study shows that it is plants with rooting depths >8 m that could access the subsurface moisture. The insight is applied to understand differential tree mortality during the 2012-2016 drought in California. The analysis shows that greater mortality trees on plutonic rocks than on sedimentary rocks. The study demonstrates, for the first time, the importance of lithology for the distribution and health of ecosystems.

54 ENVIRONMENTAL SCIENCES↗