Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Refining”

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 685 records · Page 38

A novel transfer learning framework for sorghum biomass prediction using UAV-based remote sensing data and genetic markers

Yield for biofuel crops is measured in terms of biomass, so measurements throughout the growing season are crucial in breeding programs, yet traditionally time- and labor-consuming since they involve destructive sampling. Modern remote sensing platforms, such as unmanned aerial vehicles (UAVs), can carry multiple sensors and collect numerous phenotypic traits with efficient, non-invasive field surveys. However, modeling the complex relationships between the observed phenotypic traits and biomass remains a challenging task, as the ground reference data are very limited for each genotype in the breeding experiment. In this study, a Long Short-Term Memory (LSTM) based Recurrent Neural Network (RNN) model is proposed for sorghum biomass prediction. The architecture is designed to exploit the time series remote sensing and weather data, as well as static genotypic information. As a large number of features have been derived from the remote sensing data, feature importance analysis is conducted to identify and remove redundant features. A strategy to extract representative information from high-dimensional genetic markers is proposed. To enhance generalization and minimize the need for ground reference data, transfer learning strategies are proposed for selecting the most informative training samples from the target domain. Consequently, a pre-trained model can be refined with limited training samples. Field experiments were conducted over a sorghum breeding trial planted in multiple years with more than 600 testcross hybrids. The results show that the proposed LSTM-based RNN model can achieve high accuracies for single year prediction. Further, with the proposed transfer learning strategies, a pre-trained model can be refined with limited training samples from the target domain and predict biomass with an accuracy comparable to that from a trained-from-scratch model for both multiple experiments within a given year and across multiple years.

36 MATERIALS SCIENCE↗

Complex Disorder in Type-I Clathrates: Synthesis and Structural Characterization of A8GaxSn46−x (A = Rb, Cs; 6.9 < x < 7.5)

Exploratory studies in the systems Rb–Ga–Sn and Cs–Ga–Sn yielded the cubic type-I clathrates with refined compositions Rb8GaxSn46−x and Cs8GaxSn46−x (6.9 < x < 7.5). Nearly single-phase materials with good crystallinity were obtained from stoichiometric reactions of the elements. The structures were characterized by means of single-crystal X-ray diffraction methods. Both Rb8GaxSn46−x and Cs8GaxSn46−x represents cases, where a Group 13 element randomly substitutes a Group 14 element in the structure. The extent of Ga/Sn mixing is apparently governed by the drive of the system to achieve an optimal valence electron count, and hence, Rb8GaxSn46−x and Cs8GaxSn46−x (x ≈ 8) can be regarded as Zintl phases. This notion is supported by structure refinements on a multitude of single-crystal X-ray diffraction data, which also confirm that both types of cages in the cubic type-I structure are fully occupied by Rb and Cs atoms. The open-framework, comprised of 46 nodes per formula unit, adapts to the incorporation of nearly eight Ga atoms within the matrix of Sn, whereby small, short-range distortions result. The exact nature of these effects is still unclear, as so far, the structural variations could only be modeled as both positional and occupational disorder at one of three framework sites. Since vacancies in the structures of the binary type-I clathrates A8Sn46−x☐x (A = Rb, Cs; ☐ = missing Sn atom) are also known to cause local distortions, the latter were also synthesized with the same protocols used for the synthesis of A8GaxSn46−x and structurally re-analyzed. The results from the latter studies confirm that homogeneity issues abound, and that the final structures/compositions are an intricate function of the experimental conditions.

36 MATERIALS SCIENCE↗

Mechanical Behavior of Al–Al2Cu–Si and Al–Al2Cu Eutectic Alloys

In this study, laser rapid solidification technique was used to refine the microstructure of ternary Al–Cu–Si and binary Al–Cu eutectic alloys to nanoscales. Micropillar compression testing was performed to measure the stress–strain response of the samples with characteristic microstructure in the melt pool regions. The laser-remelted Al–Al2Cu–Si ternary alloy was observed to reach the compressive strength of 1.59 GPa before failure at a strain of 28.5%, which is significantly better than the as-cast alloy with a maximum strength of 0.48 GPa at a failure strain of 4.8%. The laser-remelted Al–Cu binary alloy was observed to reach the compressive strength of 2.07 GPa before failure at a strain of 26.5%, which is significantly better than the as-cast alloy with maximum strength of 0.74 GPa at a failure strain of 3.3%. The enhanced compressive strength and improved compressive plasticity were interpreted in terms of microstructural refinement and hierarchical eutectic morphology.

