Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean 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 469 records · Page 26

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of surface evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote geographic sources of surface evaporation for precipitation, based on the implementation of three-dimensional constituent tracers of regional water vapor sources (termed water vapor tracers, WVT) in a general circulation model. The major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In the WVT approach, each tracer is associated with an evaporative source region for a prognostic three-dimensional variable that represents a partial amount of the total atmospheric water vapor. The physical processes that act on a WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be predicted within the model simulation, and can be validated against the model's prognostic water vapor. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional sources, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In most North American continental regions, the local source of precipitation is correlated with total precipitation. There is a general positive correlation between local evaporation and local precipitation, but it can be weaker because large evaporation can occur when precipitation is inhibited. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.↗

Ocean Carbon Flux, Transport, and Burial Within the Western and Eastern US Coastal Zones

This project has been to develop and apply a regional. eddy-resolving circulation and biogeochemistry model of both the western and eastern U.S. coastal regions, capable of simulating the processes that control the carbon cycle. Validation has been by statistical comparison with analyses from various satellite measurements, including those from EOS sensors, as well as from in situ measurements. Sensitivity studies were carried out to investigate how the coastal ecosystem and biogeochemical cycles respond to changes in climate, large-scale eutrophication from indus- trial pollution, and other anthropogenic induced changes. The research has been conducted in collaboration with research groups at UCLA. NASA/GSFC (Wallops), Rutgers, and SIO. Overall. the project was focused on several key modeling issues, each of which tie back into completing the primary task of developing a coastal carbon model for both the eastern and western US. coasts. Individual groups within the entire program are still collaborating to address these specific tasks. These include: implementation of the coupled circulation/biogeochemical model within the U.S. West Coast. including high-resolution, embedded subdomains for the Southern California Bight and Monterey Bay region; development of a biogeochemical model with resolved carbon, nitrogen and oxygen cycles; development of data assimilation techniques for use of satellite data sets; reconfiguration of the model domain to U.S. East Coast; development of coastal forcing fields: development of methods to compare the model against remotely sensed data; and, the test of model sensitivity to environmental conditions. Below, we present a summary of the progress made toward achieving these soak. Because this has been a multi-institutional, collaborative effort, we note the groups involved with particular activities.

McWilliams, James C.↗

WCRP Task Team for the Intercomparison of Reanalyses (TIRA): Motivation and Progress

Reanalyses have proven to be an important resource for weather and climate related research, as well as societal applications at large. Several centers have emerged to produce new atmospheric reanalyses in various forms every few years. In addition, land and ocean communities are producing disciplinary uncoupled reanalyses. Current research and development in reanalysis is directed at (1) extending the length of reanalyzed period and (2) use of coupled Earth system models for climate reanalysis. While WCRPs involvement in the reanalyses communities through its Data Advisory Council (WDAC) has been substantial, for example in organizing international conferences on reanalyses, a central team of reanalyses expertise is not in place in the WCRP structure. The differences among reanalyses and their inherent uncertainties are some of the most important questions for both users and developers of reanalyses. Therefore, a collaborative effort to systematically assess and intercompare reanalyses would be a logical progression that fills the needs of the community and contributes to the WCRP mission. The primary charge to the TIRA is to develop a reanalysis intercomparison project plan that will attain the following objectives.1)To foster understanding and estimation of uncertainties in reanalysis data by intercomparison and other means 2)To communicate new developments and best practices among the reanalyses producing centers 3)To enhance the understanding of data and assimilation issues and their impact on uncertainties, leading to improved reanalyses for climate assessment 4)To communicate the strengths and weaknesses of reanalyses, their fitness for purpose, and best practices in the use of reanalysis datasets by the scientific community. This presentation outlines the need for a task team on reanalyses, their intercomparison, the objectives of the team and progress thus far.

WDA↗

Emphasizing Spectrum Management for Sustainable Development Research and Applications in Disaster Management

