Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “persistence”

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,207 records · Page 67

Sensitivity of deep convection and cross-tropopause water transport to microphysical parameterizations in WRF

Deep convective storms can overshoot the tropopause, thus altering the composition of the stratosphere by vertically transporting tropospheric air. The transport of water vapor and ice particles into a sub-saturated environment can hydrate the stratosphere, with implications for radiative forcing and ozone chemistry. Cloud-resolved models, if employed at high spatial resolutions, are used to probe process-level questions about cross-tropopause deep convective hydration and its controls. There is considerable diversity in model representations of processes associated with water transport and transformation, and the choice of a microphysics scheme affects model skill in simulating deep convective events. This motivates our evaluation of state-of-the-art, as well as widely used standard schemes, in a high spatial- and temporal-resolution framework. Six bulk microphysics schemes were employed in a WRF-LES setup, initialized with a sounding profile representative of a tropopause-overshooting storm. Here, we used an idealized framework to isolate the effect of microphysics on the dominant processes that control the reach of deep convection and stratospheric hydration. All schemes produced the highest reaching updrafts 8-12 hours into the simulation but the strength and persistence of updrafts varied across the schemes; maximum storm heights ranged 9.1-12.6 km across the schemes. Varying microphysics produced large differences in the vertical extent and horizontal aggregation of convection, and an order of magnitude spread in above-tropopause water vapor concentrations.

54 ENVIRONMENTAL SCIENCES↗

The Shortwave Radiative Flux Response to an Injection of Sea Salt Aerosols in the Gulf of Mexico

We report marine cloud brightening (MCB) has been proposed as a potential means of geoengineering the climate, temporarily providing cooling to offset some of the effects of climate change. Marine sky brightening (MSB), involving the direct scattering of sunlight from sea salt injection into the marine boundary layer, has been proposed as an additional geoengineering method that could work in areas that are not regularly cloudy. Here, we use a regional atmospheric model to simulate MCB and MSB over the Gulf of Mexico and nearby land, a highly populated and economically important region that is not characterized by persistent marine stratocumulus cloud cover. Injection of sea salt in the Aitken mode from a region in the central Gulf of Mexico equivalent to 10.8 Tg yr -1 produces an upwards 8.4 W m -2 radiative flux change across the region at the top of the atmosphere, largely due to cloud property changes. Comparatively, a similar mass injection in the accumulation mode produces a 3.1 W m -2 radiative flux change driven primarily by direct scattering. Injection of even larger particles produces a much smaller radiative flux change. Shortwave flux changes due to clouds are largely driven by an increase in cloud droplet number concentration and an increase in cloud liquid water path (each contributing about 45% to the flux change), with a much lower contribution from cloud fraction changes (10%).

54 ENVIRONMENTAL SCIENCES↗

Cloud Phase Simulation at High Latitudes in EAMv2: Evaluation using CALIPSO Observations and Comparison with EAMv1

This study performs a comprehensive evaluation of the simulated cloud phase in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) atmosphere model version 2 (EAMv2) and version 1 (EAMv1). Enabled by the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) simulator, EAMv2 and EAMv1 predicted cloud phase is compared against the GCM-Oriented CALIPSO Cloud Product (CALIPSO-GOCCP) at high latitudes where mixed-phase clouds are prevalent. Our results indicate that the underestimation of cloud ice in simulated high-latitude mixed-phase clouds in EAMv1 has been significantly reduced in EAMv2. The increased ice clouds in the Arctic mainly result from the modification on the WBF (Wegner-Bergeron-Findeisen) process in EAMv2. The impact of the modified WBF process is moderately compensated by the low limit of cloud droplet number concentration (CDNC) in cloud microphysics and the new dCAPE_ULL trigger used in deep convection in EAMv2. Moreover, it is found that the new trigger largely contributes to the better cloud phase simulation over the Norwegian Sea and Barents Sea in the Arctic and the Southern Ocean where large errors are found in EAMv1. However, errors in simulated cloud phase in EAMv1, such as the overestimation of supercooled liquid clouds near the surface in both hemispheres and the underestimation of ice clouds over Antarctica, persist in EAMv2. In conclusion, this study highlights the impact of deep convection parameterizations, which has received little attention, on high-latitude mixed-phase clouds, and the importance of continuous improvement of cloud microphysics in climate models for accurately representing mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗

