Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reanalysis”

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 397 records · Page 22

Application of Aura OMI L2G Products Compared with NASA MERRA-2 Assimilation

The Ozone Monitoring Instrument (OMI) is one of the instruments aboard NASA's Aura satellite. It measures ozone total column and vertical profile, aerosols, clouds, and trace gases including NO2, SO2, HCHO, BrO, and OClO using absorption in the ultraviolet electromagnetic spectrum (280 - 400 nm). OMI Level-2G (L2G) products are based on the pixel-level OMI granule satellite measurements stored within global 0.25 deg. X 0.25 deg. grids, therefore they conserve all the Level 2 (L2) spatial and temporal details for 24 hours of scientific data in one file. The second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) is NASA's atmospheric reanalysis, using an upgraded version of Goddard Earth Observing System Model, version 5 (GEOS-5) data assimilation system. MERRA-2 includes aerosol data reanalysis and improved representations of stratospheric ozone, compared with its predecessor MERRA, in both instantaneous and time-averaged collections. It is found that simply comparing satellite Level-3 products might cause biases, due to lack of detailed temporal and original retrieval information. It is therefore preferable to inter-compare or implement satellite derived physical quantities directly with/to model assimilation with as high temporal and spatial resolutions as possible. This study will demonstrate utilization of OMI L2G daily aerosol and ozone products by comparing them with MERRA-2 hourly aerosol/ozone simulations, matched in both space and time aspects. Both OMI and MERRA-2 products are accessible online through NASA Goddard Earth Sciences Data Information Services Center (GES DISC, https://disc.gsfc.nasa.gov/).

OMI L2G products↗

ENSO Related Inter-Annual Lightning Variability from the Full TRMM LIS Lightning Climatology

The El Nino/Southern Oscillation (ENSO) contributes to inter-annual variability of lightning production more than any other atmospheric oscillation. This study further investigated how ENSO phase affects lightning production in the tropics and subtropics using the Tropical Rainfall Measuring Mission (TRMM) Lightning Imaging Sensor (LIS). Lightning data were averaged into mean annual warm, cold, and neutral 'years' for analysis of the different phases and compared to model reanalysis data. An examination of the regional sensitivities and preliminary analysis of three locations was conducted using model reanalysis data to determine the leading convective mechanisms in these areas and how they might respond to the ENSO phases

enso↗

File Specification for M2AMIP Products: Version Number - 1.0

The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) is an atmospheric reanalysis computed with the Goddard Earth Observing (EOS) System, Version 5.12.4 (GEOS) data assimilation system (Gelaro et al., 2017). To supplement the reanalysis, the GEOS General Circulation Model (GCM) used in MERRA-2 has been used to generate a 10-member ensemble of simulations, configured following the convention of the Atmospheric Model Intercomparison Project (AMIP; Gates et al., 1992). Each ensemble member was initialized using meteorological fields from a different date in November 1979. The AMIP simulations used the sea-surface temperature (SST) and sea-ice boundary conditions that were used in MERRA-2 (Bosilovich et al., 2016). This 10-member ensemble of AMIP simulations, denoted M2AMIP, is available for download in a group of self-describing files, which are documented in this office note. All data collections are provided on the same horizontal grid as MERRA-2. This grid has 576 points in the longitudinal direction and 361 points in the latitudinal direction, corresponding to a resolution of 0.625 degrees by 0.5 degrees. Although data collections are available at this grid, all fields are computed on a cubed-sphere grid with an approximate resolution of 50 km by 50 km and are then spatially interpolated to the latitude-longitude grid. There are no changes in the vertical grids used: variables are provided on either the native vertical grid of 72 model layers, or interpolated to 42 standard pressure levels. Unlike MERRA, no data collections are available at the vertical layer edges. More details on the grid are provided in Section 4. MERRA-2 introduced observation-based precipitation forcing for the land surface parameterization and the corresponding variable PRECTOTCORR in the MERRA-2 FLX (surface turbulent fluxes and related quantities) and LFO (land-surface forcing) collections (see Section 6; Reichle et al., 2017). While this variable is still available for M2AMIP, there was no observation-based forcing, making the value identical to the model derived precipitation, PRECTOT. Similarly, without data assimilation, the values for the analysis increments, D*DTANA, in the tendency and vertically integrated file collections are zero. The M2AMIP data are available for download online through the NASA Center for Climate Simulation (NCCS) DataPortal (https://portal.nccs.nasa.gov/datashare/gmao_m2amip/). Data are arranged in subdirectories based on ensemble member, followed by year and month. Control files that are compatible with the Grid Analysis and Display System (GrADS) are available in the ctl_daily and ctl_monthly directories for the hourly, three hourly, and monthly mean data. Control files for the monthly mean diurnal cycle can be found in the ctl_diurnal subdirectory within the directory for each individual ensemble member.

M2AMIP↗

Long-Term Ozone Variability and Trends from Reanalyses

Stratospheric ozone has a profound impact on radiation and chemistry over various spatial and temporal scales. The evolution of stratospheric ozone over the 21st century, however, is not well understood, especially in the lower stratosphere. Highly vertically resolved ozone data from satellite-borne limb sounders have proved to be invaluable for studying ozone in the middle and upper stratosphere but it was not until recently that these measurements were successfully incorporated in atmospheric reanalyses. Validation and comparison studies have demonstrated that the addition of observations from the Microwave Limb Sounder (MLS) on EOS (Earth Observing System) Aura greatly improved the quality of ozone fields in MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, Version 2) making these assimilated data sets more useful for scientific research. In this presentation we demonstrate that multidecadal lower-stratospheric ozone variability and trends can be derived from NASA's MERRA-2 reanalysis ozone. In particular, the reanalysis ozone bias-corrected using a chemistry model simulation as a transfer function agrees very well with recently reprocessed long ozonesonde records. Ozone trends in the lower stratosphere will be discussed in the context of recent findings (Ball et al., 2018) and interpreted in connection with long-term circulation changes in the lower stratosphere. Next, we show that the use of ozone data retrieved from the next generation OMPS (Ozone Mapping Profiler Suite) instruments, including the OMPS Limb Profiler, can successfully extend the reanalyses into the future allowing comprehensive monitoring of global ozone and interpretation of its evolution during the critical period of expected ozone recovery and climate change from increasing concentration of greenhouse gases.

