Comparison of Climatological Planetary Boundary Layer Depth Estimates Using the GEOS-5 AGCM

Planetary boundary layer (PBL) processes, including those influencing the PBL depth, control many aspects of weather and climate and accurate models of these processes are important for forecasting changes in the future. However, evaluation of model estimates of PBL depth are difficult because no consensus on PBL depth definition currently exists and various methods for estimating this parameter can give results that differ by hundreds of meters or more. In order to facilitate comparisons between the Goddard Earth Observation System (GEOS-5) and other modeling and observational systems, seven PBL depth estimation methods are used to produce PBL depth climatologies and are evaluated and compared here. All seven methods evaluate the same atmosphere so all differences are related solely to the definition chosen. These methods depend on the scalar diffusivity, bulk and local Richardson numbers, and the diagnosed horizontal turbulent kinetic energy (TKE). Results are aggregated by climate class in order to allow broad generalizations. The various PBL depth estimations give similar midday results with some exceptions. One method based on horizontal turbulent kinetic energy produces deeper PBL depths in the winter associated with winter storms. In warm, moist conditions, the method based on a bulk Richardson number gives results that are shallower than those given by the methods based on the scalar diffusivity. The impact of turbulence driven by radiative cooling at cloud top is most significant during the evening transition and along several regions across the oceans and methods sensitive to this cooling produce deeper PBL depths where it is most active. Additionally, Richardson number-based methods collapse better at night than methods that depend on the scalar diffusivity. This feature potentially affects tracer transport

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

Nasa GEOS Model for Sub-Seasonal to Seasonal Predictions: The Major Teleconnections, Tropical Cyclone Activity, and ENSO

Ability of the NASA GEOS seasonal prediction model system (the latest version) in predicting the large-scale teleconnections, tropical cyclone (TC) activity, and ENSO are validated. The major large-scale teleconnections include the North Atlantic Oscillation (NAO), the Arctic Oscillation (AO), and the Pacific North American (PNA) that span vast geographical area across the North Pacific/Atlantic and North America. Predictive skill of TC activity is assessed by Genesis potential index (GPI). Anomaly correlations are greater than 0.5 for winter teleconnections at up to 2 month lead and for the first four month (June- July-August-September) GPIs over the North Atlantic and the Western Pacific. Correlation for the long-range prediction of the ENSO (Nino3.4 SST) maturity reaches 0.9 at 6 month lead and 0.8 at 9 month lead

Horizontal Variability of Water and Its Relationship to Cloud Fraction near the Tropical Tropopause: Using Aircraft Observations of Water Vapor to Improve the Representation of Grid-scale Cloud Formation in GEOS-5

Large-scale models such as GEOS-5 typically calculate grid-scale fractional cloudiness through a PDF parameterization of the sub-gridscale distribution of specific humidity. The GEOS-5 moisture routine uses a simple rectangular PDF varying in height that follows a tanh profile. While below 10 km this profile is informed by moisture information from the AIRS instrument, there is relatively little empirical basis for the profile above that level. ATTREX provides an opportunity to refine the profile using estimates of the horizontal variability of measurements of water vapor, total water and ice particles from the Global Hawk aircraft at or near the tropopause. These measurements will be compared with estimates of large-scale cloud fraction from CALIPSO and lidar retrievals from the CPL on the aircraft. We will use the variability measurements to perform studies of the sensitivity of the GEOS-5 cloud-fraction to various modifications to the PDF shape and to its vertical profile

Regional Replay: A Unique Reanalysis-Based Tool for Addressing Model Error

Understanding and correcting errors in general circulation and climate models has long been part intuition and part trial and error. Efforts to diagnose the errors and provide some guidance to developers have been of some value, though such efforts, with few exceptions, have been more successful in identifying and documenting the errors in the model simulations rather than the model deficiencies that produced them. Modern atmospheric reanalyses such as MERRA-2 provide much-improved estimates of our climate system at hourly to interannual and longer time scales and have become an important tool for assessing model performance. Here we use MERRA-2 to address biases in the NASA/GMAO GEOS model by employing a "regional replay" approach developed in the GMAO. The regional replay approach constrains the model to remain close to the reanalysis over arbitrary regions and selected model variables, thus allowing us to examine how model error generated over one area is spatially translated across the globe. Several examples are given including an assessment of the global impact of errors produced over the Tibet region

The Roles of Climate Change and Climate Variability in the 2017 Atlantic Hurricane Season

