The impact of internal atmospheric variability on the North Pacific SST variability

The impact of internal atmospheric variability on North Pacific sea surface temperature (SST) variability is examined based on three coupled general circulation model simulations. The three simulations differ only in the level of atmospheric noise occuring over the ocean at the air-sea interface. The amplitude of atmospheric noise is controlled by use of the interactive ensemble technique. This technique simultaneously couples multiple realizations of a single atmospheric model to a single realization of an ocean model. The atmospheric component models all experience the same SST, but the ocean component is forced by the ensemble averaged fluxes thereby reducing the impact of internal atmospheric dynamics at the air-sea interface. The ensemble averaging is only applied at the air-sea interface so that the internal atmospheric dynamics (i.e., transients) of each atmospheric ensemble member is unaffected. This interactive ensemble technique significantly reduces the SST variance throughout the North Pacific. The reduction in SST variance is proportional to the number of ensemble members indicating that most of the variability can simply be explained as the response to atmospheric stochastic forcing. In addition, the impact of the internal atmospheric dynamics at the air-sea interface masks out much of the tropical-midlatitude SST teleconnections on interannual time scales. Once this interference is reduced (i.e., applying the interactive ensemble technique), tropical-midlatitude SST teleconnections are easily detected.     

Uncertainty in climate change projections: the role of internal variability

Uncertainty in future climate change presents a key challenge for adaptation planning. In this study, uncertainty arising from internal climate variability is investigated using a new 40-member ensemble conducted with the National Center for Atmospheric Research Community Climate System Model Version 3 (CCSM3) under the SRES A1B greenhouse gas and ozone recovery forcing scenarios during 2000–2060. The contribution of intrinsic atmospheric variability to the total uncertainty is further examined using a 10,000-year control integration of the atmospheric model component of CCSM3 under fixed boundary conditions. The global climate response is characterized in terms of air temperature, precipitation, and sea level pressure during winter and summer. The dominant source of uncertainty in the simulated climate response at middle and high latitudes is internal atmospheric variability associated with the annular modes of circulation variability. Coupled ocean-atmosphere variability plays a dominant role in the tropics, with attendant effects at higher latitudes via atmospheric teleconnections. Uncertainties in the forced response are generally larger for sea level pressure than precipitation, and smallest for air temperature. Accordingly, forced changes in air temperature can be detected earlier and with fewer ensemble members than those in atmospheric circulation and precipitation. Implications of the results for detection and attribution of observed climate change and for multi-model climate assessments are discussed. Internal variability is estimated to account for at least half of the inter-model spread in projected climate trends during 2005–2060 in the CMIP3 multi-model ensemble.     

Influence of SST biases on future climate change projections

We use a quantile-based bias correction technique and a multi-member ensemble of the atmospheric component of NCAR CCSM3 (CAM3) simulations to investigate the influence of sea surface temperature (SST) biases on future climate change projections. The simulations, which cover 1977?C1999 in the historical period and 2077?C2099 in the future (A1B) period, use the CCSM3-generated SSTs as prescribed boundary conditions. Bias correction is applied to the monthly time-series of SSTs so that the simulated changes in SST mean and variability are preserved. Our comparison of CAM3 simulations with and without SST correction shows that the SST biases affect the precipitation distribution in CAM3 over many regions by introducing errors in atmospheric moisture content and upper-level (lower-level) divergence (convergence). Also, bias correction leads to significantly different precipitation and surface temperature changes over many oceanic and terrestrial regions (predominantly in the tropics) in response to the future anthropogenic increases in greenhouse forcing. The differences in the precipitation response from SST bias correction occur both in the mean and the percent change, and are independent of the ocean?Catmosphere coupling. Many of these differences are comparable to or larger than the spread of future precipitation changes across the CMIP3 ensemble. Such biases can affect the simulated terrestrial feedbacks and thermohaline circulations in coupled climate model integrations through changes in the hydrological cycle and ocean salinity. Moreover, biases in CCSM3-generated SSTs are generally similar to the biases in CMIP3 ensemble mean SSTs, suggesting that other GCMs may display a similar sensitivity of projected climate change to SST errors. These results help to quantify the influence of climate model biases on the simulated climate change, and therefore should inform the effort to further develop approaches for reliable climate change projection.     

