Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “block structure”

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 1,675 records · Page 93

Algorithms for parallel flow solvers on message passing architectures

The purpose of this project has been to identify and test suitable technologies for implementation of fluid flow solvers -- possibly coupled with structures and heat equation solvers -- on MIMD parallel computers. In the course of this investigation much attention has been paid to efficient domain decomposition strategies for ADI-type algorithms. Multi-partitioning derives its efficiency from the assignment of several blocks of grid points to each processor in the parallel computer. A coarse-grain parallelism is obtained, and a near-perfect load balance results. In uni-partitioning every processor receives responsibility for exactly one block of grid points instead of several. This necessitates fine-grain pipelined program execution in order to obtain a reasonable load balance. Although fine-grain parallelism is less desirable on many systems, especially high-latency networks of workstations, uni-partition methods are still in wide use in production codes for flow problems. Consequently, it remains important to achieve good efficiency with this technique that has essentially been superseded by multi-partitioning for parallel ADI-type algorithms. Another reason for the concentration on improving the performance of pipeline methods is their applicability in other types of flow solver kernels with stronger implied data dependence. Analytical expressions can be derived for the size of the dynamic load imbalance incurred in traditional pipelines. From these it can be determined what is the optimal first-processor retardation that leads to the shortest total completion time for the pipeline process. Theoretical predictions of pipeline performance with and without optimization match experimental observations on the iPSC/860 very well. Analysis of pipeline performance also highlights the effect of uncareful grid partitioning in flow solvers that employ pipeline algorithms. If grid blocks at boundaries are not at least as large in the wall-normal direction as those immediately adjacent to them, then the first processor in the pipeline will receive a computational load that is less than that of subsequent processors, magnifying the pipeline slowdown effect. Extra compensation is needed for grid boundary effects, even if all grid blocks are equally sized.

Vanderwijngaart, Rob F.↗

Highly Stable Fe/Co-TPY-MIL-88(NH 2 ) Metal–Organic Framework (MOF) in Enzymatic Cascade Reactions for Chemiluminescence-Based Detection of Extracellular Vesicles

Metal-Organic Frameworks (MOFs) can deliver many advantages when acting as enzyme mimics to assist with signal amplification in molecular detection: they have abundant active catalytic sites per unit volume of the material; their structures and elemental compositions are highly tunable, and their high specific surface area and porous property can assist with target separation and enrichment. In the present work, we have demonstrated that, by adding the pore partition agent, 2,4,6-tris(4-pyridyl)pyridine (TPY) during synthesis of the bimetallic Fe/Co-MIL-88(NH 2 ) MOF to block the open metal sites, a highly porous MOF of Fe/Co-TPY-MIL-88(NH 2 ) can be produced. This material also exhibits high stability in basic solutions and biofluids and possesses high peroxidase-mimicking activity, which can be utilized to produce long-lasting chemiluminescence (CL) from luminol and H 2 O 2 . Moreover, acting as the peroxidase-mimic, the Fe/Co-TPY-MIL-88(NH 2 ) MOF can form the enzymatic cascade with glucose oxidase (GOx) for biomarker detection. Further, when applied to detect extracellular vesicles (EVs), the MOF material and GOx are brought to the proximity on the EVs through two surface proteins, which triggers the enzyme cascade to produce high CL from glucose and luminol. EVs within the concentration range of 5 × 10 5 - 4 ×10 7 particles/mL can be detected, and the method can be used to analyze EV contents in human serum without sample preparation and EV purification. Overall, our work demonstrates that the high versatility and tunability of the MOF structures could bring in significant benefits to biosensing and enable ultrasensitive detection of biomarkers with judicious material designs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conformation-specific inhibitors of activated Ras GTPases reveal limited Ras dependency of patient-derived cancer organoids

