Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Quality”

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

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

At least 1,243 records · Page 69

Grand challenges in the design, manufacture, and operation of future wind turbine systems

Abstract. Wind energy is foundational for achieving 100 % renewable electricity production, and significant innovation is required as the grid expands and accommodates hybrid plant systems, energy-intensive products such as fuels, and a transitioning transportation sector. The sizable investments required for wind power plant development and integration make the financial and operational risks of change very high in all applications but especially offshore. Dependence on a high level of modeling and simulation accuracy to mitigate risk and ensure operational performance is essential. Therefore, the modeling chain from the large-scale inflow down to the material microstructure, and all the steps in between, needs to predict how the wind turbine system will respond and perform to allow innovative solutions to enter commercial application. Critical unknowns in the design, manufacturing, and operability of future turbine and plant systems are articulated, and recommendations for research action are laid out. This article focuses on the many unknowns that affect the ability to push the frontiers in the design of turbine and plant systems. Modern turbine rotors operate through the entire atmospheric boundary layer, outside the bounds of historic design assumptions, which requires reassessing design processes and approaches. Traditional aerodynamics and aeroelastic modeling approaches are pressing against the limits of applicability for the size and flexibility of future architectures and flow physics fundamentals. Offshore wind turbines have additional motion and hydrodynamic load drivers that are formidable modeling challenges. Uncertainty in turbine wakes complicates structural loading and energy production estimates, both around a single plant and for downstream plants, which requires innovation in plant operations and flow control to achieve full energy capture and load alleviation potential. Opportunities in co-design can bring controls upstream into design optimization if captured in design-level models of the physical phenomena. It is a research challenge to integrate improved materials into the manufacture of ever-larger components while maintaining quality and reducing cost. High-performance computing used in high-fidelity, physics-resolving simulations offer opportunities to improve design tools through artificial intelligence and machine learning, but even the high-fidelity tools are yet to be fully validated. Finally, key actions needed to continue the progress of wind energy technology toward even lower cost and greater functionality are recommended.

17 WIND ENERGY↗

Evaluation of Modeling NO2 Concentrations Driven by Satellite-Derived and Bottom-Up Emission Inventories Using In-Situ Measurements Over China

Chemical transport models together with emission inventories are widely used to simulate NO2 concentrations over China, but validation of the simulations with in situ measurements has been extremely limited. Here we use ground measurements obtained from the air quality monitoring network recently developed by the Ministry of Environmental Protection of China to validate modeling surface NO2 concentrations from the CHIMERE regional chemical transport model driven by the satellite-derived DECSO and the bottom-up MIX emission inventories. We applied a correction factor to the observations to account for the interferences of other oxidized nitrogen compounds (NOz), based on the modeled ratio of NO2 to NOz. The model accurately reproduces the spatial variability in NO2 from in situ measurements, with a spatial correlation coefficient of over 0.7 for simulations based on both inventories. A negative and positive bias is found for the simulation with the DECSO (slopeD0.74 and 0.64 for the daily mean and daytime only) and the MIX (slopeD1.3 and 1.1) inventories, respectively, suggesting an underestimation and overestimation of NOx emissions from corresponding inventories. The bias between observed and modeled concentrations is reduced, with the slope dropping from 1.3 to 1.0 when the spatial distribution of NOx emissions in the DECSO inventory is applied as the spatial proxy for the MIX inventory, which suggests an improvement of the distribution of emissions between urban and suburban or rural areas in the DECSO inventory compared to that used in the bottom-up inventory. A rough estimate indicates that the observed concentrations, from sites predominantly placed in the populated urban areas, may be 10-40% higher than the corresponding model grid cell mean. This reduces the estimate of the negative bias of the DECSO-based simulation to the range of -30 to 0% on average and more firmly establishes that the MIX inventory is biased high over major cities. The performance of the model is comparable over seasons, with a slightly worse spatial correlation in summer due to the difficulties in resolving the more active NOx photochemistry and larger concentration gradients in summer by the model. In addition, the model well captures the daytime diurnal cycle but shows more significant disagreement between simulations and measurements during nighttime, which likely produces a positive model bias of about 15% in the daily mean concentrations. This is most likely related to the uncertainty in vertical mixing in the model at night.