Future Projections of the El Niño—Southern Oscillation and Tropical Pacific Mean State in CMIP6

Abstract The El Niño—Southern Oscillation (ENSO) is an important mode of tropical Pacific atmosphere‐ocean variability that drives teleconnections with weather and climate globally. However, prior studies using state‐of‐the‐art climate models lack consensus regarding future ENSO projections and are often impacted by tropical Pacific sea‐surface temperature (SST) biases. We used 173 simulations from 29 climate models participating in the Coupled Model Intercomparison Project, version 6 (CMIP6) to analyze model biases and future ENSO projections. We analyzed two ENSO indices, namely the ENSO Longitude Index (ELI), which measures zonal shifts in tropical Pacific deep convection and accounts for changes in background SST, and the Niño 3.4 index, which measures SST anomalies in the central‐eastern equatorial Pacific. We found that the warm eastern tropical‐subtropical Pacific SST bias typical of previous generations of climate models persists into many of the CMIP6 models. Future projections of ENSO shift toward more El Niño‐like conditions based on ELI in 48% of simulations and 55% of models, in association with a future weakening of the zonal equatorial Pacific SST gradient. On the other hand, none of the models project a significant shift toward La Niña‐like conditions. The standard deviation of the Niño 3.4 index indicates a lack of consensus on whether an increase or decrease in ENSO variability is expected in the future. Finally, we found a possible relationship between historical SST and low‐level cloud cover biases in the ENSO region and future changes in ELI; however, this result may be impacted by limitations in data availability.

54 ENVIRONMENTAL SCIENCES↗

Vertical Structure of Clouds and Precipitation During Arctic Cold–Air Outbreaks and Warm–Air Intrusions: Observations From COMBLE

The Arctic is marked by deep intrusions of warm, moist air, alternating with outbreaks of cold air down to lower latitudes. The typical vertical structure of clouds and precipitation during these two synoptic weather extremes is examined at a coastal site at 69°N in Norway. The Norwegian Sea is a corridor for warm-air intrusions (WAIs) and frequently witnesses cold-air outbreaks (CAOs). This study uses data from profiling radar, lidar, and microwave radiometer, radiosondes and other probes that were collected during the CAOs in the Marine Boundary Layer Experiment (COMBLE) between 1 December 2019 and 31 May 2020. Marine CAOs are defined in terms of thermal instability relative to the sea surface temperature, and WAIs in terms of equivalent potential temperature stratification between the surface and 850 hPa. Cloud structures in CAOs are convective, driven by strong surface heat fluxes over a long fetch of open water, with cloud tops rarely exceeding 6 km above sea level. The mostly open-cellular convection produces intermittent moderately-heavy precipitation at the observational site, notwithstanding the low precipitable water vapor (PWV). In contrast, WAIs are marked by high values of PWV and integrated vapor transport. Here, WAI clouds are synoptically driven, stratiform, with cloud tops often exceeding 5 km, sometimes layered, and generally producing persistent precipitation that can be heavier than in CAOs.

54 ENVIRONMENTAL SCIENCES↗

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗

The Impact of Freeze‐Thaw History on Soil Carbon Response to Experimental Freeze‐Thaw Cycles

Abstract Freeze‐thaw is a disturbance process in cold regions where permafrost soils are becoming vulnerable to temperature fluctuations above 0°C. Freeze‐thaw alters soil physical and biogeochemical properties with implications for carbon persistence and emissions in Arctic landscapes. We examined whether different freeze‐thaw histories in two soil systems led to contrasting biogeochemical responses under a laboratory‐controlled freeze‐thaw incubation. We investigated controls on carbon composition through Fourier‐transform ion cyclotron resonance mass spectrometry (FT‐ICR‐MS) to identify nominal carbon oxidation states and relative abundances of aliphatic‐type carbon molecules in both surface and subsurface soils. Soil cores (∼60 cm‐depth) were sampled from two sites in Alaskan permafrost landscapes with different in situ freeze‐thaw characteristics: Healy (>40 freeze‐thaw cycles annually) and Toolik (<15 freeze‐thaw cycles annually). FT‐ICR‐MS was coupled with in situ temperature data and soil properties (i.e., soil texture, mineralogy) to assess (a) differences in soil organic matter composition associated with previous freeze‐thaw history and (b) sensitivity to experimental freeze‐thaw in the extracted cores. Control (freeze‐only) samples showed greater carbon oxidation in Healy soils compared with Toolik, even in lower mineral horizons where freeze‐thaw history was comparable across both sites. Healy showed the most loss of carbon compounds following experimental freeze‐thaw in the lower mineral depths, including a decrease in aliphatics. Toolik soils responded more slowly to freeze‐thaw as shown by intermediary carbon oxidation distributed across multiple carbon compound classes. Variations in the response of permafrost carbon chemistry to freeze‐thaw is an important factor for predicting changes in soil function as permafrost thaws in high northern latitudes.