The small GTPases H, K, and NRAS are molecular switches indispensable for proper regulation of cellular proliferation and growth. Several mutations in the genes encoding members of this protein family are associated with cancer and result in aberrant activation of signaling processes caused by a deregulated recruitment of downstream effector proteins. In this study, we engineered variants of the Ras-binding domain (RBD) of the C-Raf proto-oncogene, Ser/Thr kinase (CRAF). These variants bound with high affinity with the effector-binding site of Ras in an active conformation. Structural characterization disclosed how the newly identified RBD mutations cooperate and thereby enhance affinity with the effector-binding site in Ras compared with WT RBD. The engineered RBD variants closely mimicked the interaction mode of naturally occurring Ras effectors and acted as dominant-negative affinity reagents that block Ras signal transduction. We report experiments with cancer cells showed that expression of these RBD variants inhibits Ras signaling, reducing cell growth and inducing apoptosis. Using these optimized RBD variants, we stratified patient-derived colorectal cancer organoids with known Ras mutational status according to their response to Ras inhibition. These results revealed that the presence of Ras mutations was insufficient to predict sensitivity to Ras inhibition, suggesting that not all of these tumors required Ras signaling for proliferation. In summary, by engineering the Ras/Raf interface of the CRAF-RBD, we identified potent and selective inhibitors of Ras in its active conformation that outcompete binding of Ras-signaling effectors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stress field rotation or block rotation: An example from the Lake Mead fault system

The Coulomb criterion, as applied by Anderson (1951), has been widely used as the basis for inferring paleostresses from in situ fault slip data, assuming that faults are optimally oriented relative to the tectonic stress direction. Consequently if stress direction is fixed during deformation so must be the faults. Freund (1974) has shown that faults, when arranged in sets, must generally rotate as they slip. Nur et al., (1986) showed how sufficiently large rotations require the development of new sets of faults which are more favorably oriented to the principal direction of stress. This leads to the appearance of multiple fault sets in which older faults are offset by younger ones, both having the same sense of slip. Consequently correct paleostress analysis must include the possible effect of fault and material rotation, in addition to stress field rotation. The combined effects of stress field rotation and material rotation were investigated in the Lake Meade Fault System (LMFS) especially in the Hoover Dam area. Fault inversion results imply an apparent 60 degrees clockwise (CW) rotation of the stress field since mid-Miocene time. In contrast structural data from the rest of the Great Basin suggest only a 30 degrees CW stress field rotation. By incorporating paleomagnetic and seismic evidence, the 30 degrees discrepancy can be neatly resolved. Based on paleomagnetic declination anomalies, it is inferred that slip on NW trending right lateral faults caused a local 30 degrees counter-clockwise (CCW) rotation of blocks and faults in the Lake Mead area. Consequently the inferred 60 degrees CW rotation of the stress field in the LMFS consists of an actual 30 degrees CW rotation of the stress field (as for the entire Great Basin) plus a local 30 degrees CCW material rotation of the LMFS fault blocks.

Ron, Hagai↗

Optical Design of the STAR-X Telescope

Top-level science goals of the Survey and Time-domain Astrophysical Research eXplorer (STAR-X) include: investigations of most violent explosions in the universe, study of growth of black holes across cosmic time and mass scale, and measure how structure formation heats majority of baryons in the universe. To meet these goals, the field-of-view of the telescope should be about 1 square-degree, the angular resolution should be 5 arc-seconds or below across large part of the field-of-view. The on-axis effective area at 1 KeV should be about 2,000 sq cm. Payload cost and launch considerations limit the outer diameter, focal length, and mass to 1.3 meters, 5 meters, and 250 kilograms, respectively. Telescope design is based on a segmented meta-shell approach we have developed at Goddard Space Flight Center for the STAR-X telescope. The telescope shells are divided into 30-degree segments. Individual telescopes and meta-shells are nested inside each other to meet the effective area requirements in 0.5 - 6.0 KeV range. We consider Wolter-Schwarzschild, and Modified-Wolter-Schwarzschild telescope designs as basic building blocks of the nested STAR-X telescope. These designs offer an excellent resolution over a large field of views. Nested telescopes are vulnerable to stray light problems. We have designed a multi-component baffle system to eliminate direct and single-reflection light paths inside the telescopes. Large number of internal and external baffle vane structures are required to prevent stray rays from reaching the focal plane. We have developed a simple ray-trace based tool to determine the dimensions and locations of the baffles. In this paper, we present the results of our trade studies, baffle design studies, and optical performance analyses of the STAR-X telescope.

Saha, Timo T.↗

Design of the STAR-X Telescope

