Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational data assimilation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Climatology of the Finite-Volume CCM and Sensitivity Studies

A new model has been developed in the Data Assimilation Office at NASA Goddard Space Flight Center. It is intended for data assimilation and general circulation modeling, including representations of radiatively active trace gases such as stratospheric ozone. The model is based on the flux-form semi-Lagrangian' dynamical core, developed by Lin and Rood. In collaboration with the National Center for Atmospheric Research, the physical parametrization package from the Community Climate model, Version 3, has been included. This paper will discuss the climatology of a multi-year integration of the model, with observed boundary conditions (sea-surface temperature, sea ice, and incoming solar radiation). The standard model version has a resolution of 2.5 degrees longitude, 2 degrees latitude and has 55 levels between the surface and 80km. The mean climate and interannual variability of this model will be discussed, with emphasis on the middle atmosphere. The sensitivity to horizontal resolution will be examined, as will the impact of increasing and decreasing the vertical resolution. Changes which differ significantly from the interannual variations of the control run will be highlighted. Impacts of using long-term mean boundary conditions or values which vary from year to year will be discussed. Changes in the tropospheric circulation induced by lowering the upper boundary of the model will be discussed. The emphasis of the paper is on the climatology of the middle atmosphere and its impacts on the troposphere, with particular attention given to the climatologically important region near the tropopause.

Pawson, Steven↗
NASA NTRS2002

TOMS UV Algorithm: Problems and Enhancements

Satellite instruments provide global maps of surface ultraviolet (UV) irradiance by combining backscattered radiance measurements with radiative transfer models. The models are limited by uncertainties in input parameters of the atmosphere and the surface. We evaluate the effects of possible enhancements of the current Total Ozone Mapping Spectrometer (TOMS) surface UV irradiance algorithm focusing on effects of diurnal variation of cloudiness and improved treatment of snow/ice. The emphasis is on comparison between the results of the current (version 1) TOMS UV algorithm and each of the changes proposed. We evaluate different approaches for improved treatment of pixel average cloud attenuation, with and without snow/ice on the ground. In addition to treating clouds based only on the measurements at the local time of the TOMS observations, the results from other satellites and weather assimilation models can be used to estimate attenuation of the incident UV irradiance throughout the day. A new method is proposed to obtain a more realistic treatment of snow covered terrain. The method is based on a statistical relation between UV reflectivity and snow depth. The new method reduced the bias between the TOMS UV estimations and ground-based UV measurements for snow periods. The improved (version 2) algorithm will be applied to re-process the existing TOMS UV data record (since 1978) and to the future satellite sensors (e.g., Quik/TOMS, GOME, OMI on EOS/Aura and Triana/EPIC).

Krotkov, Nickolay↗
NASA NTRS2004

How Consistent are Recent Variations in the Tropical Energy and Water Cycle Resolved by Satellite Measurements?

One notable aspect of Earth's climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. One question is 'How consistent, physically, are these diverse remotely-sensed data sets'? The answer is of crucial importance to understanding climate processes, improving physical models, and improving remote sensing algorithms. In this work we will evaluate two recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project 'FD' radiative flux profiles are available from mid-1983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and TOVS (TIROS Operational Vertical Sounder)thermodynamic profiles. Full and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NASA/Goddard Earth Observing System (GEOS-1) assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. We also employ precipitation measurements from the Global Precipitation Climatology Project (GPCP) and the Tropical Rainfall Measuring Mission (TRMM). Finally, ocean evaporation estimates from the Special Sensor Microwave Imager (SSM/I) are considered as well as derived evaporation from the NCAR/NCEP Reanalysis. Additional information is included in the original extended abstract.

Robertson, F. R.↗
DOE OSTI2017

Development of Geomechanical Screening Tools to Identify Risk: An Experimental and Modeling Approach for Secure CO 2 Storage Final Technical Report