NO2 modelling↗

Scale-Resolving Simulations of Low-Pressure Turbine Cascades with Wall Roughness Using A Spectral-Element Method

The accurate prediction of wall-roughness effects in turbomachinery is becoming critical as turbine designers address airfoil surface quality and degradation concerns arising from the shift to advanced ceramic matrix composite (CMC) or additively-manufactured airfoils operating in higher temperature environments. In this paper, a recently developed computational capability for accurate and efficient scale-resolving simulations of turbomachinery is extended to analyze the boundary- layer separation and transition characteristics in a rough-wall low-pressure turbine (LPT) cascade. The computational capability is based on an entropy-stable discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy, and is implemented in an efficient manner for a modern high performance computer architecture. Results from the scale-resolving simulations of both smooth and rough airfoil cascades are presented and compared to previous experiments and numerical simulations. The results show that the suction surface boundary layer undergoes laminar separation, transition, and turbulent reattachment for the smooth airfoil cascade, while in the presence of roughness the separation and transition behavior of the suction surface boundary layer is substantially modified. The differences between the smooth and rough airfoil cascades are then highlighted by a detailed analysis of their respective turbulent flow fields.

Spectral-Element↗

An Investigation of Rotorcraft Stability-Phase Margin Requirements in Hover

A cooperative study was performed to investigate the handling quality effects from reduced flight control system stability margins, and the trade-offs with higher disturbance rejection bandwidth (DRB). The piloted simulation study, perform on the NASA-Ames Vertical Motion Simulator, included three classes of rotorcraft in four configurations: a utility-class helicopter; a medium-lift helicopter evaluated with and without an external slung load; and a large (heavy-lift) civil tiltrotor aircraft. This large aircraft also allowed an initial assessment of ADS-33 handling quality requirements for an aircraft of this size. Ten experimental test pilots representing the U.S. Army, Marine Corps, NASA, rotorcraft industry, and the German Aerospace Center (DLR), evaluated the four aircraft configurations, for a range of flight control stability-margins and turbulence levels, while primarily performing the ADS-33 Hover and Lateral Reposition MTEs. Pilot comments and aircraft-task performance data were analyzed. The preliminary stability margin results suggest higher DRB and less phase margin cases are preferred as the aircraft increases in size. Extra care will need to be taken to assess the influence of variability when nominal flight control gains start with reduced margins. Phase margins as low as 20-23 degrees resulted in low disturbance-response damping ratios, objectionable oscillations, PIO tendencies, and a perception of an incipient handling qualities cliff. Pilot comments on the disturbance response of the aircraft correlated well to the DRB guidelines provided in the ADS-33 Test Guide. The A D-3S3 mid-term response-to-control damping ratio metrics can be measured and applied to the disturbance-response damping ratio. An initial assessment of LCTR yaw bandwidth shows the current Level 1 boundary needs to be relaxed to help account for a large pilot off-set from the c.g. Future efforts should continue to investigate the applicability/refinement of the current ADS-33 requirements to large vehicles, like an LCTR.

Blanken, Chris L.↗

How Climate and Data Quality Impact Photovoltaic Performance Loss Rate Estimations