NASA's spaceborne Earth and Heliospheric Observatories and airborne sensors provide a plethora of measurements. These measurements are used in science research to understand the climatology of our home planet and the solar fluxes and cycle of the only star in our solar system 'Sun' which is critical driver for the retention of life on Earth. Specifically, these measurements help us to understand the water and energy cycle, the carbon cycle, weather and climate, atmospheric chemistry, solar variability, and solid Earth and interior to feed into sophisticated mathematical models to analyze and predict the Earth's behavior as an integrated system. The main thrust of this research is on improving the prediction capability in the areas of weather, long term climate and solid Earth processes, and further help the humanity and future generations in terms of societal benefits in managing natural disasters, sustainability issues and many more. This work is further linked with our contributions in the Global Earth Observing System of Systems (GEOSS) Specifically, the data and knowledge resulting from the Earth observing systems and analytical models of the Earth can be made available for assimilation into decision support systems to serve society for disaster management. Through partnerships with national and international agencies and organizations, NASA's Science Mission Directorate's, Applied Sciences Program contributes to benchmarking practical uses of observations and predictions from Earth science remote sensing systems research. The objective is to establish innovative solutions using Earth observations and science information to provide decision support that can be adapted in applications of national and international priority. We along with the international community will continue this critical field of investigation by using our existing and future sensors from space, airborne and insitue environment. In our quest to expanding our knowledge, there will be a need for deploying additional sensors to obtain high spatial, temporal and spectral resolution measurements. These sensors operate in multiple spectral band ranging from UV, visible, infrared, microwave and radio frequency ranges. Of a particular concern is the microwave frequency bands which play a key role in land, ocean, moisture sensing. This is because of a growing commercial demand in the area of high speed broadband communication all over the world, the electronic manufacturers are looking into high frequency microwave spectral bands. This may present a risk to the remote sensing sensors because of additional sources of noise that can impair the highly sensitive passive remote sensing instruments.

Ambrose, Stephen↗

A Coupled Evaluation of Operational MODIS and Model Aerosol Products for Maritime Environments Using Sun Photometry: Evaluation of the Fine and Coarse Mode

Although satellite retrievals and data assimilation have progressed to where there is a good skill for monitoring maritime Aerosol Optical Depth (AOD), there remains uncertainty in achieving further degrees of freedom, such as distinguishing fine and coarse mode dominated species in maritime environments (e.g., coarse mode sea salt and dust versus fine mode terrestrial anthropogenic emissions, biomass burning, and maritime secondary production). For the years 2016 through 2019, we performed an analysis of 550 nm total AOD550, fine mode AOD (FAOD550; also known as FM AOD in the literature), coarse mode AOD (CAOD550), and fine mode fraction (η550) between Moderate Resolution Spectral Imaging Radiometer (MODIS) V6.1 MOD/MYD04 dark target aerosol retrievals and the International Cooperative for Aerosol Prediction (ICAP) core four multi-model consensus (C4C) of analyses/short term forecasts that assimilate total MODIS AOD550. Differences were adjudicated by the global shipboard Maritime Aerosol Network (MAN) and selected island AERONET sun photometer observations with the application of the spectral deconvolution algorithm (SDA). Through a series of conditional and regional analyses, we found divergence included regions of terrestrial influence and latitudinal dependencies in the remote oceans. Notably, MODIS and the C4C and its members, while having good correlations overall, have a persistent +0.04 to +0.02 biases relative to MAN and AERONET for typical AOD550 values (84th% < 0.28), with the C4C underestimating significant events thereafter. Second, high biases in AOD550 are largely associated with the attribution of the fine mode in satellites and models alike. Thus, both MODIS and C4C members are systematically overestimating AOD550 and FAOD550 but perform better in characterizing the CAOD550. Third, for MODIS, findings are consistent with previous reports of a high bias in the retrieved Ångström Exponent, and we diagnosed both the optical model and cloud masking as likely causal factors for the AOD550 and FAOD550 high bias, whereas for the C4C, it is likely from secondary overproduction and perhaps numerical diffusion. Fourth, while there is no wind-speed-dependent bias for surface winds <12 m s−1, the C4C and MODIS AOD550s also overestimate CAOD550 and FAOD550, respectively, for wind speeds above 12 m/s. Finally, sampling bias inherent in MAN, as well as other circumstantial evidence, suggests biases in MODIS are likely MODIS and the C4C products have their own strengths and challenges for a given climate application and discuss needed research. even larger than what was diagnosed here. We conclude with a discussion on how

Aerosols↗

Production of long-term global water vapor and liquid water data set using ultra-fast methods to assimilate multi-satellite and radiosonde observations

