ALBERT

Description: Common Land Model (CLM) and Land Surface Process (LSP) model simulations are compared to measured values for a 13-day dry-down period with a rapidly decreasing near-surface water table for a marsh wetland community in Florida. LSP was able to provide reasonable estimates without any modifications to the model physics. To obtain reasonable simulations using CLM, the baseline TOPMODEL baseflow generation and the bottom drainage mechanisms were not employed and the lower layers were allowed to remain saturated. In addition, several of CLM’s default wetland vegetation parameters were replaced with grassland parameters. Even after these modifications, CLM underestimated soil water storage. However, both model-simulated soil temperatures showed very good agreement as compared to measured temperatures, capturing both the soil warming during the study period and the diurnal fluctuations. Modeled surface energy fluxes also agreed well with measured values. LSP’s inability to consistently capture latent heat fluxes appears to be linked to its canopy resistance scaling functions. Other minor issues were that CLM’s rooting depth greatly exceeded observed depths and that CLM did not move water in the vadose zone from lower to upper layers during the nighttime as observed in the measurements. Overall, these results suggest that LSP can be applied to characterize a marsh dry down, but that minor modifications could greatly improve results. CLM demonstrated considerable potential, but requires some changes to model physics and default parameters prior to application to wetlands at a subgrid scale.

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: 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.

Description: An observed trend in the Southern Hemisphere annular mode (SAM) during recent decades has involved an intensification of the polar vortex. The source of this trend is a matter of scientific debate with stratospheric ozone losses, greenhouse gas increases, and natural variability all being possible contenders. Because it is difficult to separate the contribution of various external forcings to the observed trend, a state-of-the-art global coupled model is utilized here. Ensembles of twentieth-century simulations forced with the observed time series of greenhouse gases, tropospheric and stratospheric ozone, sulfate aerosols, volcanic aerosols, solar variability, and various combinations of these are used to examine the annular mode trends in comparison to observations, in an attempt to isolate the response of the climate system to each individual forcing. It is found that ozone changes are the biggest contributor to the observed summertime intensification of the southern polar vortex in the second half of the twentieth century, with increases of greenhouse gases also being a necessary factor in the reproduction of the observed trends at the surface. Although stratospheric ozone losses are expected to stabilize and eventually recover to preindustrial levels over the course of the twenty-first century, these results show that increasing greenhouse gases will continue to intensify the polar vortex throughout the twenty-first century, but that radiative forcing will cause widespread temperature increases over the entire Southern Hemisphere.

Description: The response of the North Atlantic thermohaline circulation to idealized climate forcing of 1% per year compound increase in CO2 is examined in three configurations of the Community Climate System Model version 3 that differ in their component model resolutions. The strength of the Atlantic overturning circulation declines at a rate of 22%–26% of the corresponding control experiment maximum overturning per century in response to the increase in CO2. The mean meridional overturning and its variability on decadal time scales in the control experiments, the rate of decrease in the transient forcing experiments, and the rate of recovery in periods of CO2 stabilization all increase with increasing component model resolution. By examining the changes in ocean surface forcing with increasing CO2 in the framework of the water-mass transformation function, we show that the decline in the overturning is driven by decreasing density of the subpolar North Atlantic due to increasing surface heat fluxes. While there is an intensification of the hydrologic cycle in response to increasing CO2, the net effect of changes in surface freshwater fluxes on those density classes that are involved in deep-water formation is to increase their density; that is, changes in surface freshwater fluxes act to maintain a stronger overturning circulation. The differences in the control experiment overturning strength and the response to increasing CO2 are well predicted by the corresponding differences in the water-mass transformation rate. Reduction of meridional heat transport and enhancement of meridional salt transport from mid- to high latitudes with increasing CO2 also act to strengthen the overturning circulation. Analysis of the trends in an ideal age tracer provides a direct measure of changes in ocean ventilation time scale in response to increasing CO2. In the subpolar North Atlantic south of the Greenland–Scotland ridge system, there is a significant increase in subsurface ages as open-ocean deep convection is diminished and ventilation switches to a predominance of overflow waters. In middle and low latitudes there is a decrease in age within and just below the thermocline in response to a decrease in the upwelling of old deep waters. However, when considering ventilation within isopycnal layers, age increases for layers in and below the thermocline due to the deepening of isopycnals in response to global warming.