ozone↗

Development of the GEOS-MIT gcm Atmosphere-Ocean Model for Coupled Data Assimilation

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MIT-GCM (General Circulation Model)) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (approximately 1/8 degrees) configuration of the GEOS atmosphere to a latitude-longitude-cap-1080 (approximately 1/12 degrees) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

The Development of the New GEOS-MITgcm Atmosphere-Ocean Model for Coupled Data Assimilation System

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8 deg) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12 deg) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

A New Atmosphere-Ocean Model for Studying Air-Sea Interactions and Coupled Data Assimilation

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8°) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12°) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

The Development of the New GEOS-MITgcm Atmosphere-Ocean Model for Coupled Data Assimilation System

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8°) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12°) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

The Role of Atmospheric Teleconnections and Local Forcings in Predicting Greenland Ice Sheet Surface Mass Loss

In recent decades, the Arctic climate has experienced substantial climactic change, including significant decreases in both sea ice extent and Greenland Ice Sheet (GrIS) surface mass balance. These trends are overlain by substantial interannual variability in atmospheric circulation driven by large-scale atmospheric teleconnection patterns. In addition, there is evidence to suggest that the removal of Arctic sea ice can alter local atmospheric circulation through increased air temperature, clouds, and water vapor, which may contribute to increased surface melting on the GrIS. Here, we seek to characterize how these processes are linked to Greenland Ice Sheet surface mass loss and constrain how the representation of these forcings can impact the prediction of meltwater runoff within the NASA Goddard Earth Observing System Model (GEOS) seasonal-to-subseasonal forecasting system (S2S v2.1). To do this, we use a combination of the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) reanalysis product, retrospective seasonal forecasts from the GEOS S2S v2.1, and independent GEOS simulations. Results from MERRA-2 reanalysis indicate that the negative phase of the North Atlantic Oscillation (NAO) results in warm surface air temperatures and reduced precipitation across Greenland, both of which act to enhance summer ice surface mass losses. When compared with MERRA-2, retrospective forecasts from the GEOS S2S v2.1 system effectively reproduce the pattern of summer GrIS surface mass loss and demonstrate reasonable skill in predicting the magnitude of meltwater runoff at leads of 1 to 3 months. However, during periods with a strong negative NAO, ice sheet surface mass balance is substantially underestimated. This pattern is also associated with an underprediction of the Greenland Blocking Index height and over prediction of sea ice extent, suggesting that both local and non-local forcings may play a role in the reduced prediction skill during these periods. Using both retrospective forecasts and independent simulations, we characterize the relative importance of local and non-local mechanisms in driving summer GrIS

Andrews, Lauren C.↗

Comparing Flood-Producing Catastrophic Hurricanes in GEOS: A Study of Juan 1985 and Harvey 2017