36 MATERIALS SCIENCE↗

Melt-Pool Dynamics and Microstructure of Mg Alloy WE43 under Laser Powder Bed Fusion Additive Manufacturing Conditions

Magnesium-based alloy WE43 is a state-of-the-art bioresorbable metallic implant material. There is a need for implants with both complex geometries to match the mechanical properties of bone and refined microstructure for controlled resorption. Additive manufacturing (AM) using laser powder bed fusion (LPBF) presents a viable fabrication method for implant applications, as it offers near-net-shape geometrical control, allows for geometry customization based on an individual patient, and fast cooling rates to achieve a refined microstructure. In this study, the laser–alloy interaction is investigated over a range of LPBF-relevant processing conditions to reveal melt-pool dynamics, pore formation, and the microstructure of laser-melted WE43. In situ X-ray imaging reveals distinct laser-induced vapor depression morphology regimes, with minimal pore formation at laser-scan speeds greater than 500 mm/s. Optical and electron microscopy of cross-sectioned laser tracks reveal three distinct microstructural regimes that can be controlled by adjusting laser-scan parameters: columnar, dendritic, and banded microstructures. These regimes are consistent with those predicted by the analytic solidification theory for conduction-mode welding, but not for keyhole-mode tracks. The results provide insight into the fundamental laser–material interactions of the WE43 alloy under AM-processing conditions and are critical for the successful implementation of LPBF-produced WE43 parts in biomedical applications.

36 MATERIALS SCIENCE↗

Transient Microstructure Evolutions and Local Properties of Dual-Phase 980 MPa Grade Steel Via Friction Stir Spot Processing

Friction stir processing is a novel solid-state process to modify microstructures and their properties by intense, localized plastic deformation. However, little research has been reported for microstructure evolutions of advanced high-strength steels during the process. The present work focuses on the study of transient microstructure changes and local mechanical properties for friction stir spot processed dual-phase (DP) 980 MPa grade steel (DP980) under different peak temperatures. A pinless silicon nitride ceramic tool was used to produce relatively simple material deformation and flow near the tool. Friction stir spot processed steel samples were characterized by optical and electron microscopies. Furthermore, Vickers microhardness and nano-indentation measurements were used to study local mechanical properties for correlation with microstructures. A swallow layer of refined grains (<0.6 µm) was obtained with a low peak temperature (under 400 °C), whereas higher peak temperatures (>Ac 1 ) led to a change in grain size with different microstructures (fine-grained DP or martensite). Electron back-scattered diffraction characterizations revealed a large deformation in the as-received microstructures (mixture of ferrite and tempered martensite) induced by friction stir spot processing, leading to recrystallization and grain refinement around the stirred zone. Also, nano-indentation measurements showed a higher hardness than the hardness of the as-received DP980. Friction stir processing with different process conditions effectively changed microstructures and local mechanical properties.

36 MATERIALS SCIENCE↗

Structure Analysis of Natural Wangdaodeite—LiNbO3-Type FeTiO3

This paper reports the first structure refinement of natural wangdaodeite, LiNbO3-type FeTiO3 from the Ries impact structure. Wangdaodeite occurs together with recrystallized ilmenite clasts in shock melt veins which have experienced peak shock pressures of between 17 and 22 GPa. Comparison of natural and synthetic wangdaodeite points toward a correlation between the distortion of ferrate- and titanate-polyhedra and the c/a ratio of the unit cell. The Raman spectrum of wangdaodeite is calculated based on the refined structure. Comparison to the reported spectrum of the type-material shows that the Raman peak at 738–740 cm−1 is indicative for this phase, whereas other features in type-wangdaodeite are tentatively assigned to disordered ilmenite.

58 GEOSCIENCES↗

Fabricating Graphene Oxide/h-BN Metal Insulator Semiconductor Diodes by Nanosecond Laser Irradiation

