Drought Predictability and Prediction in a Changing Climate: Assessing Current Predictive Knowledge and Capabilities, User Requirements and Research Priorities

Drought is fundamentally the result of an extended period of reduced precipitation lasting anywhere from a few weeks to decades and even longer. As such, addressing drought predictability and prediction in a changing climate requires foremost that we make progress on the ability to predict precipitation anomalies on subseasonal and longer time scales. From the perspective of the users of drought forecasts and information, drought is however most directly viewed through its impacts (e.g., on soil moisture, streamflow, crop yields). As such, the question of the predictability of drought must extend to those quantities as well. In order to make progress on these issues, the WCRP drought information group (DIG), with the support of WCRP, the Catalan Institute of Climate Sciences, the La Caixa Foundation, the National Aeronautics and Space Administration, the National Oceanic and Atmospheric Administration, and the National Science Foundation, has organized a workshop to focus on: 1. User requirements for drought prediction information on sub-seasonal to centennial time scales 2. Current understanding of the mechanisms and predictability of drought on sub-seasonal to centennial time scales 3. Current drought prediction/projection capabilities on sub-seasonal to centennial time scales 4. Advancing regional drought prediction capabilities for variables and scales most relevant to user needs on sub-seasonal to centennial time scales. This introductory talk provides an overview of these goals, and outlines the occurrence and mechanisms of drought world-wide

An analysis of tropopause pressure and total ozone correlations

A study of the relationship between total ozone and tropopause pressure was carried out using Total Ozone Mapping Spectrometer (TOMS) data and National Meteorological Center (NMC) global analyses. The medium scales generally show correlations greater than 0.6 throughout the middle latitudes of both hemispheres with some regions exceeding 0.8. The areas of highest correlations seem to be associated with the storm track regions of both the Northern and Southern Hemispheres. A detailed spectral analysis is performed for the medium scales on five pairs of time series of area averaged tropopause pressure and total ozone. In middle latitudes, total ozone and tropopause pressure exhibit generally similar distributions in the power spectrum. In the subtropics and tropics the power in ozone drops off more rapidly with increasing frequency than the power in tropopause pressure. Only in the Northern Hemisphere middle latitudes does one find a clear association between increased power in ozone and tropopause pressure and maxima in the coherency spectrum. Results for large scales are more complicated, showing generally positive correlations at middle latitudes

Designing an Optimal Ensemble Strategy for GMAO S2S Forecast System

The NASA Global Modeling and Assimilation Office (GMAO) Sub-seasonal to Seasonal (S2S) prediction system is being readied for a major upgrade. An important factor in successful extended range forecasting is the definition of the ensemble. Our overall strategy is to run a relatively large ensemble of about 40 members up to 3 months (focusing on the sub-seasonal forecast problem), after which we sub-sample the ensemble, and continue the forecast with about 10 members (up to 12 months). Here we present the results of our testing of various ways to generate the initial perturbations and the validation of a stratified sampling approach for choosing the members of the smaller ensemble. For the initialization of the ensemble we propose a combination of lagged and burst initial conditions. To generate perturbations for the burst ensemble members we used scaled differences of pairs of analysis states (chosen randomly from the corresponding season) separated by 1-10 days. We consider perturbing separately the atmosphere and the ocean, or both. By varying the separation times between the analysis states, we are able to produce perturbations that resemble well-known modes of variability. Focusing on the ENSO SST indices, we found that all types of perturbations are important for the ensemble spread with, however, considerable differences in the timing of the impacts on spread for the atmospheric and oceanic perturbations.Our initial (larger) ensemble size was determined so as to maximize the skill of predicting some of the leading modes of boreal winter atmospheric modes (namely the NAO, PNA and AO). Since it is not feasible for us to run with the larger ensemble beyond about 3 months, we employ a stratified sampling procedure that identifies the emerging directions of error growth to subset the ensemble. By comparing the results from the stratified ensemble with that of the randomly sampled ensemble of the same size, we find that the former provides substantially better estimates the mean of the original large ensemble

Designing an Optimal Ensemble Strategy for GMAO S2S Forecast System

GMAO Sub/Seasonal prediction system (S2S) is being readied for a major upgrade to GEOS-S2S Version 3. An important factor in successful extended range forecast is the definition of an ensemble For initialization of the ensemble we propose a combination of lagged and burst initial conditions. We plan to run a relatively large ensemble of 40 members for sub-seasonal forecast (up to 3 months), at which point we sub-sample the ensemble, and continue the forecast with 10 members (up to 12 months). Here we present the results of the extensive testing of various ways to generate the perturbations to the initial conditions and the validation of the stratified sampling strategy we chose.To generate perturbations for the burst ensemble members we used scaled differences of pairs of analysis states separated by 1-10 days, randomly chosen from a corresponding season. We considered perturbing separately only the atmospheric fields or only the ocean or both of the forecast initial conditions. Considering varying separation times between the analysis states, we were able to produce perturbations sampling various modes of variability. Focusing on the ENSO SST indices, we found that all types of perturbations are important for the ensemble spread.Our ensemble size for sub-seasonal forecasts was determined as to maximize the skill of predicting some of the leading modes of boreal winter atmospheric modes, NAO, PNA and AO. It is not feasible to run equally large ensemble for seasonal forecasts. Using a stratified sampling procedure we can identify the emerging directions of error growth. By comparing the stratified ensemble with randomly sampled ensemble of the same size, we were able to show that the former better estimates the mean of the original large ensemble

Phase Locking of the Boreal Summer Atmospheric Response to Dry Land Surface Anomalies in the Northern Hemisphere