Hurricanes Juan (1985) and Harvey (2017) bear remarkable similarities because of the exceptional amount of accumulated precipitation over coastal regions in the Gulf, consequent to the peculiar and rare situation in which a hurricane stalls for several days on a coastal area, displaying minimal progress. This is particularly dangerous because: a) terrestrial drag affect only part of the circulation so that weakening is not as fast as when there is rapid motion inland; b) a part of the circulation remains over water and as such low-level moisture flux convergence is maintained; c) precipitation lasts for a very long time over the same area; d) moisture recycling from flooded areas further contributes to maintaining flood conditions. From the forecast perspective, hurricanes with track singularities such as stalling points, loops, or multiple landfalls over the same region are exceptionally difficult to predict because they result from small variations of opposing mechanisms and are therefore intrinsically unstable. Intensity fluctuations and track singularities are strongly coupled. The representation of fluxes undergoing differential drag (ocean vs. land) inside the same circulation is critically dependent upon resolution. Recent modeling development efforts within the NASA Global Modeling and Assimilation Office (GMAO) have brought the horizontal resolution of the global cubed-sphere Goddard Earth Observing System (GEOS) atmospheric model down to 12 km. This has led to innovative findings in the field of tropical cyclone (TC) modeling within a global framework. The GEOS atmospheric model and data assimilation system (DAS), used to produce the MERRA-2 (Modern Era Retrospective analysis for Research and Applications) reanalysis, has a "Replay" capability that adds a forcing term to the model equations to constrain it to follow a specified trajectory. This capability is used to produce global, downscaled versions of the 50 km MERRA-2 reanalysis at a 12.5 km resolution. The full DAS and replay modes both use the Incremental Analysis Update or IAU for applying an analysis increment to the model's background state four times each day at 00Z, 06Z, 12Z and 18Z. In a full DAS mode, a six-hour free run of the model is made (Predictor Step), writing states that will be used as backgrounds or first-guesses for the analysis centered on 00Z. The Analysis system then produces its best estimate of the 00Z state by blending backgrounds with observations. The difference between this best estimate and the model's original first guess at 00Z is treated as an analysis increment and added uniformly as a forcing term to the model during a second, now forced, run of the 12Z to 03Z window, which produces the assimilated results. Replay mode is very similar, except only the 00Z background is needed from the predictor segment, which is thus half as long, and instead of doing an analysis, a previously analyzed state at 00Z is read in. The increment is then computed and applied to the corrector step in the same way as in the original DAS. When downscaling horizontal spatial filters can be applied to the "analysis" increment used to force the Corrector step, so as to reproduce the full meteorological state or selectively choose or discard certain forcings. In this presentation, full global hindcast simulations of Juan and operational forecasts of Harvey are compared.

Flooding↗

Estimation of Seasonal Snow Water Equivalent Using Landsat Observations

This work presents a methodology for estimating seasonal snow water equivalent (SWE) from the use of remotely sensed Visible and Near Infrared observations from the Landsat mission. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data over the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlations with in situ SWE observations of greater than 0.95. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (32 years), and accuracy provide a unique dataset for investigating snow processes. In particular, this presentation illustrates how the reanalysis dataset was used to provide climatology of the seasonal snowfall accumulation rates, distributions, and variability over the last three decades.

SWE↗

Effects of Resolution and Spectral Nudging in Simulation the Effects of Wintertime Atmospheric River Landfalls in the Western US

Landfalling atmospheric rivers (ARs) play a crucial role in the climate of the US Pacific coast region as they are frequently related with heavy precipitation and flash flooding events. Thus, the capability of climate models to accurately simulate AR landfalls and their key hydrologic effects is an important practical concern for WUS, from flood forecasting to future water resources projections. In order to examine the effects of model configuration, including the resolution and spectral nudging, in simulating the climatology of key weather events in the conterminous US, a NASA team has performed a hindcast experiment using the GEOS5 global and the NU-WRF regional models for Nov 1999 - Oct 2010. This study examines the skill of these hindcasts, with different models and their configurations, in simulating key footprints of landfalling ARs in the WUS region. Using an AR-landfall chronology based on the vertically-integrated water vapor flux calculated from the MERRA2 reanalysis, we have analyzed the observed and simulated precipitation and temperature anomalies associated with wintertime AR landfalls along the US Pacific coast. Model skill is measured using metrics including regional means, a skill score based on correlations and mean-square errors, and Taylor diagrams in four WUS Bukovsky regions. Results show that the AR-related anomalies of precipitation is more reliable than of surface temperatures. Model skill also varies according to regions. The AR temperature anomalies are well simulated in most of the WUS region except PNW. For precipitation, simulations with finer spatial resolution tend to generate larger spatial variability and agree better with the PRISM data in most regions. Such a resolution dependence of spatial variability is not found for temperatures; e.g., the MERRA2 reanalysis often outperforms, with similar spatial variability and higher pattern correlations with the PRISM data, finer-resolution NU-WRF runs in simulating temperature variations within subregions. Results from this study will be summarized to assist future (regional) climate experiments for climate change impact assessments and developing adaptationmitigation strategies, the key elements of the National Climate Assessment.

