Search

Atmospheric reanalyses are widely used for a variety of scientific endeavors in the Arctic and Antarctic. Reanalyses are used as boundary conditions for a regional and process-based models, for climate model validation, and for diagnostic analysis of physical processes, weather and climatic events. However, reanalyses are typically global and often do not account for specific, regional considerations, such as for polar regions. In this work, we provide a brief evaluation of a prototype for a new GMAO reanalysis, which incorporates higher spatial resolution, an updated approach for data assimilation, and a revised atmospheric model. We identify differences in the representation of the Arctic atmosphere in comparison to recent reanalyses. Furthermore, we provide a forum for Arctic scientists to consider the future improvements for reanalyses, and seek feedback for the following questions: 1) What are important performance factors to consider in evaluating new reanalyses? 2) What physical processes should be incorporated into new reanalyses? 3) What spatio-temporal scales should be considered?

Recently NASA's Global Modeling and Assimilation Office (GMAO) has developed a new Subseasonal to Seasonal Prediction system Version 3 (GEOS-S2S-3). This upgrade replaces the GEOS-S2S-2 which is NASA's current contribution to the North American Multi-Model Experiment seasonal prediction project (Kirtman et al., 2014). The main improvements for our S2S-3 system include 1) a higher resolution MOM5 (Griffies et al., 2005) ocean model (now 0.25o x 0.25o x 50 layers), 2) an improved atmospheric/ocean interface layer (Akella and Suarez, 2018), and 3) assimilation of a long-track satellite salinity into the ocean model (Hackert et al, 2019). Atmospheric forcing is provided by the NASA MERRA-2 reanalysis (Gelaro et al., 2017). Initialization for the ocean relies on the GMAO ocean reanalysis system which assimilates all available in situ temperature and salinity, satellite sea surface salinity, and sea level using the Local Ensemble Transform Kalman Filter (LETKF) implementation of (Penny et al., 2013) on a 5 day assimilation cycle with 20 fixed ensemble members.In this presentation, we will authenticate our new S2S-3 ocean reanalysis using standard GODAE validation metrics. For example, we will compare gridded fields of mean and standard deviation of the ocean reanalysis versus observed fields. We will show correlation/RMS of model versus observations and temperature and salinity mean profiles for the various basins and latitude bands. Basin-scale volume transports, such as the Atlantic Meridional Overturning Circulation and the Indonesian Throughflow will be validated. Equatorial ocean waves will be compared by decomposing sea level into Kelvin and Rossby components. For each of these metrics, we plan to validate the results and then compare our new S2S-3 against the current production version, S2S-2. Finally, we will compare 9-month seasonal forecasts initialized from these two systems for the tropical Pacific NINO3.4 region over the period 1981-present.

The NASA/Goddard Global Modeling and Assimilation Office (GMAO) released Version 2 of the Subseasonal to Seasonal (GEOS-S2S) forecast system in the fall of 2017, and it has been producing near-real time subseasonal to seasonal forecasts and a weakly coupled atmosphere-ocean data assimilation record since then. A new version of the coupled modeling and analysis system (Version 3) was released by the GMAO at the end of 2019. The new version runs at higher oceanic resolution than the previous (approximately 1/2 degree for the atmosphere, 1/4 degree for the ocean), and includes interactive earth system model components not typically present in seasonal prediction systems (two moment cloud microphysics for aerosol indirect effect and an interactive aerosol model). The weakly coupled atmosphere-ocean data assimilation system now includes assimilation of sea surface salinity, that has been shown to result in improved ocean mixed layer simulation and ENSO prediction skill.