Rooney, Erin C.↗

Characterizing Performance of Freshwater Wetland Methane Models Across Time Scales at FLUXNET-CH 4 Sites Using Wavelet Analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

58 GEOSCIENCES↗

Limitations of Separate Cloud and Rain Categories in Parameterizing Collision-Coalescence for Bulk Microphysics Schemes

Warm rain collision-coalescence has been persistently difficult to parameterize in bulk microphysics schemes. We use a flexible bulk microphysics scheme with bin scheme process parameterizations, called AMP, to investigate reasons for the difficulty. AMP is configured in a variety of ways to mimic bulk schemes and is compared to simulations with the bin scheme upon which AMP is built. We find that an important limitation in traditional bulk schemes is the use of separate cloud and rain categories. When the drop size distribution is instead represented by a continuous distribution, the simulation of cloud-to-rain conversion is substantially improved. We also find large sensitivity to the threshold size to distinguish cloud and rain in traditional schemes; substantial improvement is found by decreasing the threshold from 40 to 25 μm. Neither the use of an assumed functional form for the size distribution nor the choice of predicted distribution moments has a large impact on the ability of AMP to simulate rain production. When predicting four total moments of the liquid drop size distribution, either with a traditional two-category, two-moment scheme with a reduced size threshold, or a four-moment single-category scheme, errors in the evolution of mass and the cloud size distribution are similar, but the single-category scheme has a substantially better representation of the rain size distribution. Optimal moment combinations for the single-category approach are investigated and appear to be linked more to the information content they provide for constraining the size distributions than to their correlation with collision-coalescence rates.

54 ENVIRONMENTAL SCIENCES↗

Leakage From Coexisting Geologic Forcing and Injection‐Induced Pressurization: A Semi‐Analytical Solution for Multilayered Aquifers With Multiple Wells

Abstract Abnormal fluid pressures (above or below hydrostatic pressure) can develop and persist in sedimentary basins. The common occurrence of abnormal pressures may cause challenges for project permitting of geological carbon sequestration (GCS), particularly in reservoirs with pre‐injection overpressure. The leaky wells that may exist in some sedimentary basins can provide flow paths between deep brine aquifers and shallower freshwater aquifers. Pre‐injection relative overpressures can cause brine leakage through leaky wells even before any injection occurs. The tendency for flow through leaky wells is coupled with the process of pressure dissipation that occurs through aquitards. Specifically, with non‐zero permeability, aquitards can dissipate pressure over large areal extents and thereby reduce leakage rates through leaky wells. This study presents development of a semi‐analytical solution for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. The geologic forcing that causes abnormal pressures in the multilayered system can coexist with any number of injection, extraction, and leaky wells that also affect fluid pressure. The semi‐analytical model is applied to explore how leakage through leaky wells varies as functions of pressurization rate, along with aquitard and leaky well properties in an overpressured multilayered system. The results show that although injection‐induced pressures can dissipate rapidly through suitably permeable aquitards, coexisting geologic forcing may create sustained rates of brine leakage into freshwater aquifers through leaky wells. In GCS, a very low‐permeability aquitard with high capillary entry pressure to free‐phase CO 2 is desired to serve as the caprock to prevent leakage of CO 2 from the storage reservoir. Nevertheless, the results from this study show that the brine leakage impact to shallow freshwater aquifers through leaky wells might substantially decrease with increasing aquitard permeability values, as long as small aquitard permeability and high capillary entry pressure serve to prevent CO 2 leakage.