海冰模式CICE4.0与LASG/IAP气候系统模式的耦合试验

利用美国Los Alamos国家实验室发展的最新海冰模式(CICE4.0)替代了LASG/IAP气候系统模式(FGOALS_g1.1)中的海冰模式(CSIM4), 形成新的耦合模式。在此基础上, 利用新的耦合模式对20世纪中后期的全球气候进行了模拟, 来检验CICE4.0对耦合模式中海冰和海洋模拟结果的改进。结果表明CICE4.0对于FGOALS_g1.1的极地气候模拟有一定改进作用, 主要表现在:(1) 南北极海冰边缘碎冰区显著减少; (2) 南大洋海表温度和海冰的模拟明显改善, 分布特征与观测非常吻合。但是新耦合模式也存在如下不足: (1) 北大西洋海冰相对偏多, 北大西洋经圈翻转环流大大减弱, 这主要是由于北大西洋海表面温度的冷误差造成的; (2) 南北极大气环流场的模拟无明显改善。此外, 本文还比较了采用不同短波辐射方案对于耦合模拟结果的影响, 结果表明, 相对于CCSM3短波辐射方案, Delta-Eddington方案模拟的海表面温度偏冷, 海冰厚度偏厚, 北大西洋经圈翻转环流略有偏弱。     

Naturally forced multidecadal variability of the Atlantic meridional overturning circulation

The mechanisms by which natural forcing factors alone could drive simulated multidecadal variability in the Atlantic meridional overturning circulation (AMOC) are assessed in an ensemble of climate model simulations. It is shown for a new state-of-the-art general circulation model, HadGEM2-ES, that the most important of these natural forcings, in terms of the multidecadal response of the AMOC, is solar rather than volcanic forcing. AMOC strengthening occurs through a densification of the North Atlantic, driven by anomalous surface freshwater fluxes due to increased evaporation. These are related to persistent North Atlantic atmospheric circulation anomalies, driven by forced changes in the stratosphere, associated with anomalously weak solar irradiance during the late nineteenth and early twentieth centuries. Within a period of approximately 100 years the 11-year smoothed ensemble mean AMOC strengthens by 1.5 Sv and subsequently weakens by 1.9 Sv, representing respectively approximately 3 and 4 standard deviations of the 11-year smoothed control simulation. The solar-induced variability of the AMOC has various relevant climate impacts, such as a northward shift of the intertropical convergence zone, anomalous Amazonian rainfall, and a sustained increase in European temperatures. While this model has only a partial representation of the atmospheric response to solar variability, these results demonstrate the potential for solar variability to have a multidecadal impact on North Atlantic climate.     

Twentieth century North Atlantic climate change. Part II: Understanding the effect of Indian Ocean warming

Ensembles of atmospheric general circulation model (AGCM) experiments are used in an effort to understand the boreal winter Northern Hemisphere (NH) extratropical climate response to the observed warming of tropical sea surface temperatures (SSTs) over the last half of the twentieth Century. Specifically, we inquire about the origins of unusual, if not unprecedented, changes in the wintertime North Atlantic and European climate that are well described by a linear trend in most indices of the North Atlantic Oscillation (NAO). The simulated NH atmospheric response to the linear trend component of tropic-wide SST change since 1950 projects strongly onto the positive polarity of the NAO and is a hemispheric pattern distinguished by decreased (increased) Arctic (middle latitude) sea level pressure. Progressive warming of the Indian Ocean is the principal contributor to this wintertime extratropical response, as shown through additional AGCM ensembles forced with only the SST trend in that sector. The Indian Ocean influence is further established through the reproducibility of results across three different models forced with identical, idealized patterns of the observed warming. Examination of the transient atmospheric adjustment to a sudden “switch-on” of an Indian Ocean SST anomaly reveals that the North Atlantic response is not consistent with linear theory and most likely involves synoptic eddy feedbacks associated with changes in the North Atlantic storm track. The tropical SST control exerted over twentieth century regional climate underlies the importance of determining the future course of tropical SST for regional climate change and its uncertainty. Better understanding of the extratropical responses to different, plausible trajectories of the tropical oceans is key to such efforts.     