During the next decade, many programs and experiments under the Global Energy and Water Cycle Experiment (GEWEX) will utilize present day and future data sets to improve our understanding of the role of moisture in climate, and its interaction with other variables such as clouds and radiation. An important element of GEWEX will be the GEWEX Water Vapor Project (GVaP), which will eventually initiate a routine, real-time assimilation of the highest quality, global water vapor data sets including information gained from future data collection systems, both ground and space based. The comprehensive global water vapor data set being produced by METSAT Inc. uses a combination of ground-based radiosonde data, and infrared and microwave satellite retrievals. This data is needed to provide the desired foundation from which future GEWEX-related research, such as GVaP, can build. The first year of this project was designed to use a combination of the best available atmospheric moisture data including: radiosonde (balloon/acft/rocket), HIRS/MSU (TOVS) retrievals, and SSM/I retrievals, to produce a one-year, global, high resolution data set of integrated column water vapor (precipitable water) with a horizontal resolution of 1 degree, and a temporal resolution of one day. The time period of this pilot product was to be det3ermined by the availability of all the input data sets. January 1988 through December 1988 were selected. In addition, a sample of vertically integrated liquid water content (LWC) was to be produced with the same temporal and spatial parameters. This sample was to be produced over ocean areas only. Three main steps are followed to produce a merged water vapor and liquid water product. Input data from Radiosondes, TOVS, and SSMI/I is quality checked in steps one and two. Processing is done in step two to generate individual total column water vapor and liquid water data sets. The third step, and final processing task, involves merging the individual output products to produce the integrated water vapor product. A final quality control is applied to the merged data sets.

Vonderhaar, T. H.↗

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.↗

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote sources of water for precipitation, based on the implementation of passive constituent tracers of water vapor (termed water vapor tracers, WVT) in a general circulation model. In this case, the major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In this approach, each WVT is associated with an evaporative source region, and tracks the water until it precipitates from the atmosphere. By assuming that the regional water is well mixed with water from other sources, the physical processes that act on the WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be computed within the model simulation, and can be validated against the model's prognostic water vapor. Furthermore, estimates of precipitation recycling can be compared with bulk diagnostic approaches. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional tracers, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic 2 regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In general, most North American land regions showed a positive correlation between evaporation and recycling ratio (except the Southeast United States) and negative correlations of recycling ratio with precipitation and moisture transport (except the Southwestern United States). The Midwestern local source is positively correlated with local evaporation, but it is not correlated with water vapor transport. This is contrary to bulk diagnostic estimates of precipitation recycling. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.↗

Assimilation of Soil Moisture Observations Over Land Improves Analysis and Prediction of Tropical Cyclone Idai

Soil moisture conditions can impact the circulation and structure of a tropical cyclone (TC) when part or all of the circulation is over land. Dry land surface conditions may lead to faster dissipation of a TC over land, whereas very wet conditions may lead to a prolonged maintenance of its intensity. While this relationship is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tbs significantly improves modeled land surface states. Here we evaluate: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP Tb observations. We find that in the analysis with SMAP assimilation, the TC has a better-defined, more aligned vertical structure over land relative to the control run; moreover, the analyzed TC size, as measured by the wind speed radius, better matches the observed TC size. We further find significant reductions in the forecast intensity error and the forecast along-track error, measured against observations. The largest error reductions occur at lead times of 36 to 72 hours, suggesting that the land with its longer memory gains in importance as a source of predictability at this timescale. An investigation of the underlying mechanisms leading to the skill improvements from SMAP data assimilation revealed that the assimilation of SMAP leads to wetter soil moisture conditions and an increased latent heat flux in the SMAP analysis, which results in a TC with higher column-integrated total moisture content and total energy compared to the control analysis.

Jana Kolassa↗

Evaluating Nighttime CALIOP 0.532 micron Aerosol Optical Depth and Extinction Coefficient Retrievals

NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) Version 3.01 5-km nighttime 0.532 micron aerosol optical depth (AOD) datasets from 2007 are screened, averaged and evaluated at 1 deg X 1 deg resolution versus corresponding/co-incident 0.550 micron AOD derived using the US Navy Aerosol Analysis and Prediction System (NAAPS), featuring two-dimensional variational assimilation of quality-assured NASA Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) AOD. In the absence of sunlight, since passive radiometric AOD retrievals rely overwhelmingly on scattered radiances, the model represents one of the few practical global estimates available from which to attempt such a validation. Daytime comparisons, though, provide useful context. Regional-mean CALIOP vertical profiles of night/day 0.532 micron extinction coefficient are compared with 0.523/0.532 micron ground-based lidar measurements to investigate representativeness and diurnal variability. In this analysis, mean nighttime CALIOP AOD are mostly lower than daytime (0.121 vs. 0.126 for all aggregated data points, and 0.099 vs. 0.102 when averaged globally per normalised 1 deg. X 1 deg. bin), though the relationship is reversed over land and coastal regions when the data are averaged per normalised bin (0.134/0.108 vs. 0140/0.112, respectively). Offsets assessed within single bins alone approach +/- 20 %. CALIOP AOD, both day and night, are higher than NAAPS over land (0.137 vs. 0.124) and equal over water (0.082 vs. 0.083) when averaged globally per normalised bin. However, for all data points inclusive, NAAPS exceeds CALIOP over land, coast and ocean, both day and night. Again, differences assessed within single bins approach 50% in extreme cases. Correlation between CALIOP and NAAPS AOD is comparable during both day and night. Higher correlation is found nearest the equator, both as a function of sample size and relative signal magnitudes inherent at these latitudes. Root mean square deviation between CALIOP and NAAPS varies between 0.1 and 0.3 globally during both day/night. Averaging of CALIOP along-track AOD data points within a single NAAPS grid bin improves correlation and RMSD, though day/night and land/ocean biases persist and are believed systematic. Vertical profiles of extinction coefficient derived in the Caribbean compare well with ground-based lidar observations, though potentially anomalous selection of a priori lidar ratios for CALIOP retrievals is likely inducing some discrepancies. Mean effective aerosol layer top heights are stable between day and night, indicating consistent layer-identification diurnally, which is noteworthy considering the potential limiting effects of ambient solar noise during day.

aerosol optical depth↗

Use of three-cornered hat error estimates in MERRA-2 to guide an improved reanalysis-Part 1

The three-cornered hat (3CH) method estimates the uncertainties of three different co-located model or observational data sets (Anthes and Rieckh, 2018; Sjoberg et al., 2021). Rieckh et al. (2021) used the 3CH method to compare the random error statistics of different global forecast and reanalysis models, as well as radio occultation (RO) and radiosonde observations. That study showed that the MERRA-2 reanalysis, while having smaller errors in the stratosphere than its predecessor MERRA, had larger errors in the troposphere than many of the other data sets analyzed. The MERRA-2 errors were particularly large in the tropics. In a collaborative effort between UCAR’s COSMIC (Constellation Observing System for Meteorology, Ionosphere and Meteorology) program and NASA’s Global Modeling and Assimilation Office (GMAO), we carried out further 3CH error diagnostics to help isolate the causes of these larger errors and help guide the development of an improved reanalysis. This presentation summarizes random error statistics associated with MERRA-2, ECMWF’s ERA5 reanalysis, and COSMIC-2 (C2) RO observations. We compute 3CH error variance estimates of refractivity, as well as temperature and specific humidity using UCAR’s COSMIC Data Analysis and Archive Center (CDAAC) improved 1D-variational (1D-Var) retrieval (wetPf2) over 15 latitude bands from 45S to 45N. The 1D-Var retrievals of specific humidity and temperature for C2 use NCEP’s Global Forecast System (GFS) as the background. Anthes et al. (2021) showed that it gives accurate estimates of temperature and specific humidity in the tropics and subtropics, even in the challenging environment of intense Hurricane Dorian (2019). This presentation confirms the previous results that MERRA-2 has significantly larger errors in the tropics and subtropics than either C2 or ERA5. Its errors are larger between 30S and 30N compared to 30-45 N-S latitudes, and are also larger over land compared to oceans. Most of the MERRA-2 refractivity errors come from specific humidity, except over land below 3 km where temperature errors are large. These results suggest that moist convection and atmospheric boundary layer physics in MERRA-2 may be responsible for a significant part of the higher uncertainties. These results are being used to guide GMAO in developing an improved next-generation reanalysis, as shown in a companion presentation submitted to this conference (El Akkraoui et al., 2021), which extends this study and describes improvements to MERRA-2 leading to the next GMAO reanalysis.

Jeremiah Sjoberg↗