Observational paucity is a reality in the Arctic Ocean. This is especially true for near-surface variables such as temperature, moisture, heat fluxes and BL clouds. As a result, modeling has become one of the major avenues for understanding current and future Arctic trends. Reanalyses are frequently used to force global ocean circulation and sea-ice models. But in northern high latitudes, model integrations and reanalyses are known to have large uncertainties in temperature and humidity profiles, and in boundary layer cloudiness. These are common sources of error in the surface radiative budget terms.An important way to diagnose these biases spatially and temporally is by using satellite remote sensing data. However, remotely-sensed observations also have large uncertainties, especially in near-surface temperature and relative humidity profiles.In situ observational studies are important in bridging our knowledge gap in regions such as the Arctic Ocean. Here, we utilize airborne and ship observations from the ARISE, ACME-IV, and ASCOS campaigns to construct critical relative humidity (RH) profiles over the Beaufort Sea. Such profiles are used as parameterization inputs in the NASA GOES global model to derive the total water condensate in a model grid-box, which determines the cloud fraction. Currently, the critical RH profiles are derived by global AIRS data, relaying mostly on mid-latitude regions, which are not necessarily relevant to the Arctic.We derive campaign-wide mean, standard deviation, and critical RH values, for grid size of 50x50 km and altitude bins between 50 to 400 m, covering both open-ocean and sea-ice covered regions. We compare profiles over open ocean and sea-ice, and look at correlations between the observed critical RH values and water condensate (by cloud number concentration) from observations versus the modeled ones. We then input our calculated values of minimal critical RH values into a set of GEOS single column model (SCM) simulations over the ARISE and ASCOS regions and compare the differences between the predicted values of cloud liquid water path (LWP), ice water path (IWP) and surface fluxes with the observed ones under the range of input parameterizations. Finally, we discuss the implications on surface radiative budget predictions in this region.

Arctic clouds play an important role in modifying the surface energy balance. In the Arctic, clouds are thought to influence the underlying sea ice cover through changing downwelling longwave radiative fluxes to the surface and through the selective reflection of the shortwave flux in summer. Atmospheric reanalyses are generally thought to have a poor representation of cloud processes at high latitudes, although the representation of trends over the perennial Arctic sea ice pack is less well known. Here, atmospheric energy fluxes are examined at the top of the atmosphere from contemporary reanalyses in comparison to satellite measurements from the CERES-EBAF version 4.1 product. The principal reanalyses examined are the NASA MERRA-2, the ECMWF ERA5 and ERA-Interim, the JRA-55, and the regional Arctic System Reanalysis version 2. In agreement with previous observation-based studies, changes with time in the shortwave cloud radiative forcing in reanalyses are found to be negligible despite strong trends in the absorbed shortwave. Over the full satellite period, there is large disagreement in the seasonality of longwave cloud forcing trends. These trends are reduced during the CERES-EBAF observing period (2003-present). An examination of these trends with respect to sea ice cover changes in each of the reanalyses is conducted.

Reanalyses are regularly gridded, retrospective depictions of the physical earth system, which are produced through the correction of a short-term forecast to available observations. In the Arctic, reanalyses are particularly well suited to marshal the sparse observing network to provide a plausible, multivariate representation of conditions. Atmospheric reanalyses such as MERRA-2 (NASA Modern-Era Retrospective analysis for Research and Applications, version 2) and ocean reanalyses such as SODA3 (Univ. Maryland Simple Ocean Data Assimilation version 3) are widely used in Arctic research for diagnostic studies of circulation, model evaluation, and as boundary conditions for a variety of process models. Here, we provide examples that illustrate the utility of reanalyses for providing information on the spatial and temporal scales of recent, rapid changes in the Arctic. Recent trends in Arctic surface temperatures, surface melt over Greenland and Arctic glaciers, and evolving freshwater conditions in the Arctic Ocean are examples where reanalyses can provide information that cannot easily be obtained via other means. These examples provide information on the scale, magnitude, and the uncertainty of recent Arctic change and provide a context for future scenarios. We further quantify uncertainties in key reanalyses variables and approaches for addressing these issues.

Sea ice is considered an important indicator of climate change. Interest has focused on the trend towards reduced sea ice cover in the Arctic Ocean, particularly during summer months. In the Southern Hemisphere, sea ice has also received attention due to a trend towards increasing cover, in opposition to the Northern Hemisphere and to climate model projections. This all changed in the austral spring 2016, when the Southern Ocean sea ice cover dropped dramatically, and has remained consistently at below average values up to the present. Here I review some of the issues and perceived causes of the 2016 sea ice reduction. Several studies have suggested a change in surface winds, which would allow warmer waters to move upwards in the water column and towards the surface. Results from the GMAO (Global Modeling and Assimilation Office) ocean data assimilation system are consistent with this hypothesis, and indicate a significant warming of subsurface waters occurred in the winter of 2016, and has remained in place.