Influence of local and remote SST on North Atlantic tropical cyclone potential intensity

We examine the role of local and remote sea surface temperature (SST) on the tropical cyclone potential intensity in the North Atlantic using a suite of model simulations, while separating the impact of anthropogenic (external) forcing and the internal influence of Atlantic Multidecadal Variability. To enable the separation by SST region of influence we use an ensemble of global atmospheric climate model simulations forced with historical, 1856–2006 full global SSTs, and compare the results to two other simulations with historical SSTs confined to the tropical Atlantic and to the tropical Indian Ocean and Pacific. The effects of anthropogenic plus other external forcing and that of internal variability are separated by using a linear, “signal-to-noise” maximizing EOF analysis and by projecting the three model ensemble outputs onto the respective external forcing and internal variability time series. Consistent with previous results indicating a tampering influence of global tropical warming on the Atlantic hurricane potential intensity, our results show that non-local SST tends to reduce potential intensity associated with locally forced warming through changing the upper level atmospheric temperatures. Our results further indicate that the late twentieth Century increase in North Atlantic potential intensity, may not have been dominated by anthropogenic influence but rather by internal variability.     

Sensitivity of the Northern Hemisphere circulation to North Atlantic SSTs in the ARPÈGE Climate AGCM

In order to determine the response of the atmosphere to winter sea surface temperature (SST) anomalies in the North Atlantic area, we carried out ensemble runs of 20 years, forced with constant, perturbed, SST patterns using the climate version of the ARPèGE AGCM, at T42 resolution. A Monte Carlo technique was applied, in such a way that the control experiment, forced with observed climatological temperatures, and the four scenario experiments, forced with perturbed SSTs are equivalent to a length of 20 independent winters. Four anomalous winter North Atlantic sea surface temperature (SST) fields have been constructed by considering the observed SST variability in the main basins, namely the Labrador Sea and the Greenland Sea. Two patterns are of the `seesaw' type, while the two others have same the polarity in both basins. The patterns have been reinforced by a factor of 5–6 compared to presently observed multi-annual anomalies, in order to get SST anomalies which may have occurred during periods of the Little Ice Age. The differences between each of the four winter simulations with perturbed SSTs and the control run are analyzed in terms of tropospheric thickness, mean-sea-level pressure and storm activity. The `seesaw' type patterns give a weaker response in the tropospheric thickness fields than the two others. This is expected from simple considerations. In the mean circulation and synoptic activity, it appears that the Labrador Sea SST is important in determining the atmospheric response. This is probably due to enhanced temperature gradients east of New Foundland which enhances the storm activity. Received: 23 September 1998 / Accepted: 17 November 1999     

Simulated impact of vegetation on climate across the North American monsoon region in CCSM3.5

The influence of prescribed changes in vegetation on the climate of the North American monsoon region is examined using the National Center for Atmospheric Research Community Climate System Model Version 3.5 (NCAR CCSM3.5). Initial value ensemble experiments are performed in which the vegetation cover fraction over the North American monsoon region is reduced by 0.2 and the intra-annual climatic response is assessed probabilistically in each one-year ensemble experiment. Changes in the surface radiation budget include decreases in sensible and latent heat fluxes and increases in upward longwave and downward shortwave radiation fluxes, with small net changes in surface albedo. The climatic responses to reduced vegetation cover fraction include year-round increases in ground and surface air temperature, a dampened hydrologic cycle with decreased springtime evaporation, springtime and autumnal precipitation, and autumnal cloud cover, and enhanced atmospheric subsidence in late autumn. Decreased vegetation shifts the monsoon season over the Southwest United States earlier in the year. Within the North American monsoon region, the most robust vegetation feedbacks to climate are found over woody landscapes.     

