Results for “Reanalysis data” — Page 32

NASA NTRS2018

GEOS S2S-2_1 File Specification: GMAO Seasonal and Sub-Seasonal Forecast Output

The NASA GMAO seasonal (9 months) and subseasonal (45 days) forecasts are produced with the Goddard Earth Observing System (GEOS) Atmosphere-Ocean General Circulation Model and Data Assimilation System Version S2S-2_1. The new system replaces version S2S-1.0 described in Borovikov et al (2017), and includes upgrades to many components of the system. The atmospheric model includes an upgrade from a pre-MERRA-2 version running on a latitude-longitude grid at approx. 1 degree resolution to a current version running on a cubed sphere grid at approximately 1/2 degree resolution. The important developments are related to the dynamical core (Putman et al., 2011), the moist physics (''two-moment microphysics'' of Barahona et al., 2014) and the cryosphere (Cullather et al., 2014). As in the previous GMAO S2S system, the land model is that of Koster et al (2000). GMAO S2S-2_1 now includes the Goddard Chemistry Aerosol Radiation and Transport (GOCART, Colarco et al., 2010) single moment interactive aerosol model that includes predictive aerosols including dust, sea salt and several species of carbon and sulfate. The previous version of GMAO S2S specified aerosol amounts from climatology, which were used to inform the atmospheric radiation only. The ocean model includes an upgrade from MOM4 to MOM5 (Griffies 2012), and continues to be run on the tripolar grid at approximately 1/2 degree resolution in the tropics with 40 vertical levels. As in S2S-1.0, the sea ice model is from the Los Alamos Sea Ice model (CICE4, Hunke and Lipscomb 2010). The Ocean Data Assimilation System (ODAS) has been upgraded from the one described in Borovikov et al., 2017 to one that uses a modified version of the Penny, 2014 Local Ensemble Transform Kalman Filter (LETKF), and now assimilates along-track altimetry. The ODAS also does a nudging to MERRA-2 SST and sea ice boundary conditions. The atmospheric data assimilation fields used to constrain the atmosphere in the ODAS have been upgraded from MERRA to a MERRA-2 like system. The system is initialized using a MERRA-2-like atmospheric reanalysis (Gelaro et al. 2017) and the GMAO S2S-2_1 ocean analysis. Additional ensemble members for forecasts are produced with initial states at 5-day intervals, with additional members based on perturbations of the atmospheric and ocean states. Both subseasonal and seasonal forecasts are submitted to the National MultiModel Ensemble (NMME) project, and are part of the US/Canada multimodel seasonal forecasts (http://www.cpc.ncep.noaa.gov/products/NMME/). A large suite of retrospective forecasts (''hindcasts'') have been completed, and contribute to the calculation of the model's baseline climatology and drift, anomalies from which are the basis of the seasonal forecasts.

Seasonal Forecast↗

NASA NTRS2022

Automatic variable selection in ecological niche modeling: A case study using Cassin’s Sparrow (Peucaea cassinii)

MERRA/Max provides a feature selection approach to dimensionality reduction that enables direct use of global climate model outputs in ecological niche modeling. The system accomplishes this reduction through a Monte Carlo optimization in which many independent MaxEnt runs, operating on a species occurrence file and a small set of randomly selected variables in a large collection of variables, converge on an estimate of the top contributing predictors in the larger collection. These top predictors can be viewed as potential candidates in the variable selection step of the ecological niche modeling process. MERRA/Max’s Monte Carlo algorithm operates on files stored in the underlying filesystem, making it scalable to large data sets. Its software components can run as parallel processes in a high-performance cloud computing environment to yield near real-time performance. In tests using Cassin’s Sparrow (Peucaea cassinii) as the target species, MERRA/Max selected a set of predictors from Worldclim’s Bioclim collection of 19 environmental variables that have been shown to be important determinants of the species’ bioclimatic niche. It also selected biologically and ecologically plausible predictors from a more diverse set of 86 environmental variables derived from NASA’s Modern-Era Retrospective Analysis for Research and Applications Version 2 (MERRA-2) reanalysis, an output product of the Goddard Earth Observing System Version 5 (GEOS-5) modeling system. We believe these results point to a technological approach that could expand the use global climate model outputs in ecological niche modeling, foster exploratory experimentation with otherwise difficult-to-use climate data sets, streamline the modeling process, and, eventually, enable automated bioclimatic modeling as a practical, readily accessible, low-cost, commercial cloud service.

John L. Schnase↗

NASA NTRSDate not supplied

TPSAS-NF1676L-33626-DND

NASA is conducting a series of computational experiments to quantify atmospheric effects on low noise sonic booms. In the current study, simulated cruise nearfield pressure data from NASA?s X-59 Quiet Supersonic Technology aircraft was propagated from the aircraft to the ground at four cardinal headings through five years of realistic atmospheric profiles at 30 locations across the eastern USA. Statistical design of experiments was used to select the locations where primary sonic thump carpet widths (CW) and metric levels at the ground were computed. Atmospheric profiles were taken from the Climate Forecast System Reanalysis database, which contains reanalyzed atmospheric profiles four times daily. Decision tree analyses were performed to determine relative importance of predictors for CW and metric levels. Predictors included were latitude, longitude, date, time of day, season, climate, aircraft heading, and ground elevation. Results of this study indicate the propagation resolution needed to adequately characterize the distributions of CW and ground metric data, i.e., necessary separation distance between propagation locations and total number of atmospheres through which to propagate. This resolution will be used for a follow-on study of simulated X-59 carpets across the entire US mainland.

Will Doebler↗

NASA NTRS2004

A Subsonic Aircraft Design Optimization With Neural Network and Regression Approximators

The Flight-Optimization-System (FLOPS) code encountered difficulty in analyzing a subsonic aircraft. The limitation made the design optimization problematic. The deficiencies have been alleviated through use of neural network and regression approximations. The insight gained from using the approximators is discussed in this paper. The FLOPS code is reviewed. Analysis models are developed and validated for each approximator. The regression method appears to hug the data points, while the neural network approximation follows a mean path. For an analysis cycle, the approximate model required milliseconds of central processing unit (CPU) time versus seconds by the FLOPS code. Performance of the approximators was satisfactory for aircraft analysis. A design optimization capability has been created by coupling the derived analyzers to the optimization test bed CometBoards. The approximators were efficient reanalysis tools in the aircraft design optimization. Instability encountered in the FLOPS analyzer was eliminated. The convergence characteristics were improved for the design optimization. The CPU time required to calculate the optimum solution, measured in hours with the FLOPS code was reduced to minutes with the neural network approximation and to seconds with the regression method. Generation of the approximators required the manipulation of a very large quantity of data. Design sensitivity with respect to the bounds of aircraft constraints is easily generated.

Patnaik, Surya N.↗

NASA NTRS2019

A Comparison of the GPROF Algorithm Results in the 2A-CLIM and 2A Data Products

The GPROF algorithm is used to retrieve precipitation values from all the GPM mission partner radiometers (i.e. GMI, MHS, ATMS, SSMIS, AMSR2). GPM produces two precipitation retrieval products (processing at level 2) using the GPROF algorithm. One has the type designation 2A and the other the designation 2A-CLIM. As a result understanding the difference between the two products is important for using and understanding them. Both of these GPROF products use exactly the same algorithm and a-priori databases. The only difference between the two types is the ancillary products used in the initial processing step. In the case of the 2A-CLIM products, ECMWF ERA-I is used as ancillary data. In the case of the 2A products, JMA's GANAL product is used as ancillary data. This presentation will provide a comparison between the two products using zonal mean comparisons. It will provide comparisons separately for land and ocean retrievals. It will compare the GMI differences over the entire current GPM mission life. As the only difference between the two products is the ancillary information used to create them, the comparison provides an elementary sensitivity study of the variables used from the reanalysis products and their impact on the GPROF retrieval. The purpose of the presentation is to provide users with an insight into the two GPROF products and provide them key information useful when using either of the products in scientific studies. Lastly, the presentation also discusses the reasons why only the 2A-CLIM product can be produced for the 17+ years of TRMM era data products.

GPM↗

NASA NTRS2015

Estimating Planetary Boundary Layer Heights from NOAA Profiler Network Wind Profiler Data

An algorithm was developed to estimate planetary boundary layer (PBL) heights from hourly archived wind profiler data from the NOAA Profiler Network (NPN) sites located throughout the central United States. Unlike previous studies, the present algorithm has been applied to a long record of publicly available wind profiler signal backscatter data. Under clear conditions, summertime averaged hourly time series of PBL heights compare well with Richardson-number based estimates at the few NPN stations with hourly temperature measurements. Comparisons with clear sky reanalysis based estimates show that the wind profiler PBL heights are lower by approximately 250-500 m. The geographical distribution of daily maximum PBL heights corresponds well with the expected distribution based on patterns of surface temperature and soil moisture. Wind profiler PBL heights were also estimated under mostly cloudy conditions, and are generally higher than both the Richardson number based and reanalysis PBL heights, resulting in a smaller clear-cloudy condition difference. The algorithm presented here was shown to provide a reliable summertime climatology of daytime hourly PBL heights throughout the central United States.

Wind Profiler↗

NASA NTRS1992

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

NASA NTRS2012

Electrically-Active Convection in Tropical Easterly Waves and Implications for Tropical Cyclogenesis in the Atlantic and East Pacific

In this study, we investigate the characteristics of tropical easterly wave convection and the possible implications of convective structure on tropical cyclogenesis and intensification over the Atlantic Ocean and East Pacific using data from the Tropical Rainfall Measurement Mission Microwave Imager, Precipitation Radar (PR), and Lightning Imaging Sensor as well as infrared (IR) brightness temperature data from the NASA global-merged IR brightness temperature dataset. Easterly waves were partitioned into northerly, southerly, trough, and ridge phases based on the 700-hPa meridional wind from the NCEP-NCAR reanalysis dataset. Waves were subsequently divided according to whether they did or did not develop tropical cyclones (i.e., developing and nondeveloping, respectively), and developing waves were further subdivided according to development location. Finally, composites as a function of wave phase and category were created using the various datasets. Results suggest that the convective characteristics that best distinguish developing from nondeveloping waves vary according to where developing waves spawn tropical cyclones. For waves that developed a cyclone in the Atlantic basin, coverage by IR brightness temperatures .240 K and .210 K provide the best distinction between developing and nondeveloping waves. In contrast, several variables provide a significant distinction between nondeveloping waves and waves that develop cyclones over the East Pacific as these waves near their genesis location including IR threshold coverage, lightning flash rates, and low-level (<4.5 km) PR reflectivity. Results of this study may be used to help develop thresholds to better distinguish developing from nondeveloping waves and serve as another aid for tropical cyclogenesis forecasting.

Leppert, Kenneth D., II↗

NASA NTRSDate not supplied

Simulations of X-59 Sonic Thumps and Traditional Sonic Booms Propagated Around the World for Three Atmospheric Models

Propagation simulations of sonic booms from supersonic aircraft through atmospheric data over time at fixed locations provides the opportunity to assess noise exposure statistics for different climate regions. Knowledge of climate-based differences in sonic boom noise exposure statistics is important to ensure that future civil supersonic aircraft noise certification standards are globally applicable and effective. In this presentation, simulated sonic booms from the NASA X-59 Quesst quiet supersonic aircraft and conventional supersonic aircraft were propagated through atmospheric data at 100 locations across the world using PCBoom. Noise exposure statistics are compared for propagation results from three different atmospheric databases (NOAA Global Forecast System, NOAA Climate Forecast System Version 2, and the ECMWF Reanalysis Version 5). These atmospheric models were chosen due to their global coverage, popularity, and database availability. Preliminary statistical models are fit to assess the impact of several factors including flight direction, season, ground elevation, and climate on noise exposure size and loudness. Areas with prevalence of higher noise due to their climate are identified, which could help inform future supersonic aircraft noise standards.

X-59↗

NASA NTRS2005

Estimating moisture transport over oceans using space-based observations

The moisture transport integrated over the depth of the atmosphere (0) is estimated over oceans using satellite data. The transport is the product of the precipitable water and an equivalent velocity (ue), which, by definition, is the depth-averaged wind velocity weighted by humidity. An artificial neural network is employed to construct a relation between the surface wind velocity measured by the spaceborne scatterometer and coincident ue derived using humidity and wind profiles measured by rawinsondes and produced by reanalysis of operational numerical weather prediction (NWP). On the basis of this relation, 0 fields are produced over global tropical and subtropical oceans (40_N- 40_S) at 0.25_ latitude-longitude and twice daily resolutions from August 1999 to December 2003 using surface wind vector from QuikSCAT and precipitable water from the Tropical Rain Measuring Mission. The derived ue were found to capture the major temporal variability when compared with radiosonde measurements. The average error over global oceans, when compared with NWP data, was comparable with the instrument accuracy specification of space-based scatterometers. The global distribution exhibits the known characteristics of, and reveals more detailed variability than in, previous data.

Artificial Neural Network↗

NASA NTRS2013

The Use of Red Green Blue Air Mass Imagery to Investigate the Role of Stratospheric Air in a Non-convective Wind Event

Non-convective wind events commonly occur with passing extratropical cyclones and have significant societal and economic impacts. Since non-convective winds often occur in the absence of specific phenomena such as a thunderstorm, tornado, or hurricane, the public are less likely to heed high wind warnings and continue daily activities. Thus non-convective wind events result in as many fatalities as straight line thunderstorm winds. One physical explanation for non-convective winds includes tropopause folds. Improved model representation of stratospheric air and associated non-convective wind events could improve non-convective wind forecasts and associated warnings. In recent years, satellite data assimilation has improved skill in forecasting extratropical cyclones; however errors still remain in forecasting the position and strength of extratropical cyclones as well as the tropopause folding process. The goal of this study is to determine the impact of assimilating satellite temperature and moisture retrieved profiles from hyperspectral infrared (IR) sounders (i.e. Atmospheric Infrared Sounder (AIRS), Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI)) on the model representation of the tropopause fold and an associated high wind event that impacted the Northeast United States on 09 February 2013. Model simulations using the Advanced Research Weather Research and Forecasting Model (ARW) were conducted on a 12-km grid with cycled data assimilation mimicking the operational North American Model (NAM). The results from the satellite assimilation run are compared to a control experiment (without hyperspectral IR retrievals), Modern Era-Retrospective Analysis for Research and Applications (MERRA) reanalysis, and Rapid Refresh analyses.

