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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

The Ionospheric Three-Dimensional Electron Density Variations Induced by the 21 August 2017 Total Solar Eclipse by Using Global Ionospheric Specification

Global Ionospheric Specification (GIS) is based on the Gauss–Markov Kalman filter to assimilate the slant total electron content (TEC) observed from ground-based GPS receivers and space-based radio occultation instrumentations in order to reconstruct three-dimensional (3D) ionospheric electron density structure, and it can remotely sense and monitor the weather condition in space. In this study, five minutes of high temporal resolution GIS is implemented in order to reconstruct the 3D electron density structure on the 21 August 2017 total solar eclipse and analyze the variations induced by the moon’s shadow. To obtain more information of the ionosphere, from the extend 2200 GPS stations on the continental United States, are added for assimilation. The results show the ionosphere peak height (hmF2) uplift was 30–50 km altitude in latitude 25–40°N, and that the electron density depletion at higher altitudes (400 km) has a more noticeable time delay than at low altitudes (200 km), especially in low-latitude regions.

Ionosphere↗
NASA NTRS1980

Iron and titanium distribution on the moon from orbital gamma ray spectrometry with implications for crustal evolutionary models

A set of Fe and Ti maps and regional values are obtained from the Apollo 15 and 16 orbital gamma ray data by energy band analysis. High-Ti basalts predominate the early and late stages of mare volcanism with high-Fe basaltic volcanism in the interim. The first evidence of a high-Ti-KREEP basalt association is found in the Aristarchus region. A N-S asymmetry for Fe and Ti in the east limb and farside highlands complicates the E-W asymmetry for Th but substantiates crustal inhomogeneity. The observed crustal inhomogeneity adds an additional objection to the primitive source model for crustal evolution. The high-Ti-KREEP basalt association and the general trend of decreasing mare basalt Ti with time lend support to the cumulate source model; however, this model cannot account for young, high-Ti maria. The dynamic assimilation model better accounts for chemical variations observed on the moon.

Davis, P. A., Jr.↗
NASA NTRS2006

Detection of Ice Polar Stratospheric Clouds from Assimilation of Atmospheric Infrared Sounder Data

A novel technique is presented for detection of ice polar stratospheric clouds (PSCs) that form at extremely low temperatures in the lower polar stratosphere during winter. Temperature is a major factor in determining abundance of PSCs, which in turn provide surfaces for heterogeneous chemical reactions leading to ozone loss and radiative cooling. The technique infers the presence of ice PSCs using radiances from the Atmospheric Infrared Sounder (AIRS) in the Goddard Earth Observing System version 5 (GEOS-5) data assimilation system. Brightness temperatures are computed from short-term GEOS-5 forecasts for several hundred AIRS channels, using a radiation transfer module. The differences between collocated AIRS observations and these computed values are the observed-minus-forecast (O-F) residuals in the assimilation system. Because the radiation model assumes clear-sky conditions, we hypothesize that these O-F residuals contain quantitative information about PSCs. This is confirmed using sparse data from the Polar Ozone and Aerosol Measurement (POAM) III occultation instrument. The analysis focuses on 0-F residuals for the 6.79pm AIRS moisture channel. At coincident locations, when POAM III detects ice clouds, the AIRS O-F residuals for this channel are lower than -2K. When no ice PSCs are evident in POAM III data, the AIRS 0-F residuals are larger. Given this relationship, the high spatial density of AIRS data is used to construct maps of regions where 0-F residuals are lower than -2K, as a proxy for ice PSCs. The spatial scales and spatio-temporal variations of these PSCs in the Antarctic and Arctic are discussed on the basis of these maps.

Stajner, Ivanka↗
NASA NTRS1994

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗
NASA NTRS1987

Variational four-dimensional analysis using quasi-geostrophic constraints

A variational four-dimensional analysis technique using quasi-geostrophic models as constraints is examined using gridded fields as data. The analysis method uses a standard iterative nonlinear minimization technique to find the solution to the constraining forecast model which best fits the data as measured by a predefined functional. The minimization algorithm uses the derivative of the functional with respect to each of the initial condition values. This derivative vector is found by inserting the weighted differences between the model solution and the inserted data into a backwards integrating adjoint model. The four-dimensional analysis system was examined by applying it to fields created from a primitive equations model forecast and to fields created from satellite retrievals. The results show that the technique has several interesting characteristics not found in more traditional four-dimensional assimilation techniques. These features include a close fit of the model solution to the observations throughout the analysis interval and an insensitivity to the frequency of data insertion or the amount of data. The four-dimensional analysis technique is very versatile and can be extended to more complex problems with little theoretical difficulty.