To employ graphene’s rapid conduction in 2D devices, a heterostructure with a broad bandgap dielectric that is free of traps is required. Within this paradigm, h-BN is a good candidate because of its graphene-like structure and ultrawide bandgap. We show how to make such a heterostructure by irradiating alternating layers of a-C and a-BN film with a nanosecond excimer laser, melting and zone-refining constituent layers in the process. With Raman spectroscopy and ToF-SIMS analyses, we demonstrate this localized zone-refining into phase-pure h-BN and rGO films with distinct Raman vibrational modes and SIMS profile flattening after laser irradiation. Furthermore, in comparing laser-irradiated rGO-Si MS and rGO/h-BN/Si MIS diodes, the MIS diodes exhibit an increased turn-on voltage (4.4 V) and low leakage current. The MIS diode I-V characteristics reveal direct tunneling conduction under low bias and Fowler-Nordheim tunneling in the high-voltage regime, turning the MIS diode ON with improved rectification and current flow. This study sheds light on the nonequilibrium approaches to engineering h-BN and graphene heterostructures for ultrathin field effect transistor device development.

2D materials↗

Spray Drying Enzyme-Treated Cellulose Nanofibrils

Enzyme-treated cellulose nanofibrils (CNFs) were produced via a lab-scale mass colloider using bleached kraft pulp (BKP) to evaluate their processability and power requirements during refining and spray-drying operations. To evaluate the energy efficiency in the CNF refining process, the net energy consumption, degree of polymerization (DP), and viscosity were determined. Less energy was consumed to attain a given fines level by using the endoglucanase enzymes. The DP and viscosity were also decreased using the enzymes. The morphological properties of the enzyme-pretreated spray-dried CNF powders (SDCNFs) were measured. Subsequently, the enzyme-pretreated SDCNFs were added to a PP matrix with MAPP as a coupling agent. The mixture was then compounded through a co-rotating twin-screw extruder to determine whether the enzyme treatment of the CNFs affects the mechanical properties of the composites. Compared to earlier studies on enhancing PMCs with SDCNF powders, this research investigates the use of enzyme-pretreated SDCNF powders. It was confirmed that the strength properties of PP increased by adding SDCNFs, and the strength properties were maintained after adding enzyme-pretreated SDCNFs.

59 BASIC BIOLOGICAL SCIENCES↗

A Scalable Approach to Minimize Charging Costs for Electric Bus Fleets

Incorporating battery electric buses into bus fleets faces three primary challenges: a BEB’s extended refuel time, the cost of charging, both by the consumer and the power provider, and large compute demands for planning methods. When BEBs charge, the additional demands on the grid may exceed hardware limitations, so power providers divide a consumer’s energy needs into separate meters even though doing so is expensive for both power providers and consumers. Prior work has developed a number of strategies for computing charge schedules for bus fleets; however, prior work has not worked to reduce costs by aggregating meters. Additionally, because many works use mixed integer linear programs, their compute needs make planning for commercial-sized bus fleets intractable. This work presents a multi-program approach to computing charge plans for electric bus fleets. The proposed method solves a series of subproblems where the solution to the charge problem becomes more refined with each problem, moving closer to the optimal schedule. The results demonstrate how runtimes are reduced by using intermediate subproblems to refine the bus charge solution so that the proposed method can be applied to large bus fleets of 100+ buses. Not only will we demonstrate that runtimes scale linearly with the number of buses but we will also show how the proposed method scales to large bus fleets of over 100 buses while managing the monthly cost of energy.

Mortensen, Daniel (ORCID:0000000276494452)↗

Trebiskyite, the First Titanium-Decavanadate Mineral

