Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

Aerosol Models from the AERONET Database: Application to Surface Reflectance Validation Satellite Data

Aerosols play a critical role in radiative transfer within the atmosphere, and they have a significant impact on climate change. In this paper, we propose and implement a framework for developing an aerosol model using their microphysical properties. Such microphysical properties as the size distribution, the complex refractive index, and the percentage of sphericity are derived from the global AERosol RObotic NETwork (AERONET). These measurements, however, are typically retrieved when almucantar measurement procedures are performed (i.e., early mornings and late afternoons with clear sky) and might not have a temporal correspondence to a satellite overpass time, so a valid validation of satellite-derived products cannot be carried out. To address this problem of temporal inconsistency of satellite and ground-based measurements, we developed an approach to retrieve these microphysical properties (and the corresponding aerosol model) using the optical thickness at 440 nm, τ440, and the Ångström coefficient between 440 and 870 nm, α440–870. Such aerosol models were developed for 851 AERONET sites within the last 28 years. Obtained results suggest that empirically microphysical properties can be retrieved with uncertainties of up to 23 %. An exception is the imaginary part of the refractive index ni, for which the derived uncertainties reach up to 38 %. These specific parametric models of aerosol can be used for the studies when retrieval of microphysical properties is required as well as validation of satellite-derived products over land. Specifically, we demonstrate the usefulness of the aerosol models to validate surface reflectance records over land derived from optical remote sensing sensors. We then quantify the propagation of uncertainties in the surface reflectance due to uncertainties with the aerosol model retrieval that is used as a reference from radiative transfer simulations. Results indicate that individual aerosol microphysical properties can impact uncertainties in surface reflectance retrievals between 3.5 × 10(exp -5) to 1 × 10(exp -3) (in reflectance units). The overall impact of microphysical properties combined yields an overall uncertainty in surface reflectance <0.004 (in reflectance units). That corresponds, for example, to 1 to 3 % of the retrieved surface reflectance in the red spectral band (620–670 nm) by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument. These uncertainty values are well below the specification (0.005 + 0.05ρ; ρ is the retrieved surface reflectance) used for the MODIS atmospheric correction.