CO 2 injection in the subsurface introduces additional complexity in terms of interactions of the reservoir host pore fluid with natural heterogeneity of hydro-mechanical properties of the rock. Natural heterogeneity of the reservoir rock and caprock formations includes spatial variation of hydro-mechanical properties as well as the presence of natural fractures. We developed a coupled mathematical modeling and experimental framework that takes into account the effect of rock heterogeneity on effective mechanical properties of the rock in contact with CO 2 . We performed accurate laboratory experiments to determine the changes in rock mechanical properties due to mineral dissolution in the presence of carbonic acid generated by CO 2 injection. A two-scale adaptive homogenization framework was then developed to consistently upscale petrophysical and geomechanical properties to the field scale. Further, we also developed high-fidelity, numerical solution schemes, non-linear and linear solvers and preconditioners to solve the coupled flow, reactive transport and geomechanical system. The numerical model was then used to study field scale CO 2 sequestration problems for the two selected field sites: (1) Frio formation and (2) Cranfield site. We also developed robust schemes for field data assimilation, model calibration such that the residual uncertainty at the end of the data assimilation procedures can be faithfully represented while taking into account the coupled geochemical and geomechanical processes.

42 ENGINEERING↗
NASA NTRS2018

A Study of the Carbon Cycle Using NASA Observations and the GEOS Model

The Goddard Earth Observing System (GEOS) model has been developed in the Global Modeling and Assimilation Office (GMAO) at NASA's Goddard Space Flight Center. From its roots in chemical transport and as a General Circulation Model, the GEOS model has been extended to an Earth System Model based on a modular construction using the Earth System Modeling Framework (ESMF), combining elements developed in house in the GMAO with others that are imported through collaborative research. It is used extensively for research and for product generation, both as a free-running model and as the core of the GMAO's data assimilation system. In recent years, the GMAO's modeling and assimilation efforts have been strongly supported by Piers Sellers, building on both his earlier legacy as an observationally oriented model developer and his post-astronaut career as a dynamic leader into new territory. Piers' long-standing interest in the carbon cycle and the combination of models with observations motivates this presentation, which will focus on the representation of the carbon cycle in the GEOS Earth System Model. Examples will include: (i) the progression from specified land-atmosphere surface fluxes to computations with an interactive model component (Catchment-CN), along with constraints on vegetation distributions using satellite observations; (ii) the use of high-resolution satellite observations to constrain human-generated inputs to the atmosphere; (iii) studies of the consistency of the observed atmospheric carbon dioxide concentrations with those in the model simulations. The presentation will focus on year-to-year variations in elements of the carbon cycle, specifically on how the observations can inform the representation of mechanisms in the model and lead to integrity in global carbon dioxide simulations. Further, applications of the GEOS model to the planning of new carbon-climate observations will be addressed, as an example of the work that was strongly supported by Piers in the last months of his leadership of Earth Science at NASA Goddard.

Pawson, Steven↗
NASA NTRS2007

Evaluating the Capacity of Global CO2 Flux and Atmospheric Transport Models to Incorporate New Satellite Observations

As we enter the new era of satellite remote sensing for CO2 and other carbon cyclerelated quantities, advanced modeling and analysis capabilities are required to fully capitalize on the new observations. Model estimates of CO2 surface flux and atmospheric transport are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, and ultimately for future projections of carbon-climate interactions. For application to current, planned, and future remotely sensed CO2 data, it is desirable that these models are accurate and unbiased at time scales from less than daily to multi-annual and at spatial scales from several kilometers or finer to global. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 1998 through 2006. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at lxi degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-2), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to in situ observations at sites in Northern mid latitudes and the continental tropics. The influence of key process representations is inferred. We find that the model can resolve much of the hourly to synoptic variability in the observations, although there are limits imposed by vertical resolution of boundary layer processes. The seasonal cycle and its interannual variations generally respond adequately, but discrepancies in the tropics suggest the need for a refinement of the soil moisture dependence of the respiration flux in CASA. Examples and inferences for interpretation of satellite data will be discussed. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in the terrestrial CO2 sink and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗
DOE OSTI2022

