Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model physics”

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 739 records · Page 41

Modeling the Self-organized Critical Behavior of the Plasma Sheet Reconnection Dynamics

Analyses of Polar UVI auroral image data reviewed in our other presentation at this meeting (V. Uritsky, A. Klimas) show that bright night-side high-latitude UV emissions exhibit so many of the key properties of systems in self-organized criticality (SOC) that an alternate interpretation has become virtually impossible. It is now necessary to find and model the source of this behavior. We note that the most common models of self-organized criticality are numerical sandpiles. These are, at root, models that govern the transport of some quantity from a region where it is loaded to another where it is unloaded. Transport is enabled by the excitation of a local threshold instability; it is intermittent and bursty, and it exhibits a number of scale-free statistical properties. Searching for a system in the magnetosphere that is analogous and that, in addition, is known to produce auroral signatures, we focus on the reconnection dynamics of the plasma sheet. In our previous work, a driven reconnection model has been constructed and has been under study. The transport of electromagnetic (primarily magnetic) energy carried by the Poynting flux into the reconnection region of the model has been examined. All of the analysis techniques, and more, that have been applied to the auroral image data have also been applied to this Poynting flux. Here, we report new results showing that this model also exhibits so many of the key properties of systems in self-organized criticality that an alternate interpretation is implausible. Further, we find a strong correlation between these key properties of the model and those of the auroral UV emissions. We suggest that, in general, the driven reconnection model is an important step toward a realistic plasma physical model of self-organized criticality and we conclude, more specifically, that it is also a step in the right direction toward modeling the multiscale reconnection dynamics of the magnetotail.

Klimas, Alex↗

Modeling the Self-organized Critical Behavior of Earth's Plasma Sheet Reconnection Dynamics

Analyses of Polar UVI auroral image data show that bright night-side high-latitude W emissions exhibit so many of the key properties of systems in self-organized criticality that an alternate interpretation has become virtually impossible. These analyses will be reviewed. It is now necessary to find and model the source of this behavior. We note that the most common models of self-organized criticality are numerical sandpiles. These are, at root, models that govern the transport of some quantity from a region where it is loaded to another where it is unloaded. Transport is enabled by the excitation of a local threshold instability; it is intermittent and bursty, and it exhibits a number of scale-free statistical properties. Searching for a system in the magnetosphere that is analogous and that, in addition, is known to produce auroral signatures, we focus on the reconnection dynamics of the magnetotail plasma sheet. In our previous work, a driven reconnection model has been constructed and has been under study. The transport of electromagnetic (primarily magnetic) energy carried by the Poynting flux into the reconnection region of the model has been examined. All of the analysis techniques (and more) that have been applied to the auroral image data have also been applied to this Poynting flux. New results will be presented showing that this model also exhibits so many of the key properties of systems in self-organized criticality that an alternate interpretation is implausible. A strong correlation between these key properties of the model and those of the auroral UV emissions will be demonstrated. We suggest that, in general, the driven reconnection model is an important step toward a realistic plasma physical model of self-organized criticality and we conclude, more specifically, that it is also a step in the right direction toward modeling the multiscale reconnection dynamics of the magnetotail.

Klimas, Alexander J.↗

Planck Early Results. XV. Spectral Energy Distributions and Radio Continuum Spectra of Northern Extragalactic Radio Sources

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources. based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multi frequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper. physical modelling of the synchrotron bump using multiple components. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

Aatrokoski, J.↗

Brief History of Agricultural Systems Modeling

Agricultural systems science generates knowledge that allows researchers to consider complex problems or take informed agricultural decisions. The rich history of this science exemplifies the diversity of systems and scales over which they operate and have been studied. Modeling, an essential tool in agricultural systems science, has been accomplished by scientists from a wide range of disciplines, who have contributed concepts and tools over more than six decades. As agricultural scientists now consider the next generation models, data, and knowledge products needed to meet the increasingly complex systems problems faced by society, it is important to take stock of this history and its lessons to ensure that we avoid re-invention and strive to consider all dimensions of associated challenges. To this end, we summarize here the history of agricultural systems modeling and identify lessons learned that can help guide the design and development of next generation of agricultural system tools and methods. A number of past events combined with overall technological progress in other fields have strongly contributed to the evolution of agricultural system modeling, including development of process-based bio-physical models of crops and livestock, statistical models based on historical observations, and economic optimization and simulation models at household and regional to global scales. Characteristics of agricultural systems models have varied widely depending on the systems involved, their scales, and the wide range of purposes that motivated their development and use by researchers in different disciplines. Recent trends in broader collaboration across institutions, across disciplines, and between the public and private sectors suggest that the stage is set for the major advances in agricultural systems science that are needed for the next generation of models, databases, knowledge products and decision support systems. The lessons from history should be considered to help avoid roadblocks and pitfalls as the community develops this next generation of agricultural systems models.