Past modeling simulations, supported by observational composites, indicate that during boreal summer, dry soil moisture anomalies in very different locations within the United States continental interior tend to induce the same upper-tropospheric circulation pattern: a high anomaly forms over west-central North America and a low anomaly forms to the east. The present study investigates the causes of this apparent phase locking of the upper-level circulation response and extends the investigation to other land regions in the Northern Hemisphere. The phase locking over North America is found to be induced by zonal asymmetries in the local basic state originating from North American orography. Specifically, orography-induced zonal variations of air temperature, those in the lower troposphere in particular, and surface pressure play a dominant role in placing the soil moisture-forced negative Rossby wave source (dominated by upper-level divergence anomalies) over the eastern leeside of the Western Cordillera, which subsequently produces an upper-level high anomaly over west-central North America, with the downstream anomalous circulation responses phase-locked by continuity. The zonal variations of the local climatological atmospheric circulation, manifested as a climatological high over central North America, help shape the spatial pattern of the upper-level circulation responses. Considering the rest of the Northern Hemisphere, the northern Middle East exhibits similar phase locking, also induced by local orography. The Middle Eastern phase locking, however, is not as pronounced as that over North America; North America is where soil moisture anomalies have the greatest impact on the upper-tropospheric circulation

The Climate Signal in Regional Moisture Fluxes: A Comparison of Three Global Data Assimilation Products

This study assesses the quality of estimates of climate variability in moisture flux and convergence from three assimilated data sets: two are reanalysis products generated at the Goddard Data Assimilation Office (DAO) and the National Centers for Environmental Prediction/National Centers for Atmospheric Research (NCEPJNCAR), and the third consists of the operational analyses generated at the European Center for Medium Range Forecasts (ECMWF). The regions under study (the United States Great Plains, the Indian monsoon region, and Argentina east of the Andes) are characterized by frequent low level jets (LLJs) and other interannual low level wind variations tied to the large-scale flow. While the emphasis is on the reanalysis products, the comparison with the operational product is provided to help assess the improvements gained from a fixed analysis system. All three analyses capture the main moisture flux anomalies associated with selected extreme climate (drought and flood) events during the period 1985-93. The correspondence is strongest over the Great Plains and weakest over the Indian monsoon region reflecting differences in the observational coverage. For the reanalysis products, the uncertainties in the lower tropospheric winds is by far the dominant source of the discrepancies in the moisture flux anomalies in the middle latitude regions. Only in the Indian Monsoon region, where interannual variability in the low level winds is comparatively small, does the moisture bias play a substantial role. In contrast, the comparisons with the operational product show differences in moisture which are comparable torhe differences in the wind in all three regions. Compared with the fluxes, the anomalous moisture convergences show substantially larger differences among the three products. The best agreement occurs over the Great Plains region where all three products show vertically-integrated moisture convergence during the floods and divergence during the drought with differences in magnitude of about 25%. The reanalysis products, in particular, show good agreement in depicting the different roles of the mean flow and transients during the flood and drought periods. Differences between the three products in the other two regions exceed 100% reflecting differences in the low level jets and the large scale circulation patterns. The operational product tends to have locally larger amplitude convergence fields which average out in area-mean budgets: this appears to be at least in part due to errors in the surface pressure fields and aliasing from the higher resolution of the original ECMWF fields. On average, the reanalysis products show higher coherence with each other than with the operational product in the estimates of interannual variability. This result is less clear in the Indian monsoon region where differences in the input observations appears to be an important factor. The agreement in the anomalous convergence patterns is, however, still rather poor even over relatively data dense regions such as the United States Great Plains. These differences are attributed to deficiencies in the assimilating GCM's representations of the planetary boundary layer and orography, and a global observing system incapable of resolving the highly confined low level winds associated with the climate anomalies

Warm Season Subseasonal Variability and Climate Extremes in the Northern Hemisphere: The Role of Stationary Rossby Waves

This study examines the nature of boreal summer subseasonal atmospheric variability based on the new NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) for the period 1979-2010. An analysis of the June, July and August subseasonal 250hPa v-wind anomalies shows distinct Rossby wave-like structures that appear to be guided by the mean jets. On monthly subseasonal time scales, the leading waves (the first 10 rotated empirical orthogonal functions or REOFs of the 250hPa v-wind) explain about 50% of the Northern Hemisphere vwind variability, and account for more than 30% (60%) of the precipitation (surface temperature) variability over a number of regions of the northern middle and high latitudes, including the U.S. northern Great Plains, parts of Canada, Europe, and Russia. The first REOF in particular, consists of a Rossby wave that extends across northern Eurasia where it is a dominant contributor to monthly surface temperature and precipitation variability, and played an important role in the 2003 European and 2010 Russian heat waves. While primarily subseasonal in nature, the Rossby waves can at times have a substantial seasonal mean component. This is exemplified by REOF 4 which played a major role in the development of the most intense anomalies of the U.S. 1988 drought (during June) and the 1993 flooding (during July), though differed in the latter event by also making an important contribution to the seasonal mean anomalies. A stationary wave model (SWM) is used to reproduce some of the basic features of the observed waves and provide insight into the nature of the forcing. In particular, the responses to a set of idealized forcing functions are used to map the optimal forcing patterns of the leading waves. Also, experiments to reproduce the observed waves with the SWM using MERRA-based estimates of the forcing indicate that the wave forcing is dominated by sub-monthly vorticity transients

A graphics package for meteorological data, version 1.5

A plotting package has been developed to simplify the task of plotting meteorological data. The calling sequences and examples of high level yet flexible routines which allow contouring, vectors and shading of cylindrical, polar, orthographic and Mollweide (egg) projections are given. Routines are also included for contouring pressure-latitude and pressure-longitude fields with linear or log scales in pressure (interpolation to fixed grid interval is done automatically). Also included is a fairly general line plotting routine. The present version (1.5) produces plots on WMS laser printers and uses graphics primitives from WOLFPLOT