Development of an Open-source Alloy Selection and Lifetime Assessment Tool for Structural Components in CSP

Lack of sufficient data on high temperature mechanical and corrosion behavior of structural materials is a huge barrier in the technological maturity of current and future Concentrating Solar Power (CSP) technologies. Rapid development and selection of materials cannot be achieved by expensive and time-consuming acquisition of experimental data. The goal of the proposed work is development of an open-source alloy selection and lifetime prediction tool that will integrate validated physics-based models to describe influence of temperature, alloy composition, environment and component geometry (thickness) on mechanical and corrosion behavior of Ni and Fe-based alloys employed in molten salts/sCO 2 heat exchangers. This one-year project leveraged the extensive dataset on the creep\corrosion behavior of candidate materials generated at ORNL through past projects and input from current collaborations with industrial partners. Based on previous experience and the feedback provided by industry (Brayton Energy and Echogen), three candidate materials of interest, Ni-based alloys 740H, 282 and 625 and application-specific operating conditions (max. temperature of 730 °C and stress of 150 MPa) were identified for the heat exchanger. An extensive corrosion and creep dataset was assimilated for the relevant operating conditions and was supported by detailed characterization of about 100 metallographic cross-sections. The corrosion dataset consisted of scanning electron microscopy images (secondary electron and backscatter electron), measured concentration profiles of alloying elements using energy dispersive X-ray spectroscopy (EDS), widths of denuded zones (dissolution of strengthening phases) and depths of attack in molten KCl-MgCl 2 mixtures using image analyses. The creep dataset comprised of creep rupture data and creep strain curves (for 740H and 282). Coupled thermodynamic-kinetic microstructure-based models were employed to predict the stress-corrosion induced compositional and phase evolutions in the alloy during operation under the identified operating conditions. Reduced order models were developed from advanced physics-based models and were integrated in a user-friendly alloy selection tool. The corrosion model was able to predict the time to a critical Cr concentration at the oxide/alloy interface (chemical lifetime) within ±10% (1 standard deviation) of typical statistical variation in corrosion tests and EDS measurement errors (±0.5 wt%). The initial scope of the project was limited to predict creep rupture times (Larson-Miller parameter). Based on the input provided by industry, the mechanical lifetime of the heat exchanger is governed by accumulated creep strains (2%) rather than creep rupture. To be able to predict the times to specific creep strains, a more extensive creep model development was undertaken largely beyond the initial scope of the project. The continuum damage mechanics creep model was able to predict times to 2% creep strain, t 2% with an accuracy of ±500h. Ultimately, a screening protocol for SiC was generated to demonstrate the pathway for integration of one of the currently immature materials from a commercial adoption standpoint in the current material evaluation tool. The modeling tool developed here is accessible to the science community and stakeholders and lays the foundation for methods that will enable a rapid evaluation of optimum materials for CSP applications and reliable prediction of material degradation thereby considerably reducing operational costs, improving reliability and increasing overhaul intervals. However, the complete potential of such a tool to include a wider range of materials and test conditions can only be realized with a more concentrated combined experimental-characterization-computation effort.

14 SOLAR ENERGY↗
DOE OSTI2018

Documenting the function of non-cultivated microorganisms in terrestrial ecosystems