The MJO and global warming: a study in CCSM4

The change in Madden–Julian oscillation (MJO) amplitude and variance in response to anthropogenic climate change is assessed in the 1° nominal resolution community climate system model, version 4 (CCSM4), which has a reasonable representation of the MJO characteristics both dynamically and statistically. The twentieth century CCSM4 run is compared with the warmest twenty-first century projection (representative concentration pathway 8.5, or RCP8.5). The last 20 years of each simulation are compared in their MJO characteristics, including spatial variance distributions of winds, precipitation and outgoing longwave radiation, histograms of event amplitude, phase and duration, and composite maps of phases. The RCP8.5 run exhibits increased variance in intraseasonal precipitation, larger-amplitude MJO events, stronger MJO rainfall in the central and eastern tropical Pacific, and a greater frequency of MJO occurrence for phases corresponding to enhanced rainfall in the Indian Ocean sector. These features are consistent with the concept of an increased magnitude for the hydrological cycle under greenhouse warming conditions. Conversely, the number of active MJO days decreases and fewer weak MJO events occur in the future climate state. These results motivate further study of these changes since tropical rainfall variability plays such an important role in the region’s socio-economic well being.     

Interannual variations and future change of wintertime extratropical cyclone activity over North America in CCSM3

Climatology and interannual variations of wintertime extratropical cyclone frequency in CCSM3 twentieth century simulation are compared with the NCEP/NCAR reanalysis during 1950–1999. CCSM3 can simulate the storm tracks reasonably well, although the model produces slightly less cyclones at the beginning of the Pacific and Atlantic storm tracks and weaker poleward deflection over the Pacific. As in the reanalysis, frequency of cyclones stronger than 980 hPa shows significant correlation with the Pacific/North America (PNA) teleconnection pattern over the Pacific region and with the North Atlantic Oscillation (NAO) in the Atlantic sector. Composite maps are constructed for opposite phases of El Nino-Southern Oscillation (ENSO) and the NAO and all anomalous patterns coincide with observed. One CCSM3 twenty-first century A1B scenario realization indicates there is significant increase in the extratropical cyclone frequency on the US west coast and decrease in Alaska. Meanwhile, cyclone frequency increases from the Great Lakes region to Quebec and decreases over the US east coast, suggesting a possible northward shift of the Atlantic storm tracks under the warmer climate. The cyclone frequency anomalies are closely linked to changes in seasonal mean states of the upper-troposphere zonal wind and baroclinicity in the lower troposphere. Due to lack of 6-hourly outputs, we cannot apply the cyclone-tracking algorithm to the other eight CCSM3 realizations. Based on the linkage between the mean state change and the cyclone frequency anomalies, it is likely a common feature among the other ensemble members that cyclone activity is reduced on the East Coast and in Alaska as a result of global warming.     

Global SST influence on twentieth century NAO variability

Recent studies have suggested that sea surface temperature (SST) is an important source of variability of the North Atlantic Oscillation (NAO). Here, we deal with four basic aspects contributing to this issue: (1) we investigate the characteristic time scales of this oceanic influence; (2) quantify the scale-dependent hindcast potential of the NAO during the twentieth century as derived from SST-driven atmospheric general circulation model (AGCM) ensembles; (3) the relevant oceanic regions are identified, corresponding SST indices are defined and their relationship to the NAO are evaluated by means of cross spectral analysis and (4) our results are compared with long-term coupled control experiments with different ocean models in order to ensure whether the spectral relationship between the SST regions and the NAO is an intrinsic mode of the coupled climate system, involving the deep ocean circulation, rather than an artefact of the unilateral SST forcing. The observed year-to-year NAO fluctuations are barely influenced by the SST. On the decadal time scales the major swings of the observed NAO are well reproduced by various ensembles from the middle of the twentieth century onward, including the negative state in the 1960s and part of the positive trend afterwards. A six-member ECHAM4-T42 ensemble reveals that the SST boundary condition affects 25% of total decadal-mean and interdecadal-trend NAO variability throughout the twentieth century. The most coherent NAO-related SST feature is the well-known North Atlantic tripole. Additional contributions may arise from the southern Pacific and the low-latitude Indian Ocean. The coupled climate model control runs suggest only the North Atlantic SST-NAO relationship as being a true characteristic of the coupled climate system. The coherence and phase spectra of observations and coupled simulations are in excellent agreement, confirming the robustness of this decadal-scale North Atlantic air–sea coupled mode.     