Cihan, A.↗

Skillful Multi‐Month Predictions of Ecosystem Stressors in the Surface and Subsurface Ocean

Abstract Anthropogenic carbon emissions and associated climate change are driving rapid warming, acidification, and deoxygenation in the ocean, which increasingly stress marine ecosystems. On top of long‐term trends, short term variability of marine stressors can have major implications for marine ecosystems and their management. As such, there is a growing need for predictions of marine ecosystem stressors on monthly, seasonal, and multi‐month timescales. Previous studies have demonstrated the ability to make reliable predictions of the surface ocean physical and biogeochemical state months to years in advance, but few studies have investigated forecast skill of multiple stressors simultaneously or assessed the forecast skill below the surface. Here, we use the Community Earth System Model (CESM) Seasonal to Multiyear Large Ensemble (SMYLE) along with novel observation‐based biogeochemical and physical products to quantify the predictive skill of dissolved inorganic carbon (DIC), dissolved oxygen, and temperature in the surface and subsurface ocean. CESM SMYLE demonstrates high physical and biogeochemical predictive skill multiple months in advance in key oceanic regions and frequently outperforms persistence forecasts. We find up to 10 months of skillful forecasts, with particularly high skill in the Northeast Pacific (Gulf of Alaska and California Current Large Marine Ecosystems) for temperature, surface DIC, and subsurface oxygen. Our findings suggest that dynamical marine ecosystem prediction could support actionable advice for decision making.

Environmental Sciences & Ecology↗

Linking Dissolved Organic Matter Composition to Landscape Properties in Wetlands Across the United States of America

Abstract Wetlands are integral to the global carbon cycle, serving as both a source and a sink for organic carbon. Their potential for carbon storage will likely change in the coming decades in response to higher temperatures and variable precipitation patterns. We characterized the dissolved organic carbon (DOC) and dissolved organic matter (DOM) composition from 12 different wetland sites across the USA spanning gradients in climate, landcover, sampling depth, and hydroperiod for comparison to DOM in other inland waters. Using absorption spectroscopy, parallel factor analysis modeling, and ultra‐high resolution mass spectroscopy, we identified differences in DOM sourcing and processing by geographic site. Wetland DOM composition was driven primarily by differences in landcover where forested sites contained greater aromatic and oxygenated DOM content compared to grassland/herbaceous sites which were more aliphatic and enriched in N and S molecular formulae. Furthermore, surface and porewater DOM was also influenced by properties such as soil type, organic matter content, and precipitation. Surface water DOM was relatively enriched in oxygenated higher molecular weight formulae representing HUP High O/C compounds than porewaters, whose DOM composition suggests abiotic sulfurization from dissolved inorganic sulfide. Finally, we identified a group of persistent molecular formulae (3,489) present across all sites and sampling depths (i.e., the signature of wetland DOM) that are likely important for riverine‐to‐coastal DOM transport. As anthropogenic disturbances continue to impact temperate wetlands, this study highlights drivers of DOM composition fundamental for understanding how wetland organic carbon will change, and thus its role in biogeochemical cycling.

Environmental Sciences & Ecology↗

Small-Magnitude Seismic Swarms in Central Utah (US): Interactions of Regional Tectonics, Local Structures and Hydrothermal Systems

Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.

58 GEOSCIENCES↗

Impacts of a Cascadia Subduction Zone Earthquake on Water Levels and Wetlands of the Lower Columbia River and Estuary

Abstract Subsidence after a subduction zone earthquake can cause major changes in estuarine bathymetry. Here, we quantify the impacts of earthquake‐induced subsidence on hydrodynamics and habitat distributions in a major system, the lower Columbia River Estuary, using a hydrodynamic and habitat model. Model results indicate that coseismic subsidence increases tidal range, with the smallest changes at the coast and a maximum increase of ∼10% in a region of topographic convergence. All modeled scenarios reduce intertidal habitat by 24%–25% and shifts ∼93% of estuarine wetlands to lower‐elevation habitat bands. Incorporating dynamic effects of tidal change from subsidence yields higher estimates of remaining habitat by multiples of 0–3.7, dependent on the habitat type. The persistent tidal change and chronic habitat disturbance after an earthquake poses strong challenges for estuarine management and wetland restoration planning, particularly when coupled with future sea‐level rise effects.