Description: In this study the authors diagnose the sources for the contiguous U.S. seasonal forecast skill that are related to sea surface temperature (SST) variations using a combination of dynamical and empirical methods. The dynamical methods include ensemble simulations with four atmospheric general circulation models (AGCMs) forced by observed monthly global SSTs from 1950 to 1999, and ensemble AGCM experiments forced by idealized SST anomalies. The empirical methods involve a suite of reductions of the AGCM simulations. These include uni- and multivariate regression models that encapsulate the simultaneous and one-season lag linear connections between seasonal mean tropical SST anomalies and U.S. precipitation and surface air temperature. Nearly all of the AGCM skill in U.S. precipitation and surface air temperature, arising from global SST influences, can be explained by a single degree of freedom in the tropical SST field—that associated with the linear atmospheric signal of El Niño–Southern Oscillation (ENSO). The results support previous findings regarding the preeminence of ENSO as a U.S. skill source. The diagnostic methods used here exposed another skill source that appeared to be of non-ENSO origins. In late autumn, when the AGCM simulation skill of U.S. temperatures peaked in absolute value and in spatial coverage, the majority of that originated from SST variability in the subtropical west Pacific Ocean and the South China Sea. Hindcast experiments were performed for 1950–99 that revealed most of the simulation skill of the U.S. seasonal climate to be recoverable at one-season lag. The skill attributable to the AGCMs was shown to achieve parity with that attributable to empirical models derived purely from observational data. The diagnostics promote the interpretation that only limited advances in U.S. seasonal prediction skill should be expected from methods seeking to capitalize on sea surface predictors alone, and that advances that may occur in future decades could be readily masked by inherent multidecadal fluctuations in skill of coupled ocean–atmosphere systems.

Description: The Asian summer monsoon is organized into distinct convection centers, but the mechanism for this organization is not well understood. Analysis of new satellite observations reveals that narrow mountain ranges are an important organizing agent anchoring monsoon convection centers on the windward side. The Bay of Bengal convection, in particular, features the heaviest precipitation on its eastern coast because of orographic lifting as the southwest monsoon impinges on the coastal mountains of Myanmar (also known as Burma). This is in contrast to the widely held view that this convection is centered over the open ocean as implied by coarse-resolution datasets, a view that would require an entirely different explanation for its formation. Narrow in width and modest in height (≤1 km), these mountains are hardly mentioned in conceptual depictions of the large-scale monsoon and poorly represented in global climate models. The numerical simulations of this study show that orographic rainbands are not a local phenomenon but exert far-reaching effects on the continental-scale monsoon. The realization that these overlooked geographical features are an important element of the Asian monsoon has important implications for studying the monsoon in the past, present, and future.

Description: The rarity of severe coastal hurricanes implies that empirical estimates of extreme wind speed return levels will be unreliable. Here climatology models derived from extreme value theory are estimated using data from the best-track [Hurricane Database (HURDAT)] record. The occurrence of a hurricane above a specified threshold intensity level is assumed to follow a Poisson distribution, and the distribution of the maximum wind is assumed to follow a generalized Pareto distribution. The likelihood function is the product of the generalized Pareto probabilities for each wind speed estimate. A geographic region encompassing the entire U.S. coast vulnerable to Atlantic hurricanes is of primary interest, but the Gulf Coast, Florida, and the East Coast regions are also considered. Model parameters are first estimated using a maximum likelihood (ML) procedure. Results estimate the 100-yr return level for the entire coast at 157 kt (±10 kt), but at 117 kt (±4 kt) for the East Coast region (1 kt = 0.514 m s−1). Highest wind speed return levels are noted along the Gulf Coast from Texas to Alabama. The study also examines how the extreme wind return levels change depending on climate conditions including El Niño–Southern Oscillation, the Atlantic Multidecadal Oscillation, the North Atlantic Oscillation, and global temperature. The mean 5-yr return level during La Niña (El Niño) conditions is 125 (116) kt, but is 140 (164) kt for the 100-yr return level. This indicates that La Niña years are the most active for the occurrence of strong hurricanes, but that extreme hurricanes are more likely during El Niño years. Although El Niño inhibits hurricane formation in part through wind shear, the accompanying cooler lower stratosphere appears to increase the potential intensity of hurricanes that do form. To take advantage of older, less reliable data, the models are reformulated using Bayesian methods. Gibbs sampling is used to integrate the prior over the likelihood to obtain the posterior distributions for the model parameters conditional on global temperature. Higher temperatures are conditionally associated with more strong hurricanes and higher return levels for the strongest hurricane winds. Results compare favorably with an ML approach as well as with recent modeling and observational studies. The maximum possible near-coastal wind speed is estimated to be 208 kt (183 kt) using the Bayesian (ML) approach.