Abstract Trebiskyite, Na3Mg2[TiV9O28]·22H2O, from the Pickett Corral mine, Bull Canyon, Montrose County, Colorado, USA, is the first mineral containing essential Ti as a component of a decavanadate polyanion. Crystals are thin prisms or needles up to ∼0.3 mm in length, occurring as individuals and in divergent sprays. Trebiskyite crystals are translucent yellow with vitreous luster and pale-yellow streak and are non-fluorescent under both longwave and shortwave ultraviolet illumination. The crystals exhibit elongation on [001] and are modified by several 0kl prism forms, possibly corresponding to {010}, {011}, and {012}. The Mohs hardness is ca. 2 based on scratch tests. Crystals have brittle tenacity with irregular fracture. No apparent cleavage, nor twinning, was evident, and the mineral is readily soluble in H2O. The density measured by flotation in methylene iodide-toluene mixture is 2.38(2) g/cm3. The calculated density is 2.352 g/cm3 for the empirical formula and 2.343 g/cm3 for the ideal formula. Electron microprobe analyses provided the empirical formula Na3.39Mg1.87V8.75Ti1.25O50H44.00, based on Ti + V = 10 apfu, O = 50 apfu, with the Na content having been determined from the refined crystal structure. Trebiskyite is monoclinic, space group P21/c, a = 9.478(4), b = 21.426(11), c = 11.267(5) Å, β = 114.572(7)°, V = 2080.8(13) Å3, Z = 2. The structure was refined from 6642 reflections to a final R1 = 0.0517 for reflections I > 2σ(I). The structural unit of trebiskyite, (TiV9O28)7−, is built from a Ti-substituted decavanadate polyanion, (V10O28)6−, found in other members of the pascoite family. Chemically, trebiskyite is closely related to lasalite, Na2Mg2[V10O28]·20H2O, and huemulite, Na4Mg[V10O28]·24H2O; however, trebiskyite contains a novel structural arrangement that is unique among the decavanadate minerals. The interstitial complex in trebiskyite has the ideal formula {[Na3Mg2(H2O)19](H2O)3} and consists of chains of distorted octahedrally coordinated Na+ and regular Mg[O(H2O)5] octahedra. The observed low incident bond valence to the V1 site (4.45 vu) and marked deviation in vanadyl and V–Oeq bond lengths indicate that this site is host to the essential Ti4+ identified during chemical analyses.

Mineralogy↗

Fleetite, Cu 2 RhIrSb 2 , a New Species of Platinum-Group Mineral from the Miass Placer Zone, Southern Urals, Russia