The terrestrial biosphere contains a large fraction of global carbon and nearly 70% of the organic carbon in these systems is found in soils. Global changes in atmospheric CO 2 , temperature, precipitation, and ecosystem N inputs, are expected to impact primary production and carbon inputs to soils, but it remains difficult to predict the response of soil processes to anthropogenic change. Models that predict soil C-cycling as a function of ecosystem properties do not explain well variation in soil processes. Our difficulty in predicting the response of soil processes to environmental change suggests a need for a greater understanding of the biotic mechanisms that govern the soil C-cycle. Changes in microbial community structure and function have been proposed to impact soil C-cycling. However, our ability to predict the impacts of these changes on terrestrial ecosystems is constrained by our limited understanding of mechanisms that drive microbial processes in soil systems. We have applied a newly developed approach of Microbial Food Web Mapping to chart the carbon assimilation dynamics of microbial taxa in soils. This method made it possible to track different classes of 13C-labeled substrates into the DNA of nearly every member of the soil community simultaneously over time. Microbial Food Web Mapping was used to identify and characterize microbial taxa that fill key roles in the soil carbon cycle. We mapped the assimilation of carbon into thousands of discrete soil microorganisms as a function of soil carbon content and pH. In so doing we characterized assimilatory carbon metabolism, an essential component of soil C-cycling, for a vast array of uncultivated soil microorganisms. Further, we characterized the spatial and temporal dynamics of these microorganisms across a series of sites representing variation in soil characteristics. This project revealed fundamental aspects of soil C-cycling and provided ecological and metabolic insights on diverse uncultivated soil organisms that are widespread and play major roles in the global C-cycle. The genetic capacity of microbial communities can be studied through ‘omic approaches but it remains difficult to make direct links between the genetic capacity of microorganisms and their function in the soil C-cycle. Microbial Food Web Mapping makes it possible to link gene sequences to soil C-cycle processes as they occur in soil. This approach allows us to characterize the activity of non-cultivated microorganisms in a range of terrestrial systems. This data will be used to build a base of information about the role of non-cultivated organisms in critical C-cycle processes in terrestrial ecosystems and will provide insight on the manner in which soil communities metabolize soil organic matter. The results generated by this project will improve our ability to examine the impacts of management decisions, soil history, and environmental change on the behavior of microbial communities in terrestrial ecosystems, revealing the ecological mechanisms by which microbes regulate both C mineralization and carbon retention in soils, and improving our ability to predict changes in terrestrial ecosystem processes in the face of accelerating global change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Allometric-kinetic model predictions of concentration ratios for anthropogenic radionuclides across animal classes and food selection

Protection of the environment from radiation fundamentally relies on dose assessments for non-human biota. These dose assessments use measured or predicted concentrations of radionuclides in soil or water combined with Concentration Ratios (CRs) to estimate tissue concentrations in animals and plants, yet there is a paucity of CR data relative to the vast number of potential taxa and radioactive contaminants in the environment and their taxon-specific ecosystems. Because there are many taxa each having very different behaviors and biology and many possible bioavailable radionuclides, CRs have the potential to vary by orders-of-magnitude, as often seen in published data. Given the diversity of taxa, the International Commission on Radiological Protection (ICRP) has categorized non-human biota into 12 representative organisms (ROs), while the U.S. Department of Energy (DOE) uses the non-taxon specific categories of terrestrial, riparian, and aquatic animals. The question we examine here is: do these categorizations adequately “represent” the broad diversity of animals? To explore this question, here we analyze an Allometric-Kinetic (A-K) model to calculate radionuclide-specific CRs for common animal classes which are then further subcategorized into herbivores, omnivores, carnivores, and invertebrate detritivores. Comparisons in CRs among animal classes exhibited only small differences, but there was order of magnitude differences between herbivores, carnivores, and especially detritivores, for many radionuclides of interest. These findings suggest that the ICRP ROs and the DOE categories are reasonable, but their accuracy could be improved by including sub-categories related to animal dietary ecology. Finally, comparisons of A-K model predicted CR values to published CRs show order-of-magnitude variations, providing justification for additional studies of assimilation across radionuclides and animal classes.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
NASA NTRSDate not supplied

Status and Progress of All-Sky Hyperspectral IR Radiance Assimilation in GEOS

