Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial variability”

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

Evaluation of the Sensitivity of the Amazonian Diurnal Cycle to Convective Intensity in Reanalyses

Model parameterizations of tropical deep convection are unable to reproduce the observed diurnal and spatial variability of convection in the Amazon, which contributes to climatological biases in the water cycle and energy budget. Convective intensity regimes are defined using percentiles of daily minimum 3-hourly averaged outgoing longwave radiation (OLR) from Clouds and the Earth's Radiant Energy System (CERES). This study compares the observed spatial variability of convective diurnal cycle statistics for each regime to MERRA-2 and ERA-Interim (ERA) reanalysis data sets. Composite diurnal cycle statistics are computed for daytime hours (06:00-21:00 local time) in the wet season (December-January-February). MERRA-2 matches observations more closely than ERA for domain averaged composite diurnal statistics-specifically precipitation. However, ERA reproduces mesoscale features of OLR and precipitation phase associated with topography and the propagation of the coastal squall line. Both reanalysis models are shown to underestimate extreme convection.

Itterly, Kyle F.↗

Emerging mobile lidar technology to study boundary layer winds influenced by operating turbines

The development of a microjoule-class pulsed Doppler lidar and deployment of this compact system on mobile platforms such as aircraft, ships, or trucks have opened a new opportunity to characterize the dynamics of complex mesoscale wind flows. The PickUp-based Mobile Atmospheric Sounder (PUMAS) truck-based lidar system was recently used during the American Wake Experiment (AWAKEN) to assess the general structure of boundary layer (BL) wind and turbulence around wind turbines in central Oklahoma. Wind speed profiles averaged over PUMAS transects influenced by the operating turbines (waked flow) show a 1–2 m s −1 reduction compared to mean undisturbed (free flow) wind speed profiles. Spatial variability in wind speed was observed in time–height cross-sections at different distances from turbines. The wind speeds were about 9–12 m s −1 at 6 km distance compared to 5–7 m s −1 at the transects near the turbines. The PUMAS dataset from AWAKEN demonstrated the capability of the mobile Doppler lidar system to document spatial variability in wind flows at different distances from wind turbines and obtain quantitative estimates of wind speed reduction in the waked flow. The high-frequency, simultaneous measurements of the horizontal and vertical winds provide a new approach for characterizing dynamic processes critical for wind farm wake analyses.

17 WIND ENERGY↗

Time-marching solution of incompressible Navier-Stokes equations for internal flow

Primitive variables with central differencing on a staggered grid are used in the present, factored ADI finite-difference scheme for artificial compressibility method solution of the incompressible Navier-Stokes equations, leading to a close coupling between velocity and pressure that both enhances stability and eliminates the need for artificial damping. Computational efficiency is enhanced through the use of a spatially variable, fixed Courant number-based time-step. The numerical results obtained for a driven cavity at Re of 10,000, with local cell Re as high as 100, exhibits no flow variable spatial oscillations on a 40 x 40 stretched grid solution.

Soh, W. Y.↗

A Multi-Pixel Split-Window Approach to Sea Surface Temperature Retrieval from Thermal Imagers with Relatively High Radiometric Noise: Preliminary Studies