Top-level science goals of the Survey and Time-domain Astrophysical Research eXplorer (STAR-X) include: investigations of most violent explosions in the universe, study of growth of black holes across cosmic time and mass scale, and measure how structure formation heats majority of baryons in the universe. To meet these goals, the field-of-view of the telescope should be about 1 square-degree, the angular resolution should be 5 arc-seconds or below across large part of the field-of-view. The on-axis effective area at 1 KeV should be about 2,000 sq cm. Payload cost and launch considerations limit the outer diameter, focal length, and mass to 1.3 meters, 5 meters, and 250 kilograms, respectively. Telescope design is based on a segmented meta-shell approach we have developed at Goddard Space Flight Center for the STAR-X telescope. The telescope shells are divided into 30-degree segments. Individual telescopes and meta-shells are nested inside each other to meet the effective area requirements in 0.5 - 6.0 KeV range. We consider Wolter-Schwarzschild, and Modified-Wolter-Schwarzschild telescope designs as basic building blocks of the nested STAR-X telescope. These designs offer an excellent resolution over a large field of views. Nested telescopes are vulnerable to stray light problems. We have designed a multi-component baffle system to eliminate direct and single-reflection light paths inside the telescopes. Large number of internal and external baffle vane structures are required to prevent stray rays from reaching the focal plane. We have developed a simple ray-trace based tool to determine the dimensions and locations of the baffles. In this paper, we present the results of our trade studies, baffle design studies, and optical performance analyses of the STAR-X telescope.

Saha, Timo T.↗

Improved chemical and isotopic labeling of biomembranes in Bacillus subtilis by leveraging CRISPRi inhibition of beta-ketoacyl-ACP synthase (fabF)

Assessing the structure of living microbial cell membranes is a challenging analytical goal. The cell membrane is defined by its transverse structure, an approximately 5 nm-thick selectively permeable bilayer that serves many important cellular functions. Compositionally complex, dynamic, and organized in both the transverse and lateral dimensions, understanding the cell membrane structure—and the role that structure plays in cellular function, communication, and environmental sensing is an active scientific effort. Previously, we have devised a novel isotopic labeling approach for membrane lipids to enable direct in vivo structural studies of the cell membrane in the Gram-positive bacterium, Bacillus subtilis , using small-angle neutron scattering. This was accomplished through a genetic inhibition of fatty acid (FA) degradation (Δ fadN ) and a chemical inhibition of FA biosynthesis using cerulenin, an irreversible inhibitor of type II fatty acid synthases. Here, we improve upon the previous system by introducing a dCas9/sgRNA- fabF complex that blocks transcription of the essential fabF gene when under xylose induction. This leads to greater sensitivity to cerulenin in the mutant strain (JEBS102) and more robust cell growth when supplementary FAs are introduced to the culture medium. A subtle change in FA uptake is noted when compared to the prior labeling strategy. This is seen in the gas chromatography/mass spectrometry (GC/MS) data as a higher ratio of n 16:0 to a 15:0, and manifests in an apparent increase in the membrane thickness determined via neutron scattering. This represents an improved method of isotopic labeling for the cell membrane of Bacillus subtilis; enabling improved investigations of cellular uptake and utilization of FAs, cell membrane structure and organization as a phenotypic response to metabolic and environmental changes.

59 BASIC BIOLOGICAL SCIENCES↗

Feynman rules for scalar conformal blocks

We complete the proof of “Feynman rules” for constructing M-point conformal blocks with external and internal scalars in any topology for arbitrary M in any spacetime dimension by combining the rules for the blocks (based on their Witten diagram interpretation) with the rules for the construction of conformal cross ratios (based on the OPE and “flow diagrams”). The full set of Feynman rules leads to blocks as power series of the hypergeometric type in the conformal cross ratios. We then provide a proof by recursion of the Feynman rules which relies heavily on the first Barnes lemma and the decomposition of the topology of interest in comb structures. Finally, we provide a nine-point example to illustrate the rules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The role of dissolved pyrogenic carbon from biochar in the sorption of As(V) in biogenic iron (oxyhydr)oxides