Online Simulations of Global Aerosol Distributions in the NASA GEOS-4 Model and Comparisons to Satellite and Ground-Based Aerosol Optical Depth

We have implemented a module for tropospheric aerosols (GO CART) online in the NASA Goddard Earth Observing System version 4 model and simulated global aerosol distributions for the period 2000-2006. The new online system offers several advantages over the previous offline version, providing a platform for aerosol data assimilation, aerosol-chemistry-climate interaction studies, and short-range chemical weather forecasting and climate prediction. We introduce as well a methodology for sampling model output consistently with satellite aerosol optical thickness (AOT) retrievals to facilitate model-satellite comparison. Our results are similar to the offline GOCART model and to the models participating in the AeroCom intercomparison. The simulated AOT has similar seasonal and regional variability and magnitude to Aerosol Robotic Network (AERONET), Moderate Resolution Imaging Spectroradiometer, and Multiangle Imaging Spectroradiometer observations. The model AOT and Angstrom parameter are consistently low relative to AERONET in biomass-burning-dominated regions, where emissions appear to be underestimated, consistent with the results of the offline GOCART model. In contrast, the model AOT is biased high in sulfate-dominated regions of North America and Europe. Our model-satellite comparison methodology shows that diurnal variability in aerosol loading is unimportant compared to sampling the model where the satellite has cloud-free observations, particularly in sulfate-dominated regions. Simulated sea salt burden and optical thickness are high by a factor of 2-3 relative to other models, and agreement between model and satellite over-ocean AOT is improved by reducing the model sea salt burden by a factor of 2. The best agreement in both AOT magnitude and variability occurs immediately downwind of the Saharan dust plume.

Colarco, Peter↗

Modern Era Retrospective Restrospective-Analysis for Research and Applications (MERRA) Data and Services at the GES DISC

The Modern Era Retrospective-analysis for Research and Applications (MERRA) dataset is a NASA satellite era, 30 year (1979 - present), reanalysis using the Goddard Earth Observing System Data Assimilation System, Version 5 (GEOS-5). The project, run out of NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center, provides the science and application communities with a state-of-the-art global analysis with emphasis on improved estimates of the hydrological cycle over a broad range of weather and climate time scales. MERRA products are generated as a long-term synthesis that places the NASA EOS suite of observations in a climate context. The MERRA analysis is performed at a horizontal resolution of 2/3 longitude x 1/2 latitude (540x361 global gridpoints) with observational analyses every 6 hours. The MERRA output data will include 3 dimensional state fields for every 6 hourly analysis cycle on 42 pressure levels (or 72 terrain following model coordinate levels) from the surface through the stratosphere. Several data products are specifically designed to support chemistry and stratosphere transport modeling. The 2 dimensional surface and atmospheric diagnostics (numbering 259) are being stored on the native grid at 1 hourly intervals. These include radiation and vertical integrals of the atmosphere for water and energy budget studies and also surface diagnostics where the diurnal cycle is important. The one hourly surface and near surface data product will also facilitate research on the integrated analysis of Earth system observations in the land, ocean and cryosphere. The MERRA products are archived and distributed by the Goddard Earth Sciences Data and Information Services Center (GES DISC) through its Modeling DISC Web (MDISC) portal. Multiple data access methods and services are available for MERRA data through MDISC: (1) Mirador offers a quick, comprehensive search of MERRA and all GES DISC archived data holdings, allowing searches on keywords, location names or latitude/longitude box, and date/time, with responses within a few seconds. (2) Giovanni is a GES DISC developed Web application that provides data visualization and analysis online. Giovanni features popular visualizations such as latitude-longitude maps, animations, cross sections, profiles, time series, etc. and some basic statistical analysis functions such as scatter plots and correlation coefficient maps. Users are able to download results in several different formats, including Google Earth. (3) On-the-fly parameter subsetting of data within a spatial/temporal window is provided through a simple select and click Web page. (4) MERRA data are also available via OPeNDAP, GrADS Data Server (GDS) and can be converted to netCDF on the fly.

Berrick, Stephen W.↗

Integrating Remote Sensing Observations with NASA's GEOS-5 Modeling Framework in Support of Retrospective Analyses and Seasonal Prediction of Biosphere-Atmosphere CO2 Flux