Description: In this study operational rawinsonde data are used to investigate climatological features of seasonal variations in static stability in order to understand the behavior of temperature inversion layers, that is, extremely stable layers in the lower troposphere over the Indochina Peninsula region, at the southeastern edge of the Asian continent. Static stability was evaluated from the vertical gradient in potential temperature (Δθ/Δz). Stable (Δθ/Δz 〉 10 K km−1) and unstable (Δθ/Δz 〈 1 K km−1) layers frequently appear over the Indochina Peninsula region during boreal winter. Temporal and vertical variations in stability during the boreal winter can be categorized into three characteristic types, type I: the mean height of stable layers increases from 2 to 5 km from the dry to the rainy season over inland areas of the Indochina Peninsula and southern China; type II: similar to type I, with the additional occurrence of stable layers at a height of ∼1 km, mainly over coastal areas of the Indochina Peninsula; and type III: stable layers at a height of ∼2 km, mainly over the Malay Peninsula. We did not find any significant seasonal change in the vertical distribution of stable layers over the Malay Peninsula.

Description: This study focuses on modulation of lightning and convective vertical structure in the southern Amazon as a function of the South American monsoon V index (VI). Four wet seasons (December–March 1998–2001) of Tropical Rainfall Measuring Mission (TRMM) Lightning Imaging Sensor (LIS) and Precipitation Radar (PR) data are examined together with two wet seasons (2000–01) of ground-based Brazilian Lightning Detection Network (BLDN) data. These observations are composited by VI phase (northerly or southerly) for a region of the southern Amazon and discussed relative to VI-regime environmental characteristics such as thermodynamic instability and wind shear. Relative comparisons of VI-regime convective properties reveal 1) slightly larger (20%–25%) PR pixel-mean rainfall during periods of northerly VI due to increased stratiform precipitation, 2) a factor of 2 or more increase in lightning flash density and the lightning diurnal cycle amplitude during periods of southerly VI, 3) a factor of 1.5–2 increase in the conditional probability of any PR radar reflectivity pixel exceeding 30 dBZ above the −10°C level during periods of southerly VI, and 4) an associated factor of 2 or more increase in southerly VI pixel-mean ice water path, with the ice water path being highly correlated to trends in lightning activity. During periods of southerly VI, convection occurs in an environment of increased thermodynamic instability, weak southeasterly low-level, and deep upper-tropospheric easterly wind shear. During periods of northerly VI, low-level westerly shear opposes stronger deep tropospheric easterly shear in a relatively moist environment of weaker thermodynamic instability, consistent with the occurrence of more widespread stratiform precipitation. The composite results of this study point to 1) regime differences in convective forcing that alter the prevalence of ice processes and, by inference, the vertical profile of latent heating and 2) the utility of lightning observations in delineating convective regime changes.

Description: The observed structure and seasonal evolution characteristics of the tropospheric biennial oscillation (TBO) in the warm ocean areas of the Indo-Pacific region are explored using a seasonal-sequence EOF analysis approach. The major convective activity centers associated with the TBO appear in the southeast Indian Ocean (SEIO) and western North Pacific (WNP), accompanied by anticyclonic (or cyclonic) circulation patterns with a first-baroclinic-mode structure. The convection and circulation anomalies have distinctive life cycles in the SEIO and WNP: the former have a peak phase in northern fall and the latter persist from northern winter to subsequent summer. The mechanisms of the TBO in this region are investigated with a hybrid coupled GCM. Numerical results show that air–sea interaction in the warm ocean alone can support TBO variability that has many observed characteristics. Key processes involved in the TBO include the WNP monsoon variability and associated cross-equatorial flows, convective activity over Southeast Asia/the Maritime Continent and associated anomalous Walker circulation, and ocean dynamic responses to anomalous wind stress curl in the western Pacific. The coupled model experiment demonstrates that the essential element of the TBO in this region arises from the monsoon–warm ocean interaction. A possible connection between the TBO and ENSO variability is further studied in another model that excludes the delayed oscillator dynamics. The key in causing the biennial variability of ENSO arises from teleconnections between the tropical Pacific and Indian Oceans, with three “atmospheric bridges”: 1) the north–south teleconnection that connects the WNP monsoon and the SEIO, 2) the east–west teleconnection that connects the Indian Ocean and the Pacific cold tongue, and 3) the El Niño–WNP monsoon teleconnection.