agricultural systems↗

Benchmarking and Exploring Parameter Space of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation

High intensity proton pulses strike the Spallation Neutron Source (SNS)’s mercury target to provide bright neutron beams. These strikes deposit extensive energy into the mercury and its steel vessel. Prediction of the resultant loading on the target is difficult when helium gas is intentionally injected into the mercury to reduce the loading and to mitigate the pitting damage on the vessel. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. We present a study comparing the measured target strains in the SNS target station with the simulation results of the solid mechanics simulation framework. We investigate a wide range of various physical model parameters, including the number of bubble families, bubble size distribution, viscosity, surface tension, etc. to understand their impact on simulation accuracy. Our initial findings reveal that using 8-10 bubble families in the model renders a simulation strain envelope that covers the experimental ones. Further optimization studies are planned to predict the strain response more accurately.

Lin, Lianshan↗

Development of an Enhanced Radiation Physics Toolset for Modeling Spectral and Imaging Signatures in the Warm Dense Matter Experiments

Radiative and atomic processes in plasmas play a critical role in a wide variety of high energy density laboratory plasma (HEDLP) experiments. The emission, absorption, and transport of radiation can strongly affect the overall energetics and evolution of such plasmas. In addition, radiation-based diagnostics – including imaging, spectroscopy, and absolute flux measurements – are widely used to determine key features of HEDLPs. To advance our understanding of HEDLP science, it is vital to have high-fidelity computational physics tools that have well-tested radiation physics modeling, and that are readily accessible to researchers in the HEDLP community. Simulations play an extremely important role for planning and designing the experiments, as well as for post-experiment data analysis. Prism Computational Sciences develops software that is used by National Laboratories and universities (including five members of LaserNetUS network). Prominent examples of such research efforts include z-pinch and short-pulse laser experiments designed to study the basic physics of photoionized plasmas and photoionization fronts, as well as their application to astrophysical plasmas. The main effort was dedicated to the development of non-equilibrium equation-of-state (EOS) models within the HELIOS-CR code, a hydrodynamics code with inline collisional-radiative atomic kinetics. Gas cell experiments on Z and Omega demonstrated the importance of non-equilibrium effects on atomic kinetics in photoionized plasmas. Recent proof-of-principle experiments on Omega EP confirmed the advantages of using a short-pulse laser to create an intense radiation drive, leading to additional experiments being proposed. Photoionization front experiments at LLE also emphasizes the importance of radiation and atomic physics. In both studies, HELIOS-CR simulations played a crucial role in computing non-equilibrium opacities and ionization distributions. A newly developed non-LTE EOS model will help addressing possible non-equilibrium effects, on for example specific heat, and their influence in plasma evolution. Prism also implemented support for open-source atomic data generated by the Flexible Atomic Code. This allows researchers to generate custom atomic tables and use them within the complex framework of simulation tools developed by Prism. The ability to use open-source atomic data would be extremely valuable for hydrodynamics and spectroscopic simulations that include high-Z materials, e.g., picosecond x-ray pulse generation experiments. Support for new atomic structures was fully implemented, and the data can be used by all simulation tools developed at Prism: radiation-hydrodynamics, imaging and spectroscopy, EOS and opacity. The development resulted in a significant fidelity enhancement to the simulations tools developed by Prism that are currently used in other cutting-edge experiments including: opacity measurement experiments performed to both understand the basic radiative and atomic properties of plasmas as well as provide data for more accurately modeling the internal structure of the Sun and other stars, high-intensity short-pulse laser experiments performed to develop short-wavelength light sources for use as backlighters and to investigate fast ignition concepts for inertial fusion energy; capsule implosion experiments designed to develop inertial fusion as an energy source, etc.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Metal-insulator transition in transition metal dichalcogenide heterobilayer moiré superlattices

Moiré superlattices formed in two-dimensional semiconductor heterobilayers provide a new realization of Hubbard model physics in which the number of electrons per effective atom can be tuned at will. Here we report on an exact diagonalization study of the electronic properties of half-filled narrow moiré bands in which correlation strengths are varied by changing twist angles or interaction strengths. We construct a phase diagram for the bilayer, identifying where the metal-insulator phase transition occurs, estimating the sizes of the charge gaps in the insulating phase, and commenting on the nature of the transition and the importance of subdominant interaction parameters.