Extension and Intensification of the Meso-American mid-summer drought in the twenty-first century

Recent global-scale analyses of the CMIP3 model projections for the twenty-first century indicate a strong, coherent decreased precipitation response over Central America and the Intra-America Seas region. We explore this regional response and examine the models’ skill in representing present-day climate over this region. For much of Central America, the annual cycle of precipitation is characterized by a rainy season that extends from May to October with a period of reduced precipitation in July and August called the mid-summer drought. A comparison of the climate of the twentieth century simulations (20c3m) with observations over the period 1961–1990 shows that nearly all models underestimate precipitation over Central America, due in part to an underestimation of sea surface temperatures over the tropical North Atlantic and an excessively smooth representation of regional topographical features. However, many of the models capture the mid-summer drought. Differences between the A1B scenario (2061–2090) and 20c3m (1961–1990) simulations show decreased precipitation in the future climate scenario, mostly in June and July, just before and during the onset of the mid-summer drought. We thus hypothesize that the simulated twenty-first century drying over Central America represents an early onset and intensification of the mid-summer drought. An analysis of circulation changes indicates that the westward expansion and intensification of the North Atlantic subtropical high associated with the mid-summer drought occurs earlier in the A1B simulations, along with stronger low-level easterlies. The eastern Pacific inter-tropical convergence zone is also located further southward in the scenario simulations. There are some indications that these changes could be forced by ENSO-like warming of the tropical eastern Pacific and increased land–ocean heating contrasts over the North American continent.     

An overview of decadal climate predictability in a multi-model ensemble by climate model MIROC

Decadal climate predictability is examined in hindcast experiments by a multi-model ensemble using three versions of the coupled atmosphere-ocean model MIROC. In these hindcast experiments, initial conditions are obtained from an anomaly assimilation procedure using the observed oceanic temperature and salinity with prescribed natural and anthropogenic forcings on the basis of the historical data and future emission scenarios in the Intergovernmental Panel of Climate Change. Results of the multi-model ensemble in our hindcast experiments show that predictability of surface air temperature (SAT) anomalies on decadal timescales mostly originates from externally forced variability. Although the predictable component of internally generated variability has considerably smaller SAT variance than that of externally forced variability, ocean subsurface temperature variability has predictive skills over almost a decade, particularly in the North Pacific and the North Atlantic where dominant signals associated with Pacific decadal oscillation (PDO) and the Atlantic multidecadal oscillation (AMO) are observed. Initialization enhances the predictive skills of AMO and PDO indices and slightly improves those of global mean temperature anomalies. Improvement of these predictive skills in the multi-model ensemble is higher than that in a single-model ensemble.     

Extreme midlatitude cyclones and their implications for precipitation and wind speed extremes in simulations of the Maunder Minimum versus present day conditions

Extreme midlatitude cyclone characteristics, precipitation, wind speed events, their inter-relationships, and the connection to large-scale atmospheric patterns are investigated in simulations of a prolonged cold period, known as the Maunder Minimum from 1640 to 1715 and compared with today. An ensemble of six simulations for the Maunder Minimum as well as a control simulation for perpetual 1990 conditions are carried out with a coupled atmosphere-ocean general circulation model, i.e., the Climate Community System Model (CCSM). The comparison of the simulations shows that in a climate state colder than today the occurrence of cyclones, the extreme events of precipitation and wind speed shift southward in all seasons in the North Atlantic and the North Pacific. The extremes of cyclone intensity increases significantly in winter in almost all regions, which is related to a stronger meridional temperature gradient and an increase in lower tropospheric baroclinicity. Extremes of cyclone intensity in subregions of the North Atlantic are related to extremes in precipitation and in wind speed during winter. Moreover, extremes of cyclone intensity are also connected to distinct large-scale atmospheric patterns for the different subregions, but these relationships vanish during summer. Analyzing the mean 1,000 hPa geopotential height change of the Maunder Minimum simulations compared with the control simulation, we find a similar pattern as the correlation pattern with the cyclone intensity index of the southern Europe cyclones. This illustrates that changes in the atmospheric high-frequency, i.e., the simulated southward shift of cyclones in the North Atlantic and the related increase of extreme precipitation and wind speed in particular in the Mediterranean in winter, are associated with large-scale atmospheric circulation changes.     