Derber, John C.↗
NASA NTRS2006

A Test of Sensitivity to Convective Transport in a Global Atmospheric CO2 Simulation

Two approximations to convective transport have been implemented in an offline chemistry transport model (CTM) to explore the impact on calculated atmospheric CO2 distributions. GlobalCO2 in the year 2000 is simulated using theCTM driven by assimilated meteorological fields from the NASA s Goddard Earth Observation System Data Assimilation System, Version 4 (GEOS-4). The model simulates atmospheric CO2 by adopting the same CO2 emission inventory and dynamical modules as described in Kawa et al. (convective transport scheme denoted as Conv1). Conv1 approximates the convective transport by using the bulk convective mass fluxes to redistribute trace gases. The alternate approximation, Conv2, partitions fluxes into updraft and downdraft, as well as into entrainment and detrainment, and has potential to yield a more realistic simulation of vertical redistribution through deep convection. Replacing Conv1 by Conv2 results in an overestimate of CO2 over biospheric sink regions. The largest discrepancies result in a CO2 difference of about 7.8 ppm in the July NH boreal forest, which is about 30% of the CO2 seasonality for that area. These differences are compared to those produced by emission scenario variations constrained by the framework of Intergovernmental Panel on Climate Change (IPCC) to account for possible land use change and residual terrestrial CO2 sink. It is shown that the overestimated CO2 driven by Conv2 can be offset by introducing these supplemental emissions.

Bian, H.↗
NASA NTRS2004

Alternative Approaches to Land Initialization for Seasonal Precipitation and Temperature Forecasts

The seasonal prediction system of the NASA Global Modeling and Assimilation Office is used to generate ensembles of summer forecasts utilizing realistic soil moisture initialization. To derive the realistic land states, we drive offline the system's land model with realistic meteorological forcing over the period 1979-1993 (in cooperation with the Global Land Data Assimilation System project at GSFC) and then extract the state variables' values on the chosen forecast start dates. A parallel series of forecast ensembles is performed with a random (though climatologically consistent) set of land initial conditions; by comparing the two sets of ensembles, we can isolate the impact of land initialization on forecast skill from that of the imposed SSTs. The base initialization experiment is supplemented with several forecast ensembles that use alternative initialization techniques. One ensemble addresses the impact of minimizing climate drift in the system through the scaling of the initial conditions, and another is designed to isolate the importance of the precipitation signal from that of all other signals in the antecedent offline forcing. A third ensemble includes a more realistic initialization of the atmosphere along with the land initialization. The impact of each variation on forecast skill is quantified.

Koster, Randal↗
NASA NTRS2018

Studying Diurnal Variations of Aerosols with NASA MERRA-2 Reanalysis Data

Aerosols play an important role in atmospheric dynamics, climate variations, and Earth's energy cycle by altering the radiation balance in the atmosphere through interaction with clouds, providing fertilizer for forests and canopy, and as a supply of iron to the ocean over long time periods. Studies suggest that much of the feedback between dust aerosols and dynamics is associated with diurnal and synoptic scale variability. However, the lack of sub-daily resolution of aerosols from satellite observations makes it difficult to study the diurnal characteristics, especially over tropical and subtropical regions. Investigation of this topic utilizes over 37 years of simulated global aerosol products from NASA atmospheric reanalysis, in the second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) data set, available from NASA Goddard Earth Science Data and Information Services Center (GES DISC). MERRA-2 covers the period 1980-present, and is continuing as an ongoing climate analysis. Aerosol assimilation is included throughout the period, using data from MODIS, MISR, AERONET, and AVHRR (in the pre-EOS period). The aerosols are assimilated using the MERRA-2 aerosol model, which interacts directly with radiation parameterization, and is radiatively coupled with atmospheric model dynamics in the Goddard Earth Observing System Model, Version 5 (GEOS-5). Hourly, monthly, and monthly diurnal data are available at spatial resolution of 0.5o x 0.625o (latitude x longitude). By using MERRA-2 hourly and monthly diurnal products, different aerosol diurnal variabilities are observed over North America, Africa, Asia, and Australia, that may be due to different meteorological conditions and aerosol sources. The presentation will also provide an overview of MERRA-2 data services at GES DISC, such as how to find and download data, and how to quickly visualize and analyze data online with Giovanni.