2-dimensional systems↗

A two-dimensional solar model

A two-dimensional, nearly flux constant solar atmosphere which includes a physical model of granulation is presented. If the inhomogeneous character of the solar photosphere is neglected in the construction of empirical models, the result is an underestimate of the temperature gradient in the layers below a Rosseland optical depth of 0.5 and an overestimate of the convective flux penetration into the observable layers. The two-dimensional model adequately reproduces the mean limb darkening, the magnitude and center-to-limb variation of the intensity fluctuations, and the rms vertical and horizontal convective velocities.

Nelson, G. D.↗

A model of the CO2 exchanges between biosphere and atmosphere in the tundra

A physical model of the soil thermal regime in a permafrost terrain has been developed and validated with soil temperature measurements at Barrow, Alaska. The model calculates daily soil temperatures as a function of depth and average moisture contents of the organic and mineral layers using a set of five climatic variables, i.e., air temperature, precipitation, cloudiness, wind speed, and relative humidity. The model is not only designed to study the impact of climate change on the soil temperature and moisture regime, but also to provide the input to a decomposition and net primary production model. In this context, it is well known that CO2 exchanges between the terrestrial biosphere and the atmosphere are driven by soil temperature through decomposition of soil organic matter and root respiration. However, in tundra ecosystems, net CO2 exchange is extremely sensitive to soil moisture content; therefore it is necessary to predict variations in soil moisture in order to assess the impact of climate change on carbon fluxes. To this end, the present model includes the representation of the soil moisture response to changes in climatic conditions. The results presented in the foregoing demonstrate that large errors in soil temperature and permafrost depth estimates arise from neglecting the dependence of the soil thermal regime on soil moisture contents. Permafrost terrain is an example of a situation where soil moisture and temperature are particularly interrelated: drainage conditions improve when the depth of the permafrost increases; a decrease in soil moisture content leads to a decrease in the latent heat required for the phase transition so that the heat penetrates faster and deeper, and the maximum depth of thaw increases; and as excepted, soil thermal coefficients increase with moisture.

Labgaa, Rachid R.↗

Use of Physics to Improve Solar Forecast: Physics-Informed Persistence Models for Simultaneously Forecasting GHI, DNI, and DHI

Observation-based statistical models have been widely used in forecasting solar energy; however, existing models often lack a clear relation to physics and are limited largely to global horizontal irradiance (GHI) forecasts over relatively short time horizons (< 1 hour). Incorporating physics into observation-based models, increasing forecast time horizons and developing a model system for forecasting not only GHI but also direct normal irradiance (DNI) and diffuse horizontal irradiance (DHI) remain challenging, especially under cloudy conditions because of complex cloud-radiation interactions. This work attempts to address these challenges by developing a hierarchy of four new physics-informed persistence models that can be used to simultaneously forecast GHI, DNI and DHI. The decade-long measurements (1998 to 2014) at the U.S. Department of Energy's Atmospheric Radiation Measurement (ARM)'s Southern Great Plains (SGP) Central Facility site are used to evaluate the performance of the new models. Overall, the results show that the new physics-informed forecast models generally outperform the simple and smart persistence models, and improve the forecast accuracy at lead times from 1.25 hours up to 6 hours. Further analysis reveals that the forecast error is highly related to the error and temporal variability of the assumed cloud predictor. The best model for forecasting different radiative components can be explained by the relationship between solar irradiances and cloud properties.

54 ENVIRONMENTAL SCIENCES↗

Ray-tracing image simulations of transparent objects with complex shape and inhomogeneous refractive index

Optical images of transparent three-dimensional objects can be different from a replica of the object’s cross section in the image plane, due to refraction at the surface or in the body of the object. Simulations of the object’s image are thus needed for the visualization and validation of physical models. We report ray tracing image simulations that achieved high physical fidelity, reproducing optical behaviors and image features not rendered in previous studies. We replicated brightfield microscopy images of drops with complex shapes, and images of pressure and shock waves traveling inside them. For high physical fidelity, the simulations must replicate the spatial and angular distribution of illumination rays, and both the experiment and the simulation must be designed for accurate optical modeling. The simulations are highly sensitive to the properties of the drops and can be used to diagnose and refine fluid dynamics models. The simulated images can also be optimized to extract multiple 3D properties from experimental images. Compared to specialized single-shot 3D imaging methods, this approach has the advantage that it preserves the experimental simplicity, the high resolution, and the visual interpretability characteristic to basic optical imaging. The techniques introduced here are directly applicable to optical microscopy, so they can be used in other fields, such as microfluidics and biology, to expand the type and the accuracy of three-dimensional information that can be extracted from basic optical images.

