ALBERT

Description: The current Florida State University (FSU) climate model is upgraded by coupling the National Center for Atmospheric Research (NCAR) Community Land Model Version 2 (CLM2) as its land component in order to make a better simulation of surface air temperature and precipitation on the seasonal time scale, which is important for crop model application. Climatological and seasonal simulations with the FSU climate model coupled to the CLM2 (hereafter FSUCLM) are compared to those of the control (the FSU model with the original simple land surface treatment). The current version of the FSU model is known to have a cold bias in the temperature field and a wet bias in precipitation. The implementation of FSUCLM has reduced or eliminated this bias due to reduced latent heat flux and increased sensible heat flux. The role of the land model in seasonal simulations is shown to be more important during summertime than wintertime. An additional experiment that assimilates atmospheric forcings produces improved land-model initial conditions, which in turn reduces the biases further. The impact of various deep convective parameterizations is examined as well to further assess model performance. The land scheme plays a more important role than the convective scheme in simulations of surface air temperature. However, each convective scheme shows its own advantage over different geophysical locations in precipitation simulations.

Description: A number of recent papers have argued that the mechanical energy budget of the ocean places constraints on how the thermohaline circulation is driven. These papers have been used to argue that climate models, which do not specifically account for the energy of mixing, potentially miss a very important feedback on climate change. This paper reexamines the question of what energetic arguments can teach us about the climate system and concludes that the relationship between energetics and climate is not straightforward. By analyzing the buoyancy transport equation, it is demonstrated that the large-scale transport of heat within the ocean requires an energy source of around 0.2 TW to accomplish vertical transport and around 0.4 TW (resulting from cabbeling) to accomplish horizontal transport. Within two general circulation models, this energy is almost entirely supplied by surface winds. It is also shown that there is no necessary relationship between heat transport and mechanical energy supply.

Description: A simple linearized transport model of anomalous Southern Ocean sea surface temperature (SST) is studied to determine whether it can sustain anomalies of realistic amplitudes under a physically based stochastic forcing. As noted in previous studies, eigenmodes of this system with zonal wavenumbers 2 and 3 share key propagation characteristics with the SST anomalies associated with the Antarctic Circumpolar Wave (ACW). The system is solved on a grid that follows the path of the Antarctic Circumpolar Current (ACC) and is forced by a stochastic heat flux. The forcing is white in space and time and represents the advection of the mean SST gradient by high-frequency variations in the cross-ACC velocity, due to mesoscale eddy variability. The magnitude of the stochastic forcing is determined from a global eddy-permitting ocean model. Anomalous ocean surface velocity variability (8 cm s−1) coupled to a mean cross-ACC SST gradient of 0.8°C (°latitude)−1 sustains anomalous interannual SST variability at low wavenumbers and amplitudes of the order of 1°C, consistent with those associated with the ACW. In the long-term mean, variance is broadly spread among low wavenumbers, in contrast to the dominance of one or two zonal wavenumbers in the ACW observations. It is found, however, that the model produces single dominant wavenumbers over individual periods of decades, suggesting that the apparent unimodal nature of the ACW may be an artifact of the short observational record used to infer it. Alternatively, it is shown that a nonisotropic forcing may also result in a stronger preference for particular zonal wavenumbers. It is shown that if the atmosphere at mid to high southern latitudes has an equivalent barotropic response to heating, then the resulting sea level pressure anomalies reproduce the phase relationship of the observed ACW. These results are consistent with the notion that a simple stochastically forced advection of SST anomalies can explain SST variability associated with the ACW to leading order.

Description: The characteristics and trends of observed river discharge into the Hudson, James, and Ungava Bays (HJUBs) for the period 1964–2000 are investigated. Forty-two rivers with outlets into these bays contribute on average 714 km3 yr−1 [= 0.023 Sv (1 Sv ≡ 106 m3 s−1)] of freshwater to high-latitude oceans. For the system as a whole, discharge attains an annual peak of 4.2 km3 day−1 on average in mid-June, whereas the minimum of 0.68 km3 day−1 occurs on average during the last week of March. The Nelson River contributes as much as 34% of the daily discharge for the entire system during winter but diminishes in relative importance during spring and summer. Runoff rates per contributing area are highest (lowest) on the eastern (western) shores of the Hudson and James Bays. Linear trend analyses reveal decreasing discharge over the 37-yr period in 36 out of the 42 rivers. By 2000, the total annual freshwater discharge into HJUBs diminished by 96 km3 (−13%) from its value in 1964, equivalent to a reduction of 0.003 Sv. The annual peak discharge rate associated with snowmelt has advanced by 8 days between 1964 and 2000 and has diminished by 0.036 km3 day−1 in intensity. There is a direct correlation between the timing of peak spring discharge rates and the latitude of a river’s mouth; the spring freshet varies by 5 days for each degree of latitude. Continental snowmelt induces a seasonal pulse of freshwater from HJUBs that is tracked along its path into the Labrador Current. It is suggested that the annual upper-ocean salinity minimum observed on the inner Newfoundland Shelf can be explained by freshwater pulses composed of meltwater from three successive winter seasons in the river basins draining into HJUBs. A gradual salinization of the upper ocean during summer over the period 1966–94 on the inner Newfoundland Shelf is in accord with a decadal trend of a diminishing intensity in the continental meltwater pulses.

Description: In this study the fifth-generation Pennsylvania State University–National Center for Atmospheric Research (PSU–NCAR) Mesoscale Model (MM5) version 3.5.2 was used to simulate the Indian summer monsoon during the two contrasting years of 1987 and 1988, a dry year and a wet year, respectively. Three different convection parameterization schemes of Betts–Miller–Janjic, Kain–Fritsch, and Grell were used to study the sensitivity of monsoon to cumulus effects. The model was integrated for a period of 6 months, starting from three different initial conditions of 0000 UTC on 1, 2, and 3 May of each year using the NCEP–NCAR reanalysis data as input. The 6-hourly reanalysis data were used to provide the lateral boundary conditions, and the observed weekly Reynolds sea surface temperature, linearly interpolated to 6 h, was used as the lower boundary forcing. The results show that all three cumulus schemes were able to simulate the interannual and intraseasonal variabilities in the monsoon with reasonable accuracy. However, the spatial distribution of the rainfall and its quantity were different in all the schemes. The Grell scheme underestimated the rainfall in both the years. The Kain–Fritsch scheme simulated the observed rainfall well during July and August, the peak monsoon months, of the year 1988 but overestimated the rainfall in June and September of 1988 and throughout the monsoon season of 1987. The Betts–Miller–Janjic scheme simulated less rainfall in the drought year of 1987 and overestimated the rainfall in June and July of 1988. The circulation patterns simulated by the Betts–Miller–Janjic and Kain–Fritsch schemes are comparable to the observed patterns.

Description: The canonical correlation analysis (CCA) and singular value decomposition (SVD) approaches for estimating a time series from a time-dependent vector and vice versa are investigated, and their relationship to multiple linear regression (MLR) and to regression maps is discussed. Earlier findings are reviewed and combined with new aspects to provide a systematic overview. It is shown that regression maps are proportional to canonical patterns and to singular vectors and that the estimate of a time-dependent vector from a time series does not depend on whether CCA, SVD, or component-wise regressions are used. When a time series is linearly estimated from a time-dependent vector, it is known that CCA is equivalent to MLR. It is demonstrated that an estimate for the time series based on a time expansion coefficient of the regression map that is calculated by orthogonal projection is identical to an SVD estimate, but different from the CCA and MLR estimate. The two approaches also lead to different correlations between the time series and the time expansion coefficient of its signal. The CCA–MLR and the SVD–regression map approaches are compared in an example where the January Arctic Oscillation index for the period 1948–2002 was estimated from extratropical Northern Hemispheric 850-hPa temperature. For CCA–MLR the leading principal components (PCs) of the temperature field were used as predictors, while for SVD the full field was employed. For more than seven retained PCs the skill in terms of correlations and mean squared error based on cross validation was for both approaches practically identical, but CCA–MLR showed a higher bias. For a smaller number of predictor PCs the SVD–regression map approach performed better. The discrepancy between the skill on the fitting data and on the independent data used for validation was in this example larger for the CCA–MLR approach.

Description: The implications are investigated of representing ocean gyre circulations by a diffusion term in the Stommel and Rooth box models of the thermohaline circulation (THC) in one and two hemispheres, respectively. The approach includes mostly analytical solution and study of the bifurcation structure, but also numerical integration and feedback analysis. Sufficient diffusion (gyre strength) eliminates multiple equilibria from either model, highlighting the need for accurate gyre circulation strength in general circulation models (GCMs) when considering the potential for abrupt climate change associated with THC shutdown. With diffusion, steady-state flow strength in the Rooth model depends on freshwater forcing (i.e., implied atmospheric water vapor transport) in both hemispheres, not only on that in the upwelling hemisphere, as in the nondiffusive case. With asymmetric freshwater forcing, two solutions (strong stable and weak unstable) are found with sinking in the hemisphere with stronger forcing and one solution with sinking in the other hemisphere. Under increased freshwater forcing the two solutions in the hemisphere with stronger forcing meet in a saddle-node bifurcation (if diffusion is sufficiently strong to prevent a subcritical Hopf bifurcation first), followed by flow reversal. Thus, the bifurcation structure with respect to freshwater forcing of the diffusive Rooth model of two-hemisphere THC is similar to that of the Stommel model of single-hemisphere THC, albeit with a very different dynamical interpretation. Gyre circulations stabilize high-latitude sinking in the Stommel model. In the Rooth model, gyre circulations only stabilize high-latitude sinking if the freshwater forcing is weaker in the sinking hemisphere than in the upwelling hemisphere, by an amount that increases with diffusion. The values of diffusion and freshwater forcing at which qualitative change in behavior occurs correspond to the range of the values used in and obtained with GCMs, suggesting that this analysis can provide a conceptual foundation for analyzing the stability of the interhemispheric THC, and also for the potential of the Atlantic THC to undergo abrupt change.

Description: The objective of this study is to explore, based on the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalysis data, the intraseasonal variability of the South China Sea (SCS) summer monsoon (SM) in terms of its structure and propagation, as well as interannual variations. A possible mechanism that is responsible for the origin of the 10–20-day oscillation of the SCS SM is also proposed. The 30–60-day (hereafter the 3/6 mode) and 10–20-day (hereafter the 1/2 mode) oscillations are found to be the two intraseasonal modes that control the behavior of the SCSSM activities for most of the years. Both the 3/6 and 1/2 modes are distinct, but may not always exist simultaneously in a particular year, and their contributions to the overall variations differ among different years. Thus, the interannual variability in the intraseasonal oscillation activity of the SCS SM may be categorized as follows: the 3/6 category, in which the 3/6 mode is more significant (in terms of the percentage of variance explained) than the 1/2 mode; the 1/2 category, in which the 1/2 mode is dominant; and the dual category, in which both the 3/6 and 1/2 modes are pronounced. Composite analyses of the 3/6 category cases indicate that the 30–60-day oscillation of the SCS SM exhibits a trough–ridge seesaw in which the monsoon trough and subtropical ridge exist alternatively over the SCS, with anomalous cyclones (anticyclones), along with enhanced (suppressed) convection, migrating northward from the equator to the midlatitudes. The northward-migrating 3/6-mode monsoon trough–ridge in the lower troposphere is coupled with the eastward-propagating 3/6-mode divergence–convergence in the upper troposphere. It is also found that, for the years in the dual category, the SCS SM activities are basically controlled by the 3/6 mode, but modified by the 1/2 mode. Composite results of the 1/2-mode category cases show that the 10–20-day oscillation is manifest as an anticyclone–cyclone system over the western tropical Pacific, propagating northwestward into the SCS. A close coupling also exists between the upper-level convergence (divergence) and the low-level anticyclone (cyclone). It is found that the 1/2 mode of the SCS SM mainly originates from the equatorial central Pacific, although a disturbance from the northeast of the SCS also contributes to this mode. The flow patterns from an inactive to an active period resemble those associated with a mixed Rossby–gravity wave observed in previous studies.

Description: In an earlier paper by Verbitsky and Saltzman, a vertically integrated, high-resolution, nonlinearly viscous, nonisothermal ice sheet model was presented to calculate the “present-day” equilibrium regime of the Antarctic ice sheet. Steady-state solutions for the ice topography and thermodynamics, represented by the extent of the areas of basal melting, were shown to be in good agreement with both observations and results obtained from other three-dimensional thermodynamical equations. The solution for the basal temperature field of the West Antarctic Siple Coast produced areas at the pressure melting point separated by strips of frozen-to-bed ice, the structure of which looks very similar to ice streams A–E. Since the possible response of the Siple Coast basal temperature pattern to global warming and to associated changes in the snowfall rate is not obvious, a special sensitivity study was conducted. Results of such a study suggest that increased precipitation rate and associated intensification of ice advection can effectively “shut down” West Antarctic ice streams.

Description: The mechanism by which the diurnal cycle of solar radiation modulates intraseasonal SST variability in the western Pacific warm pool is investigated using a one-dimensional mixed layer model. SSTs in the model experiments forced with hourly surface fluxes during the calm–sunny phase of intraseasonal oscillation are significantly warmer than those with daily mean surface fluxes. The difference in two experiments is explained by upper-ocean mixing processes during nighttime. Surface warming during daytime creates a shallow diurnal warm layer near the surface (0–3 m), which can be easily eroded by surface cooling during nighttime. Further cooling, however, requires a substantial amount of energy because deeper waters need to be entrained into the mixed layer. Since the shallow diurnal layer is not formed in the experiment with daily mean surface fluxes, the SST for the hourly forcing case is warmer most of the time due to the diurnally varying solar radiation. Sensitivity of the intraseasonal SST variation to the penetrative component of solar radiation is examined, showing that the diurnal cycle plays an important role in the sensitivity. Solar radiation absorbed in the upper few meters significantly influences intraseasonal SST variations through changes in amplitude of diurnal SST variation.

Description: This paper compares high-elevation surface temperatures based on the Global Historical Climate Network/Climatic Research Unit (GHCN/CRU) and snow telemetry (SNOTEL) datasets, with simultaneous free-air equivalent temperatures, interpolated from NCEP–NCAR reanalysis. Mean monthly temperature anomalies from 1982 to 1999 are examined for 60 SNOTEL and 296 GHCN/CRU sites at elevations over 500 m with relatively homogenous records. The surface/free-air temperature difference ΔT (Ts − Ta) is calculated for both the SNOTEL and GHCN/CRU datasets. Topography influences the correlation between surface and free-air temperature anomalies. Physically realistic diurnal and seasonal changes in ΔTE are illustrated. Systematic secular trends in surface temperatures, free-air temperatures, and ΔT are revealed, but the sign and magnitude of change depends on location, meaning that regional signals are weak. The Ts trends are positive for most GHCN and CRU sites, and for SNOTEL sites at night. Daytime cooling in the SNOTEL network reduces the mean daily warming trend. The Ta trends are consistently positive for both networks and are often larger than Ts. Thus mean ΔT trends are negative for both datasets. The smaller sample size in the SNOTEL dataset means that error estimates for regional signals are much wider than for the GHCN/CRU dataset. Trend difference maps identify potentially anomalous SNOTEL records. Trends show no correlation with elevation and topography. Surface trends show higher variability and account for most of the uncertainty in ΔT trends. Sensitivity of trends to time period is also discussed. Such changes in the free-air/surface temperature difference may indicate change in the energy balance of mountain areas.