A new version of the coupled modeling and analysis system used to produce near real time subseasonal to seasonal forecasts was recently released by the NASA/Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution than the previous, (approximately 1/2 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter, and now includes the assimilation of along-track sea surface height. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a series of retrospective seasonal forecasts. Results show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. Analysis of the ensemble spread shows improvements relative to the previous system, including generally better reliability. The precipitation forecast skill is comparable to previous S2S systems, and the only tradeoff is an increased "double ITCZ", which is expected as we go to higher atmospheric resolution.

An accurate historical record of evolving Arctic conditions is integral to furthering our understanding of climate processes and to providing a foundation for predicting future climate scenarios in northern high latitudes. Atmospheric reanalyses are seen as an important source of information on the recent past for the data-sparse Arctic region. An assessment of near-surface Arctic air temperatures finds significant discrepancies among the various modern reanalyses. An important point is the treatment of surface boundary conditions: specifically, the sea ice cover and sea surface temperatures (SSTs) over the Arctic Ocean. Reanalyses use different methodologies and data sources for SSTs and sea ice concentration boundary forcing. Notably, the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) and the European Centre for Medium-Range Weather Forecasts Interim Re-Analysis (ERA-Interim) both use boundary forcing derived from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) over an extended, overlapping period of time. This allows for an examination of differences between the two systems while both concurrently employ the same fractional sea ice coverage. To further understand these differences, an ensemble of AMIP-style simulations using the MERRA-2 atmospheric model - but without data assimilation - shows considerable differences in Arctic temperatures as compared to reanalyses, particularly in autumn and winter months. Results from the AMIP simulations suggest that the surface representation over sea ice used in the MERRA-2 model provides an intrinsic warm bias and obfuscates Arctic Amplification, an established feature present in observations and reanalyses. An additional ensemble of AMIP-style simulations using the MERRA-2 atmospheric model was performed using boundary conditions derived from the ERA-Interim reanalysis. An in-depth comparison of surface temperatures over the Arctic from the two reanalyses and two AMIP-style ensembles will be presented, along with an assessment of the effects of the varying Arctic temperature time series on the atmospheric general circulation and energy budget.

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

Precipitation over the Arctic Ocean has a significant impact on the basin-scale freshwater and energy budgets but is one of the most poorly constrained variables in atmospheric reanalyses. Precipitation controls the snow cover on sea ice, which impedes the exchange of energy between the ocean and atmosphere, inhibiting sea ice growth. Thus, accurate precipitation amounts are needed to inform sea ice modeling, especially for the production of thickness estimates from satellite altimetry freeboard data. However, obtaining a quantitative estimate of the precipitation distribution in the Arctic is notoriously difficult because of a number of factors, including a lack of reliable, long-term in situ observations; difficulties in remote sensing over sea ice; and model biases in temperature and moisture fields and associated uncertainty of modeled cloud microphysical processes in the polar regions. Here, we compare precipitation estimates over the Arctic Ocean from eight widely used atmospheric reanalyses over the period 2000–16 (nominally the ‘‘new Arctic’’). We find that the magnitude, frequency, and phase of precipitation vary drastically, although interannual variability is similar. Reanalysis-derived precipitation does not increase with time as expected; however, an increasing trend of higher fractions of liquid precipitation (rainfall) is found. When compared with drifting ice mass balance buoys, three reanalyses (ERA-Interim, MERRA, and NCEP R2) produce realistic magnitudes and temporal agreement with observed precipitation events, while two products [MERRA, version 2 (MERRA-2), and CFSR] show large, implausible magnitudes in precipitation events. All the reanalyses tend to produce overly frequent Arctic precipitation. Future work needs to be undertaken to determine the specific factors in reanalyses that contribute to these discrepancies in the new Arctic.