Cavitation↗

Calculation of atmospheric loss from microwave radiometric noise temperature measurements

Microwave propagation loss in the atmosphere can be inferred from microwave radiometric noise temperature measurements. The relevant equations are given and a derivation and calculation is made assuming various physical models. Comparison is made with the commonly used lumped element atmospheric model (isothermal and uniform loss) and the model with linear temperature and exponential loss distributions. The results are useful for estimating the integral inversion differences due to the model selection. This indicates that the commonly used lumped element atmospheric model is a very good approximation with judicious choice of the effective physical temperature. For the worst case comparison, the lumped element model agrees with the variable parameter model within 0.2 dB up to a propagation loss of 3 dB.

Stelzried, C.↗

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗

Relationships Between the Bulk-Skin Sea Surface Temperature Difference, Wind, and Net Air-Sea Heat Flux

The primary purpose of this project was to evaluate and improve models for the bulk-skin temperature difference to the point where they could accurately and reliably apply under a wide variety of environmental conditions. To accomplish this goal, work was conducted in three primary areas. These included production of an archive of available data sets containing measurements of the skin and bulk temperatures and associated environmental conditions, evaluation of existing skin layer models using the compiled data archive, and additional theoretical work on the development of an improved model using the data collected under diverse environmental conditions. In this work we set the basis for a new physical model of renewal type, and propose a parameterization for the temperature difference across the cool skin of the ocean in which the effects of thermal buoyancy, wind stress, and microscale breaking are all integrated by means of the appropriate renewal time scales. Ideally, we seek to obtain a model that will accurately apply under a wide variety of environmental conditions. A summary of the work in each of these areas is included in this report. A large amount of work was accomplished under the support of this grant. The grant supported the graduate studies of Sandra Castro and the preparation of her thesis which will be completed later this year. This work led to poster presentations at the 1999 American Geophysical Union Fall Meeting and 2000 IGARSS meeting. Additional work will be presented in a talk at this year's American Meteorological Society Air-Sea Interaction Meeting this May. The grant also supported Sandra Castro during a two week experiment aboard the R/P Flip (led by Dr. Andrew Jessup of the Applied Physics Laboratory) to help obtain additional shared data sets and to provide Sandra with a fundamental understanding of the physical processes needed in the models. In a related area, the funding also partially supported Dr. William Emery and Daniel Baldwin in the preparation of their publication "Accuracy of in situ sea surface temperatures used to calibrate infrared satellite measurements". The remainder of this report is drawn from these publications and presentations.

Emery, William J.↗

Stream Temperature Prediction in a Shifting Environment: Explaining the Influence of Deep Learning Architecture

Stream temperature is a fundamental control on ecosystem health. Recent efforts incorporating process guidance into deep learning models for predicting stream temperature have been shown to outperform existing statistical and physical models. This performance is in part because deep learning architectures can actively learn spatiotemporal relationships that govern how water and energy propagate through a river network. However, exploration of how spatiotemporal awareness and process guidance influence a model's generalizability under shifting environmental conditions such as climate change is limited. Here, we use Explainable Artificial Intelligence (XAI) to interrogate how differing deep learning architectures affect a model's learned spatial and temporal dependencies, and how those learned dependencies affect a model's ability to maintain high accuracy when applied to unseen environmental conditions. Using the Delaware River Basin in the northeastern United States as a test case, we compare two spatiotemporally aware process–guided deep learning models for predicting stream temperature (a recurrent graph convolution network—RGCN, and a temporal convolution graph model—Graph WaveNet). Both models achieve equally high predictive performance when testing data are well represented in the training data (test root mean squared errors of 1.64°C and 1.65°C); however, Graph WaveNet significantly outperforms RGCN in 4 out of 5 experiments where test partitions represent different types of unseen environmental conditions. XAI results show that the architecture of Graph WaveNet leads to learned spatial relationships with greater fidelity to physical processes, and that this fidelity improves the generalizability of the model when applied to shifting and/or unseen environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

HED Hydrodynamics in the Common Modeling Framework

