Search NASA⌕ Search

SEARCH · Search NASA

Results for “site balance model”

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 145 records · Page 8

Analysis of soil and species composition

Measurements were made during May to October, 1987 and June to August, 1989 over a tallgrass prairie near Manhattan, Kansas. Soil at the experimental site is predominantly Dwight silty clay loam. The prairie was burned on 16 April 1987 and on 28 April 1989 to improve the mix of grasses and forbs. The experimental area was not grazed during 1986 - 1989. A summary of results are given for soil moisture and plant growth; momentum flux and canopy aerodynamic characteristics; evapotranspiration, components of energy balance and canopy conductance; modeling canopy stomatal conductance; canopy photosynthesis, photosynthetic efficiency and water use efficiency; modeling canopy photosynthesis; the carbon dioxide budget in a temperate grassland ecosystem; and photosynthesis and stomatal conductance related to reflectance on the canopy scale.

Verma, Shashi B.↗

NASA's SnowEx Campaign: Observing Seasonal Snow in a Forested Environment

SnowEx is a multi-year airborne snow campaign with the primary goal of addressing the question: How much water is stored in Earth's terrestrial snow-covered regions? Year 1(2016-17) focused on the distribution of snow-water equivalent (SWE) and the snow energy balance in a forested environment. The year 1 primary site was Grand Mesa and the secondary site was the Senator Beck Basin, both in western Colorado, USA. Nine sensors on five aircraft made observations using a broad range of sensing techniques - active and passive microwave, and active and passive optical/infrared - to determine the sensitivity and accuracy of these potential satellite remote sensing techniques, along with models, to measure snow under a range of forest conditions. SnowEx also included an extensive range of ground truth measurements - in-situ manual samples, snow pits, ground based remote sensing measurements, and sophisticated new techniques. A detailed description of the data collected will be given and some preliminary results will be presented.

cryosphere↗

Spatial patterns of snow distribution in the sub-Arctic

Abstract. The spatial distribution of snow plays a vital role in sub-Arctic and Arctic climate, hydrology, and ecology due to its fundamental influence on the water balance, thermal regimes, vegetation, and carbon flux. However, the spatial distribution of snow is not well understood, and therefore, it is not well modeled, which can lead to substantial uncertainties in snow cover representations. To capture key hydro-ecological controls on snow spatial distribution, we carried out intensive field studies over multiple years for two small (2017–2019; ∼ 2.5 km2) sub-Arctic study sites located on the Seward Peninsula of Alaska. Using an intensive suite of field observations (> 22 000 data points), we developed simple models of the spatial distribution of snow water equivalent (SWE) using factors such as topographic characteristics, vegetation characteristics based on greenness (normalized different vegetation index, NDVI), and a simple metric for approximating winds. The most successful model was random forest, using both study sites and all years, which was able to accurately capture the complexity and variability of snow characteristics across the sites. Approximately 86 % of the SWE distribution could be accounted for, on average, by the random forest model at the study sites. Factors that impacted year-to-year snow distribution included NDVI, elevation, and a metric to represent coarse microtopography (topographic position index, TPI), while slope, wind, and fine microtopography factors were less important. The characterization of the SWE spatial distribution patterns will be used to validate and improve snow distribution modeling in the Department of Energy's Earth system model and for improved understanding of hydrology, topography, and vegetation dynamics in the sub-Arctic and Arctic regions of the globe.

54 ENVIRONMENTAL SCIENCES↗

Overview of SnowEx Year 1 Activities

SnowEx is a multi-year airborne snow campaign with the primary goal of addressing the question: How much water is stored in Earths terrestrial snow-covered regions? Year 1 (2016-17) focused on the distribution of snow-water equivalent (SWE) and the snow energy balance in a forested environment. The year 1 primary site was Grand Mesa and the secondary site was the Senator Beck Basin, both in western, Colorado, USA. Nine sensors on five aircraft made observations using a broad range of sensing techniques, active and passive microwave, and active and passive optical infrared to determine the sensitivity and accuracy of these potential satellite remote sensing techniques, along with models, to measure snow under a range of forest conditions. SnowEx also included an extensive range of ground truth measurements in-situ manual samples, snow pits, ground based remote sensing measurements, and sophisticated new techniques. A detailed description of the data collected will be given and some preliminary results will be presented.