Water contamination by arsenic (As) affects millions of people around the world, making techniques to immobilize or remove this contaminant a pressing societal need. Biochar and iron (oxyhydr)oxides [in particular, biogenic iron (oxyhydr)oxides (BIOS)] offer the possibility of stabilizing As in remediation systems. However, little is known about the potential antagonism in As sorption generated by the dissolved organic carbon (DOC) from biochar, or whether DOC affects how As(V) interacts with BIOS. For this reason, our objectives were to evaluate the i) As(V) sorption potential in BIOS when there is presence of DOC from pyrolyzed biochars at different temperatures; and ii) identify whether the presence of DOC alters the surface complexes formed by As(V) sorbed in the BIOS. We conducted As(V) sorption experiments with BIOS at circumneutral pH conditions and in the presence of DOC from sugarcane (Saccharum officinarum) straw biochar at pyrolyzed 350 (BC350) and 750 °C (BC750). The As(V) content was quantified by inductively coupled plasma mass spectrometry, and the BIOS structure and As(V) sorption mechanisms were investigated by X-ray absorption spectroscopy. In addition, the organic moieties comprising the DOC from biochars were investigated by attenuated total reflectance Fourier transform infrared spectroscopy. The addition of DOC did not change the biomineral structure or As(V) oxidation state. The presence of DOC, however, reduced by 25 % the sorption of As(V), with BC350 being responsible for the greatest reduction in As(V) sorption capacity. Structural modeling revealed As(V) predominantly formed binuclear bidentate surface complexes on BIOS. The presence of DOC did not change the binding mechanism of As(V) in BIOS, suggesting that the reduction of As(V) sorption to BIOS was due to site blocking. In conclusion, our results bring insights into the fate of As(V) in surface waters and provide a basis for understanding the competitive sorption of As(V) in environments with biochar application.

54 ENVIRONMENTAL SCIENCES↗

Shape Control of Colloidal Semiconductor Nanocrystals through Thermodynamically Driven Aggregative Growth

The optoelectronic properties of colloidal semiconductor nanocrystals (NCs) can be manipulated by changing their geometric shapes. The precise synthetic control over particle morphologies, however, has remained elusive. Conventional growth techniques rely on the kinetic assembly of atomic units, where supersaturation and precipitation processes can lead to a broad distribution of particle shapes. In this paper, we demonstrate that replacing atomic precursors with small-size nanocrystals as building blocks for larger colloids offers an easier, more predictive control over nanoparticle shape evolution. The reported growth strategy is illustrated via shape-selective syntheses of CdSe and CdS NC cubes, spheres, rods, as well as unprecedented “donut” and ring-like structures. Growth of desired particle morphologies was achieved by choosing an appropriate reaction solvent that minimizes the surface free energy through a unique shape formation. This so-called "aggregative growth mechanism" is explained using a thermodynamic model for interacting viscous colloids. We expect that this novel growth strategy will offer a viable platform for shape–selective synthesis of many inorganic colloids, providing pathways to optimizable optoelectronic materials.

36 MATERIALS SCIENCE↗

Room-Temperature Ferroelectric Epitaxial Nanowire Arrays with Photoluminescence

The development of large-scale, high-quality ferroelectric semiconductor nanowire arrays with interesting light-emitting properties can address limitations in traditional wide-bandgap ferroelectrics, thus serving as building blocks for innovative device architectures and next-generation high-density optoelectronics. Here, we investigate the optical properties of ferroelectric CsGeX 3 (X = Br, I) halide perovskite nanowires that are epitaxially grown on muscovite mica substrates by vapor phase deposition. Detailed structural characterizations reveal an incommensurate heteroepitaxial relationship with the mica substrate. Furthermore, photoluminescence that can be tuned from yellow-green to red emissions by varying the halide composition demonstrates that these nanowire networks can serve as platforms for future optoelectronic applications. In addition, the room-temperature ferroelectricity and ferroelectric domain structures of these nanowires are characterized using second harmonic generation (SHG) polarimetry. Finally, the combination of room-temperature ferroelectricity with photoluminescence in these nanowire arrays unlocks new avenues for the design of novel multifunctional materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Accumulation of Large Ion Populations with High Ion Densities and Effects Due to Space Charge in Traveling Wave-Based Structures for Lossless Ion Manipulations (SLIM) IMS-MS

