Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shallow water equations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

A Proposal for an Integrated Geophysical Strategy to "Follow the Water" on Mars

The search for subsurface water has become a primary focus of Mars exploration. Its abundance and distribution (both as ground ice and groundwater) have important implications for understanding the geologic, hydrologic, and climatic evolution of the planet; the potential origin and continued survival of life; and the accessibility of a critical in situ resource for sustaining future human explorers. For these reasons, a principal goal of the Mars science, astrobiology, and the HEDS programs is to determine the 3-D distribution and state of subsurface H2O, at a resolution sufficient to permit reaching any desired volatile target by drilling. The three targets most often discussed are: groundwater, massive deposits of near-surface ground ice (associated with the ponded discharge of the outflow channels or the relic of a former ocean), and ice-saturated frozen ground. Based on the present best estimates of mean annual surface temperature, crustal thermal conductivity, geothermal heat flow, and groundwater freezing temperature, the mean thickness of frozen ground on Mars is expected to vary from approx. = 2.5 - 5 km at the equator to approx. = 6.5 - 13 km at the poles. However, natural variations in both crustal heat flow and thermal conductivity are likely to result in significant local departures from these predicted values. The recent discovery of "young" fluvial-like features, emanating from the slopes of local scarps, raises the possibility that liquid water may also exist episodically at shallow (approx. = 100 - 500 m) depth; however, the true nature and absolute age of these features remains highly uncertain. Although the belief that Mars is water-rich is supported by a wide variety geologic evidence, our ignorance about the heterogeneous nature and thermal evolution of the planet's crust effectively precludes geomorphic or theoretical attempts to quantitatively assess the current geographic and subsurface vertical distribution of ground ice and groundwater . For this reason, any exploration activity (such as drilling) whose success is contingent on the presence of subsurface water, must be preceded by a comprehensive high-resolution geophysical survey capable of assessing whether local reservoirs of water and ice actually exist. Terrestrial experience has demonstrated that the accurate identification of such targets is likely to require the application of multiple geophysical techniques. In this abstract we propose an integrated strategy for the geophysical exploration of Mars that we believe represents the fastest, most cost-effect, and technically capable approach to identifying the state and distribution of subsurface water. Additional information is contained in the original extended abstract.

Clifford, S. M.↗
NASA NTRS1991

Simplification of the Kalman filter for meteorological data assimilation

The paper proposes a new statistical method of data assimilation that is based on a simplification of the Kalman filter equations. The forecast error covariance evolution is approximated simply by advecting the mass-error covariance field, deriving the remaining covariances geostrophically, and accounting for external model-error forcing only at the end of each forecast cycle. This greatly reduces the cost of computation of the forecast error covariance. In simulations with a linear, one-dimensional shallow-water model and data generated artificially, the performance of the simplified filter is compared with that of the Kalman filter and the optimal interpolation (OI) method. The simplified filter produces analyses that are nearly optimal, and represents a significant improvement over OI.

Dee, Dick P.↗
NASA NTRS2002

Technical Challenges of Drilling on Mars

In the last year, NASA's Mars science advisory committee (MEPAG: Mars Exploration Payload Advisory Group) has formally recommended that deep drilling be undertaken as a priority investigation to meet astrobiology and geology goals. This proposed new dimension in Mars exploration has come about for several reasons. Firstly, geophysical models of the martian subsurface environment indicate that we may well find liquid water (in the form of brines) under ground-ice at depths of several kilometers near the equator. On Earth we invariably find life forms associated with any environmental niche that supports liquid water. New data from the Mars Global Surveyor have shown that the most recent volcanism on Mars is very young so we cannot rule out contemporary volcanism -- in which case subsurface temperatures consistent with having water in its liquid phase may be found at relatively shallow depths. Secondly, in recent decades we have learned to our surprise that the Earth's subsurface (microbial) biosphere extends to depths of many kilometers and this discovery provides the basis for planning to explore the martian subsurface in search of ancient or even extant microbial life forms. We know (from Viking measurements) that all the biogenic elements (C, H, O, N, P, S) are available on Mars. What we therefore hope to learn is whether or not the evolution of life is inevitable given the necessary ingredients and, by implication, whether the Universe may be teeming with life. The feasibility of drilling deep into the surface of Mars has been the subject of increasing attention within NASA (and more recently among some of its international partners) for several years and this led to a broad-based feasibility study carried out by the Los Alamos National Laboratory and, subsequently, to the development of several hardware prototypes. This paper is intended to provide a general survey of that activity.