In the following decade(s), a set of satellite missions carrying thermal infrared (TIR) imagers with a relatively high noise equivalent differential temperature (NEdT) are expected, e.g., the high resolution TIR imagers flying on the future Thermal infraRed Imaging Satellite for High-resolution Natural resource Assessment (TRISHNA), Land Surface Temperature Monitoring (LSTM) and NASA-JPL/ASI Surface Biology and Geology Thermal (SBG) missions or the secondary payload on board the ESA Earth Explorer 10 Harmony. The instruments on board these missions are expected to be characterized by an NEdT of ⪆0.1 K. In order to reduce the impact of radiometric noise on the retrieved sea surface temperature (SST), this study investigates the possibility of applying a multi-pixel atmospheric correction based on the hypotheses that (i) the spatial variability scales of radiatively active atmospheric variables are, on average, larger than those of the SST and (ii) the effect of atmosphere is accounted for via the split window (SW) difference. Based on 32 Sentinel 3 SLSTR case studies selected in oceanic regions where SST features are mainly driven by meso to sub-mesoscale turbulence (e.g., corresponding to major western boundary currents), this study documents that the local spatial variability of the SW difference term on the scale of ≃3 × 3 km2 is comparable with the noise associated with the SW difference. Similarly, the power spectra of the SW term are shown to have, for small scales, the behavior of white noise spectra. On this basis, we suggest to average the SW term and to use it for the atmospheric correction procedure to reduce the impact of radiometric noise. In principle, this methodology can be applied on proper scales that can be dynamically defined for each pixel. The applicability of our findings to high-resolution TIR missions is discussed and an example of an application to ECOSTRESS data is reported.

sea surface temperature↗

Satellite and Surface Perspectives of Snow Extent in the Southern Appalachian Mountains

Assessing snow cover patterns in mountain regions remains a challenge for a variety of reasons. Topography (e.g., elevation, exposure, aspect, and slope) strongly influences snowfall accumulation and subsequent ablation processes, leading to pronounced spatial variability of snow cover. In-situ observations are typically limited to open areas at lower elevations (<1000 m). In this paper, we use several products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to assess snow cover extent in the Southern Appalachian Mountains (SAM). MODIS daily snow cover maps and true color imagery are analyzed after selected snow events (e.g., Gulf/Atlantic Lows, Alberta Clippers, and Northwest Upslope Flow) from 2006 to 2012 to assess the spatial patterns of snowfall across the SAM. For each event, we calculate snow cover area across the SAM using MODIS data and compare with the Interactive Multi-sensor Snow and ice mapping system (IMS) and available in-situ observations. Results indicate that Gulf/Atlantic Lows are typically responsible for greater snow extent across the entire SAM region due to intensified cyclogenesis associated with these events. Northwest Upslope Flow events result in snow cover extent that is limited to higher elevations (>1000 m) across the SAM, but also more pronounced along NW aspects. Despite some limitations related to the presence of ephemeral snow or cloud cover immediately after each event, we conclude that MODIS products are useful for assessing the spatial variability of snow cover in heavily forested mountain regions such as the SAM.

Sugg, Johnathan W.↗

Growth modeling to evaluate alternative cultivation strategies to enhance national microalgal biomass production

In this paper, we present a model-based assessment of alternative cultivation strategies on algal biomass production within the conterminous United States. Our assessment focuses on two basic cultivation strategies: 1) seasonal rotation of three representative freshwater algal strains that are well suited for warm weather, cold weather, and all-season diverse weather conditions, respectively; and 2) variation between three open pond water depths (15cm, 20cm, and 25cm). The enhanced Biomass Assessment Tool (BAT) is applied on a site-specific, hourly basis at 5,832 North American Land Data Assimilation System Phase 2 1/8° meteorological model grid cells over a 30-year period (1980-2009) to evaluate the operational strategies. Recognizing that resource management decisions may also consider scales beyond the individual site, we also conduct a regional assessment of focused on 7 representative climate zones in the conterminous United States. Results demonstrate that the spatial variability of algal productivity is largely affected by the strain-specific growth response to light and temperature that vary significantly by climate zone and latitude. Therefore, with a limited set of alternative cultivation strategies, the BAT identified the best performing combination of cultivation strategies on a site-specific basis that considerably enhances national annual biomass productivity. Specifically, the appropriate choice of seasonal strain rotation is able to significantly dampen climate-driven seasonal and spatial variability.

59 BASIC BIOLOGICAL SCIENCES↗

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗

A stochastic precipitation disaggregation scheme for GCM applications

In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Gao, Xiaogang↗