Description: Using a simple stochastic model, the authors illustrate that the occurrence of a meridional dipole in the first empirical orthogonal function (EOF) of a time-dependent zonal jet is a simple consequence of the north–south excursion of the jet center, and this geometrical fact can be understood without appealing to fluid dynamical principles. From this it follows that one ought not, perhaps, be surprised at the fact that such dipoles, commonly referred to as the Arctic Oscillation (AO) or the Northern Annular Mode (NAM), have robustly been identified in many observational studies and appear to be ubiquitous in atmospheric models across a wide range of complexity.

Description: A new and distinctly interdecadal signal in the climate of the Pacific Ocean has been uncovered by examining the coupled behavior of sea surface temperatures (SSTs) and Northern Hemisphere atmospheric circulation. This interdecadal Pacific signal (IPS) of ocean–atmosphere interaction exhibits a highly statistically significant interdecadal component yet contains little to no interannual (El Niño scale) variability common to other Pacific climate anomaly patterns. The IPS thus represents the only empirically derived, distinctly interdecadal signal of Pacific Ocean SST variability that likely also represents the true interdecadal behavior of the Pacific Ocean–atmosphere system. The residual variability of the Pacific’s leading SST pattern, after removal of the IPS, is highly correlated with El Niño anomalies. This indicates that by simply including an atmospheric component, the leading mode of Pacific SST variability has been decomposed into its interdecadal and interannual patterns. Although the interdecadal signal is unrelated to interannual El Niño variability, the interdecadal ocean–atmosphere variability still seems closely linked to tropical Pacific SSTs. Because prior abrupt changes in Pacific SSTs have been related to anomalies in a variety of physical and biotic parameters throughout the Northern Hemisphere, and because of the persistence of these changes over several decades, isolation of this interdecadal signal in the Pacific Ocean–atmosphere system has potentially important and widespread implications to climate forecasting and climate impact assessment.

Description: A 110-yr simulation is conducted using a specially designed coupled ocean–atmosphere general circulation model that only allows air–sea interaction over the Atlantic Ocean within 30°S–60°N. Since the influence from the Pacific El Niño–Southern Oscillation (ENSO) over the Atlantic is removed in this run, it provides a better view of the extratropical influences on the tropical air–sea interaction within the Atlantic sector. The model results are compared with the observations that also have their ENSO components subtracted. The model reproduces the two major anomalous patterns of the sea surface temperature (SST) in the southern subtropical Atlantic (SSA) and the northern tropical Atlantic (NTA) Ocean. The SSA pattern is phase locked to the annual cycle. Its enhancement in austral summer is associated with atmospheric disturbances from the South Atlantic during late austral spring. The extratropical atmospheric disturbances induce anomalous trade winds and surface heat fluxes in its northern flank, which generate SST anomalies in the subtropics during austral summer. The forced SST anomalies then change the local sea level pressure and winds, which in turn affect the northward shift of the atmospheric disturbance and cause further SST changes in the deep Tropics during austral fall. The NTA pattern is significant throughout a year. Like the SSA pattern, the NTA pattern in boreal winter–spring is usually associated with the heat flux change caused by extratropical atmospheric disturbances, such as the North Atlantic Oscillation. The SST anomalies then feed back with the tropical atmosphere and expand equatorward. From summer to fall, however, the NTA SST anomalies are likely to persist within the subtropics for more than one season after it is generated. Our model results suggest that this feature is associated with a local feedback between the NTA SST anomalies and the atmospheric subtropical anticyclone from late boreal summer to early winter. The significance of this potential feedback in reality needs to be further examined with more observational evidence.

Description: The impact of the El Niño–Southern Oscillation (ENSO) on the atmospheric intraseasonal variability in the North Pacific is assessed, with emphasis on how ENSO modulates midlatitude circulation anomalies associated with the Madden–Julian oscillation (MJO) in the Tropics and the westward-traveling patterns (WTP) in high latitudes. The database for this study consists of the output of a general circulation model (GCM) experiment subjected to temporally varying sea surface temperature (SST) forcing in the tropical Pacific, and observational reanalysis products. Diagnosis of the GCM experiment indicates a key region in the North Pacific over which the year-to-year variation of intraseasonal activity is sensitive to the SST conditions in the Tropics. In both the simulated and observed atmospheres, the development phase of the dominant circulation anomaly in this region is characterized by incoming wave activity from northeast Asia and the subtropical western Pacific. Southeastward dispersion from the North Pacific to North America can be found in later phases of the life cycle of the anomaly. The spatial pattern of this recurrent extratropical anomaly contains regional features that are similar to those appearing in composite charts for prominent episodes of the MJO and the WTP. Both the GCM and reanalysis data indicate that the amplitude of intraseasonal variability near the key region, as well as incoming wave activity in the western Pacific and dispersion to the western United States, are enhanced in cold ENSO events as compared to warm events. Similar modulations of the MJO-related circulation patterns in the extratropics by ENSO forcing are discernible in the model simulation. It is inferred from these findings that ENSO can influence the North Pacific intraseasonal activity through its effects on the evolution of convective anomalies in the tropical western Pacific. On the other hand, there is little modification by ENSO of the circulation features associated with the WTP. The combined effect of the MJO and WTP on the intraseasonal circulation in the North Pacific is studied. Based on multiple regression analysis, it is found that the MJO and WTP make comparable contributions to the variability in the midlatitude North Pacific. These contributions may be treated as a linear combination of the anomalies attributed to the MJO and WTP separately.

Description: This study aims to examine the relationship between the tropical Atlantic latent heat flux and convective cloud coverage over northeast Brazil (NEB) during the four months of the main rainy season (February–May). The correlation with anomalies of these data is investigated, both without lag and with a 1-month lag (the heat flux in advance). In both cases, a significant positive correlation appears in the northwestern tropical Atlantic, and a significant negative correlation is obtained for a limited area off eastern NEB. These two correlation patterns are linked to anomalies in the trade wind intensity and in the meridional position of the intertropical convergence zone (ITCZ), which relate to the latent heat flux anomalies and NEB convective coverage anomalies, respectively. The positive correlation pattern is spread over a large part of the northern tropical Atlantic, whereas the negative correlation pattern is confined off NEB. This indicates the existence of different regional mechanisms in the tropical Atlantic basin. The impact of the Atlantic heat fluxes on NEB convection is somewhat different from the classical meridional dipole related to the SST variability. The analysis of the horizontal moisture flux shows that during flood years an additional meridional inflow balances the eastward loss, and the upward velocity reinforced over NEB contributes to intensify NEB convection. The positive correlation pattern indicates that the location of the northern branch of the Pilot Research moored Array in the Tropical Atlantic (PIRATA) moorings is pertinent to monitor the ocean–atmosphere interface parameters. The negative correlation pattern off NEB provides new support for the possible extension of the PIRATA array toward the Brazilian coast. Complementary results at 1-month lag and the real-time availability of the PIRATA data confirm the potential of NEB forecasting.

Description: The Indian Ocean monsoon (IOM) exhibits considerable year-to-year variations that have previously been attributed to a number of forcing mechanisms including the El Niño–Southern Oscillation (ENSO) and Eurasian snow cover anomalies. In this study, spatial data of Eurasian spring land surface temperatures are analyzed as well as proxies for soil moisture, summer IOM precipitation, and summer IOM 850-mb zonal winds for the 1979–99 period to isolate correlated modes of variability. The results indicate the existence of a prominent mode that appears to be related to the boreal winter Arctic Oscillation (AO); this mode projects strongly on the June precipitation and 850-mb zonal wind fields in the vicinity of the IOM region. Its projection on spatial fields of temperature and proxies for soil moisture shows springtime surface warming and drying in the region to the north and west of the Indian subcontinent and cooling over the higher Eurasian latitudes during years of anomalously intense June monsoon rainfall. Such surface signatures are consistent with the negative phase winter AO. It is hypothesized that the preconditioning of the spring season surface characteristics may be associated with an AO-induced quasi-stationary tropospheric circulation anomaly: the impact of this anomaly is to displace the mid-Eastern jet poleward during AO-negative phases, resulting in anomalous surface heating and drying that persist into the later spring season and finally affect the rainfall over the IOM region in June.

Description: A principal component analysis of the combined fields of sea surface temperature (SST) and surface zonal and meridional wind reveals that the dominant mode of intraseasonal (30 to 70 day) covariability during northern winter in the tropical Eastern Hemisphere is that of the Madden–Julian oscillation (MJO). Regression calculations show that the submonthly (30-day high-pass filtered) surface wind variability is significantly modulated during the MJO. Regions of increased (decreased) submonthly surface wind variability propagate eastward, approximately in phase with the intraseasonal surface westerly (easterly) anomalies of the MJO. Because of the dependence of the surface latent heat flux on the magnitude of the total wind speed, this systematic modulation of the submonthly surface wind variability produces a significant component in the intraseasonal latent heat flux anomalies, which partially cancels the latent heat flux anomalies due to the slowly varying intraseasonal wind anomalies, particularly south of 10°S. A method is derived that demodulates the submonthly surface wind variability from the slowly varying intraseasonal wind anomalies. This method is applied to the wind forcing fields of a one-dimensional ocean model. The model response to this modified forcing produces larger intraseasonal SST anomalies than when the model is forced with the observed forcing over large areas of the southwest Pacific Ocean and southeast Indian Ocean during both phases of the MJO. This result has implications for accurate coupled modeling of the MJO. A similar calculation is applied to the surface shortwave flux, but intraseasonal modulation of submonthly surface shortwave flux variability does not appear to be important to the dynamics of the MJO.

Description: A merged land–air–sea surface temperature reconstruction analysis is developed for monthly anomalies. The reconstruction is global and spatially complete. Reconstructed anomalies damp toward zero in regions with insufficient sampling. Error estimates account for the damping associated with sparse sampling, and also for bias uncertainty in both the land and sea observations. Averages of the reconstruction are similar to simple averages of the unanalyzed data for most of the analysis period. For the nineteenth century, when sampling is most sparse and the error estimates are largest, the differences between the averaged reconstruction and the simple averages are largest. Sampling is always sparse poleward of 60° latitude, and historic reconstructions for the polar regions should be used with caution.

Description: An atmospheric stationary wave teleconnection mechanism is proposed to explain how ENSO may affect the Tibetan Plateau snow depth and thereby the south Asian monsoons. Using statistical analysis, the short available record of satellite estimates of snow depth, and ray tracing, it is shown that wintertime ENSO conditions in the central Pacific may produce stationary barotropic Rossby waves in the troposphere with a northeastward group velocity. These waves reflect off the North American jet, turning equatorward, and enter the North African–Asian jet over the eastern Atlantic Ocean. Once there, the waves move with the jet across North Africa, South Asia, the Himalayas, and China. Anomalous increases in upper-tropospheric potential vorticity and increased wintertime snowfall over the Tibetan Plateau are speculated to be associated with these Rossby waves. The increased snowfall produces a larger Tibetan Plateau snowpack, which persists through the spring and summer, and weakens the intensity of the south Asian summer monsoons.

Description: Two 9-yr runs of the NCAR Community Climate Model version 3 (CCM3) are compared in their simulations of the North American summer monsoon. In a control simulation, the Zhang–McFarlane deep convection scheme is used. For an experimental simulation, the following modifications to the scheme are implemented. The closure is based on the large-scale forcing of virtual temperature, and a relative humidity threshold on convective parcels lifted from the boundary layer is applied. The sensitivity to these modifications for simulating the North American monsoon is investigated. Model validation relies on hourly precipitation rates from surface gauges over the United States, hourly precipitation rates derived from the combination of microwave and radar measurements from NASA’s Tropical Rainfall Measuring Mission (TRMM) satellite over Mexico, and CAPE values as calculated from temperature, specific humidity, and pressure fields from the NCEP–NCAR reanalysis. Results show that the experimental run improves the timing of the monsoon onset and peak in the regions of core monsoon influence considered here, though it increases a negative bias in the peak monsoon intensity in one region of northern Mexico. Sensitivity of the diurnal cycle of precipitation to modifications in the convective scheme is highly geographically dependent. Using a combination of gauge-based rainfall rates and reanalysis-based CAPE, it is found that improvements in the simulated diurnal cycle are confined to a convective regime in which the diurnal evolution of precipitation is observed to lag that of CAPE. For another regime, in which CAPE is observed to be approximately in phase with precipitation, model phase biases increase nearly everywhere. Some of the increased phase biases in the latter regime are primarily because of application of the relative humidity threshold.

Description: A catalog containing an unprecedented amount of historical data in the southeastern part of Mexico covering almost four centuries (1502–1899) is used to construct a drought time series. The catalog records information of agricultural disasters and includes events associated with hydrometeorological phenomena or hazards whose effects were mainly felt in the agricultural sector, such as droughts. An analysis of the historical series of droughts in southeastern Mexico for the period 1502–1899 is performed. The highest drought frequency occurred around the years 1650, 1782, and 1884; no droughts were reported around 1540, between 1630 and 1640, along the largest time lapse of 1672–1714, and between 1740 and 1760. From 1760 until the end of the period of study droughts definitively occur more often than they did from ∼1550 to 1760. In addition, most droughts lasted for 1–2 yr. Analyzing the frequencies of the drought time series it is found that the most conspicuous cycles are ∼3–4 and 7 yr, although cycles of ∼12, 20, 43, and 70 yr are also evident. The relation between droughts and El Niño events indicates that 38% of droughts are associated with El Niño. Sea surface temperature changes, the Southern Oscillation index, and solar activity leave their signals in the southeastern part of Mexico, with the signs in Oaxaca clearer than in the Yucatan Peninsula. However, the dominance of some phenomena over others depends on the time scales considered.

Description: The mechanisms that drove zonal wind stress (τx) changes in the near-equatorial Pacific at the end of the extreme 1997–98 El Niño event are explored using a global atmospheric general circulation model. The analysis focuses on three features of the τx evolution between October 1997 and May 1998 that were fundamental in driving the oceanic changes at the end of this El Niño event: (i) the southward shift of near-date-line surface zonal wind stress (τx) anomalies beginning November 1997, (ii) the disappearance of the easterly τx from the eastern equatorial Pacific (EEqP) in February 1998, and (iii) the reappearance of easterly τx in the EEqP in May 1998. It is shown that these wind changes represent the deterministic response of the atmosphere to the observed sea surface temperature (SST) field, resulting from changes in the meridional structure of atmospheric convective anomalies in response to the seasonally phase-locked meridional movement of the warmest SST. The southward shift of the near-date-line τx anomalies at the end of this El Niño event was controlled by the seasonal movement of the warmest SST south of the equator, which—both directly and through its influence on the atmospheric response to changes in SST anomaly—brought the convective anomalies from being centered about the equator to being centered south of the equator. The disappearance (reappearance) of easterly EEqP τx has only been evident in extreme El Niño events and has been associated with the development (northward retreat) of an equatorial intertropical convergence zone (ITCZ). The disappearance/return of EEqP easterly τx arises in the AGCM as the deterministic response to changes in the SST field, tied principally to the changes in climatological SST (given time-invariant extreme El Niño SSTA) and not to changes in the underlying SSTA field. The disappearance (return) of EEqP easterly τx in late boreal winter (late boreal spring) is a characteristic atmospheric response to idealized extreme El Niño SST anomalies; this suggests that the distinctive termination of the 1997–98 El Niño event is that to be expected for extreme El Niño events.

Description: The multidecadal climate variability in the North Pacific region is investigated by using a 2000-yr-long integration with a coupled ocean–atmosphere general circulation model. It is shown that the multidecadal variability evolves largely independent of the variations in the tropical Pacific, so that this kind of multidecadal variability may be regarded as internal to the North Pacific. The coupled model results suggest that the multidecadal variability can be explained by the dynamical ocean response to stochastic wind stress forcing. Superimposed on the red background variability, a multidecadal mode with a period of about 40 yr is simulated by the coupled model. This mode can be understood through the concept of spatial resonance between the ocean and the atmosphere.