Briggs, Geoffrey↗
NASA NTRS1973

Wind-driven currents in a shallow lake or sea.

A technique is given for obtaining exact analytical solutions for wind-driven currents in a shallow lake or sea with a large class of bottom topographies which may or may not contain an island. Welander's second-order partial differential equation for the lake's surface displacement is used for determining the lake circulation. The technique, together with some results based on conformal mapping, is used to reduce the problem of finding solutions to Welander's equation for a closed body of water to the problem of solving an ordinary differential equation. The method is applied to an elliptically shaped lake and a circular lake containing an eccentrically located circular island.

Goldstein, M.↗
NASA NTRS1983

Spatial variation of corn canopy temperature as dependent upon soil texture and crop rooting characteristics

A soil plant atmosphere model for corn (Zea mays L.) together with the scaling theory for soil hydraulic heterogeneity are used to study the sensitivity of spatial variation of canopy temperature to field averaged soil texture and crop rooting characteristics. The soil plant atmosphere model explicitly solves a continuity equation for water flux resulting from root water uptake, changes in plant water storage and transpirational flux. Dynamical equations for root zone soil water potential and the plant water storage models the progressive drying of soil, and day time dehydration and night time hydration of the crop. The statistic of scaling parameter which describes the spatial variation of soil hydraulic conductivity and matric potential is assumed to be independent of soil texture class. The field averaged soil hydraulic characteristics are chosen to be representative of loamy sand and clay loam soils. Two rooting characteristics are chosen, one shallow and the other deep rooted. The simulation shows that the range of canopy temperatures in the clayey soil is less than 1K, but for the sandy soil the range is about 2.5 and 5.0 K, respectively, for the shallow and deep rooted crops.

Choudhury, B. J.↗
DOE OSTI2023

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗
DOE OSTI2022

Westward-Propagating Moisture Mode Over the Tropical Western Hemisphere

A westward-propagating Rossby-like wave signal is found to explain a large fraction of the intraseasonal variance in cloud brightness over the Western Hemisphere. A series of diagnostic criteria suggest that this wave is a moisture mode: its moisture anomalies dominate the distribution of moist static energy (MSE) and are in phase with the precipitation anomalies; and the thermodynamic equation obeys the weak temperature gradient approximation. Further, the wave propagates westward due to zonal moisture advection by the mean flow and is maintained by radiative heating and meridional moisture advection. These properties compare favorably with the westward propagating Rossby mode in an equatorial beta-plane model with prognostic moisture, mean meridional moisture gradient, and mean zonal wind. These results underscore the importance of water vapor in the dynamics of slowly evolving tropical systems, and the limitations of dry shallow water theory that rely on a “reduced equivalent depth” to represent moist dynamics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Abrupt Changes at the Permian/Triassic Boundary: Tempo of Events from High-Resolution Cyclostratigraphy

The Permian/Triassic (P/Tr) boundary (251.4 +/- 3 Myr) is marked by the most severe mass extinction in the geologic record. Recently, precise absolute dating has bracketed the marine extinctions and associated carbon-isotope anomaly within less than 1 Myr. We improve this resolution through high-resolution stratigraphy across the P/Tr boundary in the 331-m Gartnerkofel-1 core and nearby Reppwand outcrop section (Carnic Alps, Austria) utilizing FFT and wavelet timeseries analyses of cyclic components in down-hole core logs of density and natural gamma-ray intensity, and carbon-isotopic ratios of bulk samples. The wavelet analysis indicates continuity of deposition across the P/Tr boundary interval, and the timeseries analyses show evidence for persistent cycles in the ratio of approximately 40: 10: 4.7: 2.3 meters, correlated with Milankovitch-band orbital cycles of approximately 412: 100: 40: 20 kyr (eccentricity 1 and 2, obliquity, and precession), and giving a consistent average sedimentation rate of approximately 10 cm/1,000 yr. Milankovitch periods in delta C-13 and density in these shallow-water carbonates were most likely the result of climatically induced oscillations of sea level and climate, coupled with changes in ocean circulation and productivity, that affected sedimentation. Fluctuations in gamma radiation reflect varying input of clay minerals and the presence of shaly interbeds. Throughout the P/Tr boundary interval in the core, the 100,000-year eccentricity cycle seems to be dominant. Weaker obliquity and precession cycles are in line with the location of the Austrian section in the latest Permian, close to the Equator in the western bight of the Tethys, where obliquity and precessional effects on seasonal contrast might be subdued. Using the improved resolution provided by cycle analysis in the GK-1 core, we find that the dramatic change in the faunal record that marks the P/Tr boundary takes place over less than 6m, or less than 60,000 years. In the nearby Reppwand outcrop section, the same faunal changes occurs over only 0.8 m or about 8,000 years, close to the limit of time-resolution induced by bioturbation and reworking in these sediments. The sharp negative global carbon-isotope shift took place within less than or equal to 40,000 yr, and the isotope excursions persisted for approximately 480,000 yr into the Early Triassic. The results indicate that the severe marine faunal event that marks the P/Tr boundary was very sudden, perhaps less than the resolution window in the GK-1 core, and suggest a catastrophic cause. The wavelet-analysis approach to high-resolution cyclostratigraphy can be applied to other P/Tr boundary sections, and when combined with precise absolute dating and magnetostratigraphic methods promises a significant increase in resolution in determining the correlation and tempo of the end-Permian extinctions and related events worldwide.