A Conceptual Approach to Assimilating Remote Sensing Data to Improve Soil Moisture Profile Estimates in a Surface Flux/Hydrology Model: Overview - Part 1

Knowledge of the amount of water in the soil is of great importance to many earth science disciplines. Soil moisture is a key variable in controlling the exchange of water and energy between the land surface and the atmosphere. Thus, soil moisture information is valuable in a wide range of applications including weather and climate, runoff potential and flood control, early warning of droughts, irrigation, crop yield forecasting, soil erosion, reservoir management, geotechnical engineering, and water quality. Despite the importance of soil moisture information, widespread and continuous measurements of soil moisture are not possible today. Although many earth surface conditions can be measured from satellites, we still cannot adequately measure soil moisture from space. Research in soil moisture remote sensing began in the mid 1970s shortly after the surge in satellite development. Recent advances in remote sensing have shown that soil moisture can be measured, at least qualitatively, by several methods. Quantitative measurements of moisture in the soil surface layer have been most successful using both passive and active microwave remote sensing, although complications arise from surface roughness and vegetation type and density. Early attempts to measure soil moisture from space-borne microwave instruments were hindered by what is now considered sub-optimal wavelengths (shorter than 5 cm) and the coarse spatial resolution of the measurements. L-band frequencies between 1 and 3 GHz (10-30 cm) have been deemed optimal for detection of soil moisture in the upper few centimeters of soil. The Electronically Steered Thinned Array Radiometer (ESTAR), an aircraft-based instrument operating a 1,4 GHz, has shown great promise for soil moisture determination. Initiatives are underway to develop a similar instrument for space. Existing space-borne synthetic aperture radars (SARS) operating at C- and L-band have also shown some potential to detect surface wetness. The advantage of radar is its much higher resolution than passive microwave systems, but it is currently hampered by surface roughness effects and the lack of a good algorithm based on a single frequency and single polarization. In addition, its repeat frequency is generally low (about 40 days). In the meantime, two new radiometers offer some hope for remote sensing of soil moisture from space. The Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), launched in November 1997, possesses a 10.65 GHz channel and the Advanced Microwave Scanning Radiometer (AMSR) on both the ADEOS-11 and Earth Observing System AM-1 platforms to be launched in 1999 possesses a 6.9 GHz channel. Aside from issues about interference from vegetation, the coarse resolution of these data will provide considerable challenges pertaining to their application. The resolution of TMI is about 45 km and that of AMSR is about 70 km. These resolutions are grossly inconsistent with the scale of soil moisture processes and the spatial variability of factors that control soil moisture. Scale disparities such as these are forcing us to rethink how we assimilate data of various scales in hydrologic models. Of particular interest is how to assimilate soil moisture data by reconciling the scale disparity between what we can expect from present and future remote sensing measurements of soil moisture and modeling soil moisture processes. It is because of this disparity between the resolution of space-based sensors and the scale of data needed for capturing the spatial variability of soil moisture and related properties that remote sensing of soil moisture has not met with more widespread success. Within a single footprint of current sensors at the wavelengths optimal for this application, in most cases there is enormous heterogeneity in soil moisture created by differences in landcover, soils and topography, as well as variability in antecedent precipitation. It is difficult to interpret the meaning of 'mean' soil moisture under such conditions and even more difficult to apply such a value. Because of the non-linear relationships between near-surface soil moisture and other variables of interest, such as surface energy fluxes and runoff, mean soil moisture has little applicability at such large scales. It is for these reasons that the use of remote sensing in conjunction with a hydrologic model appears to be of benefit in capturing the complete spatial and temporal structure of soil moisture. This paper is Part I of a four-part series describing a method for intermittently assimilating remotely-sensed soil moisture information to improve performance of a distributed land surface hydrology model. The method, summarized in section II, involves the following components, each of which is detailed in the indicated section of the paper or subsequent papers in this series: Forward radiative transfer model methods (section II and Part IV); Use of a Kalman filter to assimilate remotely-sensed soil moisture estimates with the model profile (section II and Part IV); Application of a soil hydrology model to capture the continuous evolution of the soil moisture profile within and below the root zone (section III); Statistical aggregation techniques (section IV and Part II); Disaggregation techniques using a neural network approach (section IV and Part III); and Maximum likelihood and Bayesian algorithms for inversely solving for the soil moisture profile in the upper few cm (Part IV).