Different data pipelines and statistical methods are applied to photovoltaic (PV) performance datasets to quantify the performance loss rate (PLR). Since the real values of PLR are unknown, a variety of unvalidated values are reported. As such, the PV industry commonly assumes PLR based on statistically extracted ranges from the literature. However, the accuracy and uncertainty of PLR depend on several parameters including seasonality, local climatic conditions, and the response of a particular PV technology. In addition, the specific data pipeline and statistical method used affect the accuracy and uncertainty. To provide insights, a framework of (≈200 million) synthetic simulations of PV performance datasets using data from different climates is developed. Time series with known PLR and data quality are synthesized, and large parametric studies are conducted to examine the accuracy and uncertainty of different statistical approaches over the contiguous US, with an emphasis on the publicly available and “standardized” library, RdTools . In the results, it is confirmed that PLRs from RdTools are unbiased on average, but the accuracy and uncertainty of individual PLR estimates vary with climate zone, data quality, PV technology, and choice of analysis workflow. Best practices and improvement recommendations based on the findings of this study are provided.

14 SOLAR ENERGY↗

Handling Qualities Considerations in Control Allocation for Multicopters

Due to multicopters’ role in advanced air mobility, there has been an increasing interest in handling qualities. This paper investigates the effect of control allocation of multicopters on handling qualities. Typically, multicopter aircraft are trimmed with equal speed on each rotor. The force generated by each rotor can be modeled, where thrust and torque of the rotor are related to the square of the rotor speed. The dynamic equations indicate that trimming pairs of motors at different RPM results in an increase of control effectiveness in the associated axis, improving bandwidth and disturbance rejection. A tradeoff occurs where power usage is no longer optimized, but handling qualities are improved in that control axis. A blade element model for multicopter, RMAC, was used to simulate how this off-nominal control allocation would affect the dynamics of a multicopter. System identification flight tests were used to validate the trends seen in RMAC. Control systems were optimized to provide improved disturbance rejection bandwidth in the axis associated with increased rotational speed. Unmanned mission task elements were performed in flight test, showing improved tracking in the lateral axis by using off-nominal mixing to increase trim speed on the rolling motors.

Handling Qualities↗

Air quality impact of the Northern California Camp Fire of November 2018

The Northern California Camp Fire that took place in November 2018 was one of the most damaging environmental events in California history. Here, we analyze ground-based station observations of airborne particulate matter that has a diameter <2.5 µm (PM 2.5 ) across Northern California and conduct numerical simulations of the Camp Fire using the Weather Research and Forecasting model online coupled with chemistry (WRF-Chem). Simulations are evaluated against ground-based observations of PM 2.5 , black carbon, and meteorology, as well as satellite measurements, such as Tropospheric Monitoring Instrument (TROPOMI) aerosol layer height and aerosol index. The Camp Fire led to an increase in Bay Area PM 2.5 to over 50 µg m –3 for nearly 2 weeks, with localized peaks exceeding 300 µg m –3 . Using the Visible Infrared Imaging Radiometer Suite (VIIRS) high-resolution fire detection products, the simulations reproduce the magnitude and evolution of surface PM 2.5 concentrations, especially downwind of the wildfire. The overall spatial patterns of simulated aerosol plumes and their heights are comparable with the latest satellite products from TROPOMI. WRF-Chem sensitivity simulations are carried out to analyze uncertainties that arise from fire emissions, meteorological conditions, feedback of aerosol radiative effects on meteorology, and various physical parameterizations, including the planetary boundary layer model and the plume rise model. Downwind PM 2.5 concentrations are sensitive to both flaming and smoldering emissions over the fire, so the uncertainty in the satellite-derived fire emission products can directly affect the air pollution simulations downwind. Our analysis also shows the importance of land surface and boundary layer parameterization in the fire simulation, which can result in large variations in magnitude and trend of surface PM 2.5 . Inclusion of aerosol radiative feedback moderately improves PM 2.5 simulations, especially over the most polluted days. Results of this study can assist in the development of data assimilation systems as well as air quality forecasting of health exposures and economic impact studies.

54 ENVIRONMENTAL SCIENCES↗

A Geostationary Instrument Simulator for Aerosol Observing System Simulation Experiments

