Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14

GCM Simulation of the Large-scale North American Monsoon Including Water Vapor Tracer Diagnostics

The geographic sources of water for the large-scale North American monsoon in a GCM are diagnosed using passive constituent tracers of regional water'sources (Water Vapor Tracers, WVT). The NASA Data Assimilation Office Finite Volume (FV) GCM was used to produce a 10-year simulation (1984 through 1993) including observed sea surface temperature. Regional and global WVT sources were defined to delineate the surface origin of water for precipitation in and around the North American i'vionsoon. The evolution of the mean annual cycle and the interannual variations of the monsoonal circulation will be discussed. Of special concern are the relative contributions of the local source (precipitation recycling) and remote sources of water vapor to the annual cycle and the interannual variation of warm season precipitation. The relationships between soil water, surface evaporation, precipitation and precipitation recycling will be evaluated.

Bosilovich, Michael G.↗

Online Simulations and Forecasts of the Global Aerosol Distribution in the NASA GEOS-5 Model

We present an analysis of simulations of the global aerosol system in the NASA GEOS-5 transport, radiation, and chemistry model. The model includes representations of all major tropospheric aerosol species, including dust, sea salt, black carbon, particulate organic matter, and sulfates. The aerosols are run online for the period 2000 through 2005 in a simulation driven by assimilated meteorology from the NASA Goddard Data Assimilation System. Aerosol surface mass concentrations are compared with existing long-term surface measurement networks. Aerosol optical thickness is compared with ground-based AERONET sun photometry and space-based retrievals from MODIS, MISR, and OMI. Particular emphasis is placed here on consistent sampling of model and satellite aerosol optical thickness to account for diurnal variations in aerosol optical properties. Additionally, we illustrate the use of this system for providing chemical weather forecasts in support of various NASA and community field missions.

Colarco, Peter↗

DAO's Next Generation Physical-Space/Finite-Volume Data Assimilation System: Formulation and Initial Evaluation

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR, and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. In this talk we will describe the general system formulation, the adaptive quality control and general aspects of the error covariance modeling. The NASA-NCAR GCM is a completely new model which replaces the CEOs GCM used in the previous GEOS-1/2/3 Data Assimilation systems. A particular configuration of adaptive Statistical Quality Control and the Physical-space Statistical Analysis System (PSAS) are currently implemented in DAO's operational Data Assimilation System. However, the unique finite-volume formulation of the NASA-NCAR GCM, combined with the generality of the observation-space formulation of PSAS, provides for a very simple and accurate model-analysis interface. The system assimilates a variety of conventional and satellite observations. In particular, TOVS Level 1B radiances are assimilated using a 1-D variational scheme, both in clear sky and cloudy conditions. Computationally, the fvDAS runs approximately 10 times faster than the operational GEOS-Terra system. We will show that the next-generation fvDAS has much improved observation-minus-6hr forecast (O-F) statistics, as well as 5-day forecast skills. Top of the atmosphere radiation fields are in closer agreement with CERES measurements, with realistic precipitation and moisture fields. We will also show that the finite-volume formulation of the fvDAS produce assimilated fields which are more suitable for driving constituent transport models.

daSilva, A.↗

Rationalizing Burned Carbon with Carbon Monoxide Exported from South America