Description: An effort was made to search for relationships between interannual variations of population, lifetime, genesis locations, and intensity of named typhoons and numbered tropical depressions in the western North Pacific during the 1979–2002 period. To support this research task, climatological relationships of tropical cyclone characteristics were also investigated for these cyclones. Major findings of this study are summarized as follows:Climatology: Measured by the intensity scale of the Japan Meteorological Agency, three groups of tropical cyclones were identified in terms of population versus intensity: Group 1 [tropical depression (TD) + typhoon (TY)], Group 2 (strong + very strong TY), and Group 3 (catastrophic TY). This group division coincides with that formed in terms of lifetime of tropical cyclones versus intensity. Weak cyclones (Group 1) have a larger population than strong cyclones (Group 3), while the former group has shorter lifetime than the latter group. For genesis locations, the monsoon trough is established as a favorable region of tropical cyclone genesis because it provides an environment of large vorticity. Therefore, the northward latitudinal displacement of the maximum genesis frequency in the three groups of tropical cyclones follows that of the monsoon trough.Interannual variation: Any mechanism that can modulate the location and intensity of the monsoon trough affects the genesis location and frequency of tropical cyclones. In response to tropical Pacific sea surface temperature anomalies, a short wave train consisting of east–west oriented cells emanates from the Tropics and progresses along the western North Pacific rim. Population of the Group-1 tropical cyclones varies interannually in phase with the oscillation of the anomalous circulation cell northeast of Taiwan and south of Japan in this short wave train, while that of Group 3 fluctuates coherently with the tropical cell of this short wave train. Because these two anomalous circulation cells exhibit opposite polarity, the out-of-phase interannual oscillation between these two cells results in the opposite interannual variation of genesis frequency between tropical cyclones of Groups 1 and 3.

Description: Intraseasonal oscillations (ISOs) are important large-amplitude and large-scale elements of the tropical Indo-Pacific climate with time scales in the 20–60-day period range, during which time they modulate higher-frequency tropical weather. Despite their importance, the ISO is poorly simulated and predicted by numerical models. A joint diagnostic and modeling study of the ISO is conducted, concentrating on the period between the suppressed and active (referred to as the “transition”) period that is hypothesized to be the defining stage for the development of the intraseasonal mode and the component that is most poorly simulated. The diagnostic study uses data from the Tropical Ocean Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE). It is found that during the transition period, the ocean and the atmosphere undergo gradual but large-scale and high-amplitude changes, especially the moistening of the lower troposphere caused jointly by the anomalously warm sea surface temperature arising from minimal cloud and low winds during the suppressed phase and the large-scale subsidence that inhibits the formation of locally deep convection. Using a cloud classification scheme based on microwave and infrared satellite data, it is observed that midtop (cloud with a top in the middle troposphere) nonprecipitating clouds are a direct response of the low-level moisture buildup. To investigate the sensitivity of ISO simulations to the transitional phase, the European Centre for Medium-Range Weather Forecasts (ECMWF) coupled ocean–atmosphere climate model is used. The ECMWF was run serially in predictive ensemble mode (five members) for 30-day periods starting from 1 December 1992 to 30 January 1993, encompassing the ISO occurring in late December. Predictability of the active convective period of the ISO is poor when initialized before the transitional phases of the ISO. However, when initialized with the correct lower-tropospheric moisture field, predictability increases substantially, although the model convective parameterization appears to trigger convection too quickly without allowing an adequate buildup of convective available potential energy during the transition period.

Description: The climate is changing because of human activities and will continue to do so regardless of any mitigation actions. Available climate observations and information are also changing as technological advances take place. Accordingly, an overview is given of a much-needed potential climate information system that embraces a comprehensive observing system to observe and track changes and the forcings of the system as they occur, and that develops the ability to relate one to the other and understand changes and their origins. Observations need to be taken in ways that satisfy the climate monitoring principles and ensure long-term continuity, and that have the ability to discern small but persistent signals. Some benchmark observations are proposed to anchor space-based observations and trends, including a much-needed step forward in the quality of water vapor observations. Satellite observations must be calibrated and validated, with orbital decay and drift effects fully dealt with if possible, and adequate overlap to ensure continuity. The health of the monitoring system must be tracked and resources identified to address issues. Fields must be analyzed into global products and delivered to users while stakeholder needs are fully considered. Data should be appropriately archived with full and open access, along with metadata that fully describe the observing system status and environment in which it operates. Reanalysis of the records must be institutionalized along with continual assessment of impacts of new observing and analysis systems. Some products will be used to validate and improve models, as well as initialize models and predict future evolution on multiple time scales using ensembles. Attribution of changes to causes is essential, and it is vital to fully assess past changes and model performance and results in making predictions to help appraise reliability and assess impacts regionally on the environment, human activities, and sectors of the economy. In particular, a revolution in the way developing countries use and apply climate information is expected. Such a system will be invaluable and further provides a framework for setting priorities of new observations and related activities. Without the end-to-end process the investments will not deliver adequate return and our understanding will be much less than it would be otherwise.

Description: The response of the atmospheric circulation to an enhanced radiative greenhouse gas forcing in a transient integration with a coupled global climate model is investigated. The spatial patterns of the leading modes of Northern Hemisphere atmospheric variability are shown to change in response to the enhanced forcing. An earlier study showed that the spatial patterns of the leading modes in the Southern Hemisphere changed in response to the enhanced forcing. These changes were associated with changes in the propagation conditions for barotropic Rossby waves. This is, however, not the case for the Northern Hemisphere, where the propagation conditions are unchanged. Other possible mechanisms for the changes in the spatial patterns of the leading modes are discussed.

Description: The tropical oceans have long been recognized as the most important region for large-scale ocean–atmosphere interactions, giving rise to coupled climate variations on several time scales. During the Tropical Ocean Global Atmosphere (TOGA) decade, the focus of much tropical ocean research was on understanding El Niño–related processes and on development of tropical ocean models capable of simulating and predicting El Niño. These studies led to an appreciation of the vital role the ocean plays in providing the memory for predicting El Niño and thus making seasonal climate prediction feasible. With the end of TOGA and the beginning of Climate Variability and Prediction (CLIVAR), the scope of climate variability and predictability studies has expanded from the tropical Pacific and ENSO-centric basis to the global domain. In this paper the progress that has been made in tropical ocean climate studies during the early years of CLIVAR is discussed. The discussion is divided geographically into three tropical ocean basins with an emphasis on the dynamical processes that are most relevant to the coupling between the atmosphere and oceans. For the tropical Pacific, the continuing effort to improve understanding of large- and small-scale dynamics for the purpose of extending the skill of ENSO prediction is assessed. This paper then goes beyond the time and space scales of El Niño and discusses recent research activities on the fundamental issue of the processes maintaining the tropical thermocline. This includes the study of subtropical cells (STCs) and ventilated thermocline processes, which are potentially important to the understanding of the low-frequency modulation of El Niño. For the tropical Atlantic, the dominant oceanic processes that interact with regional atmospheric feedbacks are examined as well as the remote influence from both the Pacific El Niño and extratropical climate fluctuations giving rise to multiple patterns of variability distinguished by season and location. The potential impact of Atlantic thermohaline circulation on tropical Atlantic variability (TAV) is also discussed. For the tropical Indian Ocean, local and remote mechanisms governing low-frequency sea surface temperature variations are examined. After reviewing the recent rapid progress in the understanding of coupled dynamics in the region, this study focuses on the active role of ocean dynamics in a seasonally locked east–west internal mode of variability, known as the Indian Ocean dipole (IOD). Influences of the IOD on climatic conditions in Asia, Australia, East Africa, and Europe are discussed. While the attempt throughout is to give a comprehensive overview of what is known about the role of the tropical oceans in climate, the fact of the matter is that much remains to be understood and explained. The complex nature of the tropical coupled phenomena and the interaction among them argue strongly for coordinated and sustained observations, as well as additional careful modeling investigations in order to further advance the current understanding of the role of tropical oceans in climate.

Description: This paper briefly surveys areas of paleoclimate modeling notable for recent progress. New ideas, including hypotheses giving a pivotal role to sea ice, have revitalized the low-order models used to simulate the time evolution of glacial cycles through the Pleistocene, a prohibitive length of time for comprehensive general circulation models (GCMs). In a recent breakthrough, however, GCMs have succeeded in simulating the onset of glaciations. This occurs at times (most recently, 115 kyr b.p.) when high northern latitudes are cold enough to maintain a snow cover and tropical latitudes are warm, enhancing the moisture source. More generally, the improvement in models has allowed simulations of key periods such as the Last Glacial Maximum and the mid-Holocene that compare more favorably and in more detail with paleoproxy data. These models now simulate ENSO cycles, and some of them have been shown to reproduce the reduction of ENSO activity observed in the early to middle Holocene. Modeling studies have demonstrated that the reduction is a response to the altered orbital configuration at that time. An urgent challenge for paleoclimate modeling is to explain and to simulate the abrupt changes observed during glacial epochs (i.e., Dansgaard–Oescher cycles, Heinrich events, and the Younger Dryas). Efforts have begun to simulate the last millennium. Over this time the forcing due to orbital variations is less important than the radiance changes due to volcanic eruptions and variations in solar output. Simulations of these natural variations test the models relied on for future climate change projections. They provide better estimates of the internal and naturally forced variability at centennial time scales, elucidating how unusual the recent global temperature trends are.

Description: It is argued that simulations of the twentieth century performed with coupled global climate models with specified historical changes in external radiative forcing can be interpreted as climate hindcasts. A simple Bayesian method for postprocessing such simulations is described, which produces probabilistic hindcasts of interdecadal temperature changes on large spatial scales. Hindcasts produced for the last two decades of the twentieth century are shown to be skillful. The suggestion that skillful decadal forecasts can be produced on large regional scales by exploiting the response to anthropogenic forcing provides additional evidence that anthropogenic change in the composition of the atmosphere has influenced the climate. In the absence of large negative volcanic forcing on the climate system (which cannot presently be forecast), it is predicted that the global mean temperature for the decade 2000–09 will lie above the 1970–99 normal with a probability of 0.94. The global mean temperature anomaly for this decade relative to 1970–99 is predicted to be 0.35°C with a 5%–95% confidence range of 0.21°–0.48°C.

Description: In this paper, a new two-dimensional blocking index is proposed by defining a difference between the daily 500-hPa geopotenial heights at the reference latitude and its north side. The reference latitude is determined by a composite latitude-dependent 500-hPa geopotenial height of blocking events in different seasons and sectors. The new index can take account of the duration, intensity, propagation, and spatial structure of a blocking event. Using this index, the characteristics (frequency, duration, intensity, and preferred occurrence region) of the blocking action in the North Hemisphere (NH) are investigated using a 42-yr sample of blocking events from the NCEP–NCAR reanalyses. It is found that blocking events in the NH are more frequent in the Atlantic–Europe sector than in the Pacific sector in winter and spring and autumn, but more persistent in the Atlantic–Europe sector than in the Pacific sector for all seasons. Blocking events in the Pacific sector tend to have larger amplitudes than the Atlantic counterparts. In addition, it is shown that in the NH independently occurring blocking events are most frequent, but simultaneously occurring blockings are rather rare, indicating that the blocking events in the NH should be a local phenomenon. On the other hand, a comparison with the existing indices [e.g., Tibaldi and Molteni (TM) index] indicates that in summer and autumn the new index shows similar longitudinal dependency of NH blocking events as does the TM index, but it shows two distinct action centers of blocking events in the Atlantic sector in winter and spring (in which the most frequent one is situated more westward) and an eastward blocking action center in the Pacific sector in spring and autumn, compared to the TM index. In addition, it is found that the new blocking index proposed here shows relatively low blocking frequency for all seasons compared to the TM index, especially in the Atlantic sector in spring and in the Pacific sector in winter and spring, which seems to be in agreement with the result obtained by Pelly and Hoskins using the PV–θ index.

Description: The warming of the near-surface air in the Arctic region has been larger than the global mean surface warming. There is general agreement that the Arctic amplification of the surface air temperature (SAT) trend to a considerable extent is due to local effects such as the retreat of sea ice, especially during the summer months, and earlier melting of snow in the spring season. There is no doubt that these processes are important causes of the Arctic SAT trend. It is less clear, however, whether the trend may also be related to recent changes in the atmospheric midlatitude circulation. This question is the focus of the present paper. Model experiments have shown that in a warmer climate responding to, for example, a doubling of CO2, the atmospheric northward energy transport (ANET) will increase and cause polar SAT amplification. In the present study, the development of the ANET across 60°N and its linkage to the Arctic SAT have been explored using the ERA-40 reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF). It is found that during 1979–2001, the ANET has experienced an overall positive but weak trend, which was largest during the period from the mid-1980s to the mid-1990s. In addition, it is found that the Arctic SAT is sensitive to variability of the ANET across 60°N and hence to variability of the midlatitude circulation: A large ANET is followed by warming of the Arctic where ANET leads by about 5 days. The warming is located primarily north of the Atlantic and Pacific sectors, indicating that baroclinic weather systems developing around the Icelandic and Aleutian lows are important for the energy transport. Furthermore, it is suggested here that a small, but statistically significant, part of the mean Arctic SAT trend is linked to the trend in the ANET. Another important indicator of the midlatitude circulation is the Arctic Oscillation (AO). Through the 1980s and early 1990s the AO index has shown a positive trend. However, even though a part of the SAT trend can be related to the AO in localized parts of the Arctic area, the mean Arctic SAT trend shows no significant linkage to the AO.

Description: There was a dramatic decrease in rainfall in the southwest of Australia (SWA) in the mid-1960s. A statistical method, based on the idea of analogous synoptic situations, is used to help clarify the cause of the drying. The method is designed to circumvent error in the rainfall simulated directly by a climate model, and to exploit the ability of the model to simulate large-scale fields reasonably well. The method uses relationships between patterns of various atmospheric fields with station records of rainfall to improve the simulation of the local rainfall spatial variability. The original technique was developed in a previous study. It is modified here for application to two four-member ensembles of simulations of the climate from 1870 to 1999 performed with the Parallel Climate Model (PCM). The first ensemble, called “natural,” is forced with natural variations in both volcanic activity and solar forcing. The second ensemble, called “full forcing,” also includes three types of human-induced forcing resulting from changes in greenhouse gases, ozone, and aerosols. The full-forcing runs provide a better match to observational changes in sea surface temperature in the vicinity of SWA. The observed rainfall decline is not well captured by rainfall changes simulated directly by the model in either ensemble. There is a hint that the fully forced ensemble is more realistic, but it is nothing more than a hint. The downscaling approach, on the other hand, provides a much more accurate reproduction of the day-to-day variability of rainfall in SWA than the rainfall simulated directly by the model. The downscaled ensemble mean rainfall in full forcing declines over the region with a spatial pattern that is similar to the observed decline. This contrasts with an increase of rainfall in the downscaled rainfall in the natural ensemble. These results give the clearest indication yet that anthropogenic forcing played a role in the drying of SWA. Note, however, that ambiguities remain. For example, although the observed decline fits within the range of downscaled model simulation, the ensemble mean rainfall decline is only about half of the observed estimate, the timing differs from the observations, drying did not occur in the downscaling of one of the four full-forced ensemble members, and not all potential forcing mechanisms are included in full forcing (e.g., land surface changes). Furthermore, while the observed rainfall decline was a sharp reduction in the 1960s, followed by a near-constant rainfall regime, the full-forcing ensemble suggests a more gradual rainfall decline over 40 yr from 1960.