Crosson, William L.↗

Grassland productivity responds unexpectedly to dynamic light and soil water environments induced by photovoltaic arrays

Abstract Agrivoltaic (AV) systems are designed to coproduce photovoltaic (PV) energy on lands simultaneously supporting agriculture (food/forage production). PV infrastructure in agroecosystems alters resources critical for plant growth, and water‐limited agroecosystems such as grasslands are likely to be particularly sensitive to the unique spatial and temporal patterns of incident sunlight and soil water inherent within AV systems. However, the impact of resource alteration on forage production, the primary ecosystem service from managed grasslands, is poorly resolved. Here, we evaluated seasonal patterns of soil moisture (SM) and diurnal variation in incident sunlight (photosynthetic photon flux density [PPFD]) in a single‐axis‐tracking AV system established in a formerly managed semiarid C 3 grassland in Colorado. Our goals were to (1) quantify dynamic patterns of PPFD and SM within a 1.2 MW PV array in a perennial grassland, and (2) determine how aboveground net primary production (ANPP) and photosynthetic parameters responded to the resource patterns created by the PV array. We hypothesized that spatial variability in ANPP would be strongly related to SM patterns, typical of most grasslands. We measured significant reductions in ANPP directly beneath PV panels, where SM and PPFD were both low. However, in locations with significantly increased SM from the shedding and redistribution of precipitation by PV panels, ANPP was not increased. Instead, ANPP was greatest in locations where plants were shaded in the afternoon but received high levels of PPFD in the morning hours, when air temperatures and vapor pressure deficits were relatively low. Thus, contrary to expectations, we found relatively weak relationships between SM and ANPP despite significant spatial variability in both. Further, there was little evidence that light‐saturated photosynthesis ( A sat ) and quantum yield of CO 2 assimilation (ϕCO 2 ) differed for plants growing directly beneath (lowest PPFD) versus between (highest PPFD) PV panels. Overall, the AV system established in this semiarid managed grassland did not alter patterns of ANPP in ways predictable from past studies of controls of ANPP in open grasslands. However, our results suggest that the diurnal timing of low versus high periods of PPFD incident on plants is an important determinant of productivity patterns in grasslands.

14 SOLAR ENERGY↗

Macro-physical Properties of Shallow Cumulus from Integrated ARM Observations (Final Report)