Description: The sea ice simulation of the Community Climate System Model version 3 (CCSM3) T42-gx1 and T85-gx1 control simulations is presented and the influence of the parameterized sea ice thickness distribution (ITD) on polar climate conditions is examined. This includes an analysis of the change in mean climate conditions and simulated sea ice feedbacks when an ITD is included. It is found that including a representation of the subgrid-scale ITD results in larger ice growth rates and thicker sea ice. These larger growth rates represent a higher heat loss from the ocean ice column to the atmosphere, resulting in warmer surface conditions. Ocean circulation, most notably in the Southern Hemisphere, is also modified by the ITD because of the influence of enhanced high-latitude ice formation on the ocean buoyancy flux and resulting deep water formation. Changes in atmospheric circulation also result, again most notably in the Southern Hemisphere. There are indications that the ITD also modifies simulated sea ice–related feedbacks. In regions of similar ice thickness, the surface albedo changes at 2XCO2 conditions are larger when an ITD is included, suggesting an enhanced surface albedo feedback. The presence of an ITD also modifies the ice thickness–ice strength relationship and the ice thickness–ice growth rate relationship, both of which represent negative feedbacks on ice thickness. The net influence of the ITD on polar climate sensitivity and variability results from the interaction of these and other complex feedback processes.

Description: Simulations of the El Niño–Southern Oscillation (ENSO) phenomenon and tropical Atlantic climate variability in the newest version of the Community Climate System Model [version 3 (CCSM3)] are examined in comparison with observations and previous versions of the model. The analyses are based upon multicentury control integrations of CCSM3 at two different horizontal resolutions (T42 and T85) under present-day CO2 concentrations. Complementary uncoupled integrations with the atmosphere and ocean component models forced by observed time-varying boundary conditions allow an assessment of the impact of air–sea coupling upon the simulated characteristics of ENSO and tropical Atlantic variability. The amplitude and zonal extent of equatorial Pacific sea surface temperature variability associated with ENSO is well simulated in CCSM3 at both resolutions and represents an improvement relative to previous versions of the model. However, the period of ENSO remains too short (2–2.5 yr in CCSM3 compared to 2.5–8 yr in observations), and the sea surface temperature, wind stress, precipitation, and thermocline depth responses are too narrowly confined about the equator. The latter shortcoming is partially overcome in the atmosphere-only and ocean-only simulations, indicating that coupling between the two model components is a contributing cause. The relationships among sea surface temperature, thermocline depth, and zonal wind stress anomalies are consistent with the delayed/recharge oscillator paradigms for ENSO. We speculate that the overly narrow meridional scale of CCSM3's ENSO simulation may contribute to its excessively high frequency. The amplitude and spatial pattern of the extratropical atmospheric circulation response to ENSO is generally well simulated in the T85 version of CCSM3, with realistic impacts upon surface air temperature and precipitation; the simulation is not as good at T42. CCSM3's simulation of interannual climate variability in the tropical Atlantic sector, including variability intrinsic to the basin and that associated with the remote influence of ENSO, exhibits similarities and differences with observations. Specifically, the observed counterpart of El Niño in the equatorial Atlantic is absent from the coupled model at both horizontal resolutions (as it was in earlier versions of the coupled model), but there are realistic (although weaker than observed) SST anomalies in the northern and southern tropical Atlantic that affect the position of the local intertropical convergence zone, and the remote influence of ENSO is similar in strength to observations, although the spatial pattern is somewhat different.