Berndt, E. B.↗

NASA NTRS2013

The Impact of the Assimilation of Hyperspectral Infrared Retrieved Profiles on Advanced Weather and Research Model Simulations of a Non-Convective Wind Event

Non-convective wind events commonly occur with passing extratropical cyclones and have significant societal and economic impacts. Since non-convective winds often occur in the absence of specific phenomena such as a thunderstorm, tornado, or hurricane, the public are less likely to heed high wind warnings and continue daily activities. Thus non-convective wind events result in as many fatalities as straight line thunderstorm winds. One physical explanation for non-convective winds includes tropopause folds. Improved model representation of stratospheric air and associated non-convective wind events could improve non-convective wind forecasts and associated warnings. In recent years, satellite data assimilation has improved skill in forecasting extratropical cyclones; however errors still remain in forecasting the position and strength of extratropical cyclones as well as the tropopause folding process. The goal of this study is to determine the impact of assimilating satellite temperature and moisture retrieved profiles from hyperspectral infrared (IR) sounders (i.e. Atmospheric Infrared Sounder (AIRS), Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI)) on the model representation of the tropopause fold and an associated high wind event that impacted the Northeast United States on 09 February 2013. Model simulations using the Advanced Research Weather Research and Forecasting Model (ARW) were conducted on a 12-km grid with cycled data assimilation mimicking the operational North American Model (NAM). The results from the satellite assimilation run are compared to a control experiment (without hyperspectral IR retrievals), Modern Era-Retrospective Analysis for Research and Applications (MERRA) reanalysis, and Rapid Refresh analyses.