The absorption of nearly half of all anthropogenic carbon dioxide (CO2) emissions by terrestrial and marine ecosystems has played a critical role in mitigating climate change. However, a persistent lack of understanding of these uptake processes impedes attribution, and thereby introduces large uncertainties into projections of the trajectory of carbon and climate in coming decades. A better understanding and quantification of the seasonal and interannual variability of carbon flux between the atmosphere and ocean and terrestrial biospheres is critically needed to improve predictions from Earth system models. We propose to bring together a diverse suite of remote sensing observations, mechanistic models, and NASA's Goddard Earth Observing System, version 5 (GEOS-5) model to address two objective 1) quantify the processes controlling the temporal variability of biosphere-atmosphere CO2 flux on seasonal and interannual timescales from 1989 to 2018 and 2) evaluating the predictability of the components of carbon flux on seasonal timescales. Because GEOS-5 has been developed to incorporate a variety of Earth system observations and includes a complete and physically consistent depiction of the atmosphere, land, and ocean carbon cycles, it represents a unique and ideal synthesis framework for this effort. By expanding the realism of processes that could be included in GEOS-5 seasonal forecasts, this effort supports NASA's goals toward an integrated Earth system modeling and prediction framework.We propose to integrate a number of observations into the GEOS-5 modeling system to refine retrospective estimates of both terrestrial and marine carbon flux and to test hypotheses regarding the effects of climate drivers on terrestrial carbon flux. This will include assimilating several satellite ocean color products in the NASA Ocean Biogeochemistry Model (NOBM) to improve the representation of ocean biology. We will also utilize remote sensing datasets to increase the temporal resolution of global land-use change estimates and fire emissions to better constrain the temporal variability of the terrestrial carbon cycle. This information will be incorporated into two NASA terrestrial biosphere models Catchment-CN, a dynamic vegetation model that is integrated within GEOS-5 and the Ecosystem Demography (ED) model, a core modeling resource supporting the Global Ecosystem Dynamics Investigation (GEDI) mission to map biomass. Fluxes will be evaluated through 1) comparison to top-down flux estimates inferred from surface CO2 observations and 2) comparison against available aircraft and satellite observations.By combining predictions of net ecosystem exchange and air-sea CO2 flux from land and ocean models with statistical predictions of land use change, fire, and fossil fuel emissions, we propose to create the first dynamical CO2 forecast on seasonal timescales. Because of the experimental nature of this endeavor, an important project objective will be to evaluate the limits to the predictability of different carbon cycle processes over lead times of 1-12 months using different satellite and surface data streams.

Ott, Lesley↗

Volatiles in High-K Lunar Basalts

Chlorine is an unusual isotopic system, being essentially unfractionated ((delta)Cl-37 approximately 0 per mille ) between bulk terrestrial samples and chondritic meteorites and yet showing large variations in lunar (approximately -4 to +81 per mille), martian, and vestan (HED) samples. Among lunar samples, the volatile-bearing mineral apatite (Ca5(PO4)3[F,Cl,OH]) has been studied for volatiles in K-, REE-, and P (KREEP), very high potassium (VHK), low-Ti and high-Ti basalts, as well as samples from the lunar highlands. These studies revealed a positive correlation between in-situ (delta)Cl-37 measurements and bulk incompatible trace elements (ITEs) and ratios. Such trends were interpreted to originate from Cl isotopic fractionation during the degassing of metal chlorides during or shortly after the differentiation of the Moon via a magma ocean. In this study, we investigate the volatile inventories of a group of samples for which new-era volatile data have yet to be reported - the high-K (greater than 2000 ppm bulk K2O), high-Ti, trace element-rich mare basalts. We used isotope imaging on the Cameca NanoSIMS 50L at JSC to obtain the Cl isotopic composition [((Cl-37/(35)Clsample/C-37l/(35)Clstandard)-1)×1000, to get a value in per thousand (per mille)] which ranges from approximately -2.7 +/- 2 per mille to +16.1 +/- 2 per mille (2sigma), as well as volatile abundances (F & Cl) of apatite in samples 10017, 10024 & 10049. Simply following prior models, as lunar rocks with high bulk-rock abundances of ITEs we might expect the high-K, high-Ti basalts to contain apatite characterized by heavily fractionated (delta)Cl-37 values, i.e., Cl obtained from mixing between unfractionated mantle Cl (approximately 0 per mille) and the urKREEP reservoir (possibly fractionated to greater than +25 per mille.). However, the data obtained for the studied samples do not conform to either the early degassing or mixing models. Existing petrogentic models for the origin of the high-K, high-Ti basalts do not include urKREEP assimilation into their LMO cumulate sources. Therefore, Cl in these basalts either originated from source region heterogeneity or through assimilation or metasomatism by volatile and incompatible trace element rich materials. The new data presented here could provide evidence for the existence of region(s) in the lunar interior that are ITE-enriched and contain Cl that does not share isotopic affinities with lunar urKREEP, possibly representing the composition of the purported 'neuKREEP'.