Description: Interannual rainfall extremes over southwest Western Australia (SWWA) are examined using observations, reanalysis data, and a long-term natural integration of the global coupled climate system. The authors reveal a characteristic dipole pattern of Indian Ocean sea surface temperature (SST) anomalies during extreme rainfall years, remarkably consistent between the reanalysis fields and the coupled climate model but different from most previous definitions of SST dipoles in the region. In particular, the dipole exhibits peak amplitudes in the eastern Indian Ocean adjacent to the west coast of Australia. During dry years, anomalously cool waters appear in the tropical/subtropical eastern Indian Ocean, adjacent to a region of unusually warm water in the subtropics off SWWA. This dipole of anomalous SST seesaws in sign between dry and wet years and appears to occur in phase with a large-scale reorganization of winds over the tropical/subtropical Indian Ocean. The wind field alters SST via anomalous Ekman transport in the tropical Indian Ocean and via anomalous air–sea heat fluxes in the subtropics. The winds also change the large-scale advection of moisture onto the SWWA coast. At the basin scale, the anomalous wind field can be interpreted as an acceleration (deceleration) of the Indian Ocean climatological mean anticyclone during dry (wet) years. In addition, dry (wet) years see a strengthening (weakening) and coinciding southward (northward) shift of the subpolar westerlies, which results in a similar southward (northward) shift of the rain-bearing fronts associated with the subpolar front. A link is also noted between extreme rainfall years and the Indian Ocean Dipole (IOD). Namely, in some years the IOD acts to reinforce the eastern tropical pole of SST described above, and to strengthen wind anomalies along the northern flank of the Indian Ocean anticyclone. In this manner, both tropical and extratropical processes in the Indian Ocean generate SST and wind anomalies off SWWA, which lead to moisture transport and rainfall extremes in the region. An analysis of the seasonal evolution of the climate extremes reveals a progressive amplification of anomalies in SST and atmospheric circulation toward a wintertime maximum, coinciding with the season of highest SWWA rainfall. The anomalies in SST can appear as early as the summertime months, however, which may have important implications for predictability of SWWA rainfall extremes.

Description: The authors perform an observational study of the relation between stratospheric final warmings (SFWs) and the boreal extratropical circulation. SFW events are found to provide a strong organizing influence upon the large-scale circulation of the stratosphere and troposphere during the period of spring onset. In contrast to the climatological seasonal cycle, SFW events noticeably sharpen the annual weakening of high-latitude circumpolar westerlies in both the stratosphere and troposphere. A coherent pattern of significant westerly (easterly) zonal wind anomalies is observed to extend from the stratosphere to the earth’s surface at high latitudes prior to (after) SFW events, coinciding with the polar vortex breakdown. This evolution is associated with a bidirectional dynamical coupling of the stratosphere–troposphere system in which tropospheric low-frequency waves induce annular stratospheric circulation anomalies, which in turn, are followed by annular tropospheric circulation anomalies. The regional tropospheric manifestation of SFW events consists of a North Atlantic Oscillation (NAO)-like phase transition in the near-surface geopotential height field, with height rises over polar latitudes and height falls over the northeast North Atlantic. This lower-tropospheric change pattern is distinct from the climatological seasonal cycle, which closely follows seasonal trends in thermal forcing at the lower boundary. Although broadly similar, the tropospheric anomaly patterns identified in the study do not precisely correspond to the canonical northern annular mode (NAM) and NAO patterns as the primary anomaly centers are retracted northward toward the pole. The results here imply that (i) high-latitude climate may be particularly sensitive to long-term trends in the annual cycle of the stratospheric polar vortex and (ii) improvements in the understanding and simulation of SFW events may benefit medium-range forecasts of spring onset in the extratropics.

Description: Since the mid-1960s, rapid regional summer warming has occurred on the east coast of the northern Antarctic Peninsula, with near-surface temperatures increasing by more than 2°C. This warming has contributed significantly to the collapse of the northern sections of the Larsen Ice Shelf. Coincident with this warming, the summer Southern Hemisphere Annular Mode (SAM) has exhibited a marked trend, suggested by modeling studies to be predominantly a response to anthropogenic forcing, resulting in increased westerlies across the northern peninsula. Observations and reanalysis data are utilized to demonstrate that the changing SAM has played a key role in driving this local summer warming. It is proposed that the stronger summer westerly winds reduce the blocking effect of the Antarctic Peninsula and lead to a higher frequency of air masses being advected eastward over the orographic barrier of the northern Antarctic Peninsula. When this occurs, a combination of a climatological temperature gradient across the barrier and the formation of a föhn wind on the lee side typically results in a summer near-surface temperature sensitivity to the SAM that is 3 times greater on the eastern side of the peninsula than on the west. SAM variability is also shown to play a less important role in determining summer temperatures at stations west of the barrier in the northern peninsula (∼62°S), both at the surface and throughout the troposphere. This is in contrast to a station farther south (∼65°S) where the SAM exerts little influence.

Description: The relative roles played by the remote El Niño–Southern Oscillation (ENSO) forcing and the local air–sea interactions in the tropical Atlantic are investigated using an intermediate coupled model (ICM) of the tropical Atlantic. The oceanic component of the ICM consists of a six-baroclinic mode ocean model and a simple mixed layer model that has been validated from observations. The atmospheric component is a global atmospheric general circulation model developed at the University of California, Los Angeles (UCLA). In a forced context, the ICM realistically simulates both the sea surface temperature anomaly (SSTA) variability in the equatorial band, and the relaxation of the Atlantic northeast trade winds and the intensification of the equatorial westerlies in boreal spring that usually follows an El Niño event. The results of coupled experiments with or without Pacific ENSO forcing and with or without explicit air–sea interactions in the equatorial Atlantic indicate that the background energy in the equatorial Atlantic is provided by ENSO. However, the time scale of the variability and the magnitude of some peculiar events cannot be explained solely by ENSO remote forcing. It is demonstrated that the peak of SSTA variability in the 1–3-yr band as observed in the equatorial Atlantic is due to the local air–sea interactions and is not a linear response to ENSO. Seasonal phase locking in boreal summer is also the result of the local coupling. The analysis of the intrinsic sustainable modes indicates that the Atlantic El Niño is qualitatively a noise-driven stable system. Such a system can produce coherent interdecadal variability that is not forced by the Pacific or extraequatorial variability. It is shown that when a simple slab mixed layer model is embedded into the system to simulate the northern tropical Atlantic (NTA) SST variability, the warming over NTA following El Niño events have characteristics (location and peak phase) that depend on air–sea interaction in the equatorial Atlantic. In the model, the interaction between the equatorial mode and NTA can produce a dipolelike structure of the SSTA variability that evolves at a decadal time scale. The results herein illustrate the complexity of the tropical Atlantic ocean–atmosphere system, whose predictability jointly depends on ENSO and the connections between the Atlantic modes of variability.

Description: The Commonwealth Scientific and Industrial Research Organisation (CSIRO) Mark 2 global coupled climatic model has been used to generate a 10 000-yr simulation of “present” climate. The resultant dataset has been used to investigate a number of aspects of extremes associated with annual mean rainfall. Multimillennial time series of normalized rainfall amounts for selected points are used to highlight secular variability, spatial variations, and the differences between pluvial and drought conditions. Global distributions are also presented for selected rainfall characteristics, including the frequency of occurrence of specified rainfall anomalies with annual durations, the frequency of occurrence of 5-yr sequences of specified rainfall anomalies, and the maximum and minimum normalized rainfall amounts attained in the simulation. Such features cannot be obtained from observations because of their limited duration. A case study is also made of a megadrought over the southwestern United States, together with an analysis of the associated causal mechanisms. Given the exclusion of all external forcing from the model, it is concluded that the extreme annual mean rainfall extremes presented in the paper are attributable to stochastic events.

Description: Centered composite analysis is described and applied to gain a better understanding of the initial phases of the Madden–Julian oscillation (MJO). Centered composite analysis identifies the dates and central locations of key events. The elements of the composite means are centered on these central locations before averages are calculated. In this way much of the spatial fuzziness, which is inherent in traditional composite analysis, is removed. The results for the MJO, based on MJO-filtered outgoing longwave radiation for the reference data and 40-yr ECMWF Re-Analysis (ERA-40) and NCEP–NCAR reanalysis products for the composites, show highly significant composites of unfiltered data for not only zero lag, but also lags back to 20 days before the target events. These composites identify propagating patterns of surface pressure, upper- and lower-troposphere zonal winds, surface temperature, and 850-hPa specific humidity associated with MJO convective events in the Indian Ocean. The propagation characteristics of important features, especially surface pressure, differ substantially for MJO convective anomalies centered over the Indian or western Pacific Oceans. This suggests that distinctly different mechanisms may be dominant in these two regions, and that many earlier analyses may be mixing properties of the two.

Description: A new coupled general circulation climate model developed at the Met Office's Hadley Centre is presented, and aspects of its performance in climate simulations run for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4) documented with reference to previous models. The Hadley Centre Global Environmental Model version 1 (HadGEM1) is built around a new atmospheric dynamical core; uses higher resolution than the previous Hadley Centre model, HadCM3; and contains several improvements in its formulation including interactive atmospheric aerosols (sulphate, black carbon, biomass burning, and sea salt) plus their direct and indirect effects. The ocean component also has higher resolution and incorporates a sea ice component more advanced than HadCM3 in terms of both dynamics and thermodynamics. HadGEM1 thus permits experiments including some interactive processes not feasible with HadCM3. The simulation of present-day mean climate in HadGEM1 is significantly better overall in comparison to HadCM3, although some deficiencies exist in the simulation of tropical climate and El Niño variability. We quantify the overall improvement using a quasi-objective climate index encompassing a range of atmospheric, oceanic, and sea ice variables. It arises partly from higher resolution but also from greater fidelity in modeling dynamical and physical processes, for example, in the representation of clouds and sea ice. HadGEM1 has a similar effective climate sensitivity (2.8 K) to a CO2 doubling as HadCM3 (3.1 K), although there are significant regional differences in their response patterns, especially in the Tropics. HadGEM1 is anticipated to be used as the basis both for higher-resolution and higher-complexity Earth System studies in the near future.

Description: To improve simulations of the Arctic climate and to quantify climate model errors, four regional climate models [the Arctic Regional Climate System Model (ARCSYM), the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS), the High-Resolution Limited-Area Model (HIRHAM), and the Rossby Center Atmospheric Model (RCA)] have simulated the annual Surface Heat Budget of the Arctic Ocean (SHEBA) under the Arctic Regional Climate Model Intercomparison Project (ARCMIP). The same lateral boundary and ocean surface boundary conditions (i.e., ice concentration and surface temperature) drive all of the models. This study evaluated modeled surface heat fluxes and cloud fields during May 1998, a month that included the onset of the surface icemelt. In general, observations agreed with simulated surface pressure and near-surface air properties. Simulation errors due to surface fluxes and cloud effects biased the net simulated surface heat flux, which in turn affected the timing of the simulated icemelt. Modeled cloud geometry and precipitation suggest that the RCA model produced the most accurate cloud scheme, followed by the HIRHAM model. Evaluation of a relationship between cloud water paths and radiation showed that a radiative transfer scheme in ARCSYM was closely matched with the observation when liquid clouds were dominant. Clouds and radiation are of course closely linked, and an additional comparison of the radiative transfer codes for ARCSYM and COAMPS was performed for clear-sky conditions, thereby excluding cloud effects. Overall, the schemes for radiative transfer in ARCSYM and for cloud microphysics in RCA potentially have some advantages for modeling the springtime Arctic.