Rampino, M. R.↗
DOE OSTI2022

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2017

Extratropical Weather Systems on Mars: Radiatively-Active Water Ice Effects

Extratropical, large-scale weather disturbances, namely transient, synoptic-period,baroclinic barotropic eddies - or - low- (high-) pressure cyclones (anticyclones), are components fundamental to global circulation patterns for rapidly rotating, differentially heated, shallow atmospheres such as Earth and Mars. Such "wave-like" disturbances that arise via (geophysical) fluid shear instability develop, mature and decay, and travel west-to-east in the middle and high latitudes within terrestrial-like planetary atmospheres. These disturbances serve as critical agents in the transport of heat and momentum between low and high latitudes of the planet. Moreover, they transport trace species within the atmosphere (e.g., water vapor/ice, other aerosols (dust), chemical species, etc). Between early autumn through early spring, middle and high latitudes on Mars exhibit strong equator-to-pole mean temperature contrasts (i.e., "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that such strong baroclinicity supports vigorous, large-scale eastward traveling weather systems [Banfield et al., 2004; Barnes et al., 1993]. A good example of traveling weather systems, frontal wave activity and sequestered dust activity from MGS/MOC image analyses is provided in Figure 1 (cf. Wang et al. [2005]). Utilizing an upgraded and evolving version of the NASA Ames Research Center (ARC) Mars global climate model, investigated here are key dynamical and physical aspects of simulated northern hemisphere (NH) large-scale extratropica lweather systems,with and without radiatively-active water ice clouds. Mars Climate Model:

cyclones↗
NASA NTRS2016

Traveling Weather Disturbances in Mars Southern Extratropics: Sway of the Great Impact Basins

As on Earth, between late autumn and early spring on Mars middle and high latitudes within its atmosphere support strong mean thermal contrasts between the equator and poles (i.e. "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that this strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e. transient synoptic-period waves). Within a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, such large-scale, extratropical weather disturbances are critical components of the global circulation. These wave-like disturbances act as agents in the transport of heat and momentum, and moreover generalized tracer quantities (e.g., atmospheric dust, water vapor and water-ice clouds) between low and high latitudes of the planet. The character of large-scale, traveling extratropical synoptic-period disturbances in Mars' southern hemisphere during late winter through early spring is investigated using a high-resolution Mars global climate model (Mars GCM). This global circulation model imposes interactively lifted (and radiatively active) dust based on a threshold value of the instantaneous surface stress. Compared to observations, the model exhibits a reasonable "dust cycle" (i.e. globally averaged, a more dusty atmosphere during southern spring and summer occurs). In contrast to their northern-hemisphere counterparts, southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are far less intense synoptically. Influences of the zonally asymmetric (i.e. east-west varying) topography on southern large-scale weather disturbances are examined. Simulations that adapt Mars' full topography compared to simulations that utilize synthetic topographies emulating essential large-scale features of the southern middle latitudes indicate that Mars' transient barotropic/baroclinic eddies are significantly influenced by the great impact basins of this hemisphere (e.g., Argyre and Hellas). In addition, the occurrence of a southern storm zone in late winter and early spring is keyed particularly to the western hemisphere via orographic influences arising from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate fundamental differences amongst such simulations and these are described.

Weather disturbances↗
NASA NTRS1996

Assimilation of Geosat Altimetric Data in a Nonlinear Shallow-Water Model of the Indian Ocean by Adjoint Approach. Part II: Some Validation and Interpretation of the Assimilated Results: Some Validation and Interpretation of the Assimilated Results - Part 2

This paper examines the results of assimilating Geosat sea level variations relative to the November 1986-November 1988 mean reference, in a nonlinear reduced-gravity model of the Indian Ocean, Data have been assimilated during one year starting in November 1986 with the objective of optimizing the initial conditions and the yearly averaged reference surface. The thermocline slope simulated by the model with or without assimilation is validated by comparison with the signal, which can be derived from expandable bathythermograph measurements performed in the Indian Ocean at that time. The topography simulated with assimilation on November 1986 is in very good agreement with the hydrographic data. The slopes corresponding to the South Equatorial Current and to the South Equatorial Countercurrent are better reproduced with assimilation than without during the first nine months. The whole circulation of the cyclonic gyre south of the equator is then strongly intensified by assimilation. Another assimilation experiment is run over the following year starting in November 1987. The difference between the two yearly mean surfaces simulated with assimilation is in excellent agreement with Geosat. In the southeastern Indian Ocean, the correction to the yearly mean dynamic topography due to assimilation over the second year is negatively correlated to the one the year before. This correction is also in agreement with hydrographic data. It is likely that the signal corrected by assimilation is not only due to wind error, because simulations driven by various wind forcings present the same features over the two years. Model simulations run with a prescribed throughflow transport anomaly indicate that assimilation is rather correcting in the interior of the model domain for inadequate boundary conditions with the Pacific.

Greiner, Eric↗
NASA NTRS1995

The development of concentration gradients in a suspension of chemotactic bacteria

When a suspension of bacterial cells of the species Bacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e., they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development when N and C depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

NASA Discipline Cell Biology↗
NASA NTRS1994

Suboptimal schemes for atmospheric data assimilation based on the Kalman filter

This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Todling, Ricardo↗
NASA NTRS2014

Modeling Mars Cyclogenesis and Frontal Waves: Seasonal Variations and Implications on Dust Activity

Between late autumn through early spring,middle and high latitudes onMars exhibit strong equator-to-polemean temperature contrasts (i.e., "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that such strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e., transient synoptic period waves) [1, 2]. For a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, these large-scale, extratropical weather disturbances are critical components of the global circulation. The wave-like disturbances serve as agents in the transport of heat and momentum between low and high latitudes of the planet. Through cyclonic/anticyclonic winds, intense shear deformations, contractions-dilatations in temperature and density, and sharp perturbations amongst atmospheric tracers (i.e., dust, volatiles (e.g., water vapor) and condensates (e.g., water-ice cloud particles)), Mars' extratropical weather systems have significant sub-synoptic scale ramifications by supporting atmospheric frontal waves (Fig. 1).

Mars Global Surveyor (MGS)↗
NASA NTRS2004

Habitability of the Shallow Subsurface on Mars: Clues from the Meteorites

The properties that define habitability are commonly understood to include the following: - Presence of water. - Temperature range allowing some or all of the water to be liquid. - A suitable physical volume or space permitting metabolism and growth. - Presence of organic compounds or the building blocks to make them. - Presence of an energy source suitable for utilization by living organisms. Interpretations of Mars Viking, Surveyor, and Odyssey orbital images have built a strong case that Mars had surface water during its past geological history. Neutron spectrometer data from Mars Odyssey show that poleward of about 60 degrees North and 60 degrees south, significant hydrogen, likely as ice or permafrost, is present in at least the upper meter or so of the Martian regolith and crust and that similar high hydrogen areas exist, even near the equator. Here we present a summary of independent data from the Mars meteorites showing that liquid water was present for at least some of the time in the upper few meters or tens of meters as early as 3.9 billion years (Ga), and was present at intervals and at various locations throughout most of Mars history.

David S McKay↗
NASA NTRS2011

Improving Simulated Soil Moisture Fields Through Assimilation of AMSR-E Soil Moisture Retrievals with an Ensemble Kalman Filter and a Mass Conservation Constraint

Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of soil moisture in a footprint area, and can be used to reduce model bias (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce model bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate the actual value of Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals to improve the mean of simulated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed.

Li, Bailing↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