In the near future, there will be several new instruments measuring atmospheric composition from geostationary orbit over North America, East Asia, and Europe. This constellation of satellites will provide high resolution, time resolved measurements of trace gases and aerosols for monitoring air quality and tracking pollution sources. This paper describes a detailed, fast, and accurate (less than 1.0% uncertainty) method for calculating synthetic top of the atmosphere (TOA) radiances from a global simulation with a mesoscale free running model, the GEOS-5 Nature Run, for remote sensing instruments in geostationary orbit that measure in the ultraviolet-visible spectral range (UV-Vis). Generating these synthetic observations is the first step of an Observing System Simulation Experiment (OSSE), a framework for evaluating the impact of a new observation or algorithm. This paper provides details of the model sampling, aerosol and cloud optical properties, surface reflectance modeling, Rayleigh scattering calculations, and a discussion of the uncertainties of the simulated TOA radiance. An application for the simulated TOA radiance observations is demonstrated in the manuscript. Simulated TEMPO (Tropospheric Emissions: Monitoring of Pollution) and GOES-R (Geostationary Operational Environmental Satellites) observations were used to show how observations from the two instruments could be combined to facilitate aerosol type discrimination. The results demonstrate the viability of a detailed instrument simulator for radiance measurements in the UV-Vis that is capable of accurately simulating high resolution, time-resolved measurements with reasonable computational efficiency.

Castellanos, Patricia↗

Piloted Simulation Assessment of a High-Speed Civil Transport Configuration

An assessment of a proposed configuration of a high-speed civil transport was conducted by using NASA and industry research pilots. The assessment was conducted to evaluate operational aspects of the configuration from a pilot's perspective, with the primary goal being to identify potential deficiencies in the configuration. The configuration was evaluated within and at the limits of the design operating envelope to determine the suitability of the configuration to maneuver in a typical mission as well as in emergency or envelope-limit conditions. The Cooper-Harper rating scale was used to evaluate the flying qualities of the configuration. A summary flying qualities metric was also calculated. The assessment was performed in the Langley six-degree-of-freedom Visual Motion Simulator. The effect of a restricted cockpit field-of-view due to obstruction by the vehicle nose was not included in this study. Tasks include landings, takeoffs, climbs, descents, overspeeds, coordinated turns, and recoveries from envelope limit excursions. Emergencies included engine failures, loss of stability augmentation, engine inlet unstarts, and emergency descents. Minimum control speeds and takeoff decision, rotation, and safety speeds were also determined.

Jackson, E. Bruce↗

Rotorcraft Research at the NASA Vertical Motion Simulator

In the 1970 s the role of the military helicopter evolved to encompass more demanding missions including low-level nap-of-the-earth flight and operation in severely degraded visual environments. The Vertical Motion Simulator (VMS) at the NASA Ames Research Center was built to provide a high-fidelity simulation capability to research new rotorcraft concepts and technologies that could satisfy these mission requirements. The VMS combines a high-fidelity large amplitude motion system with an adaptable simulation environment including interchangeable and configurable cockpits. In almost 30 years of operation, rotorcraft research on the VMS has contributed significantly to the knowledge-base on rotorcraft performance, handling qualities, flight control, and guidance and displays. These contributions have directly benefited current rotorcraft programs and flight safety. The high fidelity motion system in the VMS was also used to research simulation fidelity. This research provided a fundamental understanding of pilot cueing modalities and their effect on simulation fidelity.

Aponso, Bimal Lalith↗

Characterizing the Vertical Profile of Aerosol Particle Extinction and Linear Depolarization over Southeast Asia and the Maritime Continent: The 2007-2009 View from CALIOP