Non-stationarity of the signal and noise characteristics of seasonal precipitation anomalies

In order to improve seasonal-to-interannual precipitation forecasts and their application by decision makers, there is a clear need to understand when, where, and to what extent seasonal precipitation anomalies are driven by potentially predictable surface–atmosphere interactions versus to chaotic interannual atmospheric dynamics. Using a simple Monte Carlo approach, interannual variability and linear trends in the SST-forced signal and potential predictability of boreal winter precipitation anomalies is examined in an ensemble of twentieth century AGCM simulations. Signal and potential predictability are shown to be non-stationary over more than 80% of the globe, while chaotic noise is shown to be stationary over most of the globe. Correlation analysis with respect to magnitudes of the four leading modes of global SST variability suggests that interannual variability and trends in signal and potential predictability over 35% of the globe is associated with ENSO-related SST variability; signal and potential predictability are not significantly associated with SST modes characterized by a global SST trend, North Atlantic SST variability, and North Pacific SST variability, respectively. Results suggest that mechanisms other than SST variability contribute to the non-stationarity of signal and noise characteristics of hydroclimatic variability over mid- and high-latitude regions.     

Regional water budgets and hydroclimatic trend variations in Xinjiang from 1951 to 2000

Xinjiang is located in arid northwestern China where water cycle has accelerated due to increased precipitation and temperature. However, the regional water budget characteristics vary due to the complex topography and spatial heterogeneities of hydroclimatology. This study uses atmospheric forcing constrained by observation from 90 meteorological stations in Xinjiang as input for the optimized Community Land Model version 3.5 (CLM 3.5) to investigate Xinjiang’s regional water budgets from 1951 to 2000 between the northern and southern part divided by the Tianshan Mountain. Results show that precipitation, evapotranspiration and runoff increased in Xinjiang from 1951 to 2000, particularly after the climate shift around 1987, and the net water flux (P-E) gap between North and South Xinjiang was widened. Rapid, intense wetting occurred in North Xinjiang in response to regional climate change, whereas South Xinjiang experienced relatively small changes. North and South Xinjiang exhibited opposite trends in water table depth (WTD), which became shallower in North Xinjiang, particularly after 1987. The WTD in South Xinjiang gradually became deeper. These results suggest that water resources in North Xinjiang are more sensitive to the warmer and wetter climate than South Xinjiang, and the serious water shortage in South Xinjiang did not improve during the second half of the twentieth century.     

Influence of the Atlantic Multidecadal Oscillation on the winter climate of East China

The Atlantic Multidecadal Oscillation （AMO）, the multidecadal variation of North Atlantic sea surface temperature （SST）, exhibits an oscillation with a period of 65-80 years and an amplitude of 0.4℃. Observational composite analyses reveal that the warm phase AMO is linked to warmer winters in East China, with enhanced precipitation in the north of this region and reduced precipitation in the south, on multidecadal time scales. The pattern is reversed during the cold phase AMO. Whether the AMO acts as a forcing of the multidecadal winter climate of East China is explored by investigating the atmospheric response to warm AMO SST anomalies in a large ensemble of atmospheric general circulation model （AGCM） experiments. The results from three AGCMs are consistent and suggest that the AMO warmth favors warmer winters in East China. This influence is realized through inducing negative surface air pressure anomalies in the hemispheric-wide domain extending from the midlatitude North Atlantic to midlatitude Eurasia. These negative surface anomalies favor the weakening of the Mongolian Cold High, and thus induce a weaker East Asian Winter Monsoon.     

