Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MODEL CALIBRATION”

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 703 records · Page 39

The impact of composition choices on solar evolution: age, helio- and asteroseismology, and neutrinos

ABSTRACT The Sun is the most studied and well-known star, and as such, solar fundamental parameters are often used to bridge gaps in the knowledge of other stars, when these are required for modelling. However, the two most powerful and precise independent methodologies currently available to infer the internal solar structure are in disagreement. We aim to show the potential impact of composition choices in the overall evolution of a star, using the Sun as example. To this effect, we create two Standard Solar Models and a comparison model using different combinations of metallicity and relative element abundances and compare evolutionary, helioseismic, and neutrino-related properties for each. We report differences in age for models calibrated to the same point on the HR diagram, in red giant branch, of more than 1 Gyr, and found that the current precision level of asteroseismic measurements is enough to differentiate these models, which would exhibit differences in period spacing of 1.30–2.58 per cent. Additionally, we show that the measurement of neutrino fluxes from the carbon–nitrogen–oxygen cycle with a precision of around 17 per cent, which could be achieved by the next generation of solar neutrino experiments, could help resolve the stellar abundance problem.

Capelo, Diogo↗

Quantification of trace iodine using laser-induced breakdown spectroscopy for real-time monitoring of nuclear off-gas streams

This study evaluated the potential of laser-induced breakdown spectroscopy (LIBS) for real-time monitoring of trace gas-phase iodine, which is an element of high significance in nuclear applications due to its long radioactive half-life (as iodine-129), volatility, and biological impact. In anticipation of iodine evolving into off-gas systems in molten salt reactor and nuclear fuel recycling applications, this research aimed to assess LIBS performance in flowing argon and helium matrices; optimize measurement parameters using a multichannel spectrometer; and perform calibrations to assess predictive capabilities and limits of detection (LODs). Experimental results successfully measured gas-phase iodine in flowing argon and helium; however, trace iodine was not detected in air. Optimal delay times were determined to be 10 µs for argon and 1 µs for helium, which are consistent with the expected shorter plasma lifetime in helium relative to argon. An emission line survey was provided with the 206.16, 804.37, 902.24, and 905.83 nm peaks, which were identified as the strongest emission peaks. Calibration models were successfully built in both helium and argon, achieving LODs down to 3 ppm in helium and 5 ppm in argon. The iodine emission at 905.83 nm emerged as the most robust for calibration and was subsequently applied to a time series dataset in argon. The predictive trace confirmed the feasibility of employing LIBS for continuous, online quantification of trace iodine in flowing gas systems.