The majority of hyperspectral infrared radiance observations that peak in the lower and mid-troposphere have been excluded from the assimilation in the operational GEOS by the cloud detection scheme. The significance and challenges of assimilating these cloud-affected observations have attracted the attention of the research community, encouraged by the progress in utilizing all-sky microwave radiances affected by clouds and precipitation in many data assimilation systems, including GEOS. While the Global Modeling and Assimilation Office (GMAO) has made significant progress on the assimilation of cloud-cleared infrared radiances, this parallel study marks the first effort to directly assimilate cloud-affected hyperspectral infrared radiances in the GEOS. This initial study is focused on CrIS-FSR wa:wqter-vapor channels. The capability of the GEOS model and CRTM in simulating infrared radiances affected by clouds has been evaluated, and the all-sky brightness temperature Jacobians with respect to temperature, specific humidity, and hydrometeors are examined closely in varying cloud conditions. With hydrometeor control variables in the GEOS, several important aspects of all-sky IR radiance assimilation are investigated. The symmetric cloud effects are assessed for their robustness when used as the cloud proxy in observation error modeling that incorporates the inter-channel correlations with error situation dependency on the amount of cloud. The effective radius, bias correction and quality control processes will also be adapted to accommodate modifications necessary for all-sky infrared assimilation before comprehensive four-dimensional ensemble-variational(4dEnVar) experiments are conducted to assess their impacts on the analysis and forecast performance of GEOS.

Wei Gu↗
NASA NTRS2014

NASA Team 2 Sea Ice Concentration Algorithm Retrieval Uncertainty

Satellite microwave radiometers are widely used to estimate sea ice cover properties (concentration, extent, and area) through the use of sea ice concentration (IC) algorithms. Rare are the algorithms providing associated IC uncertainty estimates. Algorithm uncertainty estimates are needed to assess accurately global and regional trends in IC (and thus extent and area), and to improve sea ice predictions on seasonal to interannual timescales using data assimilation approaches. This paper presents a method to provide relative IC uncertainty estimates using the enhanced NASA Team (NT2) IC algorithm. The proposed approach takes advantage of the NT2 calculations and solely relies on the brightness temperatures (TBs) used as input. NT2 IC and its associated relative uncertainty are obtained for both the Northern and Southern Hemispheres using the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) TB. NT2 IC relative uncertainties estimated on a footprint-by-footprint swath-by-swath basis were averaged daily over each 12.5-km grid cell of the polar stereographic grid. For both hemispheres and throughout the year, the NT2 relative uncertainty is less than 5%. In the Southern Hemisphere, it is low in the interior ice pack, and it increases in the marginal ice zone up to 5%. In the Northern Hemisphere, areas with high uncertainties are also found in the high IC area of the Central Arctic. Retrieval uncertainties are greater in areas corresponding to NT2 ice types associated with deep snow and new ice. Seasonal variations in uncertainty show larger values in summer as a result of melt conditions and greater atmospheric contributions. Our analysis also includes an evaluation of the NT2 algorithm sensitivity to AMSR-E sensor noise. There is a 60% probability that the IC does not change (to within the computed retrieval precision of 1%) due to sensor noise, and the cumulated probability shows that there is a 90% chance that the IC varies by less than +/-3%. We also examined the daily IC variability, which is dominated by sea ice drift and ice formation/melt. Daily IC variability is the highest, year round, in the MIZ (often up to 20%, locally 30%). The temporal and spatial distributions of the retrieval uncertainties and the daily IC variability is expected to be useful for algorithm intercomparisons, climate trend assessments, and possibly IC assimilation in models.

Sea Ice↗
DOE OSTI2020

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Global Geophysical Fluids Center of IERS