Vertical profiles of 0.532 μm aerosol particle extinction coefficient and linear volume depolarization ratio are described for Southeast Asia and the Maritime Continent. Quality-screened and cloud-cleared Version 3.01 Level 2 NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) 5-km Aerosol Profile datasets are analyzed from 2007 to 2009. Numerical simulations from the U.S. Naval Aerosol Analysis and Predictive System (NAAPS), featuring two-dimensional variational assimilation of NASA Moderate Resolution Imaging Spectroradiometer and Multi-angle Imaging Spectro- Radiometer quality-assured datasets, combined with regional ground-based lidar measurements, are considered for assessing CALIOP retrieval performance, identifying bias, and evaluating regional representativeness. CALIOP retrievals of aerosol particle extinction coefficient and aerosol optical depth (AOD) are high over land and low over open waters relative to NAAPS (0.412/0.312 over land for all data points inclusive, 0.310/0.235 when the per bin average is used and each is treated as single data points; 0.102/0.151 and 0.086/0.124, respectively, over ocean). Regional means, however, are very similar (0.180/0.193 for all data points and 0.155/0.159 when averaged per normalized bin), as the two factors offset one another. The land/ocean offset is investigated, and discrepancies attributed to interpretation of particle composition and a-priori assignment of the extinction-to-backscatter ratio ("lidar ratio") necessary for retrieving the extinction coefficient from CALIOP signals. Over land, NAAPS indicates more dust present than CALIOP algorithms are identifying, indicating a likely assignment of a higher lidar ratio representative of more absorptive particles. NAAPS resolvesmore smoke overwater than identified with CALIOP, indicating likely usage of a lidar ratio characteristic of less absorptive particles to be applied that biases low AOD there. Over open waters except within the Bay of Bengal, aerosol particle scattering is largely capped below 1.5 km MSL, though ground-based lidar measurements at Singapore differ slightly from this finding. Significant aerosol particle presence over land is similarly capped near 3.0 km MSL over most regions. Particle presence at low levels regionally, except over India, is dominated by relatively non-depolarizing particles. Industrial haze, sea salt droplets and fresh smoke are thus most likely present.

linear volume depolarization ratio↗

Observing system simulations using synthetic radiances and atmospheric retrievals derived for the AMSU and HIRS in a mesoscale model

The paper addresses the concept of synthetic satellite imagery as a visualization and diagnostic tool for understanding satellite sensors of the future and to detail preliminary results on the quality of soundings from the current sensors. Preliminary results are presented on the quality of soundings from the combination of the High-Resolution Infrared Radiometer Sounder and the Advanced Microwave Sounding Unit. Results are also presented on the first Observing System Simulation Experiment using this data in a mesoscale numerical prediction model.

Diak, George R.↗

Enhancing the representation of water management in global hydrological models

Abstract. This study enhances an existing global hydrological model (GHM), Xanthos, by adding a new water management module that distinguishes between the operational characteristics of irrigation, hydropower, and flood control reservoirs. We remapped reservoirs in the Global Reservoir and Dam (GRanD) database to the 0.5∘ spatial resolution in Xanthos so that a single lumped reservoir exists per grid cell, which yielded 3790 large reservoirs. We implemented unique operation rules for each reservoir type, based on their primary purposes. In particular, hydropower reservoirs have been treated as flood control reservoirs in previous GHM studies, while here, we determined the operation rules for hydropower reservoirs via optimization that maximizes long-term hydropower production. We conducted global simulations using the enhanced Xanthos and validated monthly streamflow for 91 large river basins, where high-quality observed streamflow data were available. A total of 1878 (296 hydropower, 486 irrigation, and 1096 flood control and others) out of the 3790 reservoirs are located in the 91 basins and are part of our reported results. The Kling–Gupta efficiency (KGE) value (after adding the new water management) is ≥ 0.5 and ≥ 0.0 in 39 and 81 basins, respectively. After adding the new water management module, model performance improved for 75 out of 91 basins and worsened for only 7. To measure the relative difference between explicitly representing hydropower reservoirs and representing hydropower reservoirs as flood control reservoirs (as is commonly done in other GHMs), we use the normalized root mean square error (NRMSE) and the coefficient of determination (R2). Out of the 296 hydropower reservoirs, the NRMSE is > 0.25 (i.e., considering 0.25 to represent a moderate difference) for over 44 % of the 296 reservoirs when comparing both the simulated reservoir releases and storage time series between the two simulations. We suggest that correctly representing hydropower reservoirs in GHMs could have important implications for our understanding and management of freshwater resource challenges at regional-to-global scales. This enhanced global water management modeling framework will allow the analysis of future global reservoir development and management from a coupled human–earth system perspective.