The accumulation of very large ion populations in traveling wave (TW)-based Structures for Lossless ion Manipulations (SLIM) has been studied to better understand aspects of ‘in-SLIM’ ion accumulation, and particularly its use in conjunction with ion mobility spectrometry (IMS). A linear SLIM ion path was implemented that had a ‘gate’ for blocking and accumulating ions for arbitrary time periods. Removing the gate potential caused ions to exit and the spatial distributions of accumulated ions were examined. The ion populations for a set of peptides were observed to increase with increased accumulation times until space change effects became significant, after which the peptide precursor ion populations decreased due to growing space charge-related ion activation, reactions, and losses. Ions closest to the gate were preferentially activated, and increasingly so with added storage times, TW amplitude, and the size of the total ion population. For analytical applications, we found that relatively weak TWs in the accumulation/storage region should be utilized to prevent or minimize ion losses and/or ion heating effects that can also lead to fragmentation. Our results further indicate that an accumulation region close to the SLIM entrance is advantageous for both speeding accumulation, minimizing ion heating, and avoiding ion population profiles that result in IMS peak tailing. Importantly, space charge-driven separations were observed for large populations of accumulated species and attributed to the opposing effects of space charge and the TW. In these separations, ion species form distributions or peaks, sometimes moving against the TW, and ordered in the SLIM based on their mobilities. Such separations may offer utility for ion pre-fractionation of ions and increasing the dynamic range measurements, increasing the resolving power of IMS separations by decreasing peak widths for accumulated ion populations, and other purposes benefitting from separations of extremely large ion populations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Diffusion characteristics of water molecules in a lamellar structure formed by triblock copolymers

The distribution and diffusion of water molecules are playing important roles in determining self-assembly and transport properties of polymeric systems. Small-angle neutron scattering (SANS) experiments and molecular dynamics (MD) simulation have been applied to understand the distribution of water molecules and their dynamics in the lamellar membrane formed by Pluronic L62 block copolymers. Penetration of water molecules into the polyethylene oxide (PEO) layers of the membranes has been estimated using scattering length density (SLD) profiles obtained from SANS measurements, which agree well with the molecular distribution observed from MD simulations. The water diffusion coefficient at different regions of the lamellar membrane was further investigated using MD simulation. The diffusion characteristic shows a transition from normal to anomalous diffusion as the position of the water molecule changes from the bulk to PEO and to the polypropylene oxide (PPO) layer. Here we find that water molecules within the PEO or PPO layers follow subdiffusive dynamics, which can be interpreted by the model of fractional Brownian motion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neuromorphic learning with Mott insulator NiO

Habituation and sensitization (nonassociative learning) are among the most fundamental forms of learning and memory behavior present in organisms that enable adaptation and learning in dynamic environments. Emulating such features of intelligence found in nature in the solid state can serve as inspiration for algorithmic simulations in artificial neural networks and potential use in neuromorphic computing. In this work, we demonstrate nonassociative learning with a prototypical Mott insulator, nickel oxide (NiO), under a variety of external stimuli at and above room temperature. Similar to biological species such as Aplysia, habituation and sensitization of NiO possess time-dependent plasticity relying on both strength and time interval between stimuli. A combination of experimental approaches and first-principles calculations reveals that such learning behavior of NiO results from dynamic modulation of its defect and electronic structure. An artificial neural network model inspired by such nonassociative learning is simulated to show advantages for an unsupervised clustering task in accuracy and reducing catastrophic interference, which could help mitigate the stability–plasticity dilemma. Mott insulators can therefore serve as building blocks to examine learning behavior noted in biology and inspire new learning algorithms for artificial intelligence.

42 ENGINEERING↗

Discovery and characterization of a neutralizing pan-ELR+CXC chemokine monoclonal antibody