We present several estimates cross-checking the fluxes of carbon to the atmosphere from burning, comparing models that are based on simple land-surface parameterizations and atmospheric transport dynamics. Both estimates made by NASA Ames and USP modeling techniques are quite high compared to some detailed satellite/land-use studies of emissions. The flux of carbon liberated to the atmosphere via biomass burning is important for several reasons. This flux is a fundamental statistic for the parameterization of the large-scale flux of gases controlling the reactive greenhouse gases methane and ozone. Similarly, it is central to the estimation of the translocation of nitrogen and pyrodenitrification in the tropics. Thirdly, CO2 emitted from rainforest clearing contributes directly to carbon lost from the rainforest system as it contributes to greenhouse gas forcing. While CO2 from pasturage, agriculture, etc, is considered to be reabsorbed seasonally, and so "off budget" for the carbon cycle, it must also be accounted. CO2 anomalies related to daily weather and interannual climatic variation are strong enough to perturb our scientific perception of long-term carbon storage trends. We compare fluxes deduced from land-use statistics (originally, W.M. Hao) and from satellite hot pixels (A. Setzer) with atmospheric fluxes determined by the mesoscale/continental scale models RAMS and MM5, and point to some new work with highly resolved global models (the NASA Data Assimilation Office's GEOS4). Our simulations are tied to events, so that measured tracers like CO tie the models directly to the burning and meteorology of a specific period. We point out a particular sensitivity in estimates based on CO, and indicate how analysis of CO2 along with other biomass-burning tracers may lead to an improved multi-species estimator of carbon burned.

Chatfield, R.↗

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

Classification of Ascension Island and Natal Ozonesondes Using Self-Organizing Maps

Ozone profiles from balloon-borne ozonesondes are used for development of satellite algorithms and in chemistry-climate model initialization, assimilation and evaluation. An important issue in the application of these profiles is how best to treat variations where varying photochemical and dynamical influences can cause the ozone mixing ratio in the tropospheric segments of the profile to change by of a factor of 2-3 within a day. Clustering techniques are an ideal way to approach the statistical classification of profile data and we apply self-organizing maps to tropical tropospheric SHADOZ data, hypothesizing that the data will sort according to various influences on ozone, namely anthropogenic sources like biomass burning, meteorological conditions, and stratospheric or extra-tropical intrusions. Self-organizing maps, that use a learning algorithm to reveal the most prominent features of a data set according to a specified number of clusters, have been determined for the 1998-2009 SHADOZ profiles over Ascension Island (512 profiles, 7.98 deg. S, 14.42 deg. W) and Natal, Brazil (425 profiles, 5.42degS, 35.38degW). The 2 × 2 self-organizing map, which creates 4 clusters, reveals that deviations from the average ozone in the free troposphere include both increased ozone resulting from seasonal biomass burning in Africa and locally reduced ozone brought about by convective lifting of unpolluted boundary-layer air. Expanding to a 4 × 4 self-organizing map shows how biomass burning influences the yearly cycle of tropospheric ozone at Ascension Island and captures the seasonality of ozone at both Ascension Island and Natal. Comparing Ascension Island and Natal using a 4 × 4 self-organizing map at each site reveals similarities in mid-tropospheric ozone, but shows differences in lower-tropospheric ozone due to Ascension Island being closer to African biomass burning and more affected by descent from the mean Walker circulation, with less convective activity, than Natal.

algorithms↗

Observations, Modeling, and Predictions of Solar Activity From the Deep Interior to the Corona

Solar activity is a manifestation of magnetic self-organization processes that involve complex dynamical coupling of various layers of the Sun acting over a broad range of spatial and temporal scales. Synergy of observational, theoretical, and modeling efforts is key to understanding solar activity variation, dynamics, and evolution and to developing reliable physics-based forecasts of long-term solar cycles and short-term activity manifestations, “seasons of solar activity”, such as periods of enhanced flaring and CME activity. The session welcomes observers, modelers, and theoreticians to share their results and ideas and to discuss current challenges, development of emerging fields (such as data assimilation), and the analysis of historical and modern observational data using theoretical modeling, interpretations, and predictions.

Deep Interior↗

Synthesis and Assimilation Systems - Essential Adjuncts to the Global Ocean Observing System

Ocean assimilation systems synthesize diverse in situ and satellite data streams into four-dimensional state estimates by combining the various observations with the model. Assimilation is particularly important for the ocean where subsurface observations, even today, are sparse and intermittent compared with the scales needed to represent ocean variability and where satellites only sense the surface. Developments in assimilation and in the observing system have advanced our understanding and prediction of ocean variations at mesoscale and climate scales. Use of these systems for assessing the observing system helps identify the strengths of each observation type. Results indicate that the ocean remains under-sampled and that further improvements in the observing system are needed. Prospects for future advances lie in improved models and better estimates of error statistics for both models and observations. Future developments will be increasingly towards consistent analyses across components of the Earth system. However, even today ocean synthesis and assimilation systems are providing products that are useful for many applications and should be considered an integral part of the global ocean observing and information system.

Rienecker, Michele M.↗

Attribution of Flux Partitioning Variations Between Land Surface Models over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

uncertainty↗

Attribution of Flux Partitioning Variations Between Land Surface Models Over the Continental U.S.

Accurate quantification of the terrestrial evapotranspiration (ET) components of plant transpiration (T), soil evaporation (E) and evaporation of the intercepted water (I) is necessary for improving our understanding of the links between the carbon and water cycles. Recent studies have noted that, among the modeled estimates, large disagreements exist in the relative contributions of T, E and I to the total ET. As these models are often used in data assimilation environments for incorporating and extending ET relevant remote sensing measurements, understanding the sources of inter-model differences in ET components is also necessary for improving the utilization of such remote sensing measurements. This study quantifies the contributions of two key factors explaining inter-model disagreements to the uncertainty in total ET: (1) contribution of the local partitioning and (2) regional distribution of ET. The analysis is conducted by using outputs from a suite of land surface models in the North American Land Data Assimilation System (NLDAS) configuration. For most of these models, transpiration is the dominant component of the ET partition. The results indicate that the uncertainty in local partitioning dominates the inter-model spread in modeled soil evaporation E. The inter-model differences in T are dominated by the uncertainty in the distribution of ET over the Eastern U.S. and the local partitioning uncertainty in the Western U.S. The results also indicate that uncertainty in the T estimates is the primary driver of total ET errors. Over the majority of the U.S., the contribution of the two factors of uncertainty to the overall uncertainty is non-trivial.

evapotranspiration partitioning; transpiration; so↗

Data correlation and analysis of arc tunnel and wind tunnel tests of RSI joints and gaps, phase 2. Volume 1: Technical report

Heat transfer data measured in gaps representative of those being employed for joints in the space shuttle reusable surface insulation (RSI) thermal protection systems (TPS) were assimilated, analyzed, and correlated. Several types of gap were investigated with emphasis on simple butt joints. Gap widths ranged from 0.0 to 0.76 cm and depths ranged from 1 to 6 cm. Laminar, transitional, and turbulent boundary layer flows over the gap opening were investigated. The angle between gap axis and external flow was varied between 0 and pi/2 radians. The contoured cross section gap performed significantly better than all other wide gaps and slightly better than all other narrow gap geometries. Three dimensional heating variations were observed within gaps in the absence of external flow pressure gradients. Interactions between heating within gaps and heating of adjacent top tile surfaces were observed. Gaps aligned with the flow were observed to promote boundary layer transition. Heat transfer correlation equations were obtained for many of the tests. The TPS thickness requirements with and without gaps were computed for a current shuttle entry trajectory. Experimental data employed in the study are summarized. A description of each test facility, run schedule and test conditions, model descriptive information, and heat flux data are included.

Cristensen, H. E.↗

Stratospheric Temperature Changes: Observations and Model Simulations

This paper reviews observations of stratospheric temperatures that have been made over a period of several decades. Those observed temperatures have been used to assess variations and trends in stratospheric temperatures. A wide range of observation datasets have been used, comprising measurements by radiosonde (1940s to the present), satellite (1979 - present), lidar (1979 - present) and rocketsonde (periods varying with location, but most terminating by about the mid-1990s). In addition, trends have also been assessed from meteorological analyses, based on radiosonde and/or satellite data, and products based on assimilating observations into a general circulation model. Radiosonde and satellite data indicate a cooling trend of the annual-mean lower stratosphere since about 1980. Over the period 1979-1994, the trend is 0.6K/decade. For the period prior to 1980, the radiosonde data exhibit a substantially weaker long-term cooling trend. In the northern hemisphere, the cooling trend is about 0.75K/decade in the lower stratosphere, with a reduction in the cooling in mid-stratosphere (near 35 km), and increased cooling in the upper stratosphere (approximately 2 K per decade at 50 km). Model simulations indicate that the depletion of lower stratospheric ozone is the dominant factor in the observed lower stratospheric cooling. In the middle and upper stratosphere both the well-mixed greenhouse gases (such as CO) and ozone changes contribute in an important manner to the cooling.

Ramaswamy, V.↗

Solar Magnetic Activity and Solar-Stellar Connections

Prediction of solar magnetic activity on various temporal scales is a fundamental element of space weather, which requires a wide range of theoretical and observational expertise in solar phenomena from the deep interior to the corona. Historical observations have revealed many features of cyclic variations of the solar activity; but these data are dramatically insufficient to draw a physical picture of global magnetic field evolution. New observational data, currently available from space missions and ground-based observatories, provide us with detailed information about solar dynamics and magnetism. However, because of the relatively short duration of data series and the great variety of data types and quality, it is challenging to assimilate these data in theoretical models and make reliable forecasts. The recent unexpectedly weak solar activity cycles, as well as observations of rotational and magnetic topology transitions in solar-type stars, suggest that the Sun and its magnetic dynamo are currently in a very interesting evolutionary stage. This could relate to the difficulty in getting a model of the Sun to produce solar-like – rather than anti-solar-like – differential rotation, to reproduce the rotation profile obtained from helioseismology, and to predict solar activity cycles.

magnetic↗

An initialization procedure for assimilating geostationary satellite data into numerical weather prediction models

An attempt was made to offset the limitations of GEO satellites for supplying timely initialization data for numerical weather prediction models (NWP). The NWP considered combined an isentropic representation of the free atmosphere with a sigma-coordinate model for the lower 200 mb. A flux form of the predictive equations described vertical transport interactions at the boundary of the two model domains, thereby accounting for the poor vertical temperature and wind field resolution of GEO satellite data. A variational analysis approach was employed to insert low resolution satellite-sensed temperature data at varying rates. The model vertical resolution was limited to that available from the satellite. Test simulations demonstrated that accuracy increases with the frequency of data updates, e.g., every 0.5-1 hr. The tests also showed that extensive cloud cover negates the capabilities of IR sensors and that microwave sensors will be needed for temperature estimations for 500-1000 mb levels.

Gal-Chen, T.↗

Estimating Net Primary Productivity Using Satellite and Ancillary Data

The net primary productivity (C) or annual rate of carbon accumulation per unit ground area by terrestrial plant communities is the difference of the rate of gross photosynthesis (A(sub g)) and autotrophic respiration (R) per unit ground area. Although available observations show that R is a large and variable fraction of A(sub g), viz., 0.3 to 0.7, it is generally recognized that much uncertainties exist in this fraction due to difficulties associated with the needed measurements. Additional uncertainties arise when these measurements are extrapolated to regional or global land surface using empirical equations, for example, using regression equations relating C to mean annual precipitation and air temperature. Here, a process-based approach has been taken to calculate A(sub g) and R using satellite and ancillary data. A(sub g) has been expressed as a product of radiation use efficiency, magnitude of intercepted photosynthetically active radiation (PAR), and normalized by stresses due to soil water shortage and air temperature away from the optimum range. A biophysical model has been used to determine the radiation use efficiency from the maximum rate of carbon assimilation by a leaf, foliage temperature, and the fraction of diffuse PAR incident on a canopy. All meteorological data (PAR, air temperature, precipitation, etc.) needed for the calculation are derived from satellite observations, while a land use, land cover data (based on satellite and ground measurements) have been used to assess the maximum rate of carbon assimilation by a leaf of varied cover type based on field measurements. R has been calculated as the sum of maintenance and growth components. The maintenance respiration of foliage and live fine roots at a standard temperature of different land cover has been determined from their nitrogen content using field and satellite measurements, while that of living fraction of woody stem (viz., sapwood) from the seasonal maximum leaf area index as determined from satellite observations. These maintenance respiration values were then adjusted to that corresponding to air temperature according to a prescribed non-linear variation of respiration with temperature. The growth respiration has been calculated from the difference of Ag and maintenance respiration, according to the two-compartment model. The results of calculations will be reported for 36 consecutive months (1987-1989) over large contiguous areas (ca. 10(exp 5) sq km) Of agricultural land and tropical humid evergreen forests, and compared with available field data.

Choudhury, B. J.↗

Combined assimilation of hydrography and TOPEX data into an Indian Ocean GCM using the adjoint method

To study seasonal circulation and meridional heat transport of the Indian Ocean by synthesizing dynamics with data, climatological monthly temperatures and salinities, surface heat and freshewater fluxes, and wind stresses, together with monthly ensembles of three years (93-95) of TOPEX-derived surface geostrophic velocity anomalies, are assimilated into an Indian Ocean GCM.

Hydrography Indian Ocean TOPEX meridional heat tra↗

How Accurate is Land/Ocean Moisture Transport Variability in Reanalyses?

Quantifying the global hydrological cycle and its variability across various time scales remains a challenge to the climate community. Direct measurements of evaporation (E), evapotranspiration (ET), and precipitation (P) are not feasible on a global scale, nor is the transport of water vapor over the global oceans and sparsely populated land areas. Expanding satellite data streams have enabled development of various water (and energy) flux products, complementing reanalyses and facilitating observationally constrained modeling. But the evolution of the global observing system has produced additional complications--improvements in satellite sensor resolution and accuracy have resulted in "epochs" of observational quasi-uniformity that can adversely affect reanalysis trends. In this work we focus on vertically integrated moisture flux convergence (VMFC) variations within the period 1979 - present integrated over global land. We show that VMFC in recent reanalyses (e.g. ERA-I, NASA MERRA, NOAA CFSR and JRA55) suffers from observing system changes, though differently in each product. Land Surface Models (LSMs) forced with observations-based precipitation, radiation and near-surface meteorology share closely the interannual P-ET variations of the reanalyses associated with ENSO events. (VMFC over land and P-ET estimates are equivalent quantities since atmospheric storage changes are small on these scales.) But the long-term LSM trend over the period since 1979 is approximately one-fourth that of the reanalyses. Additional reduced observation reanalyses assimilating only surface pressure and /or specifying seasurface temperature also have a much smaller trend in P-ET like the LSMs. We explore the regional manifestation of the reanalysis P-ET / VMFC problems, particularly over land. Both principal component analysis and a simple time series changepoint analysis highlight problems associated with data poor regions such as Equatorial Africa and, for one reanalysis, the Equatorial Andes region. Onset of the availability of passive microwave Special Sensor Microwave Imager (SSMI) moisture data in July 1987 and the transition from the Microwave Sounder Unit (MSU) to an advanced version (AMSU) have significant impacts on VMFC variability. Simple accounting for these errors of leading importance results in modified reanalysis VMFC estimates that agree much better with the LSM results. Regional details of the modified reanalysis VMFC and LSM P-ET are related to changes in Pacific Decadal Variability as manifest in SST changes after the late 1990s.

Robertson, F. R.↗

Correction of Excessive Precipitation Over Steep and High Mountains in a General Circulation Model

Excessive precipitation over steep and high mountains (EPSM) is a well-known problem in GCMs and meso-scale models. This problem impairs simulation and data assimilation products. Among the possible causes investigated in this study, we found that the most important one, by far, is a missing upward transport of heat out of the boundary layer due to the vertical circulations forced by the daytime upslope winds, which are forced by the heated boundary layer on subgrid-scale slopes. These upslope winds are associated with large subgrid-scale topographic variation, which is found over steep and high mountains. Without such subgridscale heat ventilation, the resolvable-scale upslope flow in the boundary layer generated by surface sensible heat flux along the mountain slopes is excessive. Such an excessive resolvablescale upslope flow combined with the high moisture content in the boundary layer results in excessive moisture transport toward mountaintops, which in turn gives rise to EPSM. Other possible causes of EPSM that we have investigated include 1) a poorly-designed horizontal moisture flux in the terrain-following coordinates, 2) the condition for cumulus convection being too easily satisfied at mountaintops, 3) the presence of conditional instability of the computational kind, and 4) the absence of blocked flow drag. These are all minor or inconsequential. We have parameterized the ventilation effects of the subgrid-scale heated-slope-induced vertical circulation (SHVC) by removing heat from the boundary layer and depositing it in layers higher up when the topographic variance exceeds a critical value. Test results using NASA/Goddard's GEOS-S GCM have shown that this largely solved the EPSM problem.

Chao, Winston C.↗