Description: A new version of the Community Atmosphere Model (CAM) has been developed and released to the climate community. CAM Version 3 (CAM3) is an atmospheric general circulation model that includes the Community Land Model (CLM3), an optional slab ocean model, and a thermodynamic sea ice model. The dynamics and physics in CAM3 have been changed substantially compared to implementations in previous versions. CAM3 includes options for Eulerian spectral, semi-Lagrangian, and finite-volume formulations of the dynamical equations. It supports coupled simulations using either finite-volume or Eulerian dynamics through an explicit set of adjustable parameters governing the model time step, cloud parameterizations, and condensation processes. The model includes major modifications to the parameterizations of moist processes, radiation processes, and aerosols. These changes have improved several aspects of the simulated climate, including more realistic tropical tropopause temperatures, boreal winter land surface temperatures, surface insolation, and clear-sky surface radiation in polar regions. The variation of cloud radiative forcing during ENSO events exhibits much better agreement with satellite observations. Despite these improvements, several systematic biases reduce the fidelity of the simulations. These biases include underestimation of tropical variability, errors in tropical oceanic surface fluxes, underestimation of implied ocean heat transport in the Southern Hemisphere, excessive surface stress in the storm tracks, and offsets in the 500-mb height field and the Aleutian low.

Description: Simulations of regional monsoon regimes, including the Indian, Australian, West African, South American, and North American monsoons, are described for the T85 version of the Community Climate System Model version 3 (CCSM3) and compared to observations and Atmospheric Model Intercomparison Project (AMIP)-type SST-forced simulations with the Community Atmospheric Model version 3 (CAM3) at T42 and T85. There are notable improvements in the regional aspects of the precipitation simulations in going to the higher-resolution T85 compared to T42 where topography is important (e.g., Ethiopian Highlands, South American Andes, and Tibetan Plateau). For the T85 coupled version of CCSM3, systematic SST errors are associated with regional precipitation errors in the monsoon regimes of South America and West Africa, though some aspects of the monsoon simulations, particularly in Asia, improve in the coupled model compared to the SST-forced simulations. There is very little realistic intraseasonal monsoon variability in the CCSM3 consistent with earlier versions of the model. Teleconnections to the tropical Pacific are well simulated for the South Asian monsoon.

Description: The largest and potentially most important ocean near-surface biases are examined in the Community Climate System Model coupled simulation of present-day conditions. They are attributed to problems in the component models of the ocean or atmosphere, or both. Tropical biases in sea surface salinity (SSS) are associated with precipitation errors, with the most striking being a band of excess rainfall across the South Pacific at about 8°S. Cooler-than-observed equatorial Pacific sea surface temperature (SST) is necessary to control a potentially catastrophic positive feedback, involving precipitation along the equator. The strength of the wind-driven gyres and interbasin exchange is in reasonable agreement with observations, despite the generally too strong near-surface winds. However, the winds drive far too much transport through Drake Passage [〉190 Sv (1 Sv ≡ 106 m3 s−1)], but with little effect on SST and SSS. Problems with the width, separation, and location of western boundary currents and their extensions create large correlated SST and SSS biases in midlatitudes. Ocean model deficiencies are suspected because similar signals are seen in uncoupled ocean solutions, but there is no evidence of serious remote impacts. The seasonal cycles of SST and winds in the equatorial Pacific are not well represented, and numerical experiments suggest that these problems are initiated by the coupling of either or both wind components. The largest mean SST biases develop along the eastern boundaries of subtropical gyres, and the overall coupled model response is found to be linear. In the South Atlantic, surface currents advect these biases across much of the tropical basin. Significant precipitation responses are found both in the northwest Indian Ocean, and locally where the net result is the loss of an identifiable Atlantic intertropical convergence zone, which can be regained by controlling the coastal temperatures and salinities. Biases off South America and Baja California are shown to significantly degrade precipitation across the Pacific, subsurface ocean properties on both sides of the equator, and the seasonal cycle of equatorial SST in the eastern Pacific. These signals extend beyond the reach of surface currents, so connections via the atmosphere and subsurface ocean are implicated. Other experimental results indicate that the local atmospheric forcing is only part of the problem along eastern boundaries, with the representation of ocean upwelling another likely contributor.