Simulation↗

Earth Global Reference Atmospheric Model (GRAM) Overview and Comparison to Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2)

Earth-GRAM (Earth-Global Reference Atmospheric Model) Overview: Provides monthly statistics at any point in the atmosphere; Monthly, geographic, altitude variation; Current Version - Earth-GRAM 2016, C++; Output Includes - pressure, density, temperature, horizontal and vertical winds, speed of sound, atmospheric constituents; Used by engineering community to create atmospheric dispersions at a rapid runtime; Not a forecast model. MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2) Background: Developed by NASA Goddard Modelling and Assimilation Office; Horizontal Resolution: 0.625 degrees by 0.5 degrees longitude-by-latitude grid (NCEP (National Centers for Environmental Prediction) reanalysis I vs. 2.5 degrees by 2.5 degrees currently used in Earth-GRAM); Vertical resolution: 72 model layers or interpolated to 42 pressure levels to 0.1 hectopascals (hPa) (NCEP reanalysis I, surface to 10hPa at 17 pressure levels).

White, Patrick W.↗

A Global Climatology of Extratropical Transition Part II: Statistical Performance of the Cyclone Phase Space

This study analyzes the differences between an objective, automated identification of tropical cyclones (TCs) that undergo extratropical transition (ET), and the designation of ET determined subjectively by human forecasters in best-track data in all basins globally. The objective identification of ET is based on the cyclone phase space (CPS), calculated from the Japanese 55-year Reanalysis (JRA-55) or the ECMWF Interim Reanalysis (ERA-Interim). The resulting classification into “ET storms” and “non-ET storms” underlies the global climatology of ET presented in Part I of this study. Here, the authors investigate how well the CPS classifications agree with those in the best-track records calculated from JRA-55 or from ERA-Interim data. According to F1 scores and Matthews correlation coefficients (MCCs), the classification of ET storms in the CPS agrees best with the best-track classification in the Western North Pacific (MCC > 0.7) and the North Atlantic (MCC > 0.5). In other basins, the correlation between the CPS classification and the best-track classification is only slightly higher than that of a random classification. The JRA-55 classification achieves higher performance scores than does the ERA-Interim classification, and the differences are statistically significant in all basins. The lower performance of ERA-Interim is mainly due to a higher false alarm rate, particularly in the Eastern North Pacific. Overall, the results show that while the CPS-based classifications are good enough to be useful for many purposes, there is almost certainly room for improvement – in the representation of the storms in reanalyses, in our objective metrics of ET, and in our scientific understanding of the ET process.

Melanie Bieli↗

Annual Cycle of Planetary Boundary Layer Heights estimated from NOAA Profiler Network Data

The Planetary Boundary Layer (PBL) height over the Unites States Great and Central Plains during 1992-2012 is examined here using a combination of Wind Profiler-derived (WP) PBL height estimates and reanalysis fields from the Modern Era Retrospective Reanalysis Version 2 (MERRA-2). The combined analysis allows process study of the reasons behind the monthly mean behavior of the observed PBL heights under clear-sky conditions as well as the PBL height variability. WP PBL height monthly mean annual cycles were grouped into general categories of behavior, each analyzed using MERRA-2 fields of sensible and latent heat flux, surface temperature, net radiation and soil moisture. In the ’canonical’ category the latent heat plays little role in the determination of the monthly mean PBL height, and it follows the annual cycle of the net radiation. In the other categories, precipitation and latent heat flux had more influence in setting the annual cycle. An analysis of variance revealed that the role of latent heat in determining the PBL height variations is large (explaining up to 40% of PBL height variability) even in the ’canonical’ category for which latent heat played no role in setting the monthly mean. In other categories the latent heat explained up to 80% of the PBL height variations. The amount of that influence is shown to be related to the variability of column soil moisture.

MERRA-2↗

Assessment of Precipitation Anomalies in California Using TRMM and MERRA Data