Description: The direct effects of aerosols on global and regional climate during boreal spring are investigated based on numerical simulations with the NASA Global Modeling and Assimilation Office finite-volume general circulation model (fvGCM) with Microphyics of Clouds with the Relaxed–Arakawa Schubert Scheme (McRAS), using aerosol forcing functions derived from the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The authors find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excite a planetary-scale teleconnection pattern in sea level pressure, temperature, and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, that is, South Asia, East Asia, and northern and western Africa. Significant atmospheric heating is found in regions with large loading of dust (over northern Africa and the Middle East) and black carbon (over Southeast Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east–west dipole anomaly with strong cooling over the Caspian Sea and warming over central and northeastern Asia, where aerosol concentrations are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection pattern driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations. The surface temperature signature associated with the aerosol-induced teleconnection bears striking resemblance to the spatial pattern of observed long-term trend in surface temperature over Eurasia. Additionally, the boreal spring wave train pattern is similar to that reported by Fukutomi et al. associated with the boreal summer precipitation seesaw between eastern and western Siberia. The results of this study raise the possibility that global aerosol forcing during boreal spring may play an important role in spawning atmospheric teleconnections that affect regional and global climates.

Description: A recent study showed that a tropical Atlantic sea surface temperature (SST) anomaly induces a significant coupled response in late winter [February–April (FMA)] in a coupled model, in which an atmospheric general circulation model is coupled to a slab mixed layer ocean model (AGCM_ML). The coupled response comprises a dipole in the geopotential height, like the North Atlantic Oscillation (NAO), and a North Atlantic tripole in the SST. The simulated NAO response developed 1 or 2 months later in the model than in observations. To determine the possible effects of Ekman heat transport on the development of the coupled response to the tropical forcing, an extended coupled model (AGCM_EML), including Ekman transport in the slab mixed layer ocean, is now used. Large ensembles of AGCM_EML experiments are performed, parallel to the previous AGCM_ML experiments, with the model forced by the same tropical Atlantic SST anomaly over the boreal winter months (September–April). The inclusion of Ekman heat transport is found to result in an earlier development of the coupled NAO–SST tripole response in the AGCM_EML, compared to that in the AGCM_ML. The mutual reinforcement between the anomalous Ekman transport and the surface heat flux causes the tropical forcing to induce an extratropical SST response in November–January (NDJ) in the AGCM_EML that is twice as strong as that in the AGCM_ML. The feedback of this stronger extratropical SST response on the atmosphere in turn drives the development of the NAO response in NDJ. In FMA, the sign of the anomalous surface heat flux is reversed in the Gulf Stream region such that it opposes the anomalous Ekman transport. The resulting equilibrium NAO response in the AGCM_EML is similar to that in the AGCM_ML, but it is reached 1–2 months sooner in the AGCM_EML. Hence, the presence of Ekman transport causes a seasonal shift in the evolution of the coupled response. The faster development of the NAO response in the AGCM_EML suggests that tropical Atlantic SST anomalies should be able to influence the NAO, in nature, on the seasonal time scale, and that efficient interactions with the extratropical ocean play a significant role in determining the coupled response.

Description: To be confident in the analyses of long-term changes in daily climate extremes, it is necessary for the data to be homogenized because of nonclimatic influences. Here a new method of homogenizing daily temperature data is presented that is capable of adjusting not only the mean of a daily temperature series but also the higher-order moments. This method uses a nonlinear model to estimate the relationship between a candidate station and a highly correlated reference station. The model is built in a homogeneous subperiod before an inhomogeneity and is then used to estimate the observations at the candidate station after the inhomogeneity using observations from the reference series. The differences between the predicted and observed values are binned according to which decile the predicted values fit in the candidate station’s observed cumulative distribution function defined using homogeneous daily temperatures before the inhomogeneity. In this way, adjustments for each decile were produced. This method is demonstrated using February daily maximum temperatures measured in Graz, Austria, and an artificial dataset with known inhomogeneities introduced. Results show that given a suitably reliable reference station, this method produces reliable adjustments to the mean, variance, and skewness.

Description: A steady decline in Arctic sea ice has been observed over recent decades. General circulation models predict further decreases under increasing greenhouse gas scenarios. Sea ice plays an important role in the climate system in that it influences ocean-to-atmosphere fluxes, surface albedo, and ocean buoyancy. The aim of this study is to isolate the climate impacts of a declining Arctic sea ice cover during the current century. The Hadley Centre Atmospheric Model (HadAM3) is forced with observed sea ice from 1980 to 2000 (obtained from satellite passive microwave radiometer data derived with the Bootstrap algorithm) and predicted sea ice reductions until 2100 under one moderate scenario and one severe scenario of ice decline, with a climatological SST field and increasing SSTs. Significant warming of the Arctic occurs during the twenty-first century (mean increase of between 1.6° and 3.9°C), with positive anomalies of up to 22°C locally. The majority of this is over ocean and limited to high latitudes, in contrast to recent observations of Northern Hemisphere warming. When a climatological SST field is used, statistically significant impacts on climate are only seen in winter, despite prescribing sea ice reductions in all months. When correspondingly increasing SSTs are incorporated, changes in climate are seen in both winter and summer, although the impacts in summer are much smaller. Alterations in atmospheric circulation and precipitation patterns are more widespread than temperature, extending down to midlatitude storm tracks. Results suggest that areas of Arctic land ice may even undergo net accumulation due to increased precipitation that results from loss of sea ice. Intensification of storm tracks implies that parts of Europe may experience higher precipitation rates.

Description: A weeklong workshop in Brazil in August 2004 provided the opportunity for 28 scientists from southern South America to examine daily rainfall observations to determine changes in both total and extreme rainfall. Twelve annual indices of daily rainfall were calculated over the period 1960 to 2000, examining changes to both the entire distribution as well as the extremes. Maps of trends in the 12 rainfall indices showed large regions of coherent change, with many stations showing statistically significant changes in some of the indices. The pattern of trends for the extremes was generally the same as that for total annual rainfall, with a change to wetter conditions in Ecuador and northern Peru and the region of southern Brazil, Paraguay, Uruguay, and northern and central Argentina. A decrease was observed in southern Peru and southern Chile, with the latter showing significant decreases in many indices. A canonical correlation analysis between each of the indices and sea surface temperatures (SSTs) revealed two large-scale patterns that have contributed to the observed trends in the rainfall indices. A coupled pattern with ENSO-like SST loadings and rainfall loadings showing similarities with the pattern of the observed trend reveals that the change to a generally more negative Southern Oscillation index (SOI) has had an important effect on regional rainfall trends. A significant decrease in many of the rainfall indices at several stations in southern Chile and Argentina can be explained by a canonical pattern reflecting a weakening of the continental trough leading to a southward shift in storm tracks. This latter signal is a change that has been seen at similar latitudes in other parts of the Southern Hemisphere. A similar analysis was carried out for eastern Brazil using gridded indices calculated from 354 stations from the Global Historical Climatology Network (GHCN) database. The observed trend toward wetter conditions in the southwest and drier conditions in the northeast could again be explained by changes in ENSO.

Description: Outputs from two ensembles of atmospheric model simulations for 1951–98 define the influence of “realistic” land surface wetness on seasonal precipitation predictability in boreal summer. The ensembles consist of one forced with observed sea surface temperatures (SSTs) and the other forced with realistic land surface wetness as well as SSTs. Predictability was determined from correlations between the time series of simulated and observed precipitation. The ratio of forced variance to total variance determined potential predictability. Predictability occurred over some land areas adjacent to tropical oceans without land wetness forcing. On the other hand, because of the chaotic nature of the atmosphere, considerable parts of the land areas of the globe did not even show potential predictability with both land wetness and SST forcings. The use of land wetness forcing enhanced predictability over semiarid regions. Such semiarid regions are generally characterized by a negative correlation between fluxes of latent heat and sensible heat from the land surface, and are “water-regulating” areas where soil moisture plays a governing role in land–atmosphere interactions. Actual seasonal prediction may be possible in these regions if slowly varying surface conditions can be estimated in advance. In contrast, some land regions (e.g., south of the Sahel, the Amazon, and Indochina) showed little predictability despite high potential predictability. These regions are mostly characterized by a positive correlation between the surface fluxes, and are “radiation-regulating” areas where the atmosphere plays a leading role. Improvements in predictability for these regions may require further improvements in model physics.

Description: Eastward-propagating interdecadal time-scale sea surface temperature (SST) winter anomalies have been shown to exist at the North Atlantic subpolar/subtropical gyre boundary. Heat flux and surface air temperature signatures of these anomalies are investigated using satellite- and ship-based SST observations and atmospheric reanalysis. Using bandpass filter analysis, retaining periods between 9 and 25 yr, a succession of coherent propagating SST anomalies is identified. The size, speed, propagation path, and decay characteristics of propagating anomalies detected during the period 1948–2002 are documented. The behavior of the propagations changes between the periods 1948–70 and 1970–2002. In the former period, SST anomalies propagated from the east coast of North America to the British Isles in ∼10 yr. The anomalies displayed a well-defined life cycle, growing in the western basin (west of 40°W) and decaying in the eastern basin. During the period 1970–2002, SST anomalies did not propagate deep into the eastern basin, but grew in the western basin and then ceased propagating. Oceanic anomalies have a comparable marked signature in surface sensible and latent heat fluxes and in surface air temperature. Winter surface heat flux anomalies act to amplify SST anomalies during the middle of their lifetimes, normally in the west-central Atlantic. At other times, heat flux anomalies are associated with decay of anomalies. Surface heat fluxes do not always act to cause propagation, and it is likely that other processes such as advection play a role in the propagation mechanism. North European winter surface air temperatures are raised or lowered by up to ±0.5°C over decadal time scales (∼1/3 of the total variation over the United Kingdom) when an SST anomaly reaches the eastern boundary. A variety of processes can cause SST variation on decadal time scales at the eastern boundary, but in the 1950s and 1960s the variability at these periods was the signature of features that had propagated across the Atlantic from the North American coast.

Description: The water resources of the western United States depend heavily on snowpack to store part of the wintertime precipitation into the drier summer months. A well-documented shift toward earlier runoff in recent decades has been attributed to 1) more precipitation falling as rain instead of snow and 2) earlier snowmelt. The present study addresses the former, documenting a regional trend toward smaller ratios of winter-total snowfall water equivalent (SFE) to winter-total precipitation (P) during the period 1949–2004. The trends toward reduced SFE are a response to warming across the region, with the most significant reductions occurring where winter wet-day minimum temperatures, averaged over the study period, were warmer than −5°C. Most SFE reductions were associated with winter wet-day temperature increases between 0° and +3°C over the study period. Warmings larger than this occurred mainly at sites where the mean temperatures were cool enough that the precipitation form was less susceptible to warming trends. The trends toward reduced SFE/P ratios were most pronounced in March regionwide and in January near the West Coast, corresponding to widespread warming in these months. While mean temperatures in March were sufficiently high to allow the warming trend to produce SFE/P declines across the study region, mean January temperatures were cooler, with the result that January SFE/P impacts were restricted to the lower elevations near the West Coast. Extending the analysis back to 1920 shows that although the trends presented here may be partially attributable to interdecadal climate variability associated with the Pacific decadal oscillation, they also appear to result from still longer-term climate shifts.

Description: The twentieth-century simulations using by 17 coupled ocean–atmosphere general circulation models (CGCMs) submitted to the Intergovernmental Panel on Climate Change’s Fourth Assessment Report (IPCC AR4) are evaluated for their skill in reproducing the observed modes of Indian Ocean (IO) climate variability. Most models successfully capture the IO’s delayed, basinwide warming response a few months after El Niño–Southern Oscillation (ENSO) peaks in the Pacific. ENSO’s oceanic teleconnection into the IO, by coastal waves through the Indonesian archipelago, is poorly simulated in these models, with significant shifts in the turning latitude of radiating Rossby waves. In observations, ENSO forces, by the atmospheric bridge mechanism, strong ocean Rossby waves that induce anomalies of SST, atmospheric convection, and tropical cyclones in a thermocline dome over the southwestern tropical IO. While the southwestern IO thermocline dome is simulated in nearly all of the models, this ocean Rossby wave response to ENSO is present only in a few of the models examined, suggesting difficulties in simulating ENSO’s teleconnection in surface wind. A majority of the models display an equatorial zonal mode of the Bjerknes feedback with spatial structures and seasonality similar to the Indian Ocean dipole (IOD) in observations. This success appears to be due to their skills in simulating the mean state of the equatorial IO. Corroborating the role of the Bjerknes feedback in the IOD, the thermocline depth, SST, precipitation, and zonal wind are mutually positively correlated in these models, as in observations. The IOD–ENSO correlation during boreal fall ranges from −0.43 to 0.74 in the different models, suggesting that ENSO is one, but not the only, trigger for the IOD.

Description: Using an optimal detection technique and climate change simulations produced with two versions of two GCMs, we have assessed the causes of twentieth-century temperature changes from global to regional scales. Our analysis is conducted in nine spatial domains: 1) the globe; 2) the Northern Hemisphere; four large regions in the Northern Hemispheric midlatitudes covering 30°–70°N including 3) Eurasia, 4) North America, 5) Northern Hemispheric land only, 6) the entire 30°–70°N belt; and three smaller regions over 7) southern Canada, 8) southern Europe, and 9) China. We find that the effect of anthropogenic forcing on climate is clearly detectable at global through regional scales. The effect of combined greenhouse gases and sulfate aerosol forcing is detectable in all nine domains in annual and seasonal mean temperatures observed during the second half of the twentieth century. The effect of greenhouse gases can also be separated from that of sulfate aerosols over this period at continental and regional scales. Uncertainty in these results is larger in the smaller spatial domains. Detection is improved when an ensemble of models is used to estimate the response to anthropogenic forcing and the underlying internal variability of the climate system. Our detection results hold after removal of North Atlantic Oscillation (NAO)-related variability in temperature observations—variability that may or may not be associated with anthropogenic forcing. They also continue to hold when our estimates of natural internal climate variability are doubled.

Description: This paper discusses a cloud-resolving modeling study concerning the impact of warm-rain microphysics on convective–radiative quasi equilibrium with fixed surface characteristics and prescribed solar input, both mimicking the mean conditions on earth. Two limits of the concentration of cloud droplets, either 100 cm−3 (referred to as “pristine”) or 1000 cm−3 (referred to as “polluted”), are considered. In addition, three formulations of the effective radius of water droplets in diluted cloudy volumes are used, corresponding to the homogeneous, intermediate, and extremely inhomogeneous mixing scenarios. The assumed concentration of cloud droplets, together with the assumed mixing scenario, affects the local value of the effective radius of cloud droplets (the first indirect aerosol effect, also known as the Twomey effect) and the transfer of cloud water into drizzle and rain, which can affect the mean cloudiness and the hydrologic cycle (the second indirect effect). The convective–radiative quasi equilibrium mimics the estimates of globally and annually averaged water and energy fluxes across the earth’s atmosphere to within less than 10 W m−2. As on earth, the model cloudiness is dominated by shallow convection. It is found that the impact of warm microphysics is dominated by the first indirect effect, whereas the second indirect effect has a smaller impact. The assumed droplet concentration and mixing scenario impact the mean “planetary” albedo and, thus, the amount of solar energy reaching the surface, with all other components of atmospheric energy and water budgets virtually the same in all simulations. The weak second indirect effect highlights the difference between the impact of cloud microphysics on a single cloud and the impact on an ensemble of clouds, with only the latter including the feedbacks between clouds and their environment. The formulation of the effective radius in the diluted cloudy volumes turns out to be of critical importance, with the amount of solar energy reaching the surface being the same in the pristine case assuming the homogeneous mixing scenario and in the polluted case with the extremely inhomogeneous mixing. This result emphasizes the essential role of poorly understood microphysical transformations within diluted convective clouds, which strongly impact the magnitude of the first indirect (Twomey) effect. Implications for future research in this area are discussed.

Description: The analyses presented here focus on the Southern Ocean as simulated in a set of global coupled climate model control experiments conducted by several international climate modeling groups. Dominated by the Antarctic Circumpolar Current (ACC), the vast Southern Ocean can influence large-scale surface climate features on various time scales. Its climatic relevance stems in part from it being the region where most of the transformation of the World Ocean’s water masses occurs. In climate change experiments that simulate greenhouse gas–induced warming, Southern Ocean air–sea heat fluxes and three-dimensional circulation patterns make it a region where much of the future oceanic heat storage takes place, though the magnitude of that heat storage is one of the larger sources of uncertainty associated with the transient climate response in such model projections. Factors such as the Southern Ocean’s wind forcing, heat, and salt budgets are linked to the structure and transport of the ACC in ways that have not been expressed clearly in the literature. These links are explored here in a coupled model context by analyzing a sizable suite of preindustrial control experiments associated with the forthcoming Intergovernmental Panel on Climate Change’s Fourth Assessment Report. A framework is developed that uses measures of coupled model simulation characteristics, primarily those related to the Southern Ocean wind forcing and water mass properties, to allow one to categorize, and to some extent predict, which models do better or worse at simulating the Southern Ocean and why. Hopefully, this framework will also lead to increased understanding of the ocean’s response to climate changes.

Description: The commonly held view of the conditions in the North Atlantic at the last glacial maximum, based on the interpretation of proxy records, is of large-scale cooling compared to today, limited deep convection, and extensive sea ice, all associated with a southward displaced and weakened overturning thermohaline circulation (THC) in the North Atlantic. Not all studies support that view; in particular, the “strength of the overturning circulation” is contentious and is a quantity that is difficult to determine even for the present day. Quasi-equilibrium simulations with coupled climate models forced by glacial boundary conditions have produced differing results, as have inferences made from proxy records. Most studies suggest the weaker circulation, some suggest little or no change, and a few suggest a stronger circulation. Here results are presented from a three-dimensional climate model, the Hadley Centre Coupled Model version 3 (HadCM3), of the coupled atmosphere–ocean–sea ice system suggesting, in a qualitative sense, that these diverging views could all have occurred at different times during the last glacial period, with different modes existing at different times. One mode might have been characterized by an active THC associated with moderate temperatures in the North Atlantic and a modest expanse of sea ice. The other mode, perhaps forced by large inputs of meltwater from the continental ice sheets into the northern North Atlantic, might have been characterized by a sluggish THC associated with very cold conditions around the North Atlantic and a large areal cover of sea ice. The authors’ model simulation of such a mode, forced by a large input of freshwater, bears several of the characteristics of the Climate: Long-range Investigation, Mapping, and Prediction (CLIMAP) Project’s reconstruction of glacial sea surface temperature and sea ice extent.

Description: The primary regimes of local atmospheric variability are examined in a 6-km regional atmospheric model of the southern third of California, an area of significant land surface heterogeneity, intense topography, and climate diversity. The model was forced by reanalysis boundary conditions over the period 1995–2003. The region is approximately the same size as a typical grid box of the current generation of general circulation models used for global climate prediction and reanalysis product generation, and so can be thought of as a laboratory for the study of climate at spatial scales smaller than those resolved by global simulations and reanalysis products. It is found that the simulated circulation during the October–March wet season, when variability is most significant, can be understood through an objective classification technique in terms of three wind regimes. The composite surface wind patterns associated with these regimes exhibit significant spatial structure within the model domain, consistent with the complex topography of the region. These regimes also correspond nearly perfectly with the simulation’s highly structured patterns of variability in hydrology and temperature, and therefore are the main contributors to the local climate variability. The regimes are approximately equally likely to occur regardless of the phase of the classical large-scale modes of atmospheric variability prevailing in the Pacific–North American sector. The high degree of spatial structure of the local regimes and their tightly associated climate impacts, as well as their ambiguous relationship with the primary modes of large-scale variability, demonstrate that the local perspective offered by the high-resolution model is necessary to understand and predict the climate variations of the region.