Fair-weather shallow cumuli (ShCu) play an important role in many climate-related processes. Irregular geometry of ShCu and their strong temporal and spatial variability make it challenging to observe ShCu holistically and to represent them correctly in climate models. To improve ShCu parameterizations, information on both vertically and horizontally resolved cloud properties is required. Commonly, the vertically resolved cloud properties are provided by zenith pointing lidar-radar observations with a very narrow field of view (FOV). Thus, these “pencil-beam” properties may not be representative of a larger surrounding area. Limited number of areal-averaged cloud properties, such as fractional sky cover (FSC), are offered typically by wide-FOV observations. The main goal of our project was to integrate advantages of the narrow-FOV (vertical structure of clouds) and wide-FOV (spatial arrangement of clouds) observations for an improved characterization of single-layer ShCu observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for an 18-yr period (2000-2017). There are four major accomplishments of our project, First, an updated operational cloud classification for days with ShCu has been suggested and evaluated through a detailed comparison with the manually curated records. Our classification extends successfully the latest ARM cloud type Value Added Product (VAP) based on the Active Remote Sensing of Clouds (ARSCL) cloud product by incorporating both cloud fraction (CF) provided by narrow-FOV ceilometer data and FSC from wide-FOV images offered by a Total Sky Imager (TSI). Moreover, our classification allows one to identify impact of instrumentation changes at the SGP site, namely the transition to KAZRARSCL with the updated cloud radar, on the identification of periods with single-layer ShCu. Second, a new approach that resolves cloud area distributions for a given region (up to 4x4 km 2 ) has been suggested and cloud equivalent diameters (CEDs) have been estimated for the first time. These estimations have been performed over a wide range of cloud sizes (about 0.01–3.5 km) with high temporal resolution (30s) using wide-FOV TSI images and cloud base height (CBH) provided by complementary narrow-FOV lidar measurements. Our simple and computationally inexpensive approach offers a previously unavailable dataset for process studies in the convective boundary layer and evaluation of ShCu parameterizations in cloud-resolving models. Third, a long-term integrated record of ShCu macrophysical properties has been developed. The developed record represents the longest available compilation of events with ShCu and includes (i) a novel visualization of the spatial variability in cloud cover both along- and across-wind directions, (ii) updated estimates of narrow-FOV CF and wide-FOV FSC, (iii) updated narrow-FOV CBH, and (iv) complementary data, such as wind speed and direction from the 915-MHz Radar Wind Profiler (RWP) data. The developed record has been used successfully to assess conventional observational estimates of cloud cover and their sensitivity to the following two factors: (i) instrument-dependent cloud detection and data merging criteria and (ii) FOV configuration. Fourth, co-variability of the ShCu macrophysical properties and environmental parameters has been analyzed for a 3-yr period (2016-2018). Our initial analysis includes diurnal changes of FSCs obtained for clouds with small, moderate and large CEDs and several environmental parameters, such as lifted condensation level (LCL) and mixed layer height (zi). Preliminary results of our analysis suggest that the horizontal extent of ShCu is controlled substantially by the sign and magnitude of difference between these two parameters (zi-LCL): the CED tends to grow with increase of this difference (zi exceeds LCL). We have initiated relationships between the ShCu and key atmospheric parameters that control both the development and evolution of ShCu using our new data product, which combines effectively the advantages of narrow-FOV data offered by zenith pointing cloud radars and lidars and wide-FOV TSI images. While the latest instrumentation at the ARM sites may address these challenging relationships in the future, we believe that the historical ARM data at the SGP site has not yet been fully utilized. Overall, our data product can be used by researchers working on a wide range of climate-related projects. These projects may include (i) a comprehensive evaluation of outputs from the Large-Eddy Simulation (LES) and single-column models for their future improvement, (ii) the representativeness of “short-period” results obtained from the previous model and observational studies and (iii) the planning of future field campaigns with focus on improved understanding of the diurnal cycle of cumulus convection.

54 ENVIRONMENTAL SCIENCES↗

Influence of overhead HVAC and aerosol control strategies on coarse mode particle dispersion and exposure in a full-scale room experiment

Coarse mode respiratory aerosols can carry viral loads over long distances and have very different dynamics than submicron particles, but experimental studies under realistic conditions remain limited. Here, to study the differential impacts on exposure under different mixing conditions, we co-released 7–10 µm particles and carbon dioxide (CO 2 )—which served as an indicator of gas and submicron particle dynamics—in a 158 m 3 room at LBNL's FLEXLAB facility with an overhead heating, ventilation, and air conditioning (HVAC) system. The room was arranged as a distanced meeting then a classroom with eight heated manikins and a researcher. Spatial variability was measured using 16 particle counters and 22–26 CO 2 sensors throughout the space. Conditions included: HVAC off or supply air at 1000-1060 m 3 h -1 at neutral, cooling, or heating temperatures; with and without 20% outdoor air; and added HVAC filtration, portable air cleaners (PACs), or a physical barrier between the speaker and occupants. We found that good mixing via neutral or cooling supply air or use of PACs under heating lowered coarse particle exposure at some locations, but increased exposure for one-quarter to two-thirds of manikins compared to poor mixing under heating. A physical barrier reduced direct transfer of coarse particles during heating, but less during cooling. High spatial variability shows that a single measurement cannot represent occupant exposure. Instantaneous air mixing assumptions overstate the effectiveness of ventilation, HVAC filtration, and upper-room germicidal ultraviolet disinfection for coarse particles, as relatively few particles reach the return grille or upper room under most conditions.