Andrews, Hunter B. [Oak Ridge National Laboratory ↗

Automatic tuning of the Community Atmospheric Model CAM5.3 by using short-term hindcasts with an improved downhill simplex optimization method

Traditional trial-and-error tuning of uncertain parameters in global atmospheric General Circulation Models (GCM) is time consuming and subjective. This study explores the feasibility of automatic optimization of GCM parameters for fast physics by using short-term hindcasts. An automatic workflow is described and applied to the Community Atmospheric Model (CAM5) to optimize several parameters in its cloud and convective parameterizations. We show that the auto-optimization leads to 10 % reduction of the overall bias in CAM5, which is already a well calibrated model, based on a pre-defined metric that includes precipitation, temperature, humidity, and longwave/shortwave cloud forcing. The computational cost of the entire optimization procedure is about equivalent to about a single 12-year atmospheric model simulation. The tuning reduces the large underestimation in the CAM5 longwave cloud forcing by decreasing the threshold relative humidity and the sedimentation velocity of ice crystals in the cloud schemes; it reduces the overestimation of precipitation by increasing the adjustment time in the convection scheme. The physical processes behind the tuned model performance for each targeted field are discussed. Limitations of the automatic tuning are described, including the slight deterioration in some targeted fields that reflect the structural errors of the model. It is pointed out that automatic tuning can be a viable supplement to process-oriented model evaluations and improvement.

54 ENVIRONMENTAL SCIENCES↗

Characterizing particle-scale equilibrium adsorption and kinetics of uranium(VI) desorption from U-contaminated sediments

Rates of contaminant U(VI) release from individual size fractions of a composite sediment from the seasonally saturated lower vadose zone of the Hanford 300-Area were examined in flow-through batch reactors to maintain quasi-constant chemical conditions. Variability in equilibrium adsorption among the various size fractions was determined in static batch reactors and analyzed using the surface complexation modeling approach. The estimated stoichiometric coefficients of U(VI) surface complexation reactions with respect to pH and carbonate concentrations varied with size fractions. This source of variability significantly increased the uncertainty in U(VI) conditional equilibrium constants over that estimated from experimental errors alone. A minimum difference between conditional equilibrium constants was established in order to evaluate statistically significant differences between sediment adsorption properties. A set of equilibrium and kinetic expressions for cation exchange, calcite dissolution, aerobic respiration, and silica dissolution were incorporated in a reaction-rate model to describe the temporal evolution of solute concentrations observed during the flow-through batch experiments. Parameters in the reaction-rate model, calibrated using experimental data for select size fractions, predicted the changes in solute concentrations for the bulk, <2 mm, sediment sample. Kinetic U(VI) desorption was well described using a multi-rate surface complexation model with an assumed lognormal distribution for the rate constants. The estimated mean and standard deviation were the same for all < 2mm size fractions, but differed in the 2-8mm size fraction. Micropore volumes in the varied size fractions were also similar as assessed using t-plots to analyze N2 desorption data. These findings provide further support for the link between microporosity and particle-scale mass transfer rates controlling kinetic U(VI) adsorption/desorption and for the utility of N2 desorption isotherms for characterizing pore networks that influence mass transfer rates.

Stoliker, Deborah L.↗

Pu(IV) quantification via visible–near-infrared absorption spectroscopy: tackling interferences using D-optimal design and partial least squares

Here, this study presents a novel analytical approach for quantifying Pu(IV) in glove box environments using fiber-optic-based visible–near-infrared absorption spectroscopy in combination with partial least squares regression (PLSR) and design of experiments. The method addresses significant challenges posed by overlapping spectral features arising from Nd(III), which is a common fission product impurity, and the speciation variability of Pu(IV) nitrato complexes in HNO 3 concentrations ranging from 2.5 to 11 M. A curated training set consisting of data from 20 samples was developed via D-optimal design to enable robust PLSR model calibration for Pu(IV) using the near-infrared band near 1050 nm. The training set was acquired from samples in cuvettes with a 1-cm path length and was used to build the PLSR model. The robustness of the model was validated with data collected using a dip probe with a 1-cm path length and varying Pu(IV) concentrations. The strong performance of the model indicates good model transfer from cuvette to dip probe and highlights the potential for in situ measurements and online monitoring of reactions in a crystallization reactor vessel. The results demonstrate that this combined spectroscopic and chemometric approach can accurately and simultaneously quantify Pu(IV) and HNO 3 , thereby offering a promising tool for real-time monitoring in process environments.

Actinide↗

Increasing Mosquito Abundance Under Global Warming

Mosquitoes are a key virus vector that poses significant health threats globally, affecting 700 million individuals and causing 1 million deaths annually. Accurately predicting mosquito abundance and dispersion remains a challenge. Complex interactions between mosquito dynamics and various environmental factors, notably hydrology, contribute to this challenge. Existing models typically focus on precipitation and temperature and often overlook further impacts of hydrological variables within mosquito modeling. In this study, we developed an artificial intelligence‐based model for mosquito dynamics, explicitly accounting for different hydrological variables, such as precipitation, soil moisture and streamflow. Using Toronto, Canada, as a case study, we identified causal relationships between changes in mosquito populations, hydrological factors, vegetation (e.g., leaf area index), and climate variables (e.g., daylight length, precipitation, and temperature). We embedded these relationships into a Long Short‐Term Memory (LSTM) Neural Network Model capable of accurately detecting mosquito dynamics across annual, seasonal, and monthly time scales. The LSTM is able to explain, on average, approximately 40% of the variance in the observed mosquito abundance data. Using the calibrated model, we predicted that the summer season mosquito abundance would increase by ∼16% and ∼19% under an intermediate greenhouse emission scenario, Shared Socioeconomic Pathway (SSP) 2–4.5, and a high greenhouse emission scenario, SSP5‐8.5, respectively. We expect that this model can serve as a valuable tool and inform science‐based decisions affecting mosquito dynamics and public health. It can also build a foundation for future risk analysis at the regional and larger scales.

54 ENVIRONMENTAL SCIENCES↗

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗

Automatic tuning of the Community Atmospheric Model (CAM5) by using short-term hindcasts with an improved downhill simplex optimization method

Traditional trial-and-error tuning of uncertain parameters in global atmospheric general circulation models (GCMs) is time consuming and subjective. This study explores the feasibility of automatic optimization of GCM parameters for fast physics by using short-term hindcasts. An automatic workflow is described and applied to the Community Atmospheric Model (CAM5) to optimize several parameters in its cloud and convective parameterizations. We show that the auto-optimization leads to 10 % reduction of the overall bias in CAM5, which is already a well-calibrated model, based on a predefined metric that includes precipitation, temperature, humidity, and longwave/shortwave cloud forcing. The computational cost of the entire optimization procedure is about equivalent to a single 12-year atmospheric model simulation. The tuning reduces the large underestimation in the CAM5 longwave cloud forcing by decreasing the threshold relative humidity and the sedimentation velocity of ice crystals in the cloud schemes; it reduces the overestimation of precipitation by increasing the adjustment time in the convection scheme. The physical processes behind the tuned model performance for each targeted field are discussed. Limitations of the automatic tuning are described, including the slight deterioration in some targeted fields that reflect the structural errors of the model. It is pointed out that automatic tuning can be a viable supplement to process-oriented model evaluations and improvement.

54 ENVIRONMENTAL SCIENCES↗

A Robotic Lighting System for Solar Illumination Simulation

This report details our development fo a computer controlled spotlight with four actuated degrees-of-freedom for pan, tilt, linear movement and beam focusing. We review the mechanical structure, kinematic model, calibration data, and control algorithm.

space station freedom ssf solar illumination orbit↗

The connection between forest degradation and urban energy demand in sub-Saharan Africa: a characterization based on high-resolution remote sensing data

Charcoal is a key energy source for urban households in sub-Saharan Africa and charcoal production is the main cause of forest degradation across the region. We used multitemporal high-resolution remote sensing optical imagery to quantify the extent and intensity of forest degradation associated with charcoal production and its impact on forest carbon stocks for the main supplying area of an African capital. This analysis documents the advance of forest degradation and quantifies its aboveground biomass removals over a seven-year period, registering that, between 2013 and 2016, the average annual area under charcoal production was 103 sq.km and the annual aboveground biomass removals reached 1081 000 (SD = 2461) Mg. Kiln densities in the study area rose to 2 kilns-ha, with an average of 90.7 Mg ha^(−1) of extracted aboveground biomass. Charcoal production was responsible for the degradation of 55.5% of the mopane woodlands in the study area between 2013 and 2019. We estimated post-disturbance recovery times using an ecosystem model calibrated for the study area. The simulations showed that recovery times could require up to 150 years for current aboveground biomass extraction rates. The results of the remote sensing analysis and the simulations of the ecosystem model corroborate the unsustainability of the present patterns of charcoal production. The detailed characterization of the spatial and temporal patterns of charcoal production was combined with household survey information to quantify the impact of the urban energy demand of the Maputo urban area on forest carbon stocks. The analysis shows that Maputo charcoal demand was responsible for the annual degradation of up to 175.3 sq.km and that the contribution of the study area to this demand fluctuated between 75% and 33% over the study period. The extent, advance pace and distance from urban centers documented in this study support the idea that forest degradation from charcoal production cannot merely be considered a peri-urban process. The intensity of the aboveground biomass (AGB) removals and its contribution to forest carbon stocks changes is significant at the national and regional levels.

F Sedano↗

SU-E-QI-06: Design and Initial Validation of a Precise Capillary Phantom to Test Perfusion Systems

Purpose: To design a precise perfusion phantom mimicking capillaries of the brain vasculature which could be used to test various perfusion protocols and algorithms which generate perfusion maps. Methods: A perfusion phantom was designed in Solidworks and built using additive manufacturing. The phantom was an overall cylindrical shape of diameter and height 20mm and containing capillaries of 200μm or 300μm which were parallel and in contact making up the inside volume where flow was allowed. We created a flow loop using a peristaltic pump and contrast agent was injected manually. Digital Subtraction Angiographic images and low contrast images with cone beam CT were acquired after the contrast was injected. These images were analyzed by our own code in LabVIEW software and Time-Density Curve, MTT and TTP was calculated. Results: Perfused area was visible in the cone beam CT images; however, individual capillaries were not distinguishable. The Time-Density Curve acquired was accurate, sensitive and repeatable. The parameters MTT, and TTP offered by the phantom were very sensitive to slight changes in the TDC shape. Conclusion: We have created a robust calibrating model for evaluation of existing perfusion data analysis systems. This approach is extremely sensitive to changes in the flow due to the high temporal resolution and could be used as a golden standard to assist developers in calibrating and testing of imaging perfusion systems and software algorithms. Supported by NIH Grant: 2R01EB002873 and an equipment grant from Toshiba Medical Systems Corporation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dynamic reduced order modelling (ROM) of chemical and mechanical processes in CO 2 -cement systems

Damaged wells pose a significant risk of leakage of reservoir fluids stored in a geological CO 2 storage site. The leaking CO 2 can react with well cement and alter its chemical, mechanical, and hydraulic properties. Recently, we have developed an experimentally-calibrated model that couples two-phase flow, reactive transport of brine, cement geochemistry, and geomechanics to predict the leakage of reservoir fluids through a fractured pathway in a cemented well (Iyer, et al. 2018, Iyer, et al. 2017, Walsh, et al. 2013, Walsh, et al. 2014, Walsh, et al. 2014)

58 GEOSCIENCES↗

Linearized Boltzmann-Langevin model for heavy quark transport in hot and dense QCD matter

In relativistic heavy-ion collisions, the production of heavy quarks at large transverse momenta is strongly suppressed compared to proton-proton collisions. In addition, an unexpectedly large azimuthal anisotropy was observed for the emission of charmed hadrons in noncentral collisions. Both observations pose challenges to the theoretical understanding of the coupling between heavy quarks and the quark-gluon plasma produced in these collisions. Transport models for the evolution of heavy quarks in a QCD medium offer the opportunity to study these effects; two of the most successful approaches are based on the linearized Boltzmann transport equation and the Langevin equation. In this work, we develop a hybrid transport model that combines the strengths of both of these approaches: Heavy quarks scatter with medium partons using matrix-elements calculated in perturbative QCD, while between these discrete hard scatterings they evolve using a Langevin equation with empirical transport coefficients to capture the nonperturbative soft part of the interaction. With the hybrid transport model coupled to a state-of-the-art event-by-event bulk evolution model based on 2+1D relativistic viscous fluid dynamics, we study the azimuthal anisotropy and nuclear modification factor of heavy quarks in Pb+Pb collisions at √s = 5.02 TeV. The parameters related to heavy-flavor transport are calibrated using a Bayesian analysis comparing them to available D -meson and B -meson data at the Large Hadron Collider. Lastly, using the calibrated model, we study the implications on heavy-flavor transport properties and predict observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

iTOUGH2 v7.1

iTOUGH2 (inverse TOUGH2) provides inverse modeling capabilities for TOUGH2, a simulator for multi-dimensional, multi-phase, multi-component, non-isothermal flow and transport in fractured porous media. iTOUGH2 performs sensitivity analyses, data-worth analyses, parameter estimation, and uncertainty propagation analyses in geosciences and reservoir engineering and other application areas. iTOUGH2 supports a number of different combinations of fluids and components (equation-of-state (EOS) modules). In addition, the optimization routines implemented in iTOUGH2 can also be used for sensitivity analysis, automatic model calibration, and uncertainty quantification of any external code that uses text-based input and output files using the PEST protocol. iTOUGH2 solves the inverse problem by minimizing a non-linear objective function of the weighted differences between model output and the corresponding observations. Multiple minimization algorithms (derivative-free, gradient-based, and second-order; local and global) are available. iTOUGH2 also performs Latin Hypercube Monte Carlo simulations for uncertainty propagation analyses. A detailed residual and error analysis is provided. This upgrade includes (a) global sensitivity analysis methods, (b) dynamic memory allocation (c) additional input features and output analyses, (d) increased forward simulation capabilities, (e) parallel execution on multicore PCs and Linux clusters, and (f) bug fixes.

Finsterle, Stefan↗

Comparing and contrasting In-Vial and full-scale systems for sparging volatile analytes

Here, in-vial sparging was demonstrated as an effective, practical alternative to a full-scale sparging system for supporting the analysis of volatile constituents. Using elemental mercury (Hg 0 ) and toluene as representative purgeable analytes, the mass removal for various sparge configurations was measured and a reduced order model was developed and validated. In the primary experiments, Hg 0 in the sparge gas was trapped on activated carbon or gold, thermally desorbed, and quantified using atomic absorption or atomic fluorescence spectroscopy. Toluene experiments using the same in-vial sparge apparatus and sparge parameters were performed to demonstrate the applicability of the reduced order model to a broad range of compounds. Toluene removal was tracked by measuring the remaining toluene in sparged aliquots using Ultraviolet-visible (UV–Vis) spectroscopy. For the sparging, flow rates varied from 25 to 75 mL/min for periods from 0 to 30 min. Sparge performance, mass removal as a function of time, and sparge gas volume were measured for both in-vial and full-scale systems. A model based on dimensionless Henry's Law coefficient, normalized sparge gas volume, and fractional extent of equilibrium matched the experimental data for both compounds and provides a practical tool for future applications. For the conditions tested in this study, the calibrated model indicated that the sparge gas in the in-vial system reached approximately 33% of its equilibrium value before exiting the water surface, while a full-scale system reached approximately 100%. The tests validated the quality, reproducibility, and predictability of sparging performance for both full scale and in-vial sparge systems. Related factors such as waste generation, worker risk, and labor were also assessed. Full scale sparge systems provide the advantage of lower detection levels due to larger sample volume, while the in-vial sparge systems provide advantages for most other factors; including automatability, reducing secondary wastes, lessening the need to clean and check the sparge apparatus, and lowering labor and costs. The data and associated reduced order model support continued development and deployment of in-vial sparge platforms as a practical option for analysis of purgeable analytes such as volatile organic compounds and volatile metals/organometallics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the spatiotemporal variability in subsurface thermal regimes across a low-relief polygonal tundra landscape

Vast carbon stocks stored in permafrost soils of Arctic tundra are under risk of release to the atmosphere under warming climate scenarios. Ice-wedge polygons in the low-gradient polygonal tundra create a complex mosaic of microtopographic features. This microtopography plays a critical role in regulating the fine-scale variability in thermal and hydrological regimes in the polygonal tundra landscape underlain by continuous permafrost. Modeling of thermal regimes of this sensitive ecosystem is essential for understanding the landscape behavior under the current as well as changing climate. Here, we present an end-to-end effort for high-resolution numerical modeling of thermal hydrology at real-world field sites, utilizing the best available data to characterize and parameterize the models. We also develop approaches to model the thermal hydrology of polygonal tundra and apply them at four study sites near Barrow, Alaska, spanning across low to transitional to high-centered polygons, representing a broad polygonal tundra landscape. A multiphase subsurface thermal hydrology model (PFLOTRAN) was developed and applied to study the thermal regimes at four sites. Using a high-resolution lidar digital elevation model (DEM), microtopographic features of the landscape were characterized and represented in the high-resolution model mesh. The best available soil data from field observations and literature were utilized to represent the complex heterogeneous subsurface in the numerical model. Simulation results demonstrate the ability of the developed modeling approach to capture &ndash; without recourse to model calibration &ndash; several aspects of the complex thermal regimes across the sites, and provide insights into the critical role of polygonal tundra microtopography in regulating the thermal dynamics of the carbon-rich permafrost soils. Moreover, areas of significant disagreement between model results and observations highlight the importance of field-based observations of soil thermal and hydraulic properties for modeling-based studies of permafrost thermal dynamics, and provide motivation and guidance for future observations that will help address model and data gaps affecting our current understanding of the system.

54 ENVIRONMENTAL SCIENCES↗

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

Vulnerability of Amazon forests to storm-driven tree mortality

Tree mortality is a key driver of forest community composition and carbon dynamics. Strong winds associated with severe convective storms are dominant natural drivers of tree mortality in the Amazon. Why forests vary with respect to their vulnerability to wind events and how the predicted increase in storm events might affect forest ecosystems within the Amazon are not well understood. We found that windthrows are common in the Amazon region extending from northwest (Peru, Colombia, Venezuela, and west Brazil) to central Brazil, with the highest occurrence of windthrows in the northwest Amazon. More frequent winds, produced by more frequent severe convective systems, in combination with well-known processes that limit the anchoring of trees in the soil, help to explain the higher vulnerability of the northwest Amazon forests to winds. Projected increases in the frequency and intensity of convective storms in the Amazon have the potential to increase wind-related tree mortality. A forest demographic model calibrated for the northwestern and the central Amazon showed that northwestern forests are more resilient to increased wind-related tree mortality than forests in the central Amazon. Our study emphasizes the importance of including wind-related tree mortality in model simulations for reliable predictions of the future of tropical forests and their effects on the Earth' system.

54 ENVIRONMENTAL SCIENCES↗