54 ENVIRONMENTAL SCIENCES↗

Drivers of Dry Day Sensitivity to Increased CO 2

Persistent precipitation deficits are among the most impactful consequences of global warming. Here we focus on changes in the annual number of dry days (NDD) and in the annual maximum length of dry spells due to a quadrupling of atmospheric CO 2 . We use atmosphere-only simulations to decompose the projected changes into additive contributions. A fast adjustment leads to a global increase in NDD despite notable regional exceptions (e.g., South Asia and Sahel). The effect of the uniform component of the surface ocean warming is model-dependent but shapes the regional distribution of the NDD response in each model. Finally, the ocean warming pattern also contributes to large uncertainties, likely through contrasting changes in large-scale circulation. Our results thus highlight the complexity of the NDD response, with policy-relevant practical implications for mitigation and adaptation strategies.

54 ENVIRONMENTAL SCIENCES↗

Observations of Ocean Surface Wave Attenuation in Sea Ice Using Seafloor Cables

Abstract The attenuation of ocean surface waves during seasonal ice cover is an important control on the evolution of Arctic coastlines. The spatial and temporal variations in this process have been challenging to resolve with conventional sampling using sparse arrays of moorings or buoys. We demonstrate a novel method for persistent observation of wave‐ice interactions using distributed acoustic sensing (DAS) along existing seafloor fiber optic telecommunications cables. DAS measurements span a 36‐km cross‐shore cable on the Beaufort Shelf from Oliktok Point, Alaska. DAS optical sensing of fiber strain‐rate provides a proxy for seafloor pressure, which we calibrate with wave buoy measurements during the ice‐free season (August 2022). We apply this calibration during the ice formation season (November 2021) to obtain unprecedented resolution of variable wave attenuation rates in new, partial ice cover. The location and strength of wave attenuation serve as proxies for ice coverage and thickness, especially during rapidly evolving events.

Smith, Madison M.↗

Physical Insights From the Multidecadal Prediction of North Atlantic Sea Surface Temperature Variability Using Explainable Neural Networks

Abstract North Atlantic sea surface temperatures (NASST), particularly in the subpolar region, are among the most predictable in the world's oceans. However, the relative importance of atmospheric and oceanic controls on their variability at multidecadal timescales remain uncertain. Neural networks (NNs) are trained to examine the relative importance of oceanic and atmospheric predictors in predicting the NASST state in the Community Earth System Model 1 (CESM1). In the presence of external forcings, oceanic predictors outperform atmospheric predictors, persistence, and random chance baselines out to 25‐year leadtimes. Layer‐wise relevance propagation is used to unveil the sources of predictability, and reveal that NNs consistently rely upon the Gulf Stream‐North Atlantic Current region for accurate predictions. Additionally, CESM1‐trained NNs successfully predict the phasing of multidecadal variability in an observational data set, suggesting consistency in physical processes driving NASST variability between CESM1 and observations.

Geology↗

Scalar Flux Profiles in the Unstable Atmospheric Surface Layer Under the Influence of Large Eddies: Implications for Eddy Covariance Flux Measurements and the Non‐Closure Problem

How convective boundary-layer (CBL) processes modify fluxes of sensible (SH) and latent (LH) heat and CO 2 (F c ) in the atmospheric surface layer (ASL) remains a recalcitrant problem. Here, large eddy simulations for the CBL show that while SH in the ASL decreases linearly with height regardless of soil moisture conditions, LH and F c decrease linearly with height over wet soils but increase with height over dry soils. This varying flux divergence/convergence is regulated by changes in asymmetric flux transport between top-down and bottom-up processes. Such flux divergence and convergence indicate that turbulent fluxes measured in the ASL underestimate and overestimate the “true” surface interfacial fluxes, respectively. While the non-closure of the surface energy balance persists across all soil moisture states, it improves over drier soils due to overestimated LH. The non-closure does not imply that F c is always underestimated; F c can be overestimated over dry soils despite the non-closure issue.

Geology↗