Description: In recent years, there has been an increase in rainfall over northeast Brazil (Nordeste), while over the sub-Saharan region there has been a drought. The correlation coefficients between the 11-yr running means of the rainfall series over the two regions are significant (at the 95% confidence level by a two-sided t test), suggesting that both trends are related. The rainfall variations over the two regions are connected to the position of the intertropical convergence zone (ITCZ) over the Atlantic. A more southward (northward) position of the ITCZ is favorable for higher than normal rainfall over Nordeste (sub-Sahara). The correlation coefficient between the position of the ITCZ over the Atlantic and the rainfall over Nordeste (sub-Sahara) is negative (positive) and highly significant, reaching values over 0.9. Thus, this study suggests that a more southward than normal location of the ITCZ in the Atlantic may be the cause for the recent increasing (decreasing) trend of rainfall over Nordeste (sub-Sahara).

Description: This study assessed the climate and trend of cyclone activity in Canada using mainly the occurrence frequency of cyclone deepening events and deepening rates, which were derived from hourly mean sea level pressure data observed at 83 Canadian stations for up to 50 years (1953–2002). Trends in the frequency of cyclone activity were estimated by logistic regression analysis, and trends of seasonal extreme cyclone intensity, by linear regression analysis. The results of trend analysis show that, among the four seasons, winter cyclone activity has shown the most significant trends. It has become significantly more frequent, more durable, and stronger in the lower Canadian Arctic, but less frequent and weaker in the south, especially along the southeast and southwest coasts. Winter cyclone deepening rates have increased in the zone around 60°N but decreased in the Great Lakes area and southern Prairies–British Columbia. However, extreme winter cyclone activity seems to have experienced a weaker increase in northwest-central Canada but a stronger decline in the Great Lakes area and in southern Prairies. The results also show more frequent summer cyclone activity with slower deepening rates on the east coast, as well as less frequent cyclone activity with faster deepening rates in the Great Lakes area in autumn. Cyclone activity in Canada was found to be closely related to the North Atlantic Oscillation (NAO), the Pacific Decadal Oscillation (PDO), and El Niño–Southern Oscillation (ENSO). Overall, cyclone activity in Canada is most closely related to the NAO. The simultaneous NAO index explains about 44% (41%) of the winter (autumn) cyclone activity variance in the east coast, 31% of winter cyclone activity variance in the 60°–70°N zone, and 17% of autumn cyclone activity variance in the Great Lakes area. Also, in several regions (e.g., the east coast, the southwest, and the 60°–70°N zone) up to 15% of the seasonal cyclone activity variance can be explained by the NAO/PDO/ENSO index one–three seasons earlier, which is useful for seasonal forecasting.

Description: Tropospheric temperature trends based on Microwave Sounding Unit (MSU) channel 2 data are susceptible to contamination from strong stratospheric cooling. Recently, Fu et al. devised a method of removing the stratospheric contamination by linearly combining data from MSU channels 2 and 4. In this study the sensitivity of the weights of the two channels in the retrieval algorithm for the tropospheric temperatures to the choice of period of record used in the analysis and to the choice of training dataset is examined. The weights derived using monthly temperature anomalies are within about 10% of those obtained by Fu et al. irrespective of the choice of analysis period or training dataset. The trend errors in the retrieved global-mean tropospheric temperatures tested using two independent radiosonde datasets are less than about 0.01 K decade−1 for all time periods of 25 yr or longer with different starting and ending years during 1958–2004. It is found that the retrievals are more robust if they are interpreted in terms of the layer-mean temperature for the entire troposphere, rather than the mean of the 850–300-hPa layer. Because large spurious jumps remain in the reanalyses, especially prior to 1979, one should be cautious when using them as training datasets and in testing the trend errors.

Description: A midlatitude coupled ocean–atmosphere model is used to investigate interactions between the atmosphere and the wind-driven ocean circulation. This model uses idealized geometry, yet rich and complicated dynamic flow regimes arise in the ocean due to the explicit simulation of geostrophic turbulence. An interdecadal mode of intrinsic ocean variability is found, and this mode projects onto existing atmospheric modes of variability, thereby controlling the time scale of the atmospheric modes. It is also shown that ocean circulation controls the time scale of the SST response to wind forcing, and that coupled feedback mechanisms thus modify variability of the atmospheric circulation. It is concluded that ocean–atmosphere coupling in the midlatitudes is unlikely to produce new modes of variability but may control the temporal behavior of modes that exist in uncoupled systems.

Description: Using the adjoint of a fully three-dimensional primitive equation ocean model in an idealized geometry, spatial variations in the sensitivity to surface boundary forcing of the meridional overturning circulation’s strength are studied. Steady-state sensitivities to diapycnal mixing, wind stress, freshwater, and heat forcing are examined. Three different, commonly used, boundary-forcing scenarios are studied, both with and without wind forcing. Almost identical circulation is achieved in each scenario, but the sensitivity patterns show major (quantitative and qualitative) differences. Sensitivities to surface forcing and diapycnal mixing are substantially larger under mixed boundary conditions, in which fluxes of freshwater and heat are supplemented by a temperature relaxation term or under flux boundary conditions, in which climatological fluxes alone drive the circulation, than under restoring boundary conditions. The sensitivity pattern to diapycnal mixing, which peaks in the Tropics is similar both with and without wind forcing. Wind does, however, increase the sensitivity to diapycnal mixing in the regions of Ekman upwelling and decreases it in the regions of Ekman downwelling. Wind stress in the Southern Oceans plays a crucial role in restoring boundary conditions, but the effect is largely absent under mixed or flux boundary conditions. The results highlight how critical a careful formulation of the surface forcing terms is to ensuring a proper response to changes in forcing in ocean models.

Description: This paper describes the mean ocean circulation and the tropical variability simulated by the Max Planck Institute for Meteorology (MPI-M) coupled atmosphere–ocean general circulation model (AOGCM). Results are presented from a version of the coupled model that served as a prototype for the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) simulations. The model does not require flux adjustment to maintain a stable climate. A control simulation with present-day greenhouse gases is analyzed, and the simulation of key oceanic features, such as sea surface temperatures (SSTs), large-scale circulation, meridional heat and freshwater transports, and sea ice are compared with observations. A parameterization that accounts for the effect of ocean currents on surface wind stress is implemented in the model. The largest impact of this parameterization is in the tropical Pacific, where the mean state is significantly improved: the strength of the trade winds and the associated equatorial upwelling weaken, and there is a reduction of the model’s equatorial cold SST bias by more than 1 K. Equatorial SST variability also becomes more realistic. The strength of the variability is reduced by about 30% in the eastern equatorial Pacific and the extension of SST variability into the warm pool is significantly reduced. The dominant El Niño–Southern Oscillation (ENSO) period shifts from 3 to 4 yr. Without the parameterization an unrealistically strong westward propagation of SST anomalies is simulated. The reasons for the changes in variability are linked to changes in both the mean state and to a reduction in atmospheric sensitivity to SST changes and oceanic sensitivity to wind anomalies.

Description: The spatial patterns, time history, and seasonality of African rainfall trends since 1950 are found to be deducible from the atmosphere’s response to the known variations of global sea surface temperatures (SSTs). The robustness of the oceanic impact is confirmed through the diagnosis of 80 separate 50-yr climate simulations across a suite of atmospheric general circulation models. Drying over the Sahel during boreal summer is shown to be a response to warming of the South Atlantic relative to North Atlantic SST, with the ensuing anomalous interhemispheric SST contrast favoring a more southern position of the Atlantic intertropical convergence zone. Southern African drying during austral summer is shown to be a response to Indian Ocean warming, with enhanced atmospheric convection over those warm waters driving subsidence drying over Africa. The ensemble of greenhouse-gas-forced experiments, conducted as part of the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, fails to simulate the pattern or amplitude of the twentieth-century African drying, indicating that the drought conditions were likely of natural origin. For the period 2000–49, the ensemble mean of the forced experiments yields a wet signal over the Sahel and a dry signal over southern Africa. These rainfall changes are physically consistent with a projected warming of the North Atlantic Ocean compared with the South Atlantic Ocean, and a further warming of the Indian Ocean. However, considerable spread exists among the individual members of the multimodel ensemble.

Description: A neural network technique is used to quantify relationships involved in cloud–radiation feedbacks based on observations from the Surface Heat Budget of the Arctic (SHEBA) project. Sensitivities of longwave cloud forcing (CFL) to cloud parameters indicate that a bimodal distribution pattern dominates the histogram of each sensitivity. Although the mean states of the relationships agree well with those derived in a previous study, they do not often exist in reality. The sensitivity of CFL to cloud cover increases as the cloudiness increases with a range of 0.1–0.9 W m−2 %−1. There is a saturation effect of liquid water path (LWP) on CFL. The highest sensitivity of CFL to LWP corresponds to clouds with low LWP, and sensitivity decreases as LWP increases. The sensitivity of CFL to cloud-base height (CBH) depends on whether the clouds are below or above an inversion layer. The relationship is negative for clouds higher than 0.8 km at the SHEBA site. The strongest positive relationship corresponds to clouds with low CBH. The dominant mode of the sensitivity of CFL to cloud-base temperature (CBT) is near zero and corresponds to warm clouds with base temperatures higher than −9°C. The low and high sensitivity regimes correspond to the summer and winter seasons, respectively, especially for LWP and CBT. Overall, the neural network technique is able to separate two distinct regimes of clouds that correspond to different sensitivities; that is, it captures the nonlinear behavior in the relationships. This study demonstrates a new method for evaluating nonlinear relationships between climate variables. It could also be used as an effective tool for evaluating feedback processes in climate models.

Description: The seasonal and interannual predictability of ENSO variability in a version of the Zebiak–Cane coupled model is examined in a perturbation experiment. Instead of assuming that the model is “perfect,” it is assumed that a set of optimal initial conditions exists for the model. These states, obtained through a nonlinear minimization of the misfit between model trajectories and the observations, initiate model forecasts that correlate well with the observations. Realistic estimates of the observational error magnitudes and covariance structures of sea surface temperatures, zonal wind stress, and thermocline depth are used to generate ensembles of perturbations around these optimal initial states, and the error growth is examined. The error growth in response to subseasonal stochastic wind forcing is presented for comparison. In general, from 1975 to 2002, the large-scale uncertainty in initial conditions leads to larger error growth than continuous stochastic forcing of the zonal wind stress fields. Forecast ensemble spread is shown to depend most on the calendar month at the end of the forecast rather than the initialization month, with the seasons of greatest spread corresponding to the seasons of greatest anomaly variance. It is also demonstrated that during years with negative (and rapidly decaying) Niño-3 SST anomalies (such as the time period following an El Niño event), there is a suppression of error growth. In years with large warm ENSO events, the ensemble spread is no larger than in moderately warm years. As a result, periods with high ENSO variance have greater potential prediction utility. In the realistic range of observational error, the ensemble spread has more sensitivity to the initial error in the thermocline depth than to the sea surface temperature or wind stress errors. The thermocline depth uncertainty is the principal reason why initial condition uncertainties are more important than wind noise for ensemble spread.

Description: Multicentury sensitivities in a realistic geometry global ocean general circulation model are analyzed using an adjoint technique. This paper takes advantage of the adjoint model’s ability to generate maps of the sensitivity of a diagnostic (i.e., the meridional overturning’s strength) to all model parameters. This property of adjoints is used to review several theories, which have been elaborated to explain the strength of the North Atlantic’s meridional overturning. This paper demonstrates the profound impact of boundary conditions in permitting or suppressing mechanisms within a realistic model of the contemporary ocean circulation. For example, the so-called Drake Passage Effect in which wind stress in the Southern Ocean acts as the main driver of the overturning’s strength, is shown to be an artifact of boundary conditions that restore the ocean’s surface temperature and salinity toward prescribed climatologies. Advective transports from the Indian and Pacific basins play an important role in setting the strength of the overturning circulation under “mixed” boundary conditions, in which a flux of freshwater is specified at the ocean’s surface. The most “realistic” regime couples an atmospheric energy and moisture balance model to the ocean. In this configuration, inspection of the global maps of sensitivity to wind stress and diapycnal mixing suggests a significant role for near-surface Ekman processes in the Tropics. Buoyancy also plays an important role in setting the overturning’s strength, through direct thermal forcing near the sites of convection, or through the advection of salinity anomalies in the Atlantic basin.

Description: The climate feedbacks in coupled ocean–atmosphere models are compared using a coordinated set of twenty-first-century climate change experiments. Water vapor is found to provide the largest positive feedback in all models and its strength is consistent with that expected from constant relative humidity changes in the water vapor mixing ratio. The feedbacks from clouds and surface albedo are also found to be positive in all models, while the only stabilizing (negative) feedback comes from the temperature response. Large intermodel differences in the lapse rate feedback are observed and shown to be associated with differing regional patterns of surface warming. Consistent with previous studies, it is found that the vertical changes in temperature and water vapor are tightly coupled in all models and, importantly, demonstrate that intermodel differences in the sum of lapse rate and water vapor feedbacks are small. In contrast, intermodel differences in cloud feedback are found to provide the largest source of uncertainty in current predictions of climate sensitivity.

Description: The equatorial Pacific is a region with strong negative feedbacks. Yet coupled general circulation models (GCMs) have exhibited a propensity to develop a significant SST bias in that region, suggesting an unrealistic sensitivity in the coupled models to small energy flux errors that inevitably occur in the individual model components. Could this “hypersensitivity” exhibited in a coupled model be due to an underestimate of the strength of the negative feedbacks in this region? With this suspicion, the feedbacks in the equatorial Pacific in nine atmospheric GCMs (AGCMs) have been quantified using the interannual variations in that region and compared with the corresponding calculations from the observations. The nine AGCMs are the NCAR Community Climate Model version 1 (CAM1), the NCAR Community Climate Model version 2 (CAM2), the NCAR Community Climate Model version 3 (CAM3), the NCAR CAM3 at T85 resolution, the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Atmospheric Model, the Hadley Centre Atmospheric Model (HadAM3), the Institut Pierre Simon Laplace (IPSL) model (LMDZ4), the Geophysical Fluid Dynamics Laboratory (GFDL) AM2p10, and the GFDL AM2p12. All the corresponding coupled runs of these nine AGCMs have an excessive cold tongue in the equatorial Pacific. The net atmospheric feedback over the equatorial Pacific in the two GFDL models is found to be comparable to the observed value. All other models are found to have a weaker negative net feedback from the atmosphere—a weaker regulating effect on the underlying SST than the real atmosphere. Except for the French (IPSL) model, a weaker negative feedback from the cloud albedo and a weaker negative feedback from the atmospheric transport are the two leading contributors to the weaker regulating effect from the atmosphere. The underestimate of the strength of the negative feedbacks by the models is apparently linked to an underestimate of the equatorial precipitation response. All models have a stronger water vapor feedback than that indicated in Earth Radiation Budget Experiment (ERBE) observations. These results confirm the suspicion that an underestimate of the regulatory effect from the atmosphere over the equatorial Pacific region is a prevalent problem. The results also suggest, however, that a weaker regulatory effect from the atmosphere is unlikely solely responsible for the hypersensitivity in all models. The need to validate the feedbacks from the ocean transport is therefore highlighted.

Description: Observations of air and ground temperatures collected between 1993 and 2004 from Emigrant Pass Geothermal Climate Observatory in northwestern Utah are analyzed to understand the relationship between these two quantities. The influence of surface air temperature (SAT), incident solar radiation, and snow cover on surface ground temperature (SGT) variations are explored. SAT variations explain 94% of the variance in SGT. Incident solar radiation is the primary variable governing the remaining variance misfit and is significantly more important during summer months than winter months. A linear relationship between the ground–air temperature difference (ΔTsgt-sat) and solar radiation exists with a trend of 1.21 K/(100 W m−2); solar radiation accounts for 1.3% of the variance in SGT. The effects of incident solar radiation also account for the 2.47-K average offset in ΔTsgt-sat. During the winter, snow cover plays a role in governing SGT variability, but exerts only a minor influence on the annual tracking of ground and air temperatures at the site, accounting for 0.5% of the variance in SGT. These observations of the tracking of SGT and SAT confirm that borehole temperature changes mimic changes in SAT at frequencies appropriate for climatic reconstructions.