CXCR1 and CXCR2 signaling play a critical role in neutrophil migration, angiogenesis, and tumorigenesis and are therefore an attractive signaling axis to target in a variety of indications. In human, a total of seven chemokines signal through these receptors and comprise the ELR + CXC chemokine family, so named because of the conserved ELRCXC N-terminal motif. To fully antagonize CXCR1 and CXCR2 signaling, an effective therapeutic should block either both receptors or all seven ligands, yet neither approach has been fully realized clinically. In this work, we describe the generation and characterization of LY3041658, a humanized monoclonal antibody that binds and neutralizes all seven human and cynomolgus monkey ELR + CXC chemokines and three of five mouse and rat ELR+CXC chemokines with high affinity. LY3041658 is able to block ELR + CXC chemokine-induced Ca 2+ mobilization, CXCR2 internalization, and chemotaxis in vitro as well as neutrophil mobilization in vivo without affecting other neutrophil functions. In addition to the in vitro and in vivo activity, we characterized the epitope and structural basis for binding in detail through alanine scanning, crystallography, and mutagenesis. Together, these data provide a robust preclinical characterization of LY3041658 for which the efficacy and safety is being evaluated in human clinical trials for neutrophilic skin diseases.

59 BASIC BIOLOGICAL SCIENCES↗

Mutagenesis of Human Cytomegalovirus Glycoprotein L Disproportionately Disrupts gH/gL/gO over gH/gL/pUL128-131

Herein, cell-free and cell-to-cell spread of herpesviruses involves a core fusion apparatus comprised of the fusion protein glycoprotein B (gB) and the regulatory factor gH/gL. The human cytomegalovirus (HCMV) gH/gL/gO and gH/gL/pUL128-131 facilitate spread in different cell types. The gO and pUL128-131 components bind distinct receptors, but how the gH/gL portions of the complexes functionally compare is not understood. We previously characterized a panel of gL mutants by transient expression and showed that many were impaired for gH/gL-gB-dependent cell-cell fusion but were still able to form gH/gL/pUL128-131 and induce receptor interference. Here, the gL mutants were engineered into the HCMV BAC clones TB40/e-BAC4 (TB), TR, and Merlin (ME), which differ in their utilization of the two complexes for entry and spread. Several of the gL mutations disproportionately impacted gH/gL/gO-dependent entry and spread over gH/gL/pUL128-131 processes. The effects of some mutants could be explained by impaired gH/gL/gO assembly, but other mutants impacted gH/gL/gO function. Soluble gH/gL/gO containing the L201 mutant failed to block HCMV infection despite unimpaired binding to PDGFRα, indicating the existence of other important gH/gL/gO receptors. Another mutant (L139) enhanced the gH/gL/gO-dependent cell-free spread of TR, suggesting a “hyperactive” gH/gL/gO. Recently published crystallography and cryo-electron microscopy studies suggest structural conservation of the gH/gL underlying gH/gL/gO and gH/gL/pUL128-131. However, our data suggest important differences in the gH/gL of the two complexes and support a model in which gH/gL/gO can provide an activation signal for gB.

59 BASIC BIOLOGICAL SCIENCES↗

Label-free affinity screening, design and synthesis of inhibitors targeting the Mycobacterium tuberculosis L-alanine dehydrogenase

The ability of Mycobacterium tuberculosis (Mtb) to persist in its host may enable an evolutionary advantage for drug resistant variants to emerge. A potential strategy to prevent persistence and gain drug efficacy is to directly target the activity of enzymes that are crucial for persistence. We present a method for expedited discovery and structure-based design of lead compounds by targeting the hypoxia-associated enzyme L-alanine dehydrogenase (AlaDH). Biochemical and structural analyses of AlaDH confirmed binding of nucleoside derivatives and showed a site adjacent to the nucleoside binding pocket that can confer specificity to putative inhibitors. Using a combination of dye-ligand affinity chromatography, enzyme kinetics and protein crystallographic studies, we show the development and validation of drug prototypes. Crystal structures of AlaDH-inhibitor complexes with variations at the N6 position of the adenyl-moiety of the inhibitor provide insight into the molecular basis for the specificity of these compounds. We describe a drug-designing pipeline that aims to block Mtb to proliferate upon re-oxygenation by specifically blocking NAD accessibility to AlaDH. The collective approach to drug discovery was further evaluated through in silico analyses providing additional insight into an efficient drug development strategy that can be further assessed with the incorporation of in vivo studies.

60 APPLIED LIFE SCIENCES↗

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