The major stratospheric sudden warming in February 2018 occurred in unison with changes throughout the troposphere in the northern hemisphere. Changes in large-scale dynamics and regional weather patterns are described in this study, along with features in major features of the Earth system, such as sea ice. The focus is on analyses using the Goddard Earth Observing System (GEOS) analyses, weather prediction, and seasonal forecasts. Specific topics to discuss include: - The morphology of the sudden warming, including planetary waves in the troposphere and stratosphere.- Relationships with weather over North America. - Warming over the Arctic and anomalous sea-ice distributions. - Predictability on the one-to-four week timescales. The predictions will be based on routine products generated using the GEOS systems, as well as targeted experiments that focus on the roles of different mechanisms (such as ice-atmosphere interactions).

In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System with a focus on the computing time and resources and actual time it takes to complete a full set of hindcasts. The goal is to come up with some solutions to allow us to run more ensemble members for the next version of the system which will be higher resolution and take many more resources.

In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System, current users and products, and methods for validation and evaluation of the system. Methods for evaluation include baseline evaluations metrics, the ability to simulate key modes of variability, and evaluation of new development areas.

A new version of the modeling and analysis system used to produce sub-seasonal to seasonal forecasts has just been released by the NASA Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution (approximately 12 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a free-running coupled simulation with the new system and results from a series of retrospective seasonal forecasts. Results from retrospective forecasts show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. The precipitation forecast skill is comparable to previous S2S systems, and the only trade off is an increased double ITCZ, which is expected as we go to higher atmospheric resolution.

The inter-relationship between subtropical western-central Pacific sea surface temperatures (STWCPSST), sea ice concentration in the Beaufort Sea (SICBS), and the North Atlantic Oscillation (NAO) are investigated for the last 37 summers and winters (1980-2016). Lag-correlation of the STWCPSST×(-1) in spring with the NAO phase and SICBS in summer increases over the last two decades, reaching r = 0.4-0.5 with significance at 5 percent, while winter has strong correlations in approximately 1985-2005. Observational analysis and the atmospheric general circulation model experiments both suggest that STWCPSST warming acts to increase the Arctic geopotential height and temperature in the following season. This atmospheric response extends to Greenland, providing favorable conditions for developing the negative phase of the NAO. SIC and surface albedo tend to decrease over the Beaufort Sea in summer, linked to the positive surface net shortwave flux. Energy balance considering radiative and turbulent fluxes reveal that available energy that can heat surface is larger over the Arctic and Greenland and smaller over the south of Greenland, in response to the STWCPSST warming in spring. XXXX Arctic & Atlantic: Positive upper-level height/T anomaly over the Arctic and Greenland, and a negative anomaly over the central-eastern Atlantic, resembling the (-) phase of the NAO. Pacific: The negative height/T anomaly over the mid-latitudes, along with the positive anomaly over the STWCP, where 1degC warming above climatology is prescribed. Discussion: It is likely that the Arctic gets warm and the NAO is in the negative phase in response to the STWCP warming. But, there are other factors (e.g., internal variability) that contribute to determination of the NAO phase: not always the negative phase of the NAO in the event of STWCP warming (e.g.: recent winters and near neutral NAO in 2017 summer).

Reanalyses in the Arctic are widely used for model evaluation and for understanding contemporary climate change. Nevertheless, differences among reanalyses in fundamental meteorological variables including surface air temperature are large. For example, the 1980-2009 mean surface air temperature for the north polar cap (70ÂdegN-90ÂdegN) among global reanalyses span a range of 2.4 K, which approximates the average warming trend from these reanalyses over the 30-year period of 2.1 K. Understanding these differences requires evaluation over the three principal surface domains of the Arctic: glaciated land, the unglaciated terrestrial surface, and sea ice/ocean. An examination is conducted of contemporary global reanalyses of the ECMWF Interim project, NASA MERRA, MERRA-2, JRA-55, and NOAA CFSR using available in situ data and assessments of the surface energy budget. Overly-simplistic representations of the Greenland Ice Sheet surface are found to be associated with local warm air temperature biases in winter. A review of progress made in the development of the MERRA-2 land-ice representation is presented. Large uncertainty is also found in temperatures over the Arctic tundra and boreal forest zone. But a key focus of temperature differences for northern high latitudes is the Arctic Ocean. Near-surface air temperature differences over the Arctic Ocean are found to be related to discrepancies in sea ice and sea surface temperature boundary data, which are severely compromised in current reanalyses. Issues with the modeled representation of sea ice cover are an additional factor in reanalysis temperature trends. Differences in the representation of the surface energy budget among the various reanalyses are also reviewed.

Dynamical ice sheet models are being used in simulations of future sea level change resulting from changing glacier mass. One of the difficulties in doing so are the input conditions obtained from earth system models. These inputs can be of coarse spatial resolution, and may not represent surface melt in a future climate. I review various methods for overcoming this with the aim of promoting discussion among modelers.

Reanalyses have become an integral tool for evaluating regional and global climate variations, and an important component of this is modifications to the energy budget. Reductions in Arctic Sea ice extent has induced an albedo feedback, causing the Arctic to warm more rapidly than anywhere else in the world, referred to as "Arctic Amplification." This has been demonstrated by observations and numerous reanalyses, including the Modern Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). However, the Arctic Amplification signal is non-existent in a ten member ensemble of the MERRA-2 Atmospheric Model Intercomparison Project (M2AMIP) simulations, using the same prescribed climate forcing, including Sea Surface Temperature (SST) and ice. An evaluation of the temperature tendency within the lower troposphere due to radiation, moisture, and dynamics as well as the surface energy budget in MERRA-2 and M2AMIP will demonstrate that despite identical prescribed SSTs and sea ice in both versions, enhanced warming in the Arctic in MERRA-2 is in response to the analysis increment tendency due to temperature observations. Furthermore, the role of boundary conditions, model biases and changes in observing systems on the Arctic Amplification signal will be assessed. Literature on the topic of Arctic Amplification demonstrates that the enhanced warming begins in the mid-1990s. Anomalously warm Arctic SST in the early time period of MERRA-2 can mute the trend in Arctic lower troposphere temperature without the constraint of observations in M2AMIP. Additionally, MERRA-2 uses three distinct datasets of SST and sea ice concentration, which also plays a role.