Description: Seasonal and interannual variation of rainfall over Mongolia was investigated using 10-day rainfall data from 92 stations during 1993–2001 and NCEP–NCAR reanalysis data from 1979 to 2001. The break in the rainy season was found in the middle of July, and the meteorological stations with a clear break period were concentrated in eastern Mongolia where the plain prevails relative to western Mongolia without regard to the difference in annual precipitation. Clear breaks in the rainy season were recognized in 5 yr among an analysis period of 9 yr. In the break period, the stationary Rossby wave trapped in the Asian jet was predominant at 200 hPa, and a barotropic ridge associated with the Rossby wave developed over Mongolia. Furthermore, interannual variation of the break also corresponded to the variation of the stationary Rossby wave. It is considered that the break of the Mongolian rainy season is caused by the stationary Rossby wave trapped in the Asian jet. The stationary Rossby wave was climatologically phase locked in seasonal evolution and, as a result, the break period was also concentrated around the middle of July.

Description: A climatology of nonfreezing drizzle is created using surface observations from 584 stations across the United States and Canada over the 15-yr period 1976–90. Drizzle falls 50–200 h a year in most locations in the eastern United States and Canada, whereas drizzle falls less than 50 h a year in the west, except for coastal Alaska and several western basins. The eastern and western halves of North America are separated by a strong gradient in drizzle frequency along roughly 100°W, as large as about an hour a year over 2 km. Forty percent of the stations have a drizzle maximum from November to January, whereas only 13% of stations have a drizzle maximum from June to August. Drizzle occurrence exhibits a seasonal migration from eastern Canada and the central portion of the Northwest Territories in summer, equatorward to most of the eastern United States and southeast Canada in early winter, to southeastern Texas and the eastern United States in late winter, and back north to eastern Canada in the spring. The diurnal hourly frequency of drizzle across the United States and Canada increases sharply from 0900 to 1200 UTC, followed by a steady decline from 1300 to 2300 UTC. Diurnal drizzle frequency is at a maximum in the early morning, in agreement with other studies. Drizzle occurs during a wide range of atmospheric conditions at the surface. Drizzle has occurred at sea level pressures below 960 hPa and above 1040 hPa. Most drizzle, however, occurs at higher than normal sea level pressure, with more than 64% occurring at a sea level pressure of 1015 hPa or higher. A third of all drizzle falls when the winds are from the northeast quadrant (360°–89°), suggesting that continental drizzle events tend to be found poleward of surface warm fronts and equatorward of cold-sector surface anticyclones. Two-thirds of all drizzle occurs with wind speeds of 2.0–6.9 m s−1, with 7.6% in calm wind and 5% at wind speeds ⩾ 10 m s−1. Most drizzle (61%) occurs with visibilities between 1.5 and 5.0 km, with only about 20% occurring at visibilities less than 1.5 km.

Description: A global chemical transport model of the atmosphere [the Model for Ozone and Related Tracers, version 2 (MOZART-2)] driven by prescribed surface emissions and by meteorological fields provided by the ECHAM5/Max Planck Institute Ocean Model (MPI-OM-1) coupled atmosphere–ocean model is used to assess how expected climate changes (2100 versus 2000 periods) should affect the chemical composition of the troposphere. Calculations suggest that ozone changes resulting from climate change only are negative in a large fraction of the troposphere because of enhanced photochemical destruction by water vapor. In the Tropics, increased lightning activity should lead to larger ozone concentrations. The magnitude of the climate-induced ozone changes in the troposphere remains smaller than the changes produced by enhanced anthropogenic emissions when the Special Report on Emission Scenarios (SRES) A2P is adopted to describe the future evolution of these emissions. Predictions depend strongly on future trends in atmospheric methane levels, which are not well established. Changes in the emissions of NOx by bacteria in soils and of nonmethane hydrocarbons by vegetation associated with climate change could have a significant impact on future ozone levels.

Description: Extratropical and tropical transient storm tracks are investigated from the perspective of feature tracking in the ECHAM5 coupled climate model for the current and a future climate scenario. The atmosphere-only part of the model, forced by observed boundary conditions, produces results that agree well with analyses from the 40-yr ECMWF Re-Analysis (ERA-40), including the distribution of storms as a function of maximum intensity. This provides the authors with confidence in the use of the model for the climate change experiments. The statistical distribution of storm intensities is virtually preserved under climate change using the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A1B scenario until the end of this century. There are no indications in this study of more intense storms in the future climate, either in the Tropics or extratropics, but rather a minor reduction in the number of weaker storms. However, significant changes occur on a regional basis in the location and intensity of storm tracks. There is a clear poleward shift in the Southern Hemisphere with consequences of reduced precipitation for several areas, including southern Australia. Changes in the Northern Hemisphere are less distinct, but there are also indications of a poleward shift, a weakening of the Mediterranean storm track, and a strengthening of the storm track north of the British Isles. The tropical storm tracks undergo considerable changes including a weakening in the Atlantic sector and a strengthening and equatorward shift in the eastern Pacific. It is suggested that some of the changes, in particular the tropical ones, are due to an SST warming maximum in the eastern Pacific. The shift in the extratropical storm tracks is shown to be associated with changes in the zonal SST gradient in particular for the Southern Hemisphere.

Description: The individual impacts of sea surface temperature (SST) anomalies in the deep tropical eastern–central Pacific (DTEP) and Indo-western–central Pacific (IWP) on the evolution of the observed global atmospheric circulation during the 1997–2003 period have been investigated using a new general circulation model. Ensemble integrations were conducted with monthly varying SST conditions being prescribed separately in the DTEP sector, the IWP sector, and throughout the World Ocean. During the 1998–2002 subperiod, when prolonged La Niña conditions occurred in DTEP and the SST in IWP was above normal, the simulated midlatitude atmospheric responses to SST forcing in the DTEP and IWP sectors reinforced each other. The anomalous geopotential height ridges at 200 mb in the extratropics of both hemispheres exhibited a distinct zonal symmetry. This circulation change was accompanied by extensive dry and warm anomalies in many regions, including North America. During the 1997–98 and 2002–03 El Niño events, the SST conditions in both DTEP and IWP were above normal, and considerable cancellations were simulated between the midlatitude responses to the oceanic forcing from these two sectors. The above findings are contrasted with those for the 1953–58 and 1972–77 periods, which were characterized by analogous SST developments in DTEP, but by cold conditions in IWP. It is concluded that a warm anomaly in IWP and a cold anomaly in DTEP constitute the optimal SST configuration for generating zonally elongated ridges in the midlatitudes. Local diagnoses indicate that the imposed SST anomaly alters the strength of the zonal flow in certain longitudinal sectors, which influences the behavior of synoptic-scale transient eddies farther downstream. The modified eddy momentum transports in the regions of eddy activity in turn feed back on the local mean flow, thus contributing to its zonal elongation. These results are consistent with the inferences drawn from zonal mean analyses, which accentuate the role of the eddy-induced circulation on the meridional plane.

Description: The ability of coupled GCMs to correctly simulate the climatology and a prominent mode of variability of the West African monsoon is evaluated, and the results are used to make informed decisions about which models may be producing more reliable projections of future climate in this region. The integrations were made available by the Program for Climate Model Diagnosis and Intercomparison for the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. The evaluation emphasizes the circulation characteristics that support the precipitation climatology, and the physical processes of a “rainfall dipole” variability mode that is often associated with dry conditions in the Sahel when SSTs in the Gulf of Guinea are anomalously warm. Based on the quality of their twentieth-century simulations over West Africa in summer, three GCMs are chosen for analysis of the twenty-first century integrations under various assumptions about future greenhouse gas increases. Each of these models behaves differently in the twenty-first-century simulations. One model simulates severe drying across the Sahel in the later part of the twenty-first century, while another projects quite wet conditions throughout the twenty-first century. In the third model, warming in the Gulf of Guinea leads to more modest drying in the Sahel due to a doubling of the number of anomalously dry years by the end of the century. An evaluation of the physical processes that cause these climate changes, in the context of the understanding about how the system works in the twentieth century, suggests that the third model provides the most reasonable projection of the twenty-first-century climate.

Description: Paleoevidence indicates that generally wetter conditions existed in the Sahara during the mid-Holocene. Climate modeling studies addressing this issue generally agree that mid-Holocene values of the earth’s orbital parameters favored an enhanced North African summer monsoon but also suggest that land surface and vegetation feedbacks must have been important factors. Attempts to reproduce the “green” mid-Holocene Sahara in model studies with interactive vegetation may be interpreted to indicate that the problem is highly sensitive to the atmospheric dynamics of each model employed. In other work, dynamical mechanisms have been hypothesized to affect monsoon poleward extent, particularly ventilation, by import of low-moist static energy air to the continent. Here, interactive vegetation and the ventilation mechanism are studied in an intermediate complexity atmospheric model coupled to simple land and vegetation components. Interactive vegetation is found to be effective at enhancing the precipitation and vegetation amount in regions where the monsoon has advanced because of changes in orbital parameters or ventilation yet not very effective in moving the monsoon boundary if ventilation is strong. The poleward extent of the mid-Holocene monsoon and the steppe boundary are primarily controlled by the strength of ventilation in the atmospheric model. Within this boundary, the largest changes in monsoon precipitation and vegetation occur when interactive vegetation and reduced ventilation act simultaneously, as these greatly reinforce each other.

Description: An EOF analysis is used to intercompare the response of ENSO-like variability to CO2 doubling in results from 15 coupled climate models assembled for the Intergovernmental Panel on Climate Change Fourth Assessment Report. Under preindustrial conditions, 12 of the 15 models exhibit ENSO amplitudes comparable to or exceeding that observed in the second half of the twentieth century. Under CO2 doubling, three of the models exhibit statistically significant (p 〈 0.1) increases in ENSO amplitude, and five exhibit significant decreases. The overall amplitude changes are not strongly related to the magnitude or pattern of surface warming. It is, however, found that ENSO amplitude decreases (increases) in models having a narrow (wide) ENSO zonal wind stress response and ENSO amplitude comparable to or greater than observed. The models exhibit a mean fractional decrease in ENSO period of about 5%. Although many factors can influence the ENSO period, it is suggested that this may be related to a comparable increase in equatorial wave speed through an associated speedup of delayed-oscillator feedback. Changes in leading EOF, characterized in many of the models by a relative increase in the amplitude of SST variations in the central Pacific, are in most cases consistent with effects of anomalous zonal and vertical advection resulting from warming-induced changes in SST.

Description: A recently developed variance decomposition approach is applied to study the causes of the predictability of New Zealand seasonal mean rainfall. In terms of predictability, the Southern Oscillation is identified as being the most important cause of variability for both the winter and summer New Zealand rainfall, especially for the North Island. Indian Ocean sea surface temperature variability and the Southern Hemisphere annular mode are the second most important causes of variability for winter and summer rainfall, respectively. Based on this study, a statistical prediction scheme has been developed. May Niño-3 (5°N–5°S, 150°–90°W) SSTs and March–May (MAM) central Indian Ocean SSTs are identified as being the most important predictors for the winter rainfall, while September–November (SON) Niño-3 SSTs, November local New Zealand SSTs, and the SON Southern Hemisphere annular mode index are the most important predictors for the summer rainfall. The predictive skill, in term of the percentage explained variance for the verification period (1993–2000) is nearly 20%, which is considerably higher than that achieved previously.

Description: A Bayesian approach is applied to the observed global surface air temperature (SAT) changes using multimodel ensembles (MMEs) of the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) simulations and single-model ensembles (SMEs) with the ECHO-G coupled climate model. A Bayesian decision method is used as a tool for classifying observations into given scenarios (or hypotheses). The prior probability of the scenarios, which represents a degree of subjective belief in the scenarios, is changed into the posterior probability through the likelihood where observations enter, and the posterior is used as a decision function. In the identical prior case the Bayes factor (or likelihood ratio) becomes a decision function and provides observational evidence for each scenario against a predefined reference scenario. Four scenarios are used to explain observed SAT changes: “CTL” (control or no change), “Nat” (natural forcing induced change), “GHG” (greenhouse gas–induced change), and “All” (natural plus anthropogenic forcing–induced change). Observed and simulated global mean SATs are decomposed into temporal components of overall mean, linear trend, and decadal variabilities through Legendre series expansions, coefficients of which are used as detection variables. Parameters (means and covariance matrices) needed to define the four scenarios are estimated from SMEs or MMEs. Taking the CTL scenario as reference one, application results for global mean SAT changes for the whole twentieth century (1900–99) show “decisive” evidence (logarithm of Bayes factor 〉5) for the All scenario only. While “strong” evidence (log of Bayes factor 〉2.5) for both the Nat and All scenarios are found in SAT changes for the first half (1900–49), there is decisive evidence for the All scenario for SAT changes in the second half (1950–99), supporting previous results. It is demonstrated that the Bayesian decision results for global mean SATs are largely insensitive to both intermodel uncertainties and prior probabilities.

Description: An ensemble of general circulation model (GCM) integrations forced by observed sea surface temperature (SST) represents the climate response to SST forcing as well as internal variability or “noise.” Signal-to-noise analysis is used to identify the most reproducible GCM patterns of African summer precipitation related to the SST forcing. Two of these potentially predictable components are associated with the precipitation of the Guinea Coast and Sahel regions and correlate well with observations. The GCM predictable component associated with rainfall in the Sahel region reproduces observed temporal variability on both interannual and decadal time scales, though with reduced amplitude.

Description: Cloud characteristics at two sites on the North Slope of Alaska separated by ∼100 km have been examined for the warmer months of 2001–03 using data collected from microwave radiometers, ceilometers, rotating shadowband radiometers, and pyranometers. Clouds at the inland site, Atqasuk, were found to have approximately 26% greater optical depths than those at the coastal site, Barrow, and the ratio of measured irradiance to clear-sky irradiance was nearly 20% larger at Barrow under cloudy conditions. It is hypothesized that a significant factor contributing to these differences is the upward fluxes of heat and water vapor over the wet tundra and lakes. Support for this hypothesis is found from the behavior of the liquid water paths for low clouds, which tend to be higher at Atqasuk than at Barrow for onshore winds but not for offshore ones, from differences in sensible heat fluxes, which are small but significant over the tundra but are nearly zero over the ocean adjacent to Barrow, and from the mixing ratios, which are significantly higher at Atqasuk than at Barrow. Results from a simple model further indicate that latent heat fluxes over the tundra and lakes can account for a significant fraction of the differences in the estimated boundary layer water content between Barrow and Atqasuk.

Description: Global climate models forced by sea surface temperature are standard tools in seasonal climate prediction and in projection of future climate change caused by anthropogenic emissions of greenhouse gases. Assessing the ability of these models to reproduce observed atmospheric circulation given the lower boundary conditions, and thus its ability to predict climate, has been a recurrent concern. Several assessments have shown that the performance of models is seasonally dependent, but there has always been the assumption that, for a given season, the model skill is constant throughout the period being analyzed. Here, it is demonstrated that there are periods when these models perform well and periods when they do not capture observed climate variability. The variations of the model performance have temporal scales and spatial patterns consistent with decadal/interdecadal climate variability. These results suggest that there are unmodeled climate processes that affect seasonal climate prediction as well as scenarios of climate change, particularly regional climate change projections. The reliability of these scenarios may depend on the time slice of the model output being analyzed. Therefore, more comprehensive model assessment should include a verification of the long-term stability of their performance.