aerosol↗
NASA NTRS2016

Maps Suggest Transport and Source Processes of PM2.5 at 1 km x 1 km for the Whole San Joaquin Valley, Winter 2011 (Generalizations from DISCOVER-AQ)

We present interpreted data analysis using MAIAC (Multiangle implementation of Atmospheric Correction) retrievals and appropriate RAPid Update Cycle (RAP) meteorology to map respirable aerosol (PM2.5) for the period January and February, 2011. The San Joaquin Valley is one of the unhealthiest regions in the USA for PM2.5 and related morbidity. The methodology evaluated can be used for the entire moderate-resolution imaging spectrometer (MODIS, VIIRS) data record. Other difficult areas of the West: Riverside, CA, Salt Lake City, UT, and Doa Ana County, NM share similar difficulties and solutions. The maps of boundary layer depth for 1116 hr local time from RAP allows us to interpret aerosol optical thickness as a concentration of particles in a nearly well-mixed box capped by clean air. That mixing is demonstrated by DISCOVER-AQ data and afternoon samples from the airborne measurements, P3B (on-board) and B200 (HSRL2 lidar). This data and the PM2.5 gathered at the deployment sites allowed us to estimate and then evaluate consistency and daily variation of the AOT to PM2.5 relationship. Mixed-effects modeling allowed a refinement of that relation from day to day; RAP mixed layers explained the success of previous mixed-effects modeling. Compositional, size-distribution, and MODIS angle-of-regard effects seem to describe the need for residual daily correction beyond ML depth.We report on an extension method to the entire San Joaquin Valley for all days with MODIS imagery using the permanent PM2.5 stations, evaluated for representativeness. Resulting map movies show distinct sources, particularly Interstate-5 (at approx. 1km x 1km resolution) and the broader Bakersfield area. Accompanying winds suggest transport effects and variable pathways of pollution cleanout. Such estimates should allow morbiditymortality studies. They should be also useful for actual model assimilations, where composition and sources are uncertain. We conclude with a description of new work to extend these insights to similar regions, e.g. interior valleys of California, the Po Valley, the Mediterranean litoral, and the Ganges Plain.This work show generalizable use of remote sensing, a major goal of DISCOVER-AQ, Deriving Information on Surface Conditions from COlumn and VERtically Resolved Observations Relevant to Air Quality.

Chatfield, R.↗
NASA NTRS2001

Highly Siderophile Element Abundances in Martian Meteorites

Critical evaluation of new and literature data for highly siderophile elements (HSE) in Martian (SNC) meteorites allows several first order conclusions to be drawn. (i) Re concentrations in SNC meteorites are nearly constant (within a factor of two) and do not correlate with rock type. Exceptions to this rule are Chassigny and Dar al Gani (DaG) 476, both of which are inferred to have experienced terrestrial Re contamination. (ii) Fractionations between Rh and Pd are small. Excluding Shergotty, the Rh/Pd ratio of the SNC suite is 0.22\pm0.05. (iii) Os and Ir contents vary by about four orders of magnitude; and positive correlations with MgO, Cr, and Ni suggest that these variations are not controlled by sulfide fractionation. A possible exception is the orthopyroxenite ALH84001, whose HSE's (including Ni, which is compatible in opx) are very low. (iv) Zagami, Shergotty, and Nakhla have nearly identical HSE signatures. Shergotty and Zagami have experienced assimilation-fractional crystallization (AFC) and have "crustal" Sr and Nd isotopic signatures. Conversely, the Nakhla parent was a small degree partial melt of a depleted mantle that interacted little with the Martian crust. These observations suggest that "evolved" HSE signatures can be produced by either fractional crystallization or small degrees of partial melting. (v) Chassigny and other mafic SNC's have HSE signatures that are very distinct from those of Nakhla-Zagami-Shergotty. The HSE elemental ratios of mafic SNC's approach chondritic, implying that the Martian mantle has nearly chondritic relative abundances of the HSE's. (vi) This chondritic HSE signature is observed in SNC's of various ages, suggesting that this is an ancient feature that has not evolved over time. (vii) No correlation is observed between HSE's and signatures of crustal contamination (e.g., Sr isotopes), indicating that the HSE signatures of the SNC suite are not derived from the crust. (vii) The Ru/Pd for the SNC suite ratio is about 0.5 chondritic. We attribute this depletion of Ru to the presence of spinel in the SNC source regions. (viii) Assignment of this depletion of Ru to spinel seems to also imply that spinel is not the cause of Os and Ir fractionations, requiring some other mafic phase to account for the observed variations. -----.