Berndt, Emily↗

NASA NTRS2014

Impact of the Assimilation of Hyperspectral Infrared Retrieved Profiles on Advanced Weather and Research Model Simulations of a Non-Convective Wind Event

Non-convective wind events commonly occur with passing extratropical cyclones and have significant societal and economic impacts. Since non-convective winds often occur in the absence of specific phenomena such as a thunderstorm, tornado, or hurricane, the public are less likely to heed high wind warnings and continue daily activities. Thus non-convective wind events result in as many fatalities as straight line thunderstorm winds. One physical explanation for non-convective winds includes tropopause folds. Improved model representation of stratospheric air and associated non-convective wind events could improve non-convective wind forecasts and associated warnings. In recent years, satellite data assimilation has improved skill in forecasting extratropical cyclones; however errors still remain in forecasting the position and strength of extratropical cyclones as well as the tropopause folding process. The goal of this study is to determine the impact of assimilating satellite temperature and moisture retrieved profiles from hyperspectral infrared (IR) sounders (i.e. Atmospheric Infrared Sounder (AIRS), Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI)) on the model representation of the tropopause fold and an associated high wind event that impacted the Northeast United States on 09 February 2013. Model simulations using the Advanced Research Weather Research and Forecasting Model (ARW) were conducted on a 12-km grid with cycled data assimilation mimicking the operational North American Model (NAM). The results from the satellite assimilation run are compared to a control experiment (without hyperspectral IR retrievals), Modern Era-Retrospective Analysis for Research and Applications (MERRA) reanalysis, and Rapid Refresh analyses.