Description: The Community Climate System Model version 3 (CCSM3) has recently been developed and released to the climate community. CCSM3 is a coupled climate model with components representing the atmosphere, ocean, sea ice, and land surface connected by a flux coupler. CCSM3 is designed to produce realistic simulations over a wide range of spatial resolutions, enabling inexpensive simulations lasting several millennia or detailed studies of continental-scale dynamics, variability, and climate change. This paper will show results from the configuration used for climate-change simulations with a T85 grid for the atmosphere and land and a grid with approximately 1° resolution for the ocean and sea ice. The new system incorporates several significant improvements in the physical parameterizations. The enhancements in the model physics are designed to reduce or eliminate several systematic biases in the mean climate produced by previous editions of CCSM. These include new treatments of cloud processes, aerosol radiative forcing, land–atmosphere fluxes, ocean mixed layer processes, and sea ice dynamics. There are significant improvements in the sea ice thickness, polar radiation budgets, tropical sea surface temperatures, and cloud radiative effects. CCSM3 can produce stable climate simulations of millennial duration without ad hoc adjustments to the fluxes exchanged among the component models. Nonetheless, there are still systematic biases in the ocean–atmosphere fluxes in coastal regions west of continents, the spectrum of ENSO variability, the spatial distribution of precipitation in the tropical oceans, and continental precipitation and surface air temperatures. Work is under way to extend CCSM to a more accurate and comprehensive model of the earth's climate system.

Description: The simulation of Arctic sea ice and surface winds changes significantly when Community Climate System Model version 3 (CCSM3) resolution is increased from T42 (∼2.8°) to T85 (∼1.4°). At T42 resolution, Arctic sea ice is too thick off the Siberian coast and too thin along the Canadian coast. Both of these biases are reduced at T85 resolution. The most prominent surface wind difference is the erroneous North Polar summer anticyclone, present at T42 but absent at T85. An offline sea ice model is used to study the effect of the surface winds on sea ice thickness. In this model, the surface wind stress is prescribed alternately from reanalysis and the T42 and T85 simulations. In the offline model, CCSM3 surface wind biases have a dramatic effect on sea ice distribution: with reanalysis surface winds annual-mean ice thickness is greatest along the Canadian coast, but with CCSM3 winds thickness is greater on the Siberian side. A significant difference between the two CCSM3-forced offline simulations is the thickness of the ice along the Canadian archipelago, where the T85 winds produce thicker ice than their T42 counterparts. Seasonal forcing experiments, with CCSM3 winds during spring and summer and reanalysis winds in fall and winter, relate the Canadian thickness difference to spring and summer surface wind differences. These experiments also show that the ice buildup on the Siberian coast at both resolutions is related to the fall and winter surface winds. The Arctic atmospheric circulation is examined further through comparisons of the winter sea level pressure (SLP) and eddy geopotential height. At both resolutions the simulated Beaufort high is quite weak, weaker at higher resolution. Eddy heights show that the wintertime Beaufort high in reanalysis has a barotropic vertical structure. In contrast, high CCSM3 SLP in Arctic winter is found in association with cold lower-tropospheric temperatures and a baroclinic vertical structure. In reanalysis, the summertime Arctic surface circulation is dominated by a polar cyclone, which is accompanied by surface inflow and a deep Ferrel cell north of the traditional polar cell. The Arctic Ferrel cell is accompanied by a northward flux of zonal momentum and a polar lobe of the zonal-mean jet. These features do not appear in the CCSM3 simulations at either resolution.

Description: The annual cycle of precipitation and the interannual variability of the North American hydroclimate during summer months are analyzed in coupled simulations of the twentieth-century climate. The state-of-the-art general circulation models, participating in the Fourth Assessment Report for the Intergovernmental Panel on Climate Change (IPCC), included in the present study are the U.S. Community Climate System Model version 3 (CCSM3), the Parallel Climate Model (PCM), the Goddard Institute for Space Studies model version EH (GISS-EH), and the Geophysical Fluid Dynamics Laboratory Coupled Model version 2.1 (GFDL-CM2.1); the Met Office’s Third Hadley Centre Coupled Ocean–Atmosphere GCM (UKMO-HadCM3); and the Japanese Model for Interdisciplinary Research on Climate version 3.2 [MIROC3.2(hires)]. Datasets with proven high quality such as NCEP’s North American Regional Reanalysis (NARR), and the Climate Prediction Center (CPC) U.S.–Mexico precipitation analysis are used as targets for simulations. Climatological precipitation is not easily simulated. While models capture winter precipitation very well over the U.S. northwest, they encounter failure over the U.S. southeast in the same season. Summer precipitation over the central United States and Mexico is also a great challenge for models, particularly the timing. In general the UKMO-HadCM3 is closest to the observations. The models’ potential in simulating interannual hydroclimate variability over North America during the warm season is varied and limited to the central United States. Models like PCM, and in particular UKMO-HadCM3, exhibit reasonably well the observed distribution and relative importance of remote and local contributions to precipitation variability over the region (i.e., convergence of remote moisture fluxes dominate over local evapotranspiration). However, in models like CCSM3 and GFDL-CM2.1 local contributions dominate over remote ones, in contrast with warm-season observations. In the other extreme are models like GISS-EH and MIROC3.2(hires) that prioritize the remote influence of moisture fluxes and neglect the local influence of land surface processes to the regional precipitation variability.

Description: The consistency between observed changes in Subantarctic Mode Water (SAMW) properties at 32°S in the Indian Ocean and model simulations is explored using the Third Hadley Centre Coupled Ocean–Atmosphere GCM (HadCM3). Hydrographic data collected in 2002 show that the water mass is warmer and saltier on isopycnals than in 1987, in contrast to the isopycnal freshening observed between 1962 and 1987. The response of HadCM3 under a range of forcing scenarios is explored and the observed freshening is only seen in experiments that include greenhouse gas forcing; however, there is no subsequent return to more saline conditions in 2002. The response of the model to greenhouse gas forcing is dominated by a persistent freshening trend, the simulated water mass variability agrees well with that suggested by the limited observations. Comparing model isopycnal changes from the forced experiments with a control run shows that the changes from the 1960s to 2002 are best explained by internal variability. This is in contrast to earlier work, which attributed the observed isopycnal freshening to anthropogenic forcing. Although the model shows that at present an anthropogenic climate change signal is not detectable in SAMW, the model water mass freshens on isopycnals during the twenty-first century under increased greenhouse gas forcing. This is consistent with recent heat content observations, which suggest that the salting is unlikely to persist. In HadCM3, this freshening is due to an increasing surface heat flux and Ekman heat and freshwater flux into the water mass formation region. This paper emphasizes the importance of higher-frequency observations of SAMW if detection and attribution statements are to be made.

Description: This study investigates the North Atlantic Oscillation (NAO) on an intraseasonal time scale. The authors investigate the question of how the characteristics of NAO events are influenced by the choice of its definitions using daily NCEP–NCAR reanalysis data spanning 51 boreal winters. Four different NAO indexes are used in this study, including one station/gridpoint–based index and three pattern-based indexes. It is found that the NAO events obtained using pattern–based indexes are quite similar to each other, while some notable differences are observed when the NAO is defined using the station/gridpoint–based index (NAO1). The characteristics of the pattern-based NAO are found to be more antisymmetric for its two phases, including its time-averaged spatial structures, its lifetime distributions, and time-evolving spatial structures. The NAO1, on the other hand, reveals some asymmetric characteristics between the two phases. Emphasis is placed on comparing the characteristics of the NAO events obtained using the NAO1 index and one of the pattern-based indices, that is, NAO2. The time-averaged spatial structures for the NAO2 expand across more of the polar region than the NAO1. The positive NAO1 shows a wave train signal over the Pacific–North American region during the setup phase, while the negative NAO1 is found to develop more locally over northern Europe and the North Atlantic. The wave activity flux for the NAO2 is primarily in the zonal direction while for the NAO1, on the other hand, it is mostly concentrated over the North Atlantic with a pronounced southward component. The barotropic vorticity equation is used to examine the physical mechanisms that drive the life cycle of the NAO.

Description: Coupled climate models generally have a small residual radiative flux at the top of the atmosphere. In the Met Office climate model, Hadley Centre Global Environmental Model version 1 (HadGEM1), it is incoming (heating the planet) and reduces over a 350-yr period from 0.4 to 0.1 W m−2. The process of the adjustment in HadGEM1 is examined and is shown to be linked to excessive heat gain. In the tropical and South Atlantic, cold, fresh Antarctic Intermediate Water is replaced by anomalously warm, salty intermediate water. The loss of Antarctic Intermediate Water in the South Atlantic is related to a weak Agulhas retroflection. The erosion is enhanced in the tropical Atlantic by strong upwelling. The warm, salty anomalies are advected northward outcropping in the North Atlantic subpolar gyre. In the outcrop zone, sea surface temperature and salinity are increased, which lead to an increase in global mean surface temperature and a reduction in the sea ice area. This adjusts the top of the atmosphere balance via increased outgoing longwave radiation and is partly offset by a decrease in outgoing shortwave radiation. The increased surface salinity triggers convection in the Labrador Sea and leads to a strong flushing of the thermohaline circulation. These results demonstrate that adjustment time scales for coupled climate models can be in excess of 350 yr. The potential implications of the adjustment time scale of climate models need to be considered when planning scenario and sensitivity experiments, as model drifts can be nonlinear.

Description: This study considers the probability distribution of sea surface wind speeds, which have historically been modeled using the Weibull distribution. First, non-Weibull structure in the observed sea surface wind speeds (from SeaWinds observations) is characterized using relative entropy, a natural information theoretic measure of the difference between probability distributions. Second, empirical models of the probability distribution of sea surface wind speeds, parameterized in terms of the parameters of the vector wind probability distribution, are developed. It is shown that Gaussian fluctuations in the vector wind cannot account for the observed features of the sea surface wind speed distribution, even if anisotropy in the fluctuations is accounted for. Four different non-Gaussian models of the vector wind distribution are then considered: the bi-Gaussian, the centered gamma, the Gram–Charlier, and the constrained maximum entropy. It is shown that so long as the relationship between the skewness and kurtosis of the along-mean sea surface wind component characteristic of observations is accounted for in the modeled probability distribution, then all four vector wind distributions are able to simulate the observed mean, standard deviation, and skewness of the sea surface wind speeds with an accuracy much higher than is possible if non-Gaussian structure in the vector winds is neglected. The constrained maximum entropy distribution is found to lead to the best simulation of the wind speed probability distribution. The significance of these results for the parameterization of air/sea fluxes in general circulation models is discussed.

Description: A 44-yr climatology of nonconvective wind events (NCWEs) for the Great Lakes region has been created using hourly wind data for 38 first-order weather stations during the months of November through April. The data were analyzed in terms of the two National Weather Service (NWS) criteria for a high-wind watch or warning: sustained winds of at least 18 m s−1 for at least 1 h or a wind gust of at least 26 m s−1 for any duration. The results indicate a pronounced southwest quadrant directional preference for nonconvective high winds in this region. Between 70% and 76% of all occurrences that satisfied the NWS criteria for NCWEs were associated with wind directions from 180° through 270°. Within the southwest quadrant, the west-southwest direction is preferred, with 14%–35% of all NCWEs coming from this particular compass heading. This directional preference is borne out in five out of six stations with high occurrences of cold-season NCWEs (Buffalo, New York; Dayton, Ohio; Lansing, Michigan; Moline, Illinois; Springfield, Illinois). Given the geographic spread of these stations, a nontopographic cause for the directional preference of cold-season NCWEs is indicated. The connection between NCWEs and low pressure systems found in this climatology and in case studies suggests that midlatitude cyclone dynamics may be a possible cause of the directional preference.

Description: Using satellite cloud data from the Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) and collocated precipitation rates from the Advanced Microwave Scanning Radiometer (AMSR), it is shown that rain rate is closely related to the amount of very thick high cloud, which is a better proxy for precipitation than outgoing longwave radiation (OLR). It is also shown that thin high cloud, which has a positive net radiative effect on the top-of-atmosphere (TOA) energy balance, is nearly twice as abundant in the west Pacific compared to the east Pacific. For a given rain rate, anvil cloud is also more abundant in the west Pacific. The ensemble of all high clouds in the east Pacific induces considerably more TOA radiative cooling compared to the west Pacific, primarily because of more high, thin cloud in the west Pacific. High clouds are also systematically colder in the west Pacific by about 5 K. The authors examine whether the anvil cloud temperature is better predicted by low-level equivalent potential temperature (ΘE), or by the peak in upper-level convergence associated with radiative cooling in clear skies. The temperature in the upper troposphere where ΘE is the same as that at the lifting condensation level (LCL) seems to influence the temperatures of the coldest, thickest clouds, but has no simple relation to anvil cloud. It is shown instead that a linear relationship exists between the median anvil cloud-top temperature and the temperature at the peak in clear-sky convergence. The radiatively driven clear-sky convergence profiles are thus consistent with the warmer anvil clouds in the EP versus the WP.

Description: Two new high-resolution sea surface temperature (SST) analysis products have been developed using optimum interpolation (OI). The analyses have a spatial grid resolution of 0.25° and a temporal resolution of 1 day. One product uses the Advanced Very High Resolution Radiometer (AVHRR) infrared satellite SST data. The other uses AVHRR and Advanced Microwave Scanning Radiometer (AMSR) on the NASA Earth Observing System satellite SST data. Both products also use in situ data from ships and buoys and include a large-scale adjustment of satellite biases with respect to the in situ data. Because of AMSR’s near-all-weather coverage, there is an increase in OI signal variance when AMSR is added to AVHRR. Thus, two products are needed to avoid an analysis variance jump when AMSR became available in June 2002. For both products, the results show improved spatial and temporal resolution compared to previous weekly 1° OI analyses. The AVHRR-only product uses Pathfinder AVHRR data (currently available from January 1985 to December 2005) and operational AVHRR data for 2006 onward. Pathfinder AVHRR was chosen over operational AVHRR, when available, because Pathfinder agrees better with the in situ data. The AMSR–AVHRR product begins with the start of AMSR data in June 2002. In this product, the primary AVHRR contribution is in regions near land where AMSR is not available. However, in cloud-free regions, use of both infrared and microwave instruments can reduce systematic biases because their error characteristics are independent.

Description: The Community Atmosphere Model version 3 (CAM3) is the latest generation of a long lineage of general circulation models produced by a collaboration between the National Center for Atmospheric Research (NCAR) and the scientific research community. Many aspects of the hydrological cycle have been changed relative to earlier versions of the model. It is the goal of this paper to document some aspects of the tropical variability of clouds and the hydrologic cycle in CAM3 on time scales shorter than 30 days and to discuss the differences compared to the observed atmosphere and earlier model versions, with a focus on cloud-top brightness temperature, precipitation, and cloud liquid water path. The transient behavior of the model in response to changes in resolution to various numerical methods used to solve the equations for atmospheric dynamics and transport and to the underlying lower boundary condition of sea surface temperature and surface fluxes has been explored. The ratio of stratiform to convective rainfall is much too low in CAM3, compared to observational estimates. It is much higher in CAM3 (10%) than the Community Climate Model version 3 (CCM3; order 1%–2%) but is still a factor of 4–5 too low compared to observational estimates. Some aspects of the model transients are sensitive to resolution. Higher-resolution versions of CAM3 show too much variability (both in amplitude and spatial extent) in brightness temperature on time scales of 2–10 days compared to observational estimates. Precipitation variance is underestimated on time scales from a few hours to 10 days, compared to observations over ocean, although again the biases are reduced compared to previous generations of the model. The diurnal cycle over tropical landmasses is somewhat too large, and there is not enough precipitation during evening hours. The model tends to produce maxima in precipitation and liquid water path that are a few hours earlier than that seen in the observations over both oceans and land.