The 2017 hurricane season was extremely active with six major hurricanes, the third most on record. The sea-surface temperatures (SSTs) over the eastern Main Development Region (EMDR), where many tropical cyclones (TCs) developed during active months of August/September, were approximately 0.96 degrees Centigrade above the 1901-2017 average (warmest on record): about 0.42 degrees Centigrade from a long-term upward trend and the rest (around 80 percent) attributed to the Atlantic Meridional Mode (AMM). The contribution to the SST from the North Atlantic Oscillation over the EMDR was a weak warming, while that from ENSO was negligible. Nevertheless, ENSO, the NAO, and the AMM all contributed to favorable wind shear conditions, while the AMM also produced enhanced atmospheric instability. Compared with the strong hurricane years of 2005-2010, the ocean heat content (OHC) during 2017 was larger across the tropics, with higher SST anomalies over the EMDR and Caribbean Sea. On the other hand, the dynamical/thermodynamical atmospheric conditions, while favorable for enhanced TC activity, were less prominent than in 2005-2010 across the tropics. The results suggest that unusually warm SST in the EMDR together with the long fetch of the resulting storms in the presence of record-breaking OHC were key factors in driving the strong TC activity in 2017

Improved Boundary Layer Depth Retrievals from MPLNET

Continuous lidar observations of the planetary boundary layer (PBL) depth have been made at the Micropulse Lidar Network (MPLNET) site in Greenbelt, MD since April 2001. However, because of issues with the operational PBL depth algorithm, the data is not reliable for determining seasonal and diurnal trends. Therefore, an improved PBL depth algorithm has been developed which uses a combination of the wavelet technique and image processing. The new algorithm is less susceptible to contamination by clouds and residual layers, and in general, produces lower PBL depths. A 2010 comparison shows the operational algorithm overestimates the daily mean PBL depth when compared to the improved algorithm (1.85 and 1.07 km, respectively). The improved MPLNET PBL depths are validated using radiosonde comparisons which suggests the algorithm performs well to determine the depth of a fully developed PBL. A comparison with the Goddard Earth Observing System-version 5 (GEOS-5) model suggests that the model may underestimate the maximum daytime PBL depth by 410 m during the spring and summer. The best agreement between MPLNET and GEOS-5 occurred during the fall and they diered the most in the winter

The Impact of SST-Forced and Unforced Teleconnections on 2015/16 El Nino Winter Precipitation over the Western United States

The factors impacting western U.S. winter precipitation during the 2015/16 El Nio are investigated using the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data, and simulations with the Goddard Earth Observing System version 5 (GEOS-5) atmospheric general circulation model forced with specified sea surface temperatures (SSTs). Results reveal that the simulated response to the tropical Pacific SST associated with the 2015/16 El Nio was to produce wetter than normal conditions over much of the west coast including California - a result at odds with the negative precipitation anomalies observed over much of the Southwestern U.S. It is shown that two factors acted to partly counter the canonical ENSO response in that region. First, a potentially predictable but modest response to the unusually strong and persistent warm SST in the northeastern Pacific decreased precipitation in the Southwestern U.S. by increasing sea level pressure, driving anticyclonic circulation and atmospheric descent, and reducing moisture transport into that region. Second, large-scale unforced (by SST) components of atmospheric variability (consisting of the leading modes of unpredictable intra-ensemble variability) resembling the positive phase of the North Atlantic Oscillation and Arctic Oscillation are found to be an important contributor to the drying over the western U.S. While a statistical reconstruction of the precipitation from our simulations that account for internal atmospheric variability does much to close the gap between the ensemble mean and observed precipitation in the Southwestern U.S., some differences remain, indicating that model error is also playing a role

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 GrI

Large-Scale Controls on Atlantic Tropical Cyclone Activity on Seasonal Time Scales

Interannual variations in seasonal tropical cyclone (TC) activity (e.g., genesis frequency and location, track pattern, and landfall) over the Atlantic are explored by employing observationally-constrained simulations with the NASA Goddard Earth Observing System version (GEOS-5) atmospheric general circulation model. The climate modes investigated are El Nino-Southern Oscillation (ENSO), the North Atlantic Oscillation (NAO), and the Atlantic Meridional Mode (AMM). The results show that the NAO and AMM can strongly modify and even oppose the well- known ENSO impacts, like in 2005, when a strong positive AMM (associated with warm SSTs and a negative SLP anomaly over the western tropical Atlantic), led to a very active TC season with enhanced TC genesis over the Caribbean Sea and a number of landfalls over North America, under a neutral ENSO condition. On the other end, the weak TC activity during 2013 (characterized by weak negative Nio index) appears caused by a NAO-induced positive SLP anomaly with enhanced vertical wind shear over the tropical North Atlantic. During 2010, the combined impact of the three modes produced positive SST anomalies across the entire low-latitudinal Atlantic and a weaker subtropical high, leading to more early recurvers and thus fewer landfalls despite enhanced TC genesis. The study provides evidence that TC number and track are very sensitive to the relative phases and intensities of these three modes, and not just to ENSO alone. Examination of seasonal predictability reveals that predictive skill of the three modes is limited over tropics to sub-tropics, with the AMM having the highest predictability over the North Atlantic, followed by ENSO and NAO