Jones, J. H.↗
NASA NTRS1979

The origin of KREEP

KREEP is a lunar material having very high concentrations of incompatible elements; its name is an acronym for the incompatibles K, rare-earth elements (REE), and P. Although a few pristine (endogenously igneous) KREEPy samples were returned from the Apollo 15 and 17 sites, most KREEPy samples are polymict breccias. Most models of KREEP petrogenesis have been based on partial melting of a variety of sources. Such models fail to explain the veritable absence of variations in incompatible element patterns over the sampled portion of the moon. We have defined a KREEP component based on the average composition of Apollo 14 breccias having extremely high concentrations of incompatible elements. Normalization of accurate incompatible data for KREEPy samples from the Apollo 12, 14, 15, 16, and 17 sites to this component virtually always shows no resolvable fractionation (e.g., <10% variation in the La/Lu ratio), whereas partial melting models typically produce larger fractionations (±20–25% in La/Lu) from a factor of 2 difference in degree of partial melting. Required is a single major source that could provide KREEP to widely separated locations on the nearside of the moon. The anorthositic crust of the moon is commonly attributed to the flotation of plagioclase on a deep, moon-wide magma ocean. Fractional crystallization of this magma ocean would have produced large enrichments of incompatibles in a residual liquid. No other plausible major source of incompatibles has been proposed. We borrow the German prefix ur—meaning primeval and designate this residual liquid ‘urKREEP.’ We propose that all KREEPy rocks originated by dilution of urKREEP with crustal or mantle materials during assimilation, or zone-refining (pristine samples), or impact-induced brecciation (breccias and melt rocks). The formation of urKREEP cannot be dated precisely. Correction of breccia Rb-Sr model ages for Rb loss or gain during the early intense bombardments yields ages that cluster in the range 4.4–4.5 Gy. This implies that crystallization of the magma ocean was essentially complete at this time and is in general agreement with U-Pb evidence indicating crustal formation at 4.4 Gy. Assuming that the moon had the composition of an H-group chondrite depleted in Fe-Ni and FeS and that half the incompatibles fractionated into materials other than urKREEP, the thickness of a moon-wide urKREEP layer was <2 km. Thorium concentrations determined by gamma ray spectroscopy indicate that about 4% of the incompatibles in an H chondritic moon are now in the outermost kilometer.

Paul H. Warren↗
NASA NTRS1999

Monitoring Global Geophysical Fluids by Space Geodesy

Since its establishment on 1/1/1998 by the International Earth Rotation Service, the Coordinating Center for Monitoring Global Geophysical Fluids (MGGF) and its seven Special Bureaus have engaged in an effort to support and facilitate the understanding of the geophysical fluids in global geodynamics research. Mass transports in the atmosphere-hydrosphere-solid Earth-core system (the "global geophysical fluids") will cause the following geodynamic effects on a broad time scale: (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, with ever increasing precision/accuracy and temporal/spatial resolution. Each of the seven Special Bureaus within MGGF is responsible for calculations related to a specific Earth component or aspect -- Atmosphere, Ocean, Hydrology, Ocean Tides, Mantle, Core, and Gravity/Geocenter. Angular momenta and torques, gravitational coefficients, and geocenter shift will be computed for geophysical fluids based on global observational data, and from state-of-the-art models, some of which assimilate such data. The computed quantities, algorithm and data formats are standardized. The results are archived and made available to the scientific research community. This paper reports the status of the MGGF activities and current results.

Chao, Benjamin F.↗
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.↗
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.↗
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↗