Jean-Claude Roger↗

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution (PMAD), electric propulsion subsystem (EPS), and a primary heat rejection system. Specific mass, or αe (kg / kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass estimates for individual components. To inform technology maturation planning, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium (HALEU) reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively cooled heat pipe reactor concepts. Each concept requires specific geometries and working fluids to reach the performance goals of PCS interface conditions (temperature, pressure, flow rate) and system mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support technology down-selection along with related technology development planning. The reactor and shield αe are a function of several PCS design choices, and reactor scaling with these parameters must be considered to enable an informed decision on reactor geometry and working fluid combination.

Nuclear Electric Propulsion↗

Optical depth to reionization from perturbative 21 cm clustering

The optical depth τ is the least well determined parameter in the standard model of cosmology, and one whose precise value is important for both understanding reionization and for inferring fundamental physics from cosmological measurements. Here, we forecast how well future epoch of reionization experiments could constraint τ using a symmetries-based bias expansion that highlights the special role played by anisotropies in the power spectrum on large scales. Given a parametric model for the ionization evolution inspired by the physical behavior of more detailed reionization simulations, we find that future 21 cm experiments could place tight constraints on the timing and duration of reionization and hence constraints on τ that are competitive with proposed, space-based CMB missions provided they can measure k ≈ 0.1 h Mpc -1 with a clean foreground wedge across redshifts spanning the most active periods of reionization, corresponding to ionization fractions 0.2 ≲ x ≲ 0.8. Significantly improving upon existing CMB-based measurements with next-generation 21 cm surveys would require substantially longer observations (~ 5 years) than standard $\mathscr{O}$(1000 hour) integration times. Precise measurements of smaller scales will not improve constraints on τ until a better understanding of the astrophysics of reionization is achieved. In the presence of noise and foregrounds even future 21 cm experiments will struggle to constrain τ if the ionization evolution deviates significantly from simple parametric forms.

79 ASTRONOMY AND ASTROPHYSICS↗

The Fast Evolution of SN 2010bh Associated with XRF 100316D

The first observational evidence of a connection between supernovae (SNe) and gamma-ray bursts (GRBs) was found about a decade ago. Since then, only half a dozen spectroscopically confirmed associations have been discovered and XRF 1003160 associated with the type-Ic SN 20lObh is among the latest. Aims. We constrain the progenitor radius, the host-galaxy extinction, and the physical parameters of the explosion of XRF l00316D/SN 20lObh at z 0.059. We study the SN brightness and colours in the context of GRB-SNe. Methods. We began observations with the Gamma-Ray burst Optical and Near-infrared Detector (GROND) 12 hours after the GRB trigger and continued until 80 days after the burst. GROND provided excellent photometric data in six filter bands covering a wavelength range from approximately 350 to 1800 nm, significantly expanding the pre-existing data set for this event. Combining GROND and Swift data, the early broad-band spectral energy distribution (SED) is modelled with a blackbody and afterglow component attenuated by dust and gas absorption. The temperature and radius evolution of the thermal component are analysed and combined with earlier measurements available from the literature. Templates of SN 1998bw are fitted to the SN itself to directly compare the lightcurve properties. Finally, a two-component parametrized model is fitted to the quasi-bolometric light curve. which delivers physical parameters of the explosion. Results. The best-fit models to the broad-band SEDs imply moderate reddening along the line of sight through the host galaxy (A(sub v.host = 1.2 +/- 0.1 mag). Furthermore, the parameters of the blackbody component reveal a cooling envelope at an apparent initial radius of 7 x 10(exp 11) cm, which is compatible with a dense wind surrounding a Wolf-Rayet star. A multicolor comparison shows that SN 2010bh is 60 - 70% as bright as SN 1998bw. Reaching maximum brightness at 8 - 9 days after the burst in the blue bands, SN 20lObh proves to be the most rapidly evolving GRB-SNe to date. Modelling of the quasi-bolometric light curve yields M(sub Ni) = 0.21 +/- 0.03 solar M and M(sub ej) = 2.6 +/- 0.2 solar M, typical of values within the GRB-SN population. The kinetic energy is E(sub k) = (2.4 +/- 0.7) x 10(exp 52) erg, which is making this SN the second most energetic GRB-SN after SN 1998bw. Conclusions. This supernova has one of the earliest peaks ever recorded and thereafter fades more rapidly than other GRB-SNe, hypernovae, or typical type-Ic SNe. This could be explained by a thin envelope expanding at very high velocities, which is therefore unable to retain the gamma-rays that would prolong the duration of the SN event.

Olivares E., F.↗

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Estimating the Uncertainty in Passive-Microwave Rain-Retrievals

Current passive-microwave rain-retrieval methods are largely based on databases built off-line using cloud models. The vertical distribution of hydrometeors within the cloud has a large impact on upwelling brightness temperatures ([31,[5]). Thus, a forward radiative transfer model can predict off-line the radiance associated with different rain scenarios. To estimate the rain from measured brightness temperatures, one simply looks for the rain scenario whose associated radiances are closest to the measurements. To understand the uncertainties in this process, we first study the dependence of the simulated brightness temperatures on different hydrometeor size distribution (DSD) models. We then analyze the marginal and joint distributions of the radiances observed by the Tropical Rainfall Measuring Mission satellite and of those in the databases used in the TRMM rain retrievals. We finally calculate the covariances of the rain profiles and brightness temperatures in the TRMM passive-microwave database and derive a simple parametric model for the conditional uncertainty, given measured radiances. These results are used to characterize the uncertainty inherent in the passive-microwave retrieval.

Coppens, Dorothee↗

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides for development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the "open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media "Darcy" relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations, Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the 'open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media 'Darcy' relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations. Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

Synergy of tensile strength and high cycle fatigue properties in a novel additively manufactured Al-Ni-Ti-Zr alloy with a heterogeneous microstructure

Alloy design strategies in additive manufacturing (AM) to achieve grain refinement and terminal eutectic solidification have been introduced to engineer Al alloys having microstructural hierarchy and heterogeneity. Such alloy design strategies enable crack-free builds with an expanded AM processing window and pushed the strength limit in Al alloys. However, fatigue performance of Al alloys made by AM is restricted by the presence of process induced defects and its stochasticity. In this work, tensile and high cycle fatigue (HCF) behavior of a novel Al-Ni-Ti-Zr alloy with a heterogeneous microstructure is studied in the as-built condition, supplemented by detailed microstructural and mechanical characterization. Excellent strength-ductility synergy of 342 MPa and 16% failure strain achieved in the alloy was associated with the microstructural attributes that pertain to the novel alloy. Additionally, the alloy showed excellent HCF performance with a fatigue endurance limit to ultimate tensile strength ratio of 0.29 in flexural fatigue mode. The study revealed the existence of multiple crack retardation mechanisms and favorable crack propagation pathways through the fine-grained regions which enabled good fatigue performance to the alloy. Further, a probabilistic model has been used to estimate the fatigue life of the alloy as a function of the stochastic microstructure by utilizing the statistical distribution of pores, solid-state inclusions, and grains in the AM Al alloy. Finally, the model parametric trends are consistent with the experimental observations.

36 MATERIALS SCIENCE↗

A copula-based Bayesian method for probabilistic solar power forecasting

With increased penetration of solar energy sources, solar power forecasting has become more crucial and challenging. This paper proposes a copula-based Bayesian approach to improve probabilistic solar power forecasting by capturing the joint distribution between solar power and ambient temperature. A prior forecast distribution is first obtained using different underlying point forecasting models. Parametric and empirical copulas of solar power and temperature are then developed to update the prior distribution to the posterior forecast distribution. A public solar power database is used to demonstrate effectiveness of the proposed method. Numerical results show that the copula-based Bayesian method outperforms the forecasting method that directly uses temperature as a feature. The Bayesian method is also compared with persistent models and show improved performance. Furthermore, this article includes supplementary material (data and code) for reproducibility.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modeling personnel turnover in the parametric organization

A model is developed for simulating the dynamics of a newly formed organization, credible during all phases of organizational development. The model development process is broken down into the activities of determining the tasks required for parametric cost analysis (PCA), determining the skills required for each PCA task, determining the skills available in the applicant marketplace, determining the structure of the model, implementing the model, and testing it. The model, parameterized by the likelihood of job function transition, has demonstrated by the capability to represent the transition of personnel across functional boundaries within a parametric organization using a linear dynamical system, and the ability to predict required staffing profiles to meet functional needs at the desired time. The model can be extended by revisions of the state and transition structure to provide refinements in functional definition for the parametric and extended organization.

Dean, Edwin B.↗

DSO+T: Valuation Methodology and Economic Metrics (DSO+T Study: Volume 4)

This report summarizes a rigorous valuation analysis methodology used by the Distribution System Operator with Transactive (DSO+T) study to estimate the financial benefits and costs of adopting Transactive Energy coordination of distributed energy resources for key stakeholders (for example distribution system operators and customers). This was achieved by modeling the value exchanges between stakeholders and determining the annualized costs and revenues experienced by stakeholders, enabling the evaluation of overall impact on stakeholder’s annualized cash flow. Extensive work was conducted developing methods to estimate the operating costs of distribution system operators at a level of granularity that would allow the financial impact of implementing a transactive energy approach to be estimated. This work included developing parametric models for labor and software costs, distribution system capital and maintenance costs, growth rates, and factors to determine annualized costs of capital investments. Simulation results were used to calculate wholesale energy costs and revenues from retail sales. Valuation analysis methods were also developed for other stakeholders including customers, the transmission system operator, independent system operator, and generators. The resulting capability allows a complete mapping of the flow of financial value between stakeholders that can be directly integrated with the results of demand flexibility simulations. Example results are provided for a business-as-usual case and compared to a transactive energy case as well as to actual cost data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Practical Aspects of the Frequency Domain Approach for Aircraft System Identification

Practical aspects of the frequency-domain approach for aircraft system identification are explained and demonstrated. Topics related to experiment design, flight data analysis, and dynamic modeling are included. For demonstration purposes, simulated time series data and simulated flight data from an F-16 nonlinear simulation with realistic noise are used. This approach enables detailed evaluations of the techniques and results, because the true characteristics of the data and aircraft dynamics are known for the simulated data. Analytical techniques and practical considerations are examined for the finite Fourier transform, nonparametric frequency response estimation, parametric modeling in the frequency domain, experiment design for frequency-domain modeling, data analysis and modeling in the frequency domain, and real-time calculations. Flight data from a subscale jet transport aircraft are used to demonstrate some of the techniques and technical issues.

Morelli, Eugene A.↗

Pairing properties of the $t$–$t$'–$t$"–$J$ model

We study the pairing properties of the two-dimensional t–t'–t"–J model, where t' and t" are second and third neighbor hoppings, at a doping level x ≈ 0.1. Recent studies of the t–t'–J model find strong pairing for t' > 0, associated with electron doping, but an absence of pairing for t'<0, associated with hole doping. This is in contrast to the cuprates, where the highest transition temperatures appear for hole doping. Model parametrizations for the cuprates estimate a t" comparable to t', which, in principle, might fix this discrepancy. However, we find that it does not; we observe a suppression of pairing for the hole-doped system (t' < 0,t" > 0), while for the electron-doped system (t' > 0,t" < 0) d-wave pairing is robust. As a result, extended hoppings appear to be insufficient to make the one-band t–t'–J model capable of describing the pairing in the hole-doped system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗

Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks

Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.

54 ENVIRONMENTAL SCIENCES↗

Modeling Fire Severity in Black Spruce Stands in the Alaskan Boreal Forest Using Spectral and Non-Spectral Geospatial Data

Biomass burning in the Alaskan interior is already a major disturbance and source of carbon emissions, and is likely to increase in response to the warming and drying predicted for the future climate. In addition to quantifying changes to the spatial and temporal patterns of burned areas, observing variations in severity is the key to studying the impact of changes to the fire regime on carbon cycling, energy budgets, and post-fire succession. Remote sensing indices of fire severity have not consistently been well-correlated with in situ observations of important severity characteristics in Alaskan black spruce stands, including depth of burning of the surface organic layer. The incorporation of ancillary data such as in situ observations and GIS layers with spectral data from Landsat TM/ETM+ greatly improved efforts to map the reduction of the organic layer in burned black spruce stands. Using a regression tree approach, the R2 of the organic layer depth reduction models was 0.60 and 0.55 (pb0.01) for relative and absolute depth reduction, respectively. All of the independent variables used by the regression tree to estimate burn depth can be obtained independently of field observations. Implementation of a gradient boosting algorithm improved the R2 to 0.80 and 0.79 (pb0.01) for absolute and relative organic layer depth reduction, respectively. Independent variables used in the regression tree model of burn depth included topographic position, remote sensing indices related to soil and vegetation characteristics, timing of the fire event, and meteorological data. Post-fire organic layer depth characteristics are determined for a large (N200,000 ha) fire to identify areas that are potentially vulnerable to a shift in post-fire succession. This application showed that 12% of this fire event experienced fire severe enough to support a change in post-fire succession. We conclude that non-parametric models and ancillary data are useful in the modeling of the surface organic layer fire depth. Because quantitative differences in post-fire surface characteristics do not directly influence spectral properties, these modeling techniques provide better information than the use of remote sensing data alone.

Barrett, K.↗