Within the Common Modeling Framework (CMF), the HED Hydrodynamics Project has implemented a common methodology, or "Authority", for modeling high energy density (HED) experiments by extending development of the Eulerian Applications Project code Authority, "EAPA," to include relevant physics models, customization tools and templates, and in situ and post-processing capabilities. Additionally, several models for small-scale HED experiments were successfully added to the CMF and simulations reproduce the results of preexisting, benchmarked input decks. We expect use of the CMF to improve our effectiveness in developing predictive capabilities for experiments ranging from small-scale planar single-interface, single-shock to multi-interface, multi-shock configurations and up through multi-shell inertial confinement fusion (ICF) implosions. Additionally, this will strengthen our ability to design future targeted experiments. Using the CMF infrastructure and Authorities for HED modeling has three key impacts. First, it enables a version-controlled, pedigreed and archived base model for each experiment with tools that allow easy setup, execution and analysis of simulations to assess and/or expose sensitivities to various parameters. Second, having a common base set of inputs and physics definitions helps identify systematic differences between models in the codes and experimental data to support code validation, and build confidence in our predictive capabilities. Finally, designed to work in concert with multiple programs and projects including PEM, IC, OES, V&V and DSW, the CMF supports improved collaboration and integration through shared model definitions, data, simulations and results, enhancing the workflow required to underwrite the knowledge base and capabilities ultimately required for stockpile stewardship

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-energy-density-physics measurements in implosions using Bayesian inference

Convergent high-energy-density (HED) experimental platforms are used to study matter under some of the most extreme conditions that can be produced on Earth, comparable to the interior of stars. There are many challenges in using these systems for fundamental measurements currently being addressed by new analysis methods, such as the combination of a reduced physics model and Bayesian inference, allowing a self-consistent inference of physical quantities with a robust error analysis. These methods in combination with simple (as compared to inertial confinement fusion implosions) implosion platforms, which can be modified to show sensitivity to different physical mechanisms of interest, are used to study the physical properties of matter under extreme conditions. This work discusses a subset of implosion targets for studying opacity effects, electron-ion equilibration, and thermal conductivity and, as an example, a system consisting of a thick-shelled, gas-filled laser-direct-drive implosion is used to show how a reduced model and Bayesian inference can help inform experimental design decisions such as diagnostic choice. It is shown that for this system that a combination of neutron and x-ray self-emission diagnostics is critical for constraining the details of the thermodynamic states in the system and that the conductivity exponent in a Spitzer like framework can be constrained to the 30% level in deuterium at gigabar conditions. Furthermore, this process can be applied to many HED systems to make underlying model assumptions explicit and facilitate experimental design and analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulating water dynamics related to pedogenesis across space and time: Implications for four-dimensional digital soil mapping

Digital soil mapping (DSM) relies on machine-learning and geostatistics to represent soil property observations across space. DSM techniques are powerful but often empirical, being limited to the quality and density of point samples. Water dynamics are closely related to soil variability, and the physics that govern water movement are well known. Hydrological properties can hence be simulated by physical models through space and time, unveiling key characteristics about soils. We propose the use of hydrologic models to map soils across the surface (2D), depth (1D), and time (1D)–which provides a 4D approach to digital soil mapping (4DSM). The Distributed Hydrology Soil Vegetation Model (DHSVM) was applied to a watershed currently under pasture. Moisture sensors and wells were installed at different depths in the watershed on summit, sideslope and toeslope positions to validate the model. DHSVM simulations of soil moisture distribution and depth to saturation were performed during the hydrological year (October 2008-September 2009). Clusters of similar pixels based on soil moisture values were determined using Dynamic Time Warping (DTW) to align temporal data and K-means. Clustering was performed both seasonally and for the entire year. Temporal patterns simulated by DHSVM matched measurements given by moisture sensors and wells. Seasonal clusters differed from the annual cluster. Distinct clusters were observed for each season and with depth, showing that spatiotemporal soil variability is lost when statically assessing soils. Spatiotemporal clusters corroborated field observations of fragipan occurrence not explicitly spatially mapped by Soil Survey Geographic Database (SSURGO). If a connection can be made between water and soils, static and dynamic soil variability can be predicted using physically based hydrologic models. Hydrologic models can benefit soil mapping by enabling reliable 4D simulation of water dynamics, which are fundamental to soil variability and soil classification and directly relate to biological, physical and chemical soil processes not captured by typical soil sampling protocols.

54 ENVIRONMENTAL SCIENCES↗