Fleetite, Cu 2 RhIrSb 2 , a new species of platinum-group mineral (PGM), was discovered intergrown with an Os–Ir–Ru alloy in the Miass Placer Zone (Au–PGE), southern Urals, Russia. A single grain 50 μm across was found. Osmium, ruthenium, and iridium are the main associated minerals; also present are Pt–Fe alloys, laurite, Sb-rich irarsite, Rh-rich tolovkite, kashinite, anduoite, ferronickelplatinum, heazlewoodite, PGE-bearing pentlandite and digenite, as well as micrometric inclusions of forsterite (Fo 93.7 ), chromite–magnesiochromite, and Mg-rich edenite. In reflected light, fleetite is light gray; it is opaque, isotropic, non-pleochroic, and non-bireflectant. We report reflectance values measured in air. A mean of seven point-analyses (wavelength-dispersive spectrometry) gave Cu 13.93, Ni 8.60, Fe 0.10, Ir 28.07, Rh 7.91, Ru 1.96, Sb 39.28, total 99.85 wt.%, corresponding to (Cu 1.41 Ni 0.58 Fe 0.01 ) Σ2.00 (Rh 0.49 Ni 0.36 Ru 0.12 ) Σ0.97 Ir 0.95 Sb 2.08 on the basis of six atoms per formula unit, taking into account the structural results. The calculated density is 10.83 g/cm 3 . Single-crystal X-ray studies show that fleetite is cubic, space group Fdm (#227), a = 11.6682(8) Å, V = 1588.59(19) Å 3 , and Z = 16. A least-squares refinement of X-ray powder-diffraction data gave a = 11.6575(5) Å and V = 1584.22(19) Å 3 . The strongest five reflections in the powder pattern [d in Å(I)(hkl)] are: 6.70(75)(111), 4.13(100)(220), 3.52(30)(311), 2.380(50)(422), 2.064(40)(440). Results of synchrotron micro-Laue diffraction experiments are consistent [a = 11.66(2) Å]. The crystal structure of fleetite was solved and refined to R = 0.0340 based upon 153 reflections with F o > 4σ(F o ). It is isotypic with Pd 11 Bi 2 Se 2 and best described as intermetallic, with all metal atoms in 12-fold coordination. Fleetite and other late exotic phases were formed by reaction of the associated alloy phases with a fluid phase enriched in Sb, As, and S in circulation in the cooling ophiolite source-rock. The mineral is named after Michael E. Fleet (1938–2017) in recognition of his significant contributions to the Earth Sciences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying the impacts of marine aerosols over the southeast Atlantic Ocean using a chemical transport model: implications for aerosol–cloud interactions

The southeast Atlantic region, characterized by persistent stratocumulus clouds, has one of the highest uncertainties in aerosol radiative forcing and significant variability across climate models. In this study, we analyze the seasonally varying role of marine aerosol sources and identify key uncertainties in aerosol composition at cloud-relevant altitudes over the southeast Atlantic using the GEOS-Chem chemical transport model. We evaluate simulated aerosol optical depth (AOD) and speciated aerosol concentrations against those collected from ground observations and aircraft campaigns such as LASIC, ORACLES, and CLARIFY, conducted during 2017. The model consistently underestimates AOD relative to AERONET, particularly at remote locations like Ascension Island. However, when compared with aerosol mass concentrations from aircraft campaigns during the biomass burning period, it performs adequately at cloud-relevant altitudes, with a normalized mean bias (NMB) between -3.5 % (CLARIFY) and -7.5 % (ORACLES). At these altitudes, in the model, organic aerosols (63 %) dominate during the biomass burning period, while sulfate (41 %) prevails during austral summer, when dimethylsulfide (DMS) emissions peak in the model. Our findings indicate that marine sulfate can account for up to 69 % of total sulfate during the high-DMS period. Sensitivity analyses indicate that refining DMS emissions and oxidation chemistry may increase sulfate aerosol produced from marine sources, highlighting that there remains large uncertainty as to the role of DMS emissions in the marine boundary layer. Additionally, we find marine primary organic aerosol emissions may substantially increase total organic aerosol concentrations, particularly during austral summer. This study underscores the imperative need to refine marine emissions and their chemical transformations, as aerosols from marine sources are a major component of total aerosols at cloud-relevant altitudes and may impact uncertainties in aerosol radiative forcing over the southeast Atlantic.

54 ENVIRONMENTAL SCIENCES↗

Version 4 CALIPSO Imaging Infrared Radiometer ice and liquid water cloud microphysical properties – Part I: The retrieval algorithms

Following the release of the version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version (version 4; V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. The IIR Level 2 algorithm has been modified in the V4 data release to improve the accuracy of the retrievals in clouds of very small (close to 0) and very large (close to 1) effective emissivities. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear-sky brightness temperatures over oceans has been updated and tuned for the simulations using Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data to match IIR observations in clear-sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5 km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ±0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also improved retrievals in ice clouds having large emissivity by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3 and increases the number of retrievals. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 and 10.6 µm and at 12.05 and 08.65 µm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice habit models from the recent “TAMUice2016” database, namely the single-hexagonal-column model and the eight-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 and 10.6 µm and effective diameter as derived from in situ measurements at tropical and midlatitudes during the Tropical Composition, Cloud, and Climate Coupling (TC4) and Small Particles in Cirrus Science and Operations Plan (SPARTICUS) field experiments.

54 ENVIRONMENTAL SCIENCES↗

Disturbance triggers non-linear microbe–environment feedbacks

Abstract. Conceptual frameworks linking microbial community membership, properties, and processes with the environment and emergent function have been proposed but remain untested. Here we refine and test a recent conceptual framework using hyporheic zone sediments exposed to wetting–drying transitions. Our refined framework includes relationships between cumulative properties of a microbial community (e.g., microbial membership, community assembly properties, and biogeochemical rates), environmental features (e.g., organic matter thermodynamics), and emergent ecosystem function. Our primary aim was to evaluate the hypothesized relationships that comprise the conceptual framework and contrast outcomes from the whole and putatively active bacterial and archaeal communities. Throughout the system we found threshold-like responses to the duration of desiccation. Membership of the putatively active community – but not the whole bacterial and archaeal community – responded due to enhanced deterministic selection (an emergent community property). Concurrently, the thermodynamic properties of organic matter (OM) became less favorable for oxidation (an environmental component), and respiration decreased (a microbial process). While these responses were step functions of desiccation, we found that in deterministically assembled active communities, respiration was lower and thermodynamic properties of OM were less favorable. Placing the results in context of our conceptual framework points to previously unrecognized internal feedbacks that are initiated by disturbance and mediated by thermodynamics and that cause the impacts of disturbance to be dependent on the history of disturbance.

59 BASIC BIOLOGICAL SCIENCES↗

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