Description: The Community Land Model version 3 (CLM3) Dynamic Global Vegetation Model (CLM–DGVM) is used diagnostically to identify land and atmospheric model biases that lead to biases in the simulated vegetation. The CLM–DGVM driven with observed atmospheric data (offline simulation) underestimates global forest cover, overestimates grasslands, and underestimates global net primary production. These results are consistent with earlier findings that the soils in CLM3 are too dry. In the offline simulation an increase in simulated transpiration by changing this variable's soil moisture dependence and by decreasing canopy-intercepted precipitation results in better global plant biogeography and global net primary production. When CLM–DGVM is coupled to the Community Atmosphere Model version 3 (CAM3), the same modifications do not improve simulated vegetation in the eastern United States and Amazonia where the most serious vegetation biases appear. The dry bias in eastern U.S. precipitation is so severe that the simulated vegetation is insensitive to changes in the hydrologic cycle. In Amazonia, strong coupling among soil moisture, vegetation, evapotranspiration, and precipitation produces a highly complex hydrologic cycle in which small perturbations in precipitation are accentuated by vegetation. These interactions in Amazonia lead to a dramatic precipitation decrease and a collapse of the forest. These results suggest that the accurate parameterization of convection poses a complex and challenging scientific issue for climate models that include dynamic vegetation. The results also emphasize the difficulties that may arise when coupling any two highly nonlinear systems that have only been tested uncoupled.

Description: New features that may affect the behavior of the upper ocean in the Community Climate System Model version 3 (CCSM3) are described. In particular, the addition of an idealized diurnal cycle of solar forcing where the daily mean solar radiation received in each daily coupling interval is distributed over 12 daylight hours is evaluated. The motivation for this simple diurnal cycle is to improve the behavior of the upper ocean, relative to the constant forcing over each day of previous CCSM versions. Both 1- and 3-h coupling intervals are also considered as possible alternatives that explicitly resolve the diurnal cycle of solar forcing. The most prominent and robust effects of all these diurnal cycles are found in the tropical oceans, especially in the Pacific. Here, the mean equatorial sea surface temperature (SST) is warmed by as much as 1°C, in better agreement with observations, and the mean boundary layer depth is reduced. Simple rectification of the diurnal cycle explains about half of the shallowing, but less than 0.1°C of the warming. The atmospheric response to prescribed warm SST anomalies of about 1°C displays a very different heat flux signature. The implication, yet to be verified, is that large-scale air–sea coupling is a prime mechanism for amplifying the rectified, daily averaged SST signals seen by the atmosphere. Although the use of upper-layer temperature for SST in CCSM3 underestimates the diurnal cycle of SST, many of the essential characteristics of diurnal cycling within the equatorial ocean are reproduced, including boundary layer depth, currents, and the parameterized vertical heat and momentum fluxes associated with deep-cycle turbulence. The conclusion is that the implementation of an idealized diurnal cycle of solar forcing may make more frequent ocean coupling and its computational complications unnecessary as improvements to the air–sea coupling in CCSM3 continue. A caveat here is that more frequent ocean coupling tends to reduce the long-term cooling trends typical of CCSM3 by heating already too warm ocean depths, but longer integrations are needed to determine robust features. A clear result is that the absence of diurnal solar forcing of the ocean has several undesirable consequences in CCSM3, including too large ENSO variability, much too cold Pacific equatorial SST, and no deep-cycle turbulence.

Description: Because land surface emissivity (ɛ) has not been reliably measured, global climate model (GCM) land surface schemes conventionally set this parameter as simply constant, for example, 1 as in the National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Prediction (NCEP) model, and 0.96 for bare soil as in the National Center for Atmospheric Research (NCAR) Community Land Model version 2 (CLM2). This is the so-called constant-emissivity assumption. Accurate broadband emissivity data are needed as model inputs to better simulate the land surface climate. It is demonstrated in this paper that the assumption of the constant emissivity induces errors in modeling the surface energy budget, especially over large arid and semiarid areas where ɛ is far smaller than unity. One feasible solution to this problem is to apply the satellite-based broadband emissivity into land surface models. The Moderate Resolution Imaging Spectroradiometer (MODIS) instrument has routinely measured spectral emissivities (ɛλ) in six thermal infrared bands. The empirical regression equations have been developed in this study to convert these spectral emissivities to broadband emissivity (ɛ) required by land surface models. The observed emissivity data show strong seasonality and land-cover dependence. Specifically, emissivity depends on surface-cover type, soil moisture content, soil organic composition, vegetation density, and structure. For example, broadband ɛ is usually around 0.96–0.98 for densely vegetated areas [(leaf area index) LAI 〉 2], but it can be lower than 0.90 for bare soils (e.g., desert). To examine the impact of variable surface broadband emissivity, sensitivity studies were conducted using offline CLM2 and coupled NCAR Community Atmosphere Models, CAM2–CLM2. These sensitivity studies illustrate that large impacts of surface ɛ occur over deserts, with changes up to 1°–2°C in ground temperature, surface skin temperature, and 2-m surface air temperature, as well as evident changes in sensible and latent heat fluxes.