13 HYDRO ENERGY↗

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Use of a Scale Model in the Design of Modifications to the NASA Glenn Icing Research Tunnel

Major modifications were made in 1999 to the 6- by 9-Foot (1.8- by 2.7-m) Icing Research tunnel (IRT) at the NASA Glenn Research Center, including replacement of its heat exchanger and associated ducts and turning vanes, and the addition of fan outlet guide vanes (OGV's). A one-tenth scale model of the IRT (designated as the SMIRT) was constructed with and without these modifications and tested to increase confidence in obtaining expected improvements in flow quality around the tunnel loop. The SMIRT is itself an aerodynamic test facility whose flow patterns without modifications have been shown to be accurate, scaled representations of those measured in the IRT prior to the 1999 upgrade program. In addition, tests in the SMIRT equipped with simulated OGV's indicated that these devices in the IRT might reduce flow distortions immediately downstream of the fan by two thirds. Flow quality parameters measured in the SMIRT were projected to the full-size modified IRT, and quantitative estimates of improvements in flow quality were given prior to construction. In this paper, the results of extensive flow quality studies conducted in the SMIRT are documented. Samples of these are then compared with equivalent measurements made in the full-scale IRT, both before and after its configuration was upgraded. Airspeed, turbulence intensity, and flow angularity distributions are presented for cross sections downstream of the drive fan, both upstream and downstream of the replacement flat heat exchanger, in the stilling chamber, in the test section, and in the wakes of the new comer turning vanes with their unique expanding and contracting designs. Lessons learned from these scale-model studies are discussed.

Canacci, Victor A.↗

Importance of inlet boundary conditions for numerical simulation of combustor flows

Fluid dynamic computer codes for the mathematical simulation of problems in gas turbine engine combustion systems are required as design and diagnostic tools. To eventually achieve a performance standard with these codes of more than qualitative accuracy it is desirable to use benchmark experiments for validation studies. Typical of the fluid dynamic computer codes being developed for combustor simulations is the TEACH (Teaching Elliptic Axisymmetric Characteristics Heuristically) solution procedure. It is difficult to find suitable experiments which satisfy the present definition of benchmark quality. For the majority of the available experiments there is a lack of information concerning the boundary conditions. A standard TEACH-type numerical technique is applied to a number of test-case experiments. It is found that numerical simulations of gas turbine combustor-relevant flows can be sensitive to the plane at which the calculations start and the spatial distributions of inlet quantities for swirling flows.

Sturgess, G. J.↗

The Influence of Observation Errors on Analysis Error and Forecast Skill Investigated with an Observing System Simulation Experiment

The Global Modeling and Assimilation Office (GMAO) observing system simulation experiment (OSSE) framework is used to explore the response of analysis error and forecast skill to observation quality. In an OSSE, synthetic observations may be created that have much smaller error than real observations, and precisely quantified error may be applied to these synthetic observations. Three experiments are performed in which synthetic observations with magnitudes of applied observation error that vary from zero to twice the estimated realistic error are ingested into the Goddard Earth Observing System Model (GEOS-5) with Gridpoint Statistical Interpolation (GSI) data assimilation for a one-month period representing July. The analysis increment and observation innovation are strongly impacted by observation error, with much larger variances for increased observation error. The analysis quality is degraded by increased observation error, but the change in root-mean-square error of the analysis state is small relative to the total analysis error. Surprisingly, in the 120 hour forecast increased observation error only yields a slight decline in forecast skill in the extratropics, and no discernable degradation of forecast skill in the tropics.

Observation Error↗

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↗