Berndt, E. B.↗

NASA NTRS2017

Aerosol EnKF at GMAO

In the GEOS near real-time system, as well as in MERRA-2 which is the latest reanalysis produced at NASAs Global Modeling and Assimilation Office(GMAO), the assimilation of aerosol observations is performed by means of a so-called analysis splitting method. In line with the transition of the GEOS meteorological data assimilation system to a hybrid Ensemble-Variational formulation, we are updating the aerosol component of our assimilation system to an ensemble square root filter(EnSRF; Whitaker and Hamill (2002)) type of scheme.We present a summary of our preliminary results of the assimilation of column integrated aerosol observations (Aerosol Optical Depth; AOD) using an Ensemble Square Root Filters (EnSRF) scheme and the ensemble members produced routinely by the meteorological assimilation.

MERRA-2 EnSRF↗

NASA NTRSDate not supplied

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s Global Earth Observing System

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy, Jr.↗

NASA NTRSDate not supplied

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗

NASA NTRSDate not supplied

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and increases in number from approximately 6 thousand RO profiles per day during 2019 to approximately 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2022. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗

NASA NTRS2015

Global Ocean Evaporation: How Well Can We Estimate Interannual to Decadal Variability?

Evaporation from the world's oceans constitutes the largest component of the global water balance. It is important not only as the ultimate source of moisture that is tied to the radiative processes determining Earth's energy balance but also to freshwater availability over land, governing habitability of the planet. Here we focus on variability of ocean evaporation on scales from interannual to decadal by appealing to three sources of data: the new MERRA-2 (Modern-Era Retrospective analysis for Research and Applications -2); climate models run with historical sea-surface temperatures, ice and atmospheric constituents (so-called AMIP experiments); and state-of-the-art satellite retrievals from the Seaflux and HOAPS (Hamburg Ocean-Atmosphere Parameters and Fluxes from Satellite) projects. Each of these sources has distinct advantages as well as drawbacks. MERRA-2, like other reanalyses, synthesizes evaporation estimates consistent with observationally constrained physical and dynamical models-but data stream discontinuities are a major problem for interpreting multi-decadal records. The climate models used in data assimilation can also be run with lesser constraints such as with SSTs and sea-ice (i.e. AMIPs) or with additional, minimal observations of surface pressure and marine observations that have longer and less fragmentary observational records. We use the new ERA-20C reanalysis produced by ECMWF embodying the latter methodology. Still, the model physics biases in climate models and the lack of a predicted surface energy balance are of concern. Satellite retrievals and comparisons to ship-based measurements offer the most observationally-based estimates, but sensor inter-calibration, algorithm retrieval assumptions, and short records are dominant issues. Our strategy depends on maximizing the advantages of these combined records. The primary diagnostic tool used here is an analysis of bulk aerodynamic computations produced by these sources and uses a first-order Taylor series analysis of wind speed, SST, near-surface stability and relative humidity variations around climatology to gauge the importance of these components. We find that the MERRA-2 evaporation record is strongly influenced by the availability of wind speed and humidity from passive microwave imagers beginning in the late 1980s as well as by the SST record. The trend over the period 1980 to present is nearly 10%. AMIP or the ERA-20C trends are much smaller. We find that ENSO-related signals involving both wind speed and thermodynamic variability remain the primary signal in the latter and are confirmed by satellite retrievals. We present uncertainty estimates based on the various data sources and discuss the implications for GEWEX water and energy budget science challenges.

Robertson, Franklin R.↗