Barnes, Jessica J.↗

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗

A Reference Ocean Surface Emission and Backscatter Model from Microwaves to Infrared

Satellite observations are vital for the initialization of Numerical Weather Prediction models, and very important for climate monitoring and prediction, as well as other applications such as hydrology and flood awareness prediction. Knowledge of radiative contributions from the Earth's surface is needed to sound the lower troposphere from space. The lack of a reference quality ocean emission and backscatter model is a major gap in our ability to provide absolute calibration of the satellite based observing system. Uncertainty in emissivity models is not well characterized and different models are used for different spectral bands, for active and passive instruments. An International Space Science Institute (ISSI) team was put together [4] to address these issues. The objectives of the team are to provide a reference model as a community software (i.e., documented and freely available code), that is maintained and supported, has traceable uncertainty estimations, and that enables new science from microwaves to infrared with bidirectional reflectance distribution function (BRDF) capability. We will present the model and its various components, discussing the choices between various parameterizations, building on the LOCEAN model of [2]. The model predictions will be evaluated at various frequencies, including comparisons to radiometric observations by SMAP, AMSR2 and GMI (e.g., [5]). We will discuss early model evaluation in the infrared and for active microwave sensors. Areas of ongoing research include improving the foam parametrization (coverage and emissivity) to provide consistent performances across frequencies, building on [1], and the azimuthal dependence of the active and passive signals. The model will be used to generate training data for fast models e.g., Fastem, [3], that are used in operational data assimilation and climate re-analysis.

Emmanuel Dinnat↗

Evolution of the 1997-98 El Nino and La Nina: the Role of Altimetry in ENSO Studies and Prediction

El Nino and the Southern Oscillation (ENSO) are known to be a quasi-periodic coupled ocean-atmosphere phenomenon occurring roughly every 3-7 years. Its largest oceanic manifestations are found in the equatorial Pacific, but it is of global climatic significance. During the flight of the TOPEX]Poseidon altimeter, the equatorial Pacific was in the midst of an unusual period of prolonged ENSO warming in the mid-1990's. In 1997, the tropical Pacific witnessed the development of a major El Nino event, rivaling in strength the 1982-83 El Nino. The intensity of this El Nino, and its rapid intensification early in the calendar year, caught the scientific community by surprise. None of the forecast models, statistical or dynamical, captured these features of the 1997 El Nino, although those models that assimilated in situ ocean observations did better than most in predicting that the event would occur in 1997. The onset of the El Nino coincided with the occurrence of several westerly wind events in the western Pacific rather than to delayed oscillator physics. However, western Pacific sea levels were anomalously elevated a year and a half prior to the onset, which may have helped precondition the system to a particularly strong ENSO episode. The coming of La Nina was foreshadowed by large-scale drops in sea level in the western tropical Pacific and cyclonic wind stress curl anomalies far beyond the +\- 8 degrees latitude bounds of the TAO mooring array. In past El Nino events, the evolution of the coupled system could only be studied from the perspective of in situ observations. The 1997-1998 El Nino event will be the first time a major El Nino event will have been observed from start to finish via remotely-sensed measurements of sea surface topography, sea surface temperature, sea surface winds, ocean color, and precipitation. These observations have considerable potential for advancing coupled climate forecasts on time scales from seasons to decades. This presentation will address questions along the lines of: (1) What are the scientific problems in ENSO that can benefit most from altimetry data? (2) Are there any precursors in the global ocean-atmosphere system for the 1997-98 extreme event? (3) Is there a unique role of altimetry in detecting potential precursors?(4) Has Topex/Poseidon data been adequately used to improve ENSO forecasts? Is there any firm evidence? (6) What is the status of the transition of El Nino into a La Nina?

Busalacchi, Antonio J.↗