Description: On the premise that urban heat islands are strongest in calm conditions but are largely absent in windy weather, daily minimum and maximum air temperatures for the period 1950–2000 at a worldwide selection of land stations are analyzed separately for windy and calm conditions, and the global and regional trends are compared. The trends in temperature are almost unaffected by this subsampling, indicating that urban development and other local or instrumental influences have contributed little overall to the observed warming trends. The trends of temperature averaged over the selected land stations worldwide are in close agreement with published trends based on much more complete networks, indicating that the smaller selection used here is sufficient for reliable sampling of global trends as well as interannual variations. A small tendency for windy days to have warmed more than other days in winter over Eurasia is the opposite of that expected from urbanization and is likely to be a consequence of atmospheric circulation changes.

Description: A passive acoustic method of detecting breaking waves of different scales has been developed. The method also showed promise for measuring breaking severity. Sounds were measured by a subsurface hydrophone in various wind and wave states. A video record of the surface was made simultaneously. Individual sound pulses corresponding to the many individual bubble formations during wave-breaking events typically last only a few tens of milliseconds. Each time a sound-level threshold was exceeded, the acoustic signal was captured over a brief window typical of a bubble formation pulse, registering one count. Each pulse was also analyzed to determine the likely bubble size generating the pulse. Using the time series of counts and visual observations of the video record, the sound-level threshold that detected bubble formations at a rate optimally discriminating between breaking and nonbreaking waves was determined by a classification-accuracy analysis. This diagnosis of breaking waves was found to be approximately 70%–75% accurate once the optimum threshold had been determined. The method was then used for detailed analysis of wave-breaking properties across the spectrum. When applied to real field data, a breaking probability distribution could be obtained. This is the rate of occurrence of wave-breaking events at different wave scales. With support from a separate, laboratory experiment, the estimated bubble size is argued to be dependent on the severity of wave breaking and thus to provide information on the energy loss due to the breaking at the measured spectral frequencies. A combination of the breaking probability distribution and the bubble size could lead to direct estimates of spectral distribution of wave dissipation.

Description: This paper highlights recent results obtained with the Turbulent Eddy Profiler (TEP), which was developed by the University of Massachusetts. This unique 915-MHz radar has up to 64 spatially separated receiving elements, each with an independent receiver. The calibrated raw data provided by this array could be processed using sophisticated imaging algorithms to resolve the horizontal structures within each range gate. After collecting all of the closely spaced horizontal slices, the TEP radar can produce three-dimensional images of echo power, radial velocity, and spectral width. From the radial velocity measurements, it is possible to estimate the three-dimensional wind with high horizontal and vertical resolution. Given the flexibility of the TEP system, various array configurations are possible. In the present work exploitation of the flexibility of TEP is attempted to enhance the rejection of clutter from unwanted biological targets. From statistical studies, most biological clutter results from targets outside the main imaging field of view, that is, the sidelobes and grating lobes (if they exist) of the receiving beam. Because the TEP array's minimum receiver separation exceeds the spatial Nyquist sampling requirement, substantial possibilities for grating-lobe clutter exist and are observed in actual array data. When imaging over the transmit beam volume, the receiving array main lobe is scanned over a ±12.5° region. This scanning also sweeps the grating lobes over a wide angular region, virtually guaranteeing that a biological scatterer outside of the main beam will appear somewhere in the imaged volume. This paper focuses on suppressing pointlike targets in the grating-lobe regions. With a subtle change to the standard TEP array hardware configuration, it is shown via simulations and actual experimental observations (collected in June 2003) that adaptive beamforming methods can subsequently be used to significantly suppress the effects of point targets on the wind field estimates. These pointlike targets can be birds or planes with strong reflectivity. By pointlike the authors mean its appearance is a distinct point (up to the imaging resolution) in the images. The pointlike strong reflectivity signature exploits the capability of adaptive beamforming to suppress the interference using the new array configuration. It should be noted that this same array configuration does not exhibit this beneficial effect when standard Fourier beamforming is employed.