Twentieth century Walker Circulation change: data analysis and model experiments

Recent studies indicate a weakening of the Walker Circulation during the twentieth century. Here, we present evidence from an atmospheric general circulation model (AGCM) forced by the history of observed sea surface temperature (SST) that the Walker Circulation may have intensified rather than weakened. Observed Equatorial Indo-Pacific Sector SST since 1870 exhibited a zonally asymmetric evolution: While the eastern part of the Equatorial Pacific showed only a weak warming, or even cooling in one SST dataset, the western part and the Equatorial Indian Ocean exhibited a rather strong warming. This has resulted in an increase of the SST gradient between the Maritime Continent and the eastern part of the Equatorial Pacific, one driving force of the Walker Circulation. The ensemble experiments with the AGCM, with and without time-varying external forcing, suggest that the enhancement of the SST gradient drove an anomalous atmospheric circulation, with an enhancement of both Walker and Hadley Circulation. Anomalously strong precipitation is simulated over the Indian Ocean and anomalously weak precipitation over the western Pacific, with corresponding changes in the surface wind pattern. Some sensitivity to the forcing SST, however, is noticed. The analysis of twentieth century integrations with global climate models driven with observed radiative forcing obtained from the Coupled Model Intercomparison Project (CMIP) database support the link between the SST gradient and Walker Circulation strength. Furthermore, control integrations with the CMIP models indicate the existence of strong internal variability on centennial timescales. The results suggest that a radiatively forced signal in the Walker Circulation during the twentieth century may have been too weak to be detectable.     

Forced and unforced variability of twentieth century North American droughts and pluvials

Research on the forcing of drought and pluvial events over North America is dominated by general circulation model experiments that often have operational limitations (e.g., computational expense, ability to simulate relevant processes, etc). We use a statistically based modeling approach to investigate sea surface temperature (SST) forcing of the twentieth century pluvial (1905?C1917) and drought (1932?C1939, 1948?C1957, 1998?C2002) events. A principal component (PC) analysis of Palmer Drought Severity Index (PDSI) from the North American Drought Atlas separates the drought variability into five leading modes accounting for 62% of the underlying variance. Over the full period spanning these events (1900?C2005), the first three PCs significantly correlate with SSTs in the equatorial Pacific (PC 1), North Pacific (PC 2), and North Atlantic (PC 3), with spatial patterns (as defined by the empirical orthogonal functions) consistent with our understanding of North American drought responses to SST forcing. We use a large ensemble statistical modeling approach to determine how successfully we can reproduce these drought/pluvial events using these three modes of variability. Using Pacific forcing only (PCs 1?C2), we are able to reproduce the 1948?C1957 drought and 1905?C1917 pluvial above a 95% random noise threshold in over 90% of the ensemble members; the addition of Atlantic forcing (PCs 1?C2?C3) provides only marginal improvement. For the 1998?C2002 drought, Pacific forcing reproduces the drought above noise in over 65% of the ensemble members, with the addition of Atlantic forcing increasing the number passing to over 80%. The severity of the drought, however, is underestimated in the ensemble median, suggesting this drought intensity can only be achieved through internal variability or other processes. Pacific only forcing does a poor job of reproducing the 1932?C1939 drought pattern in the ensemble median, and less than one third of ensemble members exceed the noise threshold (28%). Inclusion of Atlantic forcing improves the ensemble median drought pattern and nearly doubles the number of ensemble members passing the noise threshold (52%). Even with the inclusion of Atlantic forcing, the intensity of the simulated 1932?C1939 drought is muted, and the drought itself extends too far into the southwest and southern Great Plains. To an even greater extent than the 1998?C2002 drought, these results suggest much of the variance in the 1932?C1939 drought is dependent on processes other than SST forcing. This study highlights the importance of internal noise and non SST processes for hydroclimatic variability over North America, complementing existing research using general circulation models.