After more than a decade of moderate seasonal deviations from the expected climate, it is easy to forget that California is actually prone to instabilities in precipitation patterns that occur on various scales. Using modern satellite and reanalysis data we reassess certain aspects of the precipitation climate in California from the past three decades. California has a well-pronounced rain season that peaks in December-February. However, the 95% confidence interval around the climatological precipitation during these months imply that deviations on the order of 60% of the expected amounts are very likely during the most important period of the rain season. While these positive and negative anomalies alternate almost every year and tend to cancel each other, severe multi-year declines of precipitation in California seem to appear on decadal scales. The 1986-1994 decline of precipitation was similar to the current one that started in 2011, and is apparent in the reanalysis data. In terms of accumulated deficits of precipitation, that episode was no less severe than the current one. While El Niño (the warm phase of the El Niño Southern Oscillation, ENSO) is frequently cited as the natural forcing expected to bring a relief, our assessment is that ENSO has been driving at best only 6% of precipitation variability in California in the past three decades. It means El Niño needs to be stronger and longer, in order to have a higher likelihood of a positive impact, and the current one does not match these criteria. Using fractional risk analysis of precipitation populations during normal and dry periods, we show that the likelihood of losing the most intensive precipitation events drastically increases during the multi-year drying events. Since storms delivering up to 50% of precipitation in California are driven by atmospheric rivers making landfall, thus the importance of their suppression and blockage by persistent ridges of atmospheric pressure in the northeast Pacific.

Savtchenko, Andrey K.↗

Toward Coupled Data Assimilation in NASA’s GEOS: Developments in the Ocean Context

The Global Modeling & Assimilation Office (GMAO) at NASA GSFC produces analyses and predictions of the Earth system using various configurations of the Goddard Earth Observing System (GEOS) model and assimilation system. The current sub-seasonal-to-seasonal prediction system (GEOS-S2S) is based on a coupled atmosphere-ocean-land-ice configuration of GEOS which includes the Modular Ocean Model version 5 (MOM5) run at approximately 50-km resolution and a de-coupled OI-based ocean analysis that uses an initialization of MOM5 forced by the MERRA-2 reanalysis. GMAO will soon implement an updated GEOS-S2S system that will run at 25-km resolution and adopt aspects of the hybrid four-dimensional ensemble-variational (H4DEnVar) system already running in the production-version atmospheric analysis system, including a Local Ensemble Transform Kalman Filter (LETKF) to provide initial conditions for the oceanic state. This presentation will focus on developments to sustain the GMAO's systems on longer time horizons, where more radical transformations will be required to adapt to advanced computing environments, higher resolution and more diverse model components, and new observations for the Earth system. Results will describe progress toward a version of the GEOS coupled system that will be based around the Joint Effort for Data assimilation Integration (JEDI) framework being developed within Joint Center for Satellite Data Assimilation (JCSDA) and include an updated ocean model, MOM6. Discussion will focus specifically on the use of a Unified Forward Operator (UFO) for simulating observations and the Object Oriented Prediction System (OOPS) for providing the state estimate. These features are being developed as a multi-agency effort under the auspices of the JCSDA and are being adopted in the GMAO for all its applications of coupled data assimilation including S2S, numerical weather prediction, and reanalysis.

Mahajan, Rahul↗

Impact of Bias in the Marine Air Temperature Observation Set on Atmospheric Reanalyses

Marine air temperature observations are collected by volunteer observing ships carrying automatic meteorological sensor packages. Close proximity of the package mount to a dark colored ship's deck is known to introduce a warm bias on sunny afternoons in summer under low wind conditions. This observation set is used as an influential constraint in the MERRA-2 reanalysis, as well as several other reanalyses, in the poorly sampled surface marine boundary layer. Here we explore the consequences of bias in the marine air temperature observation set during the 36 year period 1980-2015. Marine air temperature observations have mainly been collected along major shipping routes, with the densest coverage in the North Atlantic and North Pacific. In the first part of this presentation we review the historical observation coverage and compare air temperature and contemporaneous SST measurements to identify systematic differences. The comparison suggests the presence of diurnal bias in the air temperature observations peaking in subtropics and midlatitudes latitudes in mid-afternoon local time. In the second part we explore the impact of the marine air temperature observations on MERRA-2 through analysis of the terms in the temperature forecast tendency equation, the equation which describes the contributions to changes in temperature during each analysis cycle. Examination of the temperature forecast tendency equation allows us to separate systematic errors due to such processes and radiation, mixing, and advection and allows us to identify the major source of systematic forecast tendencies. The results show that the marine air temperature measurements introduce an error that can reach 0.4K/dy which also peaks in local afternoon at subtropical latitudes. This bias varies with season and has gradually reduced from the beginning of the MERRA-2 reanalysis in 1980. In the final part we conduct numerical experiments with a version of the MERRA-2 analysis system, comparing experiments when the marine air temperature observations are included or withheld from the analysis system.

Akella, Santharam↗