snow↗

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimated Hourly Electricity Demand Profiles for Each County in the Contiguous United States

In this data descriptor, we present hourly electricity demand estimates for each county in the contiguous United States from 2016 to 2023. The demand profiles represent the sum of two components for each county. The first is the weighted average of reported hourly demand profiles for North American Electric Reliability Corporation balancing authority (BA) regions and subregions, scaled to match annual estimates of county-level retail sales and direct use of electricity and weighted by the estimated percentage of county load served by each BA region or subregion. The second is the weighted average of modeled hourly, county- and sector-level distributed photovoltaic (DPV) capacity factor profiles, scaled to match annual estimates of on-site consumption of DPV-generated electricity for each county and weighted by the percentage of consumption attributable to each sector.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Three Paradigms of Lunar Regolith Evolution

Integration of diverse datasets on the Moon may render some paradigms of lunar science either better-defended or vulnerable. We will consider three paradigms commonly used for understanding the processes of lunar regolith evolution in light of new and accumulated data. Our premise is that all data-sets should converge to a single interpretation if a concept or model is to be accepted as a paradigm. If a convergence is lacking, the paradigm needs fresh scrutiny. SteadyState: Lunar regolith evolution is currently understood in terms of comminution, agglutination, and replenishment as described by McKay and coworkers). Briefly, the model envisages continued micrometeoritic bombardment to comminute exposed soil particles to finer sizes while continued agglutination consumes finer sizes to produce larger constructional particles. Eventually, a balance between these two opposing processes achieves a steady state; soils at steady state maintain their mean grain size (M(sub z)). Episodic higher-energy impacts excavate fresh coarse material from below the soil cover, disturb the steady state, and restart the process to achieve a new steady state. It follows that the thickness of the regolith at any site would control the frequency of replenishment; indeed, the thickness of the regolith at Apollo landing sites was predicted by McKay et al. from the average M(sub z) of local soils. However, replenishment may come also from disintegrating boulders and cobbles at the lunar surface, and rates of comminution and agglutination may depend on the properties of target material. Regression between M(sub z) and I(sub s)/Fe(sup 0) (a measure of maturity or total surface exposure) of Apollo soils at different sites shows the following relations and estimated M(sub z) at a high maturity of I(sub s)/Fe(sup 0)= 100. It is possible that Apollo 12 and 15 sites have the thickest regolith and the Apollo 16 site has the thinnest. It is also possible that Apollo 12 and 15 basalts are comminuted faster than Apollo 16 highland rocks and Apollo 14 and 17 soils are products of mixed parentage. If a soil becomes continually finer as it matures until agglutination catches up, and if comminution is differential-dependent on the physical properties of the constituents, then the composition of the bulk soil has to match the composition of some "fulcrum" grain size fraction, say X Grain size fractions >X and <X will complement each other; their mass balance is the bulk soil. It appears that the 10-20-micron size fraction may be the fulcrum. In general, trace-element chemistry and IR reflectance spectra of this size fraction are closest to that of the bulk soil, regardless of maturity that is surprising. Disaggregated products of regolith breccias may also show similar relationships. If the 10-20 gm is the fulcrum (i.e., X as above) for many soil properties (e.g., major element composition, FMR, solar-wind-implanted elements), then this may be the ultimate mean grain size of lunar soils at steady state. However, different properties of soils may find steady states at different grain size fractions. The steady state of solar-wind-implanted elements, on the other hand, will climb up the grain-size scale as agglutinates transfer surface-correlated components into volume correlated components until a saturation level is reached or the rates of replenishment and implantation become equal. The same will be the case with vapor-deposited reduced metals as they too are incorporated inside constructional particles. Properties that are directly affected by soil-maturation processes will thus have different pathways of achieving steady states. Maturity, i.e., cumulative surface exposure, of lunar soils is best quantified by the amount of nanophase superparamagnetic Fe(sup 0) (np-Fe(sup 0)) normalized to Fe content (=I(sub s)/Fe(sup 0). The majority consensus (paradigm?) for the production of np-Fe(sup 0) is associated with the production of agglutinates. Because large doses of solar-wind H are implanted in all lunar soils upon exposure, any melting (e.g., during agglutinate production) triggers a chemical reduction of Fe-bearing minerals resulting in np-Fe(sup 0) production. The quantity of np-Fe(sup 0) is thus dependent on melting events, (i.e., exposure), and limited by the Fe content of the soil. All freshly produced np-Fe(sup 0) resides in agglutinitic glass, as new TEM images show. Apparently, the correction procedure developed by Lucey et al. to estimate the Fe content of the lunar surface from IR-reflectance spectra depends on accepting the above. However, the process of producing np-Fe(sup 0) may be physical rather than chemical. All np-Fe(sup 0) could be deposits from a vapor produced by micrometeoritic impact on lunar soils. If metal-O bonds in target phases are broken, O being "most volatile" will escape leaving an O-deficient vapor to facilitate the production of np-Fe(sup 0). If so, the quantity of np-Fe(sup 0) is dependent on the vaporizing events, (i.e., exposure), and limited by the efficiency of breaking metal-O bonds and the escape of 0. To the extent that strengths of metal-O bonds are dependent on the local crystal field, production of np-Fe(sup 0) may be limited by the mineral composition of target soils and not by their total Fe content. According to this model, vapor-deposited np-Fe(sup 0) should be found at any retentive sites on lunar soil grains. Indeed, TEM images show np-Fe(sup 0) on plagioclase and ilmenite. Incorporation of such pre-irradiated np-Fe(sup 0)-bearing grains into agglutinates may account for eventual increased emplacement of np-Fe(sup 0) in agglutinates. Such a paradigm shift in understanding the origin of np-Fe(sup 0) will raise questions ranging from the unquestionable use of Is/FeO as the universal maturity parameter of lunar soils to global elemental maps of the Moon from remote-sensing data. Additional information is contained in the original.

Basu, A.↗

Connecting the Radiative Influences of Aerosol upon the Mass Flux Profiles of Shallow Cumuli across the Southeast Atlantic Ocean Basin and its Boundaries (Final Report)

The Atlantic Ocean covers approximately 25% of Earth’s surface and the atmosphere above it is home to a complex array of clouds and aerosols that have important influences on regional and global weather and climate. These influences must be accurately depicted in short range, medium range, seasonal, and climate forecast models. Conditions over the tropical Atlantic are particularly complex due to continental scale plumes of dust from the Sahara Desert and smoke from agricultural burning in Africa that drift across the Atlantic Ocean basin toward the Americas. These plumes are often found meandering in the lower atmosphere above shallow tropical clouds that form above the ocean surface, presumably mingling with these clouds on occasion due to convective mixing processes. Elevated dust and smoke particles absorb incoming sunlight and substantially warm the marine atmosphere in the layer in which they are present. This warming may alter the thermal stability of the marine atmosphere and may throttle or enhance the development of clouds, change their internal structure and the rate at which they precipitate. Alternatively, it may isolate the lower atmosphere from drier layers above enabling water vapor to accumulate near the ocean surface potentially leading to the development of deeper convection. Our study was organized around the principal concept of determining the impact of African smoke and dust plumes upon cloud development at ASI and understanding how a popular shallow convection parameterization used in models responds to the presence of this aerosol. We analyzed observations collected during the US Department of Energy (DOE) Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign using the Atmospheric Radiation Measurement Program’s Mobile Facility #1 (AMF-1), which was deployed to Ascension Island (ASI) in Southeastern Atlantic for a one-year period. Ascension Island is immediately downwind from an African source of these plumes, but far enough removed to enable the lower atmosphere to have reacted to their presence. A main goal of our study was to compute radiative heating rate profiles over Ascension Island and along the trajectory from the biomass burning regions along coastal Africa to Ascension Island. To compute the radiative heating rate due to aerosols and clouds, we employed observed profiles of temperature, humidity, and clouds from LASIC alongside aerosol optical properties from the Modern-era Retrospective analysis for Research and Applications Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM). Radiative heating was also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We were successful in this effort and the resulting publication is already being cited despite its relatively short lifetime in the literature. The second part of our study involved the process-level representation of the clouds observed over the Southeastern Atlantic in models. Our initial task, upon which the balance of this portion of the study depended, was to evaluate the representativeness of the clouds observed by AMF-1. What orographic influence did Ascension Island have on the measured cloud structure? We set out to answer this question by trying to separate observations that were clearly indicative of orographic forcing from those that were representative of open ocean. The siting of AMF-1, which was 340-m above sea-level on the slope of a steep escarpment, proved demonstrably problematic for cloud process measurements. Despite considerable effort using multiple approaches including artificial intelligence (AI), we were unable to successfully compensate for the orographic forcing present at the AMF-1 site on ASI and, hence, produce process-level summaries of the convective mass flux representative of open ocean over the Southeastern Atlantic. Even though the project has officially ended, we are making a final attempt using two new types of clustering algorithm (i.e., AI) on the recommendation of a former student, who is an expert in this area. Since we have only recently begun to test these new AI schemes, the recommendations from the second part of our project outlined below are based on our experience at the time of this report.

54 ENVIRONMENTAL SCIENCES↗

No evidence of canopy-scale leaf thermoregulation to cool leaves below air temperature across a range of forest ecosystems

Understanding and predicting the relationship between leaf temperature (T leaf ) and air temperature (T air ) is essential for projecting responses to a warming climate, as studies suggest that many forests are near thermal thresholds for carbon uptake. Based on leaf measurements, the limited leaf homeothermy hypothesis argues that daytime T leaf is maintained near photosynthetic temperature optima and below damaging temperature thresholds. Specifically, leaves should cool below T air at higher temperatures (i.e., > ~25–30°C) leading to slopes <1 in T leaf /T air relationships and substantial carbon uptake when leaves are cooler than air. This hypothesis implies that climate warming will be mitigated by a compensatory leaf cooling response. A key uncertainty is understanding whether such thermoregulatory behavior occurs in natural forest canopies. We present an unprecedented set of growing season canopy-level leaf temperature (T can ) data measured with thermal imaging at multiple well-instrumented forest sites in North and Central America. Our data do not support the limited homeothermy hypothesis: canopy leaves are warmer than air during most of the day and only cool below air in mid to late afternoon, leading to T can /T air slopes >1 and hysteretic behavior. We find that the majority of ecosystem photosynthesis occurs when canopy leaves are warmer than air. Using energy balance and physiological modeling, we show that key leaf traits influence leaf-air coupling and ultimately the T can /T air relationship. Canopy structure also plays an important role in T can dynamics. Future climate warming is likely to lead to even greater T can , with attendant impacts on forest carbon cycling and mortality risk.

59 BASIC BIOLOGICAL SCIENCES↗

Spin prediction techniques

The NASA Langley Research Center has the responsibility to advance the state-of-the-art in the area of stall/spin technology. This includes the development and validation of experimental and analytical techniques for the prediction of stall/spin characteristics. As a part of this effort, a rotary balance facility located in the Langley spin tunnel was developed two and a half years ago to rapidly identify an airplane's aerodynamic characteristics in a rotational flow environment. On-line rotary balance data plots and on-line predicted steady spin modes permit the designer to develop, on site, a configuration highly resistant to spins, or one which has good spin characteristics, i.e., no spin equilibrium conditions possible with lateral-directional controls neutral, if the airplane is to be used for acrobatic maneuvers or training. The rotary balance data are also used to compute time histories of a spin's incipient, developed, and recovery phases. These spin analysis techniques and their correlation with spin tunnel model and full-scale flight results are discussed herein.

Bihrle, W., Jr.↗

Watching the Grand Ethiopian Renaissance Dam from a distance: Implications for sustainable water management of the Nile water

Increased demands for sustainable water and energy resources in densely populated basins have led to the construction of dams, which impound waters in artificial reservoirs. In many cases, scarce field data led to the development of models that underestimated the seepage losses from reservoirs and ignored the role of extensive fault networks as preferred pathways for groundwater flow. We adopt an integrated approach (remote sensing, hydrologic modeling, and field observations) to assess the magnitude and nature of seepage from such systems using the Grand Ethiopian Renaissance Dam (GERD), Africa's largest hydropower project, as a test site. The dam was constructed on the Blue Nile within steep, highly fractured, and weathered terrain in the western Ethiopian Highlands. The GERD Gravity Recovery and Climate Experiment Terrestrial Water Storage (GRACETWS), seasonal peak difference product, reveals significant mass accumulation (43 ± 5 BCM) in the reservoir and seepage in its surroundings with progressive south-southwest mass migration along mapped structures between 2019 and 2022. Seepage, but not a decrease in inflow or increase in outflow, could explain, at least in part, the observed drop in the reservoir's water level and volume following each of the three fillings. Using mass balance calculations and GRACETWS observations, we estimate significant seepage (19.8 ± 6 BCM) comparable to the reservoir's impounded waters (19.9 ± 1.2 BCM). Investigating and addressing the seepage from the GERD will ensure sustainable development and promote regional cooperation; overlooking the seepage would compromise hydrological modeling efforts on the Nile Basin and misinform ongoing negotiations on the Nile water management.

GRACE and GRACE-FO↗

Understanding Solar Photovoltaic System Performance: An Assessment of 75 Federal Photovoltaic Systems

This report presents a performance analysis of 75 photovoltaic systems based on PV system production data collected as part of a FEMP Federal PV Performance Assessment project combined with co-incident insolation, and ambient temperature to analyze how actual performance compares with a performance model. FEMP collaborated with 17 Federal agencies and sub-agencies to collect the information required to analyze the performance of each system. The systems represent a total capacity of 30,714 kW and range in size from 1 kW to 4,043 kW, with an average size of 410 kW, and were installed between 2011 and 2020. The data is analyzed for Key Performance Indicators, Availability, Performance Ratio and Energy Ratio by comparing the measured production data to model production data. The System Advisor Model (SAM) combines a description of the system (such as inverter capacity, de-rating for temperature, balance-of-system efficiency) with environmental parameters (coincident solar and temperature data) to calculate predicted performance. The performance metrics are calculated by lining up the measured production data with the model estimate on an hour-by-hour, day-by-day, or month-by-month basis (depending on the interval resolution of the production data). A report with system description, photo of the system, special assumptions made for the site, graph of measured production and model production, table of key performance indicators, and links to O&M resources that might improve performance was produced and delivered to site and agency staff with a short on-line briefing.

14 SOLAR ENERGY↗

Soil Moisture Active Passive (SMAP) Project Assessment Report for Version 4 of the L4_SM Data Product

This report provides an assessment of Version 4 of the SMAP Level 4 Surface and Root Zone Soil Moisture (L4_SM) product, released on 14 June 2018. The assessment includes comparisons of L4_SM soil moisture and temperature estimates with in situ measurements from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product, including observation-minus-forecast (O-F) brightness temperature residuals and soil moisture analysis increments.Together, the core validation site comparisons and the statistics of the assimilation diagnostics areconsidered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to upscaling errors from the point-scale to the grid-cell scale of the data product.The Version 4 L4_SM product benefits from an improved land surface modeling system and from retrospective surface meteorological forcing data that are as consistent as possible with the present-day datain terms of their climatology. Specifically, the model changes include revised parameters and parameterizations for (i) the surface energy balance, (ii) recharge from below of the model's surface excess reservoir, and (iii) the snow depletion curve. Updated ancillary inputs include improved datasets for landcover, topography, and vegetation height. The Version 4 algorithm further includes a revised approach to precipitation corrections that improves the precipitation climatology in Africa and the high-latitudes. Moreover, for system calibration the model is forced retrospectively with MERRA-2 reanalysis data, which are more consistent with the near-real time GEOS forward processing (FP) data used during the SMAP period than the retrospective GEOS data that were available for previous L4_SM versions. An analysis of the time-average surface and root zone soil moisture shows that the global pattern ofarid and humid regions is captured by the Version 4 L4_SM estimates. Owing to the changes in the landsurface modeling system, surface soil moisture is typically drier by several volumetric percent in Version 4 compared to Version 3, whereas root zone soil moisture is wetter in Version 4 in some regions and drierin others. Because of these climatological differences, the Version 3 and Version 4 products should not be combined into a single dataset for use in applications.Results from the core validation site comparisons indicate that Version 4 of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the RMSE after removal of the long-term mean difference (ubRMSE). The overall ubRMSE of the 3-hourly L4_SM dataat the 9 km scale is 0.039 m3 m-3 for surface soil moisture and 0.029 m3 m-3 for root zone soil moisture,below the 0.04 m3 m-3 requirement. The L4_SM estimates are an improvement over estimates from a model-only Nature Run version 7.2 (NRv7.2), which demonstrates the beneficial impact of the SMAP brightness temperature data. Overall, L4_SM surface and root zone soil moisture estimates are more skillful than NRv7.2 estimates, with statistically significant improvements at the 5% level for surface soil moisture R and anomaly R values. Results from comparisons of the L4_SM product to i

Reichle, Rolf H.↗

LAI Assimilation Schedules to Constrain Uncertain Cultivars and Soils in CERES-Maize

In anticipation of large-domain crop model applications where precise local configuration and calibration is not possible, we describe benefits and potential drawbacks of employing a crop pest module to achieve leaf area index (LAI) assimilation into a high performing CERES-Maize crop model configuration at a field experiment site in Perry, Iowa. Simulation experiments explore the use of MODIS satellite-derived LAI to constrain and adjust LAI to counter imprecise cultivar and soil configurations often occurring in the absence of high-quality local information. Simulations using single-day, window, and continuous LAI replacement across 5 cultivars and 2 soil calibration approaches for nine corn rotation years from 2004 to 2020 led to different yield outcomes and reverberations throughout the field environment. Evaluating variance and mean bias, results indicate minimal interventions in early vegetative and grain-filling stages were more beneficial than use of continuous LAI adjustments, as they minimized disruptions to the internal resource balances governing plant stresses and grain production. LAI adjustment was particularly helpful in constraining growth related to uncertain thermal unit requirements and leaf tip appearance rates (P1 and PHINT cultivar parameters, respectively). Findings underscore the need to assimilate additional state variables to ensure internal biophysical coherence. This approach shows promise for applications spanning wider domains with prediction time pressure where detailed configuration, more complex assimilation methods, or recalibration of crop model parameters may not be practical.

Phenology↗

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗

Bedrock weathering contributes to subsurface reactive nitrogen and nitrous oxide emissions

Atmospheric nitrous oxide contributes directly to global warming, yet models of the nitrogen cycle do not account for bedrock, the largest pool of terrestrial nitrogen, as a source of nitrous oxide. Although it is known that release rates of nitrogen from bedrock are large, there is an incomplete understanding of the connection between bedrock-hosted nitrogen and atmospheric nitrous oxide. Here, we quantify nitrogen fluxes and mass balances at a hillslope underlain by marine shale. We found that, at this site, bedrock weathering contributes 78% of the subsurface reactive nitrogen, while atmospheric sources (commonly regarded as the sole sources of reactive nitrogen in pristine environments) account for only the remaining 22%. About 56% of the total subsurface reactive nitrogen denitrifies, including 14% emitted as nitrous oxide. The remaining reactive nitrogen discharges in porewaters to a floodplain where additional denitrification probably occurs. We also found that the release of bedrock nitrogen occurs primarily within the zone of the seasonally fluctuating water table and suggest that the accumulation of nitrate in the vadose zone, often attributed to fertilization and soil leaching, may also include contributions from weathered nitrogen-rich bedrock. Our hillslope study suggests that, under oxygenated and moisture-rich conditions, weathering of deep, nitrogen-rich bedrock makes an important contribution to the nitrogen cycle.

58 GEOSCIENCES↗

Bedrock weathering rates, reactive nitrogen influxes and effluxes, and nitrous oxide emissions rates from the Pumphouse Hillslope, East River Watershed, Colorado

Atmospheric nitrous oxide contributes directly to global warming, yet models of the nitrogen cycle do not account for bedrock, the largest pool of terrestrial nitrogen, as a source of nitrous oxide. Although it is known that release rates of nitrogen from bedrock are large, there is an incomplete understanding of the connection between bedrock-hosted nitrogen and atmospheric nitrous oxide. Here, we quantify nitrogen fluxes and mass balances at a hillslope underlain by marine shale. We found that at this site bedrock weathering contributes 78% of the subsurface reactive nitrogen, while atmospheric sources (commonly regarded as the sole sources of reactive nitrogen in pristine environments) account for only the remaining 22%. About 56% of the total subsurface reactive nitrogen denitrifies, including 14% emitted as nitrous oxide. The remaining reactive nitrogen discharges in porewaters to a floodplain where additional denitrification likely occurs. We also found that the release of bedrock nitrogen occurs primarily within the zone of the seasonally fluctuating water table and suggest that the accumulation of nitrate in the vadoes zone, often attributed to fertilization and soil leaching, may also include contributions from weathered nitrogen-rich bedrock. Our hillslope study suggests that under oxygenated and moisture-rich conditions, weathering of deep, nitrogen-rich bedrock makes an important contribution to the nitrogen cycle. The data files are in Excel, which can be accessed using Microsoft Office, and consist of many data sets from the Pumphouse Hillslope PLM (Pumphouse Lower Montane) 1, 2, 3, and 4. They include soil to rock (0-10 meters) solid phase minerals and elements compositions; time- and depth-resolved pore-water chemistry and pore-gas compositions; time- and depth-resolved water table depths and water fluxes; subsurface weathering rates; nitrogen influxes and effluxes and mass balance. The attached paper is in a Word document, which can be accessed using Microsoft Office, and in pdf format.

54 ENVIRONMENTAL SCIENCES↗

Gusev-Meridiani-Type Soil Component Dissolved in Some Shock Glasses in Shergottites

Modal analysis, based on APXS, MiniTES and Mossbauer results obtained at Gusev and Meridiani sites on Mars, indicates that Martian soils consist predominantly of igneous minerals such as olivine, pyroxene and feldspar (approximately70 - 80%), with the balance consisting of alteration minerals such as sulfates, silica and chlorides]. These studies also showed that soil alteration did not occur in-situ and igneous and alteration components are derived from different sources. Below, we analyse the chemical abundance data obtained from shock glasses in shergottites using mass balance mixing models. In these models, the two main end members used are (a) host rock chemical composition and (b) the GM soils average composition as the second component. Here, we consider the S-bearing phases as indicators of added alteration phases in the shock glasses and GM soils. Although the S-bearing phase in shock glasses occurs as micron sized sulfide blebs, we showed in earlier abstracts that sulfur was originally present as sulfate in impact glass-precursor materials and was subsequently reduced to sulfide during shock melting. This conclusion is based on results obtained from S-K XANES studies, Fe/S atomic ratios in sulfide blebs and 34S/32S isotopic measurements in these sulfide blebs. Additionally, sulfur in several EET79001 Lith. A glasses is found to correlate positively with Al2O3 and CaO (and negatively with FeO and MgO), suggesting the presence of Al- and Ca- sul-fate-bearing phases. The distribution of the 87Sr/86Sr iso-topic ratios determined in Lith. A glasses (,27 &,188 and,54) indicate that Martian soil gypsum and host rock material were mixed with each other in the glass precursors. In some vugs in Lith A glass,27 detected gypsum laths. Furthermore, the Martian regolith-de-rived component (where sulfur typically occurs as sul-fate) is identified in these glasses by determining neutron produced isotopic excesses/deficits in 80Kr and 149Sm isotopes. Moreover, the suggestion that the large amount of sulfur found in,507 was sourced from pyrrhotite in the host rock, would require that excessive quantities of host rock would need to be stripped of sulfur to make this sulfide-rich glass. These results provide ample evidence that S occurred as sulfate and was added to glass precursor materials prior to impact shock.

Ross, D. K.↗