The Global Geophysical Fluids Center (GGFC) and its seven Special Bureaus (SB, for Atmosphere, Oceans, Tides, Hydrology, Mantle, Core and Gravity/Geocenter) were establishes by the International Earth Rotation Service in 1998, to support global geodynamic research. Mass transports in the geophysical fluids of the Earth system will cause observable geodynamic effects on a broad time scale.These include (1) variations in the solid Earth's rotation (in length-of-day and polar motion/nutation) via the conservation of angular momentum and effected by torques at the fluid-solid Earth interface; (2) changes in the global gravitational field according to Newton's gravitational law; and (3) motion in the center of mass of the solid Earth relative to that of the whole Earth ("geocenter") via the conservation of linear momentum. These minute signals have become observable by space geodetic techniques, primarily VLBI, SLR, GPS, and DORIS, and new exciting data will be available by space gravity, altimetry, SAR, and magnetic missions. In this sense the precise space geodetic techniques have become effective means of remote sensing of global mass transports. The GGFC and its SBs have the responsibility of supporting, facilitating, and providing services to the worldwide research community in the related research areas. We compute, analyze, compare, archive, and disseminate the time series of the angular momenta and the related torques, gravitational coefficients, and geocenter shift for all geophysical fluids, based on global observational data, and/or products from state-of-the-art models some of which assimilate such data. The computed quantities, algorithm and data formats are standardized. This paper reviews our activities, reports the status, and looks forward into the future.

Chao, Benjamin F.↗
NASA NTRS2007

Salinity Remote Sensing and the Study of the Global Water Cycle

The SMOS and AquariusISAC-D satellite missions will begin a new era to map the global sea surface salinity (SSS) field and its variability from space within the next twothree years. They will provide critical data needed to study the interactions between the ocean circulation, global water cycle and climate. Key scientific issues to address are (1) mapping large expanses of the ocean where conventional SSS data do not yet exist, (2) understanding the seasonal and interannual SSS variations and the link to precipitation, evaporation and sea-ice patterns, (3) links between SSS and variations in the oceanic overturning circulation, (4) air-sea coupling processes in the tropics that influence El Nino, and (4) closing the marine freshwater budget. There is a growing body of oceanographic evidence in the form of salinity trends that portend significant changes in the hydrologic cycle. Over the past several decades, highlatitude oceans have become fresher while the subtropical oceans have become saltier. This change is slowly spreading into the subsurface ocean layers and may be affecting the strength of the ocean's therrnohaline overturning circulation. Salinity is directly linked to the ocean dynamics through the density distribution, and provides an important signature of the global water cycle. The distribution and variation of oceanic salinity is therefore attracting increasing scientific attention due to the relationship to the global water cycle and its influence on circulation, mixing, and climate processes. The oceans dominate the water cycle by providing 86% of global surface evaporation (E) and receiving 78% of global precipitation (P). Regional differences in E-P, land runoff, and the melting or freezing of ice affect the salinity of surface water. Direct observations of E-P over the ocean have large uncertainty, with discrepancies between the various state-of-the-art precipitation analyses of a factor of two or more in many regions. Quantifying the climatic influence of the oceanic water cycle requires more accurately resolving the net air-sea water flux. Measuring global SSS trends on seasonal to interannual timescales by satellite is fundamental to this problem because the SSS trends represent detectable time-integrated signals of the variable marine hydrological cycle. Satellite measurements, coupled with an array of in situ observations, will provide global synoptic SSS fields for the first time history. These data will provide a strong constraint on climate models and data assimilation efforts, which must properly represent the freshwater budget in terms of E-P, ocean advection and surface layer mixing in order to accurately simulate the true ocean state. The SSS fields will allow us to quantify the covariability between the SSS and the strong seasonal E-P cycle in the tropics and high latitudes. Field measurement campaigns to exploit satellite and in situ measurements to close the seasonal E-P cycle over an ocean region are being considered. Lastly the satellite systems will monitor and trace the large long-lived SSS anomalies from year to year that have the potential to influence El Nino and the large scale ocean circulation.

Lagerloef, G. S. E.↗
DOE OSTI2023

Extending Parsimonious Bayesian Inference