Air mixing↗

Calculating Soil Wetness, Evapotranspiration and Carbon Cycle Processes Over Large Grid Areas Using a New Scaling Technique

Soil wetness typically shows great spatial variability over the length scales of general circulation model (GCM) grid areas (approx 100 km ), and the functions relating evapotranspiration and photosynthetic rate to local-scale (approx 1 m) soil wetness are highly non-linear. Soil respiration is also highly dependent on very small-scale variations in soil wetness. We therefore expect significant inaccuracies whenever we insert a single grid area-average soil wetness value into a function to calculate any of these rates for the grid area. For the particular case of evapotranspiration., this method - use of a grid-averaged soil wetness value - can also provoke severe oscillations in the evapotranspiration rate and soil wetness under some conditions. A method is presented whereby the probability distribution timction(pdf) for soil wetness within a grid area is represented by binning. and numerical integration of the binned pdf is performed to provide a spatially-integrated wetness stress term for the whole grid area, which then permits calculation of grid area fluxes in a single operation. The method is very accurate when 10 or more bins are used, can deal realistically with spatially variable precipitation, conserves moisture exactly and allows for precise modification of the soil wetness pdf after every time step. The method could also be applied to other ecological problems where small-scale processes must be area-integrated, or upscaled, to estimate fluxes over large areas, for example in treatments of the terrestrial carbon budget or trace gas generation.

Sellers, Piers↗

Representing lateral groundwater flow from land to river in Earth system models

Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.

54 ENVIRONMENTAL SCIENCES↗

Three dimensional simulation of spatial and temporal variability of stratospheric hydrogen chloride

Spatial and temporal variability of atmospheric HCl columns are calculated for January 1979 using a three-dimensional chemistry-transport model designed to provide the best possible representation of stratospheric transport. Large spatial and temporal variability of the HCl columns is shown to be correlated with lower stratospheric potential vorticity and thus to be of dynamical origin. Systematic longitudinal structure is correlated with planetary wave structure. These results can help place spatially and temporally isolated column and profile measurements in a regional and/or global perspective.

Kaye, Jack A.↗

Multi‐year hydroclimatic droughts and pluvials across the conterminous United States

Abstract Time series of water‐year runoff for 2,109 hydrologic units (HUs) across the conterminous United States (CONUS) for the 1900 through 2014 period were used to identify drought and pluvial (i.e., wet) periods. Characteristics of the drought and pluvial events including frequency, duration, and severity were examined and compared. Additionally, a similar analysis was performed using gridded tree‐ring reconstructions of the Palmer Drought Severity Index (PDSI) for the period 1475 through 2005 to place the drought and pluvial characteristics determined using water‐year runoff for 1900 through 2014 in the context of multi‐century climate variability. The temporal and spatial variability of droughts and pluvials determined using runoff for the 1900 through 2014 period indicated that most drought events in the CONUS occurred before about 1970, whereas most pluvial periods occurred after about 1970. This change in the frequencies of drought and pluvial events around 1970 was largely related to an increase in fall (October through December) precipitation across much of the central United States. Also, the duration and severity of droughts and pluvials identified using runoff for the 1900 through 2014 period generally were not significantly different from the drought and pluvial characteristics identified using the PDSI for the 1475 through 2005 period.

McCabe, Gregory J.↗