Improved seasonal forecasts of Arctic sea ice are important for regional stakeholders, but also for obtaining a better understanding of the Arctic climate system. An important part of the forecasts is the initial sea ice, ocean, and atmosphere initial conditions. I briefly give an overview of the initial conditions currently being used in seasonal sea ice predictions. I also identify available sources of observational data and prospects for coupled atmosphere/ocean reanalyses.

The GMAO seasonal forecast is produced from coupled model integrations that are initialized every five days, with seven additional ensemble members generated by coupled model breeding and initialized on the date closest to the beginning of the month. The main components of the AOGCM are the GEOS-5 atmospheric model, the MOM4 ocean model, and CICE sea ice model. Forecast fields were re-gridded to the passive microwave grid for averaging.

The Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) is a global atmospheric reanalysis produced by the NASA Global Modeling and Assimilation Office (GMAO). It spans the satellite observing era from 1980 to the present. The goals of MERRA-2 are to provide a regularly-gridded, homogeneous record of the global atmosphere, and to incorporate additional aspects of the climate system including trace gas constituents (stratospheric ozone), and improved land surface representation, and cryospheric processes. MERRA-2 is also the first satellite-era global reanalysis to assimilate space-based observations of aerosols and represent their interactions with other physical processes in the climate system. The inclusion of these additional components are consistent with the overall objectives of an Integrated Earth System Analysis (IESA). MERRA-2 is intended to replace the original MERRA product, and reflects recent advances in atmospheric modeling and data assimilation. Modern hyperspectral radiance and microwave observations, along with GPS-Radio Occultation and NASA ozone datasets are now assimilated in MERRA-2. Much of the structure of the data files remains the same in MERRA-2. While the original MERRA data format was HDF-EOS, the MERRA-2 supplied binary data format is now NetCDF4 (with lossy compression to save space).