Parsimonious Bayesian inference is a theoretical framework for efficient data assimilation that seeks to balance increased consistency between predictions and training data against corresponding increases in model complexity. Within this framework, over-training is understood as optimization that encodes excessive information within model parameters while only achieving small improvements between predictions and training data. This project aims to develop practical methods of limiting excess model information during optimization. One key observation is that practical heuristics for parsimonious learning in high-dimensions must balance expressivity, i.e. the ability of the model to capture diverse predictions with only a few non-zero parameters, against discoverability, i.e. the ability to train the model with gradient-based optimization and drive parameters to low information states. As such, we developed logical activation functions that are able to adaptively approximate arbitrary truth tables that define Boolean logic operations within a probabilistic framework. These functions have demonstrated the ability to learn exclusive disjunction (XOR) and conditioned disjunction (if [condition] then [result_if_true] else [result_if_false]) within a single layer of a neural network. To efficiently exploit these activation functions to drive parsimonious learning required several other advances within the domain of variational inference. The most efficient form of complexity suppression is structured sparsification, driving most model parameters to zero while achieving the structural coherence among nonzeros needed for bandwidth reduction. Such models are not only far more efficient at suppressing information-theoretic complexity, they also reduce the other forms of complexity (computations, communication, storage, and the number of dependencies needed to evaluate predictions). Aiming to support enhanced sparsification, this project examined new approaches to high-dimensional variational inference that allow us to calibrate and control parameter uncertainty during optimization. By identifying which parameters can sustain sparsifying perturbations with little impact on prediction quality, we can develop better pruning strategies by framing them as approximate Bayesian inference. These advances also open paths to mitigate concerns with deploying advanced learning methods in resource-constrained environments, such as running models on power-limited or communication-limited devices.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗
NASA NTRS2010

Using Leaf Chlorophyll to Parameterize Light-Use-Efficiency Within a Thermal-Based Carbon, Water and Energy Exchange Model

Chlorophylls absorb photosynthetically active radiation and thus function as vital pigments for photosynthesis, which makes leaf chlorophyll content (C(sub ab) useful for monitoring vegetation productivity and an important indicator of the overall plant physiological condition. This study investigates the utility of integrating remotely sensed estimates of C(sub ab) into a thermal-based Two-Source Energy Balance (TSEB) model that estimates land-surface CO2 and energy fluxes using an analytical, light-use-efficiency (LUE) based model of canopy resistance. The LUE model component computes canopy-scale carbon assimilation and transpiration fluxes and incorporates LUE modifications from a nominal (species-dependent) value (LUE(sub n)) in response to short term variations in environmental conditions, However LUE(sub n) may need adjustment on a daily timescale to accommodate changes in plant phenology, physiological condition and nutrient status. Day to day variations in LUE(sub n) were assessed for a heterogeneous corn crop field in Maryland, U,S.A. through model calibration with eddy covariance CO2 flux tower observations. The optimized daily LUE(sub n) values were then compared to estimates of C(sub ab) integrated from gridded maps of chlorophyll content weighted over the tower flux source area. The time continuous maps of daily C(sub ab) over the study field were generated by focusing in-situ measurements with retrievals generated with an integrated radiative transfer modeling tool (accurate to within +/-10%) using at-sensor radiances in green, red and near-infrared wavelengths acquired with an aircraft imaging system. The resultant daily changes in C(sub ab) within the tower flux source area generally correlated well with corresponding changes in daily calibrated LUE(sub n) derived from the tower flux data, and hourly water, energy and carbon flux estimation accuracies from TSEB were significantly improved when using C(sub ab) for delineating spatio-temporal variations in LUE(sub n). The results demonstrate the synergy between thermal infrared and shortwave reflective wavebands in producing valuable remote sensing data for operational monitoring of carbon and water fluxes.

Houlborg, Rasmus↗
DOE OSTI2021

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79–86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. Furthermore, these findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

59 BASIC BIOLOGICAL SCIENCES↗