ALBERT

Description: Advances in computer power, new forecasting challenges, and new diagnostic techniques have brought about changes in the way atmospheric development and vertical motion are diagnosed in an operational setting. Many of these changes, such as improved model skill, model resolution, and ensemble forecasting, have arguably been detrimental to the ability of forecasters to understand and respond to the evolving atmosphere. The use of nondivergent wind in place of geostrophic wind would be a step in the right direction, but the advantages of potential vorticity suggest that its widespread adoption as a diagnostic tool on the west side of the Atlantic is overdue. Ertel potential vorticity (PV), when scaled to be compatible with pseudopotential vorticity, is generally similar to pseudopotential vorticity, so forecasters accustomed to quasigeostrophic reasoning through the height tendency equation can transfer some of their intuition into the Ertel-PV framework. Indeed, many of the differences between pseudopotential vorticity and Ertel potential vorticity are consequences of the choice of definition of quasigeostrophic PV and are not fundamental to the quasigeostrophic system. Thus, at its core, PV thinking is consistent with commonly used quasigeostrophic diagnostic techniques.

Description: Synoptic and mesoscale meteorology underwent a revolution in the 1940s and 1950s with the widespread deployment of novel weather observations, such as the radiosonde network and the advent of weather radar. These observations provoked a rapid increase in our understanding of the structure and dynamics of the atmosphere by pioneering analysts such as Fred Sanders. The authors argue that we may be approaching an analogous revolution in our ability to study the structure and dynamics of atmospheric phenomena with the advent of probabilistic objective analyses. These probabilistic analyses provide not only best estimates of the state of the atmosphere (e.g., the expected value) and the uncertainty about this state (e.g., the variance), but also the relationships between all locations and all variables at that instant in time. Up until now, these relationships have been determined by sampling in time by, for example, case studies, composites, and time-series analysis. Here the authors propose a new approach, ensemble synoptic analysis, which exploits the information contained in probabilistic samples of analyses at one or more instants in time. One source of probabilistic analyses is ensemble-based state-estimation methods, such as ensemble-based Kalman filters. Analyses from such a filter may be used to study atmospheric phenomena and the relationships between fields and locations at one or more instants in time. After a brief overview of a research-based ensemble Kalman filter, illustrative examples of ensemble synoptic analysis are given for an extratropical cyclone, including relationships between the cyclone minimum sea level pressure and other synoptic features, statistically determined operators for potential-vorticity inversion, and ensemble-based sensitivity analysis.

Description: The pioneering large-scale studies of cyclone frequency, location, and intensity conducted by Fred Sanders prompt similar questions about lesser-studied anticyclone development. The results of a climatology of closed anticyclones (CAs) at 200, 500, and 850 hPa, with an emphasis on the subtropics and midlatitudes, is presented to assess the seasonally varying distribution and hemispheric differences of these features. To construct the CA climatology, a counting program was applied to twice-daily 2.5° NCEP–NCAR reanalysis 200-, 500-, and 850-hPa geopotential height fields for the period 1950–2003. Stationary CAs, defined as those CAs that were located at a particular location for consecutive time periods, were counted only once. The climatology results show that 200-hPa CAs occur preferentially during summer over subtropical continental regions, while 500-hPa CAs occur preferentially over subtropical oceans in all seasons and over subtropical continents in summer. Conversely, 850-hPa CAs occur preferentially over oceanic regions beneath upper-level midocean troughs, and are most prominent in the Northern Hemisphere, and over midlatitude continents in winter. Three case studies of objectively identified CAs that produced heal waves over the United States, Europe, and Australia in 1995, 2003, and 2004, respectively, are presented to supplement the climatological results. The case studies, examining the subset of CAs than can produce heat waves, illustrate how climatologically hot continental tropical air masses produced over arid and semiarid regions of the subtropics and lower midlatitudes can become abnormally hot in conjunction with dynamically driven upper-level ridge amplification. Subsequently, these abnormally hot air masses are advected downstream away from their source regions in conjunction with transient disturbances embedded in anomalously strong westerly jets.

Description: Oklahoma Mesonetwork data are used to illustrate important atmospheric features that are not well shown by the usual synoptic data. For example, some shifts of wind from south to north that are shown as cold fronts on synoptic charts are not cold fronts by any plausible definition. As previously discussed by Fred Sanders and others, such errors in analysis can be reduced by knowledge of the wide variety of weather phenomena that actually exists, and by more attention to temperatures at the earth's surface as revealed by conventional synoptic data. Mesoscale data for four cases reinforce previous discussions of the ephemeral nature of fronts and deficiencies in the usual analyses of cold fronts. One type of misanalyzed case involves post-cold-frontal boundary layer air that is warmer than the prefrontal air. A second type is usually nocturnal, with a rise of local temperature during disruption of an inversion and a wind shift with later cooling that accompanies advection of a climatological gradient of temperature.

Description: The advent of the polar front theory of cyclones in Norway early in the last century held that the development of fronts and air masses is central to understanding midlatitude weather phenomena. While work on fronts continues to this day, the concept of air masses has been largely forgotten, superseded by the idea of a continuum. The Norwegians placed equal emphasis on the thermodynamics of airmass formation and on the dynamical processes that moved air masses around; today, almost all the emphasis is on dynamics, with little published literature on diabatic processes acting on a large scale. In this essay, the author argues that a lack of understanding of large-scale diabatic processes leads to an incomplete picture of the atmosphere and contributes to systematic errors in medium- and long-range weather forecasts. At the same time, modern concepts centered around potential vorticity conservation and inversion lead one to a redefinition of the term "air mass" that may have some utility in conceptualizing atmospheric physics and in weather forecasting.

Description: Historically, the atmospheric sciences have tended to treat problems of weather and climate separately. The real physical system, however, is a continuum, with short-term (minutes to days) “weather” fluctuations influencing climate variations and change, and, conversely, more slowly varying aspects of the system (typical time scales of a season or longer) affecting the weather that is experienced. While this past approach has served important purposes, it is becoming increasingly apparent that in order to make progress in addressing many socially important problems, an improved understanding of the connections between weather and climate is required. This overview summarizes the progress over the last few decades in the understanding of the phenomena and mechanisms linking weather and climate variations. The principal emphasis is on developments in understanding key phenomena and processes that bridge the time scales between synoptic-scale weather variability (periods of approximately 1 week) and climate variations of a season or longer. Advances in the ability to identify synoptic features, improve physical understanding, and develop forecast skill within this time range are reviewed, focusing on a subset of major, recurrent phenomena that impact extratropical wintertime weather and climate variations over the Pacific–North American region. While progress has been impressive, research has also illuminated areas where future gains are possible. This article concludes with suggestions on near-term directions for advancing the understanding and capabilities to predict the connections between weather and climate variations.

Description: Fred Sanders' career extended over 55 yr, touching upon many of the revolutionary transformations in the field of meteorology during that period. In this paper, his contributions to the transformation of synoptic meteorology, his research into the nature of explosive cyclogenesis, and related advances in the ability to predict these storms are reviewed. In addition to this review, the current status of forecasting oceanic cyclones 4.5 days in advance is presented, illustrating the progress that has been made and the challenges that persist, especially for forecasting those extreme extratropical cyclones that are marked by surface wind speeds exceeding hurricane force. Last, Fred Sanders' participation in a forecast for the historic 1947 snowstorm (that produced snowfall amounts in the New York City area that set records at that time) is reviewed along with an attempt to use today's operational global model to simulate this storm using data that were available at the time. The study reveals the predictive limitations involved with this case based on the scarcity of upper-air data in 1947, while confirming Fred Sanders' forecasting skills when dealing with these types of major storm events, even as a young aviation forecaster at New York's LaGuardia Airport.

Description: A case study of a double dryline on 22 May 2002 is presented. Mobile, 3-mm-wavelength Doppler radars from the University of Massachusetts and the University of Wyoming (Wyoming cloud radar) were used to collect very fine resolution vertical-velocity data in the vicinity of each of the moisture gradients associated with the drylines. Very narrow (50–100 m wide) channels of strong upward vertical velocity (up to 8 m s–1) were measured in the convergence zone of the easternmost dryline, larger in magnitude than reported with previous drylines. Distinct areas of descending motion were evident to the east and west of both drylines. Radar data are interpreted in the context of other observational platforms available during the International H2O Project (IHOP-2002). a variational ground-based mobile radar data processing technique was developed and applied to pseudo-dual-Doppler data collected during a rolling range-height indicator deployment. It was found that there was a secondary (vertical) circulation normal to the easternmost moisture gradient; the circulation comprised an easterly component near-surface flow to the east, a strong upward vertical component in the convergence zone, a westerly return, flow above the convective boundary layer, and numerous regions of descending motion, the most prominent approximately 3–5 km to the east of the surface convergence zone.

Description: Fred Sanders' teaching and research contributions in the area of quasigeostrophic theory are highlighted in this paper. The application of these contributions is made to the topic of extreme cold-season precipitation events in the Saint Lawrence valley in the northeastern United States and southern Quebec. This research focuses on analyses of Saint Lawrence valley heavy precipitation events. Synoptic- and planetary-scale circulation anomaly precursors are typically identified several days prior to these events. These precursors include transient upper-level troughs, strong moisture transports into the region, and anomalously large precipitable water amounts. The physical insight of Fred Sanders' work is used in the analysis of these composite results. Further details of this insight are provided in analyses of one case of heavy precipitation.

Description: Sanders designed a barotropic tropical cyclone (TC) track prediction model for the North Atlantic TC basin that became known as the Sanders barotropic (SANBAR) model. It predicted the streamfunction of the deeplayer mean winds (tropical circulation vertically averaged from 1000 to 100 hPa) that represents the vertically averaged tropical circulations. Originally, the wind input for the operational objective analysis (OA) consisted of winds measured by radiosondes and 44 bogus winds provided by analysis at the National Hurricane Center (NHC), which corresponded to the vertically averaged flow over sparsely observed tropical, subtropical, and midlatitude oceanic regions. The model covered a fixed regional area and had a grid size of ~ 154 km. It estimated the initial storm motion solely on the basis of the large-scale flow from the OA, not taking into account the observed storm motion. During 1970, the SANBAR model became the first dynamical TC track model to be run operationally at NHC. Track forecasts of SANBAR were verified from the 1971 TC season when track model verifications began at NHC until its retirement after the 1989 Atlantic TC season. The average annual SANBAR forecast track errors were verified relative to Climatology and Persistence (CLIPER), the standard no-skill track forecast. Comparison with CLIPER determines the skill of track forecast methods. Verifications are presented for two different versions of the SANBAR model system used operationally during 1973–84 and 1985–89. In homogeneous comparisons (i.e., includes only forecasts for the same initial times) for the former period, SANBAR's track forecasts were slightly better than CLIPER at 24–48-h forecast intervals; however, from 1985 to 1989 the average SANBAR track forecast errors from 24–72 h were ~10% more skillful than homogeneous CLIPER track forecasts.

Description: In the last decade, Fred Sanders was often critical of current surface analysis techniques. This led to his promoting the use of surface potential temperatures to distinguish between fronts, baroclinic troughs, and non-frontal baroclinic zones, and to the development of a climatology of surface baroclinic zones. In this paper, criticisms of current surface analysis techniques and the usefulness of surface potential temperature analyses are discussed. Case examples are used to compare potential temperature analyses and current National Centers for Environmental Prediction analyses. The 1-yr climatology of Sanders and Hoffman is reconstructed using a composite technique. Annual and seasonal mean potential temperature analyses over the continental United States, southern Canada, northern Mexico, and adjacent coastal waters are presented. In addition, gridpoint frequencies of moderate and strong potential temperature gradients are calculated. The results of the mean potential temperature analyses show that moderate and strong surface baroclinic zones are favored along the coastlines and the slopes of the North American cordillera. Additional subsynoptic details, not found in Sanders and Hoffman, are identified. The availability of the composite results allows for the calculation of potential temperature gradient anomalies. It is shown that these anomalies can be used to identify significant frontal baroclinic zones that are associated with weak potential temperature gradients. Together the results and reviews in this paper show that surface potential temperature analyses are a valuable forecasting and analysis tool allowing analysts to distinguish and identify fronts, baroclinic troughs, and nonfrontal baroclinic zones.

Description: The nature of the different types of surface boundaries that appear in the southern plains of the United States during the convectively active season is reviewed. The following boundaries are discussed: fronts, the dryline, troughs, and outflow boundaries, The boundaries are related to their environment and to local topography. The role these boundaries might play in the initiation of convective storms is emphasized. The various types of boundary-related vertical circulations and their dynamics are discussed. In particular, quasigeostrophic and semigeostrophic dynamics, and the dynamics of solenoidal circulations, density currents, boundary layers, and gravity waves are considered. Miscellaneous topics pertinent to convective storms and their relationship to surface boundaries such as along-the-boundary variability, boundary collisions, and the role of vertical shear are also discussed. Although some cases of storm initiation along surface boundaries have been well documented using research datasets collected during comprehensive field experiments, much of what we know is based only on empirical forecasting and nowcasting experience. It is suggested that many problems relating to convective-storm formation need to be explored in detail using real datasets with new observing systems and techniques, in conjunction with numerical simulation studies, and through climatological studies.

Description: This paper begins with a review of basic surface frontogenesis concepts with an emphasis on fronts located over sloping terrain adjacent to mountain barriers and fronts located in large-scale baroclinic zones close to coastlines. The impact of cold-air damming and differential diabatic heating and cooling on frontogenesis is considered through two detailed case studies of intense surface fronts. The first case, from 17 to 18 April 2002, featured the westward passage of a cold (side-door) front across coastal eastern New England in which 15°–20°C temperature decreases were observed in less than one hour. The second case, from 28 February to 4 March 1972, featured a long-lived front that affected most of the United States from the Rockies to the Atlantic coast and was noteworthy for a 50°C temperature contrast between Kansas and southern Manitoba, Canada. In the April 2002 case most of New England was initially covered by an unusually warm, dry air mass. Dynamical anticyclogenesis over eastern Canada set the stage for a favorable pressure gradient to allow chilly marine air to approach coastal New England from the east. Diabatic cooling over the chilly (5°–8°C) waters of the Gulf of Maine allowed surface pressures to remain relatively high offshore while diabatic heating over the land (31°–33°C temperatures) enabled surface pressures to fall relative to over the ocean. The resulting higher pressures offshore resulted in an onshore cold push. Frontal intensity was likely enhanced prior to leaf out and grass green-up as virtually all of the available insolation went into sensible heating. The large-scale environment in the February–March 1972 case favored the accumulation of bitterly cold arctic air in Canada. Frontal formation occurred over northern Montana and North Dakota as the arctic air moved slowly southward in conjunction with surface pressure rises east of the Canadian Rockies. The arctic air accelerated southward subsequent to lee cyclogenesis–induced pressure falls ahead of an upstream trough that crossed the Rockies. The southward acceleration of the arctic air was also facilitated by dynamic anticyclogenesis in southern Canada beneath a poleward jet-entrance region. Frontal intensity varied diurnally in response to differential diabatic heating. Three types of cyclogenesis events were observed over the lifetime of the event: 1) low-amplitude frontal waves with no upper-level support, 2) low-amplitude frontal waves that formed in a jet-entrance region, and 3) cyclones that formed ahead of advancing upper-level troughs. All cyclones were either nondeveloping or weak developments despite extreme baroclinicity, likely the result of large atmospheric static stability in the arctic frontal zone and unfavorable alongfront stretching deformation. Significant frontal–mountain interactions were observed over the Rockies and the Appalachians.

Description: Over 50 yr have passed since the publication of Sanders' 1955 study, the first quantitative study of the structure and dynamics of a surface cold front. The purpose of this chapter is to reexamine some of the results of that study in light of modern methods of numerical weather prediction and diagnosis. A simulation with a resolution as high as 6-km horizontal grid spacing was performed with the fifth-generation-Pennsylvania State University-National Center for Atmospheric Research (PSU-NCAR) Mesoscale Model (MM5), given initial and lateral boundary conditions from the National Centers for Environmental Precipitation-National Center for Atmospheric Research (NCEP-NCAR) reanalysis project data from 17 to 18 April 1953. The MM5 produced a reasonable simulation af the front, albeit its strength was not as intense and its movement was not as fast as was analyzed by Sanders. The vertical structure of the front differed from that analyzed by Sanders in several significant ways. First, the strongest horizontal temperature gradient associated with the cold front in the simulation occurred above a surface-based inversion, not at the earth's surface. Second, the ascent plume at the leading edge of the front was deeper and more intense than that analyzed by Sanders. The reason was an elevated mixed layer that had moved over the surface cold front in the simulation, allowing a much deeper vertical circulation than was analyzed by Sanders. This structure is similar to that of Australian cold fronts with their deep, well-mixed, prefrontal surface layer. These two differences between the model simulation and the analysis by Sanders may be because upper-air data from Fort Worth, Texas, was unavailable to Sanders. Third, the elevated mixed layer also meant that isentropes along the leading edge of the front extended vertically. Fourth, the field of frontogenesis of the horizontal temperature gradient calculated from the three-dimensional wind differed in that the magnitude of the maximum of the deformation term was larger than the magnitude of the maximum of the tilting term in the simulation, in contrast to Sanders' analysis and other previously published cases. These two discrepancies may be attributable to the limited horizontal resolution of the data that Sanders used in constructing his cross section. Last, a deficiency of the model simulation was that the postfrontal surface superadiabatic layer in the model did not match the observed well-mixed boundary layer. This result raises the question of the origin of the well-mixed postfrontal boundary layer behind cold fronts. To address this question, an additional model simulation without surface fluxes was performed, producing a well-mixed, not superadiabatic, layer. This result suggests that surface fluxes were not necessary for the development of the well-mixed layer, in agreement with previous research. Analysis of this event also amplifies two research themes that Sanders returned to later in his career, First, a prefrontal wind shift occurred in both the observations and model simulation at stations in western Oklahoma. This prefrontal wind shift was caused by a lee cyclone departing the leeward slopes of the Rockies slightly equatorward of the cold front, rather than along the front as was the case farther eastward. Sanders' later research showed how the occurrence of these prefrontal wind shifts leads to the weakening of fronts. Second, this study shows the advantage of using surface potential temperature, rather than surface temperature, for determining the locations of the surface fronts on sloping terrain.

Description: Today, even with state-of-the-art observational, data assimilation, and modeling systems run routinely on supercomputers, there are often surprises in the prediction of snowstorms, especially the “big ones,” affecting coastal regions of the mid-Atlantic and northeastern United States. Little did the author know that lessons from Fred Sanders' synoptic meteorology class at the Massachusetts Institute of Technology (1967) would later (late 1980s) inspire him to pursue practical issues of predictability in the context of the development of ensemble prediction systems, strategies, and applications for providing information on the inevitable case-dependent uncertainties in forecasts. This paper is a brief qualitative and somewhat colloquial overview, based upon this author's personal involvement and experiences, intended to highlight some basic aspects of the source and nature of uncertainties in forecasts and to illustrate the sort of value added information ensembles can provide in dealing with uncertainties in predictions of East Coast snowstorms.

Description: One characteristic of Fred Sanders' research is his ability to take a topic that is believed to be well understood by the research community and show that interesting research problems still exist. Among Sanders' considerable contributions to synoptic meteorology, those concerned with surface cold fronts have been especially influential. After a brief historical review of fronts and frontal analysis, this chapter presents three stages in Sanders' career when he performed research on the structure, dynamics, and analysis of surface cold fronts. First, his 1955 paper, "An investigation of the structure and dynamics of an intense surface frontal zone," was the first study to discuss quantitatively the dynamics of a surface cold front. In the 1960s, Sanders and his students further examined the structure of cold fronts, resulting in the unpublished 1967 report to the National Science Foundation, "Frontal structure and the dynamics of frontogenesis." For a third lime in his career, Sanders published several papers (1995–2005) revisiting the structure and dynamics of cold fronts. His 1967 and 1995–2005 work raises the question of the origin and dynamics of the surface pressure trough and/or wind shift that sometimes precedes the temperature gradient (hereafter called a prefrontal trough or prefrontal wind shift, respectively). Sanders showed that the relationship between this prefrontal feature and the temperature gradient is fundamental to the strength of the front. When the wind shift is coincident with the temperature gradient, frontogenesis (strengthening of the front) results; when the wind shift lies ahead of the temperature gradient, frontolysis (weakening of the front) results. a number of proposed mechanisms for the formation of prefrontal troughs and prefrontal wind shifts exist. Consequently, much research remains to be performed on these topics.

Description: When a forecast is assessed, a single value for a verification measure is often quoted. This is of limited use, as it needs to be complemented by some idea of the uncertainty associated with the value. If this uncertainty can be quantified, it is then possible to make statistical inferences based on the value observed. There are two main types of inference: confidence intervals can be constructed for an underlying “population” value of the measure, or hypotheses can be tested regarding the underlying value. This paper will review the main ideas of confidence intervals and hypothesis tests, together with the less well known “prediction intervals,” concentrating on aspects that are often poorly understood. Comparisons will be made between different methods of constructing confidence intervals—exact, asymptotic, bootstrap, and Bayesian—and the difference between prediction intervals and confidence intervals will be explained. For hypothesis testing, multiple testing will be briefly discussed, together with connections between hypothesis testing, prediction intervals, and confidence intervals.

Description: Studies using International Satellite Cloud Climatology Project (ISCCP) data have reported decreases in cloud optical depth with increasing temperature, thereby suggesting a positive feedback in cloud optical depth as climate warms. The negative cloud optical depth and temperature relationships are questioned because ISCCP employs threshold assumptions to identify cloudy pixels that have included partly cloudy pixels. This study applies the spatial coherence technique to one month of Advanced Very High Resolution Radiometer (AVHRR) data over the Pacific Ocean to differentiate overcast pixels from the partly cloudy pixels and to reexamine the cloud optical depth–temperature relationships. For low-level marine stratus clouds studied here, retrievals from partly cloudy pixels showed 30%–50% smaller optical depths, 1°–4°C higher cloud temperatures, and slightly larger droplet effective radii, when they were compared to retrievals from the overcast pixels. Despite these biases, retrievals for the overcast and partly cloudy pixels show similar negative cloud optical depth–temperature relationships and their magnitudes agree with the ISCCP results for the midlatitude and subtropical regions. There were slightly negative droplet effective radius–temperature relationships, and considerable positive cloud liquid water content–temperature relationships indicated by aircraft measurements. However, cloud thickness decreases appear to be the main reason why cloud optical depth decreases with increasing temperature. Overall, cloud thickness thinning may explain why similar negative cloud optical depth–temperature relationships are found in both overcast and partly cloudy pixels. In addition, comparing the cloud-top temperature to the air temperature at 740 hPa indicates that cloud-top height generally rises with warming. This suggests that the cloud thinning is mainly due to the ascending of cloud base. The results presented in this study are confined to the midlatitude and subtropical Pacific and may not be applicable to the Tropics or other regions.

Description: A statistical model to analyze different time scales of the variability of extreme high sea levels is presented. This model uses a time-dependent generalized extreme value (GEV) distribution to fit monthly maxima series and is applied to a large historical tidal gauge record (San Francisco, California). The model allows the identification and estimation of the effects of several time scales—such as seasonality, interdecadal variability, and secular trends—in the location, scale, and shape parameters of the probability distribution of extreme sea levels. The inclusion of seasonal effects explains a large amount of data variability, thereby allowing a more efficient estimation of the processes involved. Significant correlation with the Southern Oscillation index and the nodal cycle, as well as an increase of about 20% for the secular variability of the scale parameter have been detected for the particular dataset analyzed. Results show that the model is adequate for a complete analysis of seasonal-to-interannual sea level extremes providing time-dependent quantiles and confidence intervals.

Description: The interannual and intraseasonal variability of West African vegetation over the period 1982–2002 is studied using the normalized difference vegetation index (NDVI) from the Advanced Very High Resolution Radiometer (AVHRR).The novel independent component analysis (ICA) technique is applied to extract the main modes of the interannual variability of the vegetation, among which two modes are worth describing. The first component (IC1) describes NDVI variability over the Sahel from August to October. A strong photosynthetic activity over the Sahel is related to above-normal convection and rainfall within the intertropical convergence zone (ITCZ) in summertime and is partly associated with colder (warmer) SST in the eastern tropical Pacific (the Mediterranean). The second component (IC2) depicts a dipole pattern between the Sahelian and Guinean regions during the northern summer followed by a southward-propagating signal from October to December. It is associated with a north–south dipole in convection and rainfall induced by variations in the latitudinal location of the ITCZ as a response to the occurrence of the tropical Atlantic dipole.The analysis of the intraseasonal variability of the Sahelian vegetation relies on the analysis of the seasonal marches and their main phenological stages. Green-up usually starts in early July and shows a very low year-to-year variability, while senescence ends by mid-November and is prone to larger interannual variability. Six types of vegetative seasonal marches are discriminated according to variations in the timing of phenological stages as well as in the greening intensity. These types appear to be strongly dependent on rainfall distribution and amount, particularly those recorded in late August. Finally, year-to-year memory effects are highlighted: NDVI recorded during the green-up phase in year j appears to be strongly related to the maximum NDVI value recorded at year j − 1.

Description: Based on the data of optimum interpolation sea surface temperature (OISST), the temporal correlations of the sea surface temperature anomaly (SSTA) in the South China Sea (SCS) are studied by using the rescaled range analysis (R/S) and detrended fluctuation analysis (DFA). The results show that the scaling exponents of SSTAs are larger than 0.8. This finding indicates that the SSTAs in the SCS exhibit persistent long-range time correlation of the fluctuations and the interval spreads over a wide period, from about 1 month to 4.5 yr (4∼235 weeks). In addition, the “degree” of the correlations depends very much on the geographic locations: near to the coastal regions, the value is small, while far from the coastline, the value is relatively larger. This means that SSTAs in the central SCS are smoother than those of the coastal regions. The persistence of SST in the SCS may be used as a “minimum skill” to assess the ocean models and to evaluate their performance.

Description: Winds at the Salt Lake City International Airport (SLC) during the April–October period from 1948 to 2003 have been observed to shift to the north (up-valley direction) between late morning and afternoon on over 70% of the days without precipitation. Lake-breeze fronts that develop as a result of the differential heating between the air over the nearby Great Salt Lake and that over the lake’s surroundings are observed at SLC only a few times each month. Fewer lake-breeze fronts are observed during late July–early September than before or after that period. Interannual fluctuations in the areal extent of the shallow Great Salt Lake contribute to year-to-year variations in the number of lake-breeze frontal passages at SLC. Data collected during the Vertical Transport and Mixing Experiment (VTMX) of October 2000 are used to examine the structure and evolution of a lake-breeze front that moved through the Salt Lake Valley on 17 October. The onset of upslope and up-valley winds occurred within the valley prior to the passage of the lake-breeze front. The lake-breeze front moved at roughly 3 m s−1 up the valley and was characterized near the surface by an abrupt increase in wind speed and dewpoint temperature over a distance of 3–4 km. Rapid vertical mixing of aerosols at the top of the 600–800-m-deep boundary layer was evident as the front passed.

Description: This study combines the experimental measurements with large-eddy simulation (LES) data of a neutral planetary boundary layer (PBL) documented by a 60-m tower instrumented with eight sonic anemometers, and a high-resolution Doppler lidar during the 1999 Cooperative Atmospheric and Surface Exchange Study (CASES-99) experiment. The target of the paper is (i) to investigate the multiscale nature of the turbulent eddies in the surface layer (SL), (ii) to explain the existence of a −1 power law in the velocity fluctuation spectra, and (iii) to investigate the different nature of turbulence in the two sublayers within the SL, which are the eddy surface layer (ESL; lower sublayer of the SL lying between the surface and about 20-m height) and the shear surface layer (SSL; lying between the ESL top and the SL top). The sonic anemometers and Doppler lidar prove to be proper instruments for LES validation, and in particular, the Doppler lidar because of its ability to map near-surface eddies.This study shows the different nature of turbulence in the ESL and the SSL in terms of organized eddies, velocity fluctuation spectra, and second-order moment statistics. If quantitative agreement is found in the SSL between the LES and the measurements, only qualitative similarity is found in the ESL due to the subgrid-scale model, indicating that the LES captures part of the physics of the ESL. The LES helps explain the origin of the −1 power-law spectral subrange evidence in the velocity fluctuation spectra computed in the SL using the CASES-99 dataset: strong interaction between the mean flow and the fluctuating vorticities is found in the SL, and turbulent production is found to be larger than turbulent energy transfer for k1z 〉 1 (k1 being the longitudinal wavenumber and z the height), which is the condition of the −1 power-law existence.

Description: The effects of natural and anthropogenic heterogeneity on a hydrological simulation are evaluated using a distributed biosphere hydrological model (DBHM) system. The DBHM embeds a biosphere model into a distributed hydrological scheme, representing both topography and vegetation in a mesoscale hydrological simulation, and the model system includes an irrigation scheme. The authors investigated the effects of two kinds of variability, precipitation variability and the variability of irrigation redistributing runoff, representing natural and anthropogenic heterogeneity, respectively, on hydrological processes. Runoff was underestimated if rainfall was placed spatially uniformly over large grid cells. Accounting for precipitation heterogeneity improved the runoff simulation. However, the negative runoff contribution, namely, the situation that mean annual precipitation is less than evapotranspiration, cannot be simulated by only considering the natural heterogeneity. This constructive model shortcoming can be eliminated by accounting for anthropogenic heterogeneity caused by irrigation water withdrawals. Irrigation leads to increased evapotranspiration and decreased runoff, and surface soil moisture in irrigated areas increases because of irrigation. Simulations performed for the Yellow River basin of China indicated streamflow decreases of 41% due to irrigation effects. The latent heat flux in the peak irrigation season [June–August (JJA)] increased 3.3 W m−2 with a decrease in the ground surface temperature of 0.1 K for the river basin. The maximum simulated increase in the latent heat flux was 43 W m−2, and the ground temperature decrease was 1.6 K in the peak irrigation season.

Description: A method for routinely verifying numerical weather prediction surface marine winds with satellite scatterometer winds is introduced. The marine surface winds from the Australian Bureau of Meteorology’s operational global and regional numerical weather prediction systems are evaluated. The model marine surface layer is described. Marine surface winds from the global and limited-area models are compared with observations, both in situ (anemometer) and remote (scatterometer). A 2-yr verification shows that wind speeds from the regional model are typically underestimated by approximately 5%, with a greater bias in the meridional direction than the zonal direction. The global model also underestimates the surface winds by around 5%–10%. A case study of a significant marine storm shows that where larger errors occur, they are due to an underestimation of the storm intensity, rather than to biases in the boundary layer parameterizations.

Description: The applicability of axisymmetric theory of angular momentum conserving circulations to the large-scale steady monsoon is studied in a general circulation model with idealized representations of continental geometry and simple physics. Results from an aquaplanet setup with localized subtropical forcing are compared with a continental case. It is found that the meridional circulation that develops is close to angular momentum conserving for cross-equatorial circulation cells, both in the aquaplanet and in the continental cases. The equator proves to be a substantial barrier to boundary layer meridional flow; flow into the summer hemisphere from the winter hemisphere tends to occur in the free troposphere rather than in the boundary layer. A theory is proposed to explain the location of the monsoon; assuming quasiequilibrium, the poleward boundary of the monsoon circulation is collocated with the maximum in subcloud moist static energy, with the monsoon rains occurring near and slightly equatorward of this maximum. The model results support this theory of monsoon location, and it is found that the subcloud moist static energy distribution is determined by a balance between surface forcing and advection by the large-scale flow.

Description: This paper concerns the calculation of the probability of exceedance of wave crest elevation. The statistics have been calculated for broadbanded, unidirectional, deep-water sea states by incorporating a fully nonlinear wave model into a spectral response surface method. This is a novel approach to the calculation of statistics and, as all of the calculations are performed in the probability domain, avoids the need for long time-domain simulations. Furthermore, in contrast to theoretical distributions, the broadbanded, fully nonlinear nature of the sea state can be considered. The results have been compared with those of fully nonlinear time-domain simulations and are shown to be in good agreement. The results indicate that in unidirectional seas the crest elevations of the largest waves can be much higher than would be predicted by linear or second-order theory. Hence, the occurrence of locally long crested sea states offers a possible explanation for the formation of freak or rogue waves.

Description: Nearly one-half of the earth’s terrestrial surface is susceptible to drought, which can have significant social, economic, and environmental impacts. Therefore, it is important to develop better descriptions and models of the processes linking the land surface and atmosphere during drought. Using data collected during the International H2O Project, the study presented here investigates the effects of variations in the environmental factors driving the latent heat flux (λE) during drought conditions at a rangeland site located in the panhandle of Oklahoma. Specifically, this study focuses on the relationships of λE with vapor pressure deficit, wind speed, net radiation, soil moisture content, and greenness fraction. While each of these environmental factors has an influence, soil moisture content is the key control on λE. The role of soil moisture in regulating λE is explained in terms of the surface resistance to water vapor transfer. The results show that λE transitioned between being water or energy limited during the course of the drought. The implications of this on the ability to understand and model drought conditions and transitions into or out of droughts are discussed.

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: In March 2003 several autonomous underwater vehicle (AUV) missions were carried out under sea ice in the western Bellingshausen Sea. Data from the upward-looking acoustic Doppler current profiler (ADCP) on the “Autosub” AUV indicate a strongly oscillating horizontal velocity of the ice due to ocean swell. Swell period, height, direction, and directional spread have been computed every 800 m from the ice edge to 10 km inward for three missions. Exponential, period-dependent attenuation of waves propagating through sea ice was observed. Mean period increased with distance from the ice edge. The wave field refracted during propagation. The directional wave spread does not seem to relate to distance from the ice edge, although higher frequencies tended to be more spread. If suitably deployed, an ordinary ADCP may be used with this technique to study both scalar and directional properties of waves in open or ice-covered water.

Description: Wind-sea generation was observed during two experiments off the coast of North Carolina. One event with offshore winds of 9–11 m s−1 directed 20° from shore normal was observed with eight directional stations recording simultaneously and spanning a fetch from 4 to 83 km. An opposing swell of 1-m height and 10-s period was also present. The wind-sea part of the wave spectrum conforms to established growth curves for significant wave height and peak period, except at inner-shelf stations where a large alongshore wind-sea component was observed. At these short fetches, the mean wave direction θm was observed to change abruptly across the wind-sea spectral peak, from alongshore at lower frequencies to downwind at higher frequencies. Waves from another event with offshore winds of 6–14 m s−1 directed 20°–30° from shore normal were observed with two instrument arrays. A significant amount of low-frequency wave energy was observed to propagate alongshore from the region where the wind was strongest. These measurements are used to assess the performance of some widely used parameterizations in wave models. The modeled transition of θm across the wind-sea spectrum is smoother than that in the observations and is reproduced very differently by different parameterizations, giving insights into the appropriate level of dissipation. Calculations with the full Boltzmann integral of quartet wave–wave interactions reveal that the discrete interaction approximation parameterization for these interactions is reasonably accurate at the peak of the wind sea but overpredicts the directional spread at high frequencies. This error is well compensated by parameterizations of the wind input source term that have a narrow directional distribution. Observations also highlight deficiencies in some parameterizations of wave dissipation processes in mixed swell–wind-sea conditions.

Description: Seven sets of 2D particle image velocimetry data obtained in the bottom boundary layer of the coastal ocean along the South Carolina and Georgia coast [at the South Atlantic Bight Synoptic Offshore Observational Network (SABSOON) site] are examined, covering the accelerating and decelerating phases of a single tidal cycle at several heights above the seabed. Additional datasets from a previous deployment are also included in the analysis. The mean velocity profiles are logarithmic, and the vertical distribution of Reynolds stresses normalized by the square of the free stream velocity collapse well for data obtained at the same elevation but at different phases of the tidal cycle. The magnitudes of 〈u′u′〉, 〈w′w′〉, and −〈u′w′〉 decrease with height above bottom in the 25–160-cm elevation range and are consistent with the magnitudes and trends observed in laboratory turbulent boundary layers. If a constant stress layer exists, it is located below 25-cm elevation. Two methods for estimating dissipation rate are compared. The first, a direct estimate, is based on the measured in-plane instantaneous velocity gradients. The second method is based on fitting the resolved part of the dissipation spectrum to the universal dissipation spectrum available in Gargett et al. Being undervalued, the direct estimates are a factor of 2–2.5 smaller than the spectrum-based estimates. Taylor microscale Reynolds numbers for the present analysis range from 24 to 665. Anisotropy is present at all resolved scales. At the transition between inertial and dissipation range the longitudinal spectra exhibit a flatter than −5/3 slope and form spectral bumps. Second-order statistics of the velocity gradients show a tendency toward isotropy with increasing Reynolds number. Dissipation exceeds production at all measurement heights, but the difference varies with elevation. Close to the bottom, the production is 40%–70% of the dissipation, but it decreases to 10%–30% for elevations greater than 80 cm.

Description: Scalar exchange between San Francisco Bay and the coastal ocean is examined using shipboard observations made across the Golden Gate Channel. The study consists of experiments during each of the following three “seasons”: winter/spring runoff (March 2002), summer upwelling (July 2003), and autumn relaxation (October 2002). Within each experiment, transects across the channel were repeated approximately every 12 min for 25 h during both spring and neap tides. Velocity was measured from a boat-mounted ADCP. Scalar concentrations were measured at the surface and from a tow-yoed SeaSciences Acrobat. Net salinity exchange rates for each season are quantified with harmonic analysis. Accuracy of the net fluxes is confirmed by comparison with independently measured values. Harmonic results are then decomposed into flux mechanisms using temporal and spatial correlations. In this study, the temporal correlation of cross-sectionally averaged salinity and velocity (tidal pumping flux) is the largest part of the dispersive flux of salinity into the bay. From the tidal pumping flux portion of the dispersive flux, it is shown that there is less exchange than was found in earlier studies. Furthermore, tidal pumping flux scales strongly with freshwater flow resulting from the density-driven movement of a tidally trapped eddy and stratification-induced increases in ebb–flood frictional phasing. Complex bathymetry makes salinity exchange scale differently with flow than would be expected from simple tidal asymmetry and gravitational circulation models.

Description: As part of a program aimed at developing a long-duration, subsurface mooring, known as Ultramoor, several modern acoustic current meters were tested. The instruments with which the authors have the most experience are the Aanderaa RCM11 and the Nortek Aquadopp, which measure currents using the Doppler shift of backscattered acoustic signals, and the Falmouth Scientific ACM, which measures changes in travel time of acoustic signals between pairs of transducers. Some results from the Doppler-based Sontek Argonaut and the travel-time-based Nobska MAVS are also reported. This paper concentrates on the fidelity of the speed measurement but also presents some results related to the accuracy of the direction measurement. Two procedures were used to compare the instruments. In one, different instruments were placed close to one another on three different deep-ocean moorings. These tests showed that the RCM11 measures consistently lower speeds than either a vector averaging current meter or a vector measuring current meter, both more traditional instruments with mechanical velocity sensors. The Aquadopp in use at the time, but since updated to address accuracy problems in low scattering environments, was biased high. A second means of testing involved comparing the appropriate velocity component of each instrument with the rate of change of pressure when they were lowered from a ship. Results from this procedure revealed no depth dependence or measurable bias in the RCM11 data, but did show biases in both the Aquadopp and Argonaut Doppler-based instruments that resulted from low signal-to-noise ratios in the clear, low scattering conditions beneath the thermocline. Improvements in the design of the latest Aquadopp have reduced this bias to a level that is not significant.

Description: The energy pathways in geostrophic turbulence are explored using a two-layer, flat-bottom, f-plane, quasigeostrophic model forced by an imposed, horizontally homogenous, baroclinically unstable mean flow and damped by bottom Ekman friction. A systematic presentation of the spectral energy fluxes, the mean flow forcing, and dissipation terms allows for a comprehensive understanding of the sources and sinks for baroclinic and barotropic energy as a function of length scale. The key new result is a robust inverse cascade of kinetic energy for both the baroclinic mode and the upper layer. This is consistent with recent observations of satellite altimeter data over the South Pacific Ocean. The well-known forward cascade of baroclinic potential and total energy was found to be very robust. Decomposing the spectral fluxes into contributions from different terms provided further insight. The inverse baroclinic kinetic energy cascade is driven mostly by an efficient interaction between the baroclinic velocity and the barotropic vorticity, the latter playing a crucial catalytic role. This cascade can be further enhanced by the baroclinic mode self-interaction, which is only present with nonuniform stratification (unequal layer depths). When model parameters are set such that modeled eddies compare favorably with observations, the inverse baroclinic kinetic energy cascade is actually much stronger than the well-known inverse cascade in the barotropic mode. The upper-layer kinetic energy cascade was found to dominate the lower-layer cascade over a wide range of parameters, suggesting that the surface cascade and time mean density stratification may be sufficient for estimating the depth-integrated cascade from ocean observations. This may find useful application in inferring the kinetic to gravitational potential energy conversion rate from satellite measurements.

Description: A statistical forecast model, referred to as the snow-cast (sCast) model, has been developed using observed October mean snow cover and sea level pressure anomalies to predict upcoming winter land surface temperatures for the extratropical Northern Hemisphere. In operational forecasts since 1999, snow cover has been used for seven winters, and sea level pressure anomalies for three winters. Presented are skill scores for these seven real-time forecasts and also for 33 winter hindcasts (1972/73–2004/05). The model demonstrates positive skill over much of the eastern United States and northern Eurasia—regions that have eluded skillful predictions among the existing major seasonal forecast centers. Comparison with three leading dynamical forecast systems shows that the statistical model produces superior skill for the same regions. Despite the increasing complexity of the dynamical models, they continue to derive their forecast skill predominantly from tropical atmosphere–ocean coupling, in particular from ENSO. Therefore, in the Northern Hemisphere extratropics, away from the influence of ENSO, the sCast model is expected to outperform the dynamical models into the foreseeable future.

Description: The AmeriFlux network continues to improve the understanding of carbon, water, and energy fluxes across temporal and spatial scales. The network includes ∼120 research sites that contribute to the understanding of processes within and among ecosystems. To improve the networks ability and confidence to synthesize data across multiple sites, the AmeriFlux quality assurance and quality control laboratory was established to reduce the within- and among-site uncertainties. This paper outlines the design of the portable eddy covariance system (PECS) and subsequent data processing procedures used for site comparisons. Because the PECS makes precision measurements of atmospheric CO2, the authors also present the results of uncertainty analyses in determining the polynomials for an infrared gas analyzer, estimating the CO2 in secondary standards, and estimating ambient CO2 in field measurements. Under field conditions, drift in the measurement of CO2 increased the uncertainty in flux measurements across 7 days by 5% and was not dependent on the magnitude or direction of the flux. The maximum relative flux measurement error for unstable conditions was 10.03 μmol CO2 m−2 s−1.

Description: Observations of the turbulent exchange between a river surface and the atmosphere in a mountainous area in southern Brazil are presented and discussed. A micrometeorological tower was installed directly above the surface of a 60-m-wide river. This paper describes the observed turbulent fluxes over 12 days of observations at this site. Eddy correlation sensible and latent heat fluxes are directed toward the river during daytime and from the river at night, and they are controlled by differences between water and air temperatures. The magnitude of the vertical fluxes between the river and the atmosphere increases during daytime with increasing temperature gradient up to a threshold, beyond which the increasing stability starts to dampen the fluxes. Water and air temperatures show very little variations across the width of the river, indicating that the measurements taken at one margin may be representative of the mean river exchange. Local scalar budgets show that daytime warming and moistening rates above the river are controlled by local transport from the riverbanks. The main vertical fluxes have a very small magnitude: 0.8 W m−2 for sensible heat and 1.1 W m−2 for latent heat. Events of very large sensible heat fluxes from the river to the atmosphere and very large latent heat fluxes from the atmosphere to the river happened on 3 days, following nights with a very deep fog layer in the valley. These events represented the passage of a warm and dry air mass down the river. A process to explain the occurrence of these large fluxes is suggested that is associated with differential fog dissipation over the valley.

Description: Following recent results showing the potential for using surface observations of temperature, water vapor mixing ratio, and winds to determine PBL profiles, this paper reports on experiments with real observations. A 1D column model with soil, surface-layer, and PBL parameterization schemes that are the same as in the Weather Research and Forecasting model is used to estimate PBL profiles with an ensemble filter. Surface observations over the southern Great Plains are assimilated during the spring and early summer period of 2003. To strictly quantify the utility of the observations for determining PBL profiles in the ensemble filter framework, only climatological information is provided for initialization and forcing. The analysis skill, measured against rawinsondes for an independent verification, is compared against climatology to quantify the influence of the observations. Sensitivity to changing parameterization schemes, and to prescribed values of observation error variance, is examined. Temporal propagation of skillful analyses is also assessed, separating the effects of good prior state estimates from the impact of assimilation at night when covariance is weak. Results show that accurate profiles of temperature, mixing ratio, and winds are estimated with the column model and ensemble filter assimilating only surface observations. Results are largely insensitive to choice of parameterization scheme and specified observation error variance. The effects of using different parameterization schemes within the column model depend on whether assimilation is included, showing the importance of evaluating models within assimilation systems. At night, skillful estimates are possible because the influence of the observations from the previous day is temporally propagated, and atmospheric dynamics in the residual layer operate on slow time scales. It is expected that these profiles will have applications for nowcasting and secondary models (e.g., plume dispersion models) that rely on accurate specification of PBL structure.

Description: Continuous time series of soil water content over a period of more than 9 months for a midlatitude sandy loam soil covered by grass are calculated with the Campbell and the van Genuchten soil hydraulic functions and the Clapp–Hornberger, Cosby et al., and Rawls–Brakensiek parameter sets. The results are compared with soil water content observed at several soil depths, and the water balance components are evaluated. The Campbell soil hydraulic functions are often used by meteorologists, whereas the van Genuchten functions are widespread among hydrologists. The simulations are performed with the “VEG3D” soil–vegetation model in stand-alone mode forced by on-site meteorological observations. The soil water content and meteorological observations were obtained within the Regional Climate Project (REKLIP) at a site in the Rhine valley in southern Germany with 10-min temporal resolution. Apart from the different soil hydraulic functions and parameter sets, the effects of different lower boundary conditions and initializations on the simulations are compared in terms of statistical quantities like mean error, bias, correlation coefficient, and least squares fit. Large differences between the various combinations are found. For the situation considered in this paper, the van Genuchten–Clapp–Hornberger, the Campbell–Cosby et al., and the van Genuchten–Rawls–Brakensiek combinations give the best overall agreement with the observations.

Description: The extremes of near-surface temperature and 24-h and 5-day mean precipitation rates are examined in simulations performed with atmospheric general circulation models (AGCMs) participating in the second phase of the Atmospheric Model Intercomparison Project (AMIP-2). The extremes are evaluated in terms of 20-yr return values of annual extremes. The model results are validated against the European Centre for Medium-Range Weather Forecasts and National Centers for Environmental Prediction reanalyses and station data. Precipitation extremes are also validated against the pentad dataset of the Global Precipitation Climatology Project, which is a blend of rain gauge observations, satellite data, and model output. On the whole, the AGCMs appear to simulate temperature extremes reasonably well. Model disagreements are larger for cold extremes than for warm extremes, particularly in wet and cloudy regions, and over sea ice and snow-covered areas. Many models exhibit an exaggerated clustering behavior for temperatures near the freezing point of water. Precipitation extremes are less reliably reproduced by the models and reanalyses. The largest disagreements are found in the Tropics where the parameterizations of deep convection affect the simulated daily precipitation extremes.

Description: Climatic and hydrologic observations and results from a terrestrial ecosystem model coupled to a regional-scale river-routing algorithm are used to document the associations between the El Niño–Southern Oscillation (ENSO) phenomenon and anomalies in climate, surface water balance, and river hydrology within the Mississippi River basin. While no ENSO signal is found in streamflow at the outlet of the basin in Vicksburg, Mississippi, significant anomalies in all water balance components are found in certain regions within the basin. ENSO is mainly associated with positive winter temperature anomalies, but hydrologic patterns vary with season, location, and ENSO phase. El Niño precipitation anomalies tend to affect evapotranspiration (ET) in the western half of the basin and runoff in the eastern half. La Niña events are associated with ET anomalies in the central portion of the basin and runoff anomalies in the southern and eastern portions of the basin. Both ENSO phases are associated with decreased snow depth. Anomalous soil moisture patterns occur at seasonal time scales and filter noisier spatial patterns of precipitation anomalies into coherent patterns with larger field significance; however, for all water budget components, there is a large amount of variability in response within a particular ENSO phase. With anomalies that are up to 4 times those of a typical event, it is clear that improved predictability of the onset and strength of an upcoming ENSO event is important for both water resource management and disaster mitigation.

Description: The extreme phases of El Niño–Southern Oscillation (ENSO) are known to dominate the interannual variability of tropical rainfall. However, the relationship between ENSO and the spatial extent of drought and excessively wet conditions is an important characteristic of the tropical climate that has received relatively less attention from researchers. Here, a standardized precipitation index is computed from monthly rainfall analyses and the temporal variability of the spatial extent of such extremes, for various levels of severity, is examined from a Tropics-wide perspective (land areas only, 30°S–30°N). Maxima in the spatial extent of both precipitation extremes are compared across multiple ENSO events that occurred during the period 1950–2003. The focus on tropical land areas is motivated by the numerous, often negative, impacts of ENSO-related precipitation variability on human populations. Results show that major peaks in the spatial extent of drought and excessively wet conditions are generally associated with extreme phases of ENSO. A remarkably robust linear relationship is documented between the spatial extent of drought in the Tropics and El Niño strength (based on Niño-3.4 sea surface temperature anomalies), with a comparatively weaker relationship for La Niña and excessive wetness. Both conditions are found to increase by about a factor of 2 between strong and weak ENSO events, and in several locations they are shown to be more likely during ENSO events than at all other times, especially for severe categories. Relatively stronger El Niño events during recent decades are associated with increased drought extent in tropical land areas with increasing surface temperatures likely acting to exacerbate these dry conditions.

Description: A currently popular idea is that El Niño–Southern Oscillation (ENSO) can be viewed as a linear deterministic system forced by noise representing processes with periods shorter than ENSO. Also, there is observational evidence to suggest that the Madden–Julian oscillation (MJO) acts to trigger and/or amplify the warm phase of ENSO in this way. The feedback of the slower process, ENSO, to higher-frequency atmospheric phenomena, of which a large part of the variability in the intraseasonal band is due to the MJO, has received little attention. This paper considers the hypothesis that the probability of an El Niño event is modified by high MJO activity and that, in turn, the MJO is regulated by ENSO activity. If this is indeed the case, then viewing ENSO as a low-frequency oscillation forced by additive stochastic noise would not present a complete picture. This paper tests the above hypothesis using a stochastically forced intermediate coupled model by allowing ENSO to directly influence the stochastic forcing. The model response to a variety of stochastic forcing types is found to be sensitive to the type of forcing applied. When the model is operated beyond its intrinsic Hopf bifurcation, its probability distribution function (PDF) is fundamentally altered when the stochastic forcing is changed from additive to multiplicative. The model integration period also influences the shape of the PDF, which is also compared to the PDF derived from observations. It is found that multiplicative stochastic forcing reproduces some measures of the observations better than the additive stochastic forcing.

Description: Set cover models are used to develop two reference station networks that can serve as near-term substitutes (as well as long-term backups) for the recently established Climate Reference Network (CRN) in the United States. The first network contains 135 stations distributed in a relatively uniform fashion in order to match the recommended spatial density for CRN. The second network contains 157 well-distributed stations that are generally not in urban areas in order to minimize the impact of future changes in land use. Both networks accurately reproduce the historical temperature and precipitation variations of the twentieth century.

Description: The impacts of the sea surface temperatures (SSTs) in the northern Gulf of California (GC) on warm-season rainfall in the Arizona–New Mexico (AZNM) and the northwestern Mexico (NWM) regions associated with the North American monsoon (NAM) are examined from two sets of seasonal simulations in which different SSTs were prescribed in the GC. The simulations reproduced important features in the low-level mesoscale circulations and upper air fields around the time of monsoon rainfall onset in AZNM such as sea-breeze-like diurnal variations in the low-level winds between the GC and the land, development of south-southeasterly winds over the GC and the western slope of the Sierra Madre Occidental after the onset of rainfall, and the strengthening of the 500-hPa high over AZNM around the onset of monsoon rainfall in AZNM. The simulated temporal variations in the upper air fields and daily rainfall, as well as the mesoscale circulation around the GC, suggest that the GC SSTs affect the water cycle around the GC mainly by altering mesoscale circulation and water vapor fluxes, but they have minimal impacts on the onset timing of monsoon rainfall in NWM and AZNM. With higher SSTs in the NGC, rainfall in NWM and AZNM increases in response to enhanced water vapor fluxes from the GC into the land. The enhanced onshore component of the low-level water vapor fluxes from the GC with higher GC SSTs results from two opposing effects: weakened sea-breeze-like circulation between the GC and the surrounding lands that tends to reduce the water vapor fluxes from the GC, and increased evaporation from the GC that tends to increase the water vapor fluxes. The simulations also suggest that the development of south-southeasterly low-level winds over the GC after monsoon rainfall onset plays an important role in enhancing rainfall as longer fetches over the GC can provide more water vapor into the low atmosphere.

Description: Previous studies show that the climatological precipitation over South America, particularly the Nordeste region, is influenced by the presence of the African continent. Here the influence of African topography and surface wetness on the Atlantic marine ITCZ (AMI) and South American precipitation are investigated. Cross-equatorial flow over the Atlantic Ocean introduced by north–south asymmetry in surface conditions over Africa shifts the AMI in the direction of the flow. African topography, for example, introduces an anomalous high over the southern Atlantic Ocean and a low to the north. This results in a northward migration of the AMI and dry conditions over the Nordeste region. The implications of this process on variability are then studied by analyzing the response of the AMI to soil moisture anomalies over tropical Africa. Northerly flow induced by equatorially asymmetric perturbations in soil moisture over northern tropical Africa shifts the AMI southward, increasing the climatological precipitation over northeastern South America. Flow associated with an equatorially symmetric perturbation in soil moisture, however, has a very weak cross-equatorial component and very weak influence on the AMI and South American precipitation. The sensitivity of the AMI to soil moisture perturbations over certain regions of Africa can possibly improve the skill of prediction.

Description: Recent studies of regime behavior in the extratropical variability have been based on a nonlinear extension to principal component analysis. Multimodality has been identified in the nonlinear principal component, and the multimodality has been interpreted as evidence for the existence of multiple circulation regimes. Here, multimodality is shown to be abundant in nonlinear principal component analysis when applied to sufficiently isotropic data even if these data are inherently unimodal. It is recommended that the nonlinear principal component analysis should not be used for detection of multimodality and regime behavior.

Description: The sensitivity of tropical atmospheric hydrologic processes to cloud microphysics is investigated using the NASA Goddard Earth Observing System (GEOS) general circulation model (GCM). Results show that a faster autoconversion rate leads to (a) enhanced deep convection in the climatological convective zones anchored to tropical land regions; (b) more warm rain, but less cloud over oceanic regions; and (c) an increased convective-to-stratiform rain ratio over the entire Tropics. Fewer clouds enhance longwave cooling and reduce shortwave heating in the upper troposphere, while more warm rain produces more condensation heating in the lower troposphere. This vertical differential heating destabilizes the tropical atmosphere, producing a positive feedback resulting in more rain and an enhanced atmospheric water cycle over the Tropics. The feedback is maintained via secondary circulations between convective tower and anvil regions (cold rain), and adjacent middle-to-low cloud (warm rain) regions. The lower cell is capped by horizontal divergence and maximum cloud detrainment near the freezing–melting (0°C) level, with rising motion (relative to the vertical mean) in the warm rain region connected to sinking motion in the cold rain region. The upper cell is found above the 0°C level, with induced subsidence in the warm rain and dry regions, coupled to forced ascent in the deep convection region. It is that warm rain plays an important role in regulating the time scales of convective cycles, and in altering the tropical large-scale circulation through radiative–dynamic interactions. Reduced cloud–radiation feedback due to a faster autoconversion rate results in intermittent but more energetic eastward propagating Madden–Julian oscillations (MJOs). Conversely, a slower autoconversion rate, with increased cloud radiation produces MJOs with more realistic westward-propagating transients embedded in eastward-propagating supercloud clusters. The implications of the present results on climate change and water cycle dynamics research are discussed.

Description: Equatorial Pacific sea surface temperature (SST) anomalies in the Center for Ocean–Land–Atmosphere Studies (COLA) interactive ensemble coupled general circulation model show near-annual variability as well as biennial El Niño–Southern Oscillation (ENSO) variability. There are two types of near-annual modes: a westward propagating mode and a stationary mode. For the westward propagating near-annual mode, warm SST anomalies are generated in the eastern equatorial Pacific in boreal spring and propagate westward in boreal summer. Consistent westward propagation is seen in precipitation, surface wind, and ocean current. For the stationary near-annual mode, warm SST anomalies develop near the date line in boreal winter and decay locally in boreal spring. Westward propagation of warm SST anomalies also appears in the developing year of the biennial ENSO mode. However, warm SST anomalies for the westward propagating near-annual mode occur about two months earlier than those for the biennial ENSO mode and are quickly replaced by cold SST anomalies, whereas warm SST anomalies for the biennial ENSO mode only experience moderate weakening. Anomalous zonal advection contributes to the generation and westward propagation of warm SST anomalies for both the westward propagating near-annual mode and the biennial ENSO mode. However, the role of mean upwelling is markedly different. The mean upwelling term contributes to the generation of warm SST anomalies for the biennial ENSO mode, but is mainly a damping term for the westward propagating near-annual mode. The development of warm SST anomalies for the stationary near-annual mode is partially due to anomalous zonal advection and upwelling, similar to the amplification of warm SST anomalies in the equatorial central Pacific for the biennial ENSO mode. The mean upwelling term is negative in the eastern equatorial Pacific for the stationary near-annual mode, which is opposite to the ENSO mode. The development of cold SST anomalies in the aftermath of warm SST anomalies for the westward propagating near-annual mode is coupled to large easterly wind anomalies, which occur between the warm and cold SST anomalies. The easterly anomalies contribute to the cold SST anomalies through anomalous zonal, meridional, and vertical advection and surface evaporation. The cold SST anomalies, in turn, enhance the easterly anomalies through a Rossby-wave-type response. The above processes are most effective during boreal spring when the mean near-surface-layer ocean temperature gradient is the largest. It is suggested that the westward propagating near-annual mode is related to air–sea interaction processes that are limited to the near-surface layers.

Description: One of the possible consequences of global warming is that there will be more days with precipitation throughout the year, and also that the level of precipitation will be higher. In this paper a detailed statistical analysis of a century of daily precipitation levels is provided for the central meteorological station in the Netherlands. This paper shows that the often-considered gamma distribution does not fit well to samples of yearly data. It is argued that its incorrect use can lead to spuriously high probabilities of extreme precipitation levels. Relying on advanced nonparametric techniques, it is first found that there are fewer rainy days in the central part of the Netherlands. Next, more rainy days involve higher precipitation levels. Most importantly, no statistically significant shift is found in the annual largest values of daily rainfall over the course of the century, which suggest that the probability of extremely high levels has not changed over time.

Description: The decadal trend behavior of the Northern Hemisphere atmospheric circulation is investigated utilizing long-term simulations with different state-of-the-art coupled general circulation models (GCMs) for present-day climate conditions (1990), reconstructions of the past 500 yr, and observations. The multimodel simulations show that strong positive winter North Atlantic Oscillation (NAO) trends are connected with the underlying sea surface temperature (SST) and exhibit an SST tripole trend pattern and a northward shift of the storm-track tail. Strong negative winter trends of the Aleutian low are associated with SST changes in the El Niño–Southern Oscillation (ENSO) region and a westward shift of the storm track in the North Pacific. The observed simultaneous appearance of strong positive NAO and negative Aleutian low trends is very unlikely to occur by chance in the unforced simulations and reconstructions. The positive winter NAO trend of the last 50 yr is not statistically different from the level of internal atmosphere–ocean variability. The unforced simulations also show a strong link between positive SST trends in the ENSO region and negative Aleutian low trends. With much larger observed SST trends in the ENSO region, this suggests that the observed negative Aleutian low trend is possibly influenced by external forcing, for example, global warming, volcanism, and/or solar activity change.

Description: Two widely used statistical approaches to reconstructing past climate histories from climate “proxy” data such as tree rings, corals, and ice cores are investigated using synthetic “pseudoproxy” data derived from a simulation of forced climate changes over the past 1200 yr. These experiments suggest that both statistical approaches should yield reliable reconstructions of the true climate history within estimated uncertainties, given estimates of the signal and noise attributes of actual proxy data networks.

Description: Modeling storm occurrences has become a vital part of hurricane prediction. In this paper, a method for simulating event occurrences using a simulated annealing algorithm is described. The method is illustrated using annual counts of hurricanes and of tropical storms in the Atlantic Ocean and Gulf of Mexico. Simulations closely match distributional properties, including possible correlations, in the historical data. For hurricanes, traditionally used Poisson and negative binomial processes also predict univariate properties well, but for tropical storms parametric methods are less successful. The authors determined that simulated annealing replicates properties of both series. Simulated annealing can be designed so that simulations mimic historical distributional properties to whatever degree is desired, including occurrence of extreme events and temporal patterning.

Description: This study presents the simulation of the Madden–Julian oscillation (MJO) in the NCAR CCM3 using a modified Zhang–McFarlane convection parameterization scheme. It is shown that, with the modified scheme, the intraseasonal (20–80 day) variability in precipitation, zonal wind, and outgoing longwave radiation (OLR) is enhanced substantially compared to the standard CCM3 simulation. Using a composite technique based on the empirical orthogonal function (EOF) analysis, the paper demonstrates that the simulated MJOs are in better agreement with the observations than the standard model in many important aspects. The amplitudes of the MJOs in 850-mb zonal wind, precipitation, and OLR are comparable to those of the observations, and the MJOs show clearly eastward propagation from the Indian Ocean to the Pacific. In contrast, the simulated MJOs in the standard CCM3 simulation are weak and have a tendency to propagate westward in the Indian Ocean. Nevertheless, there remain several deficiencies that are yet to be addressed. The time period of the MJOs is shorter, about 30 days, compared to the observed time period of 40 days. The spatial scale of the precipitation signal is smaller than observed. Examination of convective heating from both deep and shallow convection and its relationship with moisture anomalies indicates that near the mature phase of the MJO, regions of shallow convection developing ahead of the deep convection coincide with regions of positive moisture anomalies in the lower troposphere. This is consistent with the recent observations and theoretical development that shallow convection helps to precondition the atmosphere for MJO by moistening the lower troposphere. Sensitivity tests are performed on the individual changes in the modified convection scheme. They show that both change of closure and use of a relative humidity threshold for the convection trigger play important roles in improving the MJO simulation. Use of the new closure leads to the eastward propagation of the MJO and increases the intensity of the MJO signal in the wind field, while imposing a relative humidity threshold enhances the MJO variability in precipitation.

Description: A detailed view of Southern Hemisphere storm tracks is obtained based on the application of filtered variance and modern feature-tracking techniques to a wide range of 45-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) data. It has been checked that the conclusions drawn in this study are valid even if data from only the satellite era are used. The emphasis of the paper is on the winter season, but results for the four seasons are also discussed. Both upper- and lower-tropospheric fields are used. The tracking analysis focuses on systems that last longer than 2 days and are mobile (move more than 1000 km). Many of the results support previous ideas about the storm tracks, but some new insights are also obtained. In the summer there is a rather circular, strong, deep high-latitude storm track. In winter the high-latitude storm track is more asymmetric with a spiral from the Atlantic and Indian Oceans in toward Antarctica and a subtropical jet–related lower-latitude storm track over the Pacific, again tending to spiral poleward. At all times of the year, maximum storm activity in the higher-latitude storm track is in the Atlantic and Indian Ocean regions. In the winter upper troposphere, the relative importance of, and interplay between, the subtropical and subpolar storm tracks is discussed. The genesis, lysis, and growth rate of lower-tropospheric winter cyclones together lead to a vivid picture of their behavior that is summarized as a set of overlapping plates, each composed of cyclone life cycles. Systems in each plate appear to feed the genesis in the next plate through downstream development in the upper-troposphere spiral storm track. In the lee of the Andes in South America, there is cyclogenesis associated with the subtropical jet and also, poleward of this, cyclogenesis largely associated with system decay on the upslope and regeneration on the downslope. The genesis and lysis of cyclones and anticyclones have a definite spatial relationship with each other and with the Andes. At 500 hPa, their relative longitudinal positions are consistent with vortex-stretching ideas for simple flow over a large-scale mountain. Cyclonic systems near Antarctica have generally spiraled in from lower latitudes. However, cyclogenesis associated with mobile cyclones occurs around the Antarctic coast with an interesting genesis maximum over the sea ice near 150°E. The South Pacific storm track emerges clearly from the tracking as a coherent deep feature spiraling from Australia to southern South America. A feature of the summer season is the genesis of eastward-moving cyclonic systems near the tropic of Capricorn off Brazil, in the central Pacific and, to a lesser extent, off Madagascar, followed by movement along the southwest flanks of the subtropical anticyclones and contribution to the “convergence zone” cloud bands seen in these regions.

Description: The response of the climate system to natural, external forcing during the Maunder Minimum (ca. a.d. 1645–1715) is investigated using a comprehensive climate model. An ensemble of six transient simulations is produced in order to examine the relative importance of externally forced and internally generated variability. The simulated annual Northern Hemisphere and zonal-mean near-surface air temperature agree well with proxy-based reconstructions on decadal time scales. A mean cooling signal during the Maunder Minimum is masked by the internal unforced variability in some regions such as Alaska, Greenland, and northern Europe. In general, temperature exhibits a better signal-to-noise ratio than precipitation. Mean salinity changes are found in basin averages. The model also shows clear response patterns to volcanic eruptions. In particular, volcanic forcing is projected onto the winter North Atlantic Oscillation index following the eruptions. It is demonstrated that the significant spread of ensemble members is possible even on multidecadal time scales, which has an important implication in coordinating comparisons between model simulations and regional reconstructions.

Description: The authors demonstrate through atmospheric general circulation model (the Community Climate Model version 3.10) simulations of the 1997/98 El Niño that the observed “remote” (i.e., outside the Pacific) tropical land and ocean surface warming appearing a few months after the peak of the El Niño event is causally linked to the Tropics-wide warming of the troposphere resulting from increased atmospheric heating in the Pacific, with the latter acting as a conduit for the former. Unlike surface temperature, the surface flux behavior in the remote Tropics in response to El Niño is complex, with sizable spatial variation and compensation between individual flux components; this complexity suggests a more fundamental control (i.e., tropospheric temperature) for the remote tropical surface warming. Over the remote oceans, latent heat flux acting through boundary layer humidity variations is the important regulator linking the surface warming in the model simulations to the tropospheric warming over the remote tropical oceans. Idealized 1997/98 El Niño simulations using an intermediate tropical circulation model (the Quasi-Equilibrium Tropical Circulation Model) in which individual surface fluxes are directly manipulated confirms this result. The findings over the remote ocean are consistent with the “tropospheric temperature mechanism” previously proposed for the tropical ENSO teleconnection, with equatorial planetary waves propagating tropospheric temperature anomalies from the eastern Pacific to the remote Tropics and moist convective processes mediating the troposphere-to-remote-surface connection. The latter effectively requires the boundary layer moist static energy to vary in concert with the free tropospheric moist static energy. Over the remote land regions, idealized model simulations suggest that sensible heat flux regulates the warming response to El Niño, though the underlying mechanism has not yet been fully determined.

Description: To investigate relationships between large-scale circulation and regional-scale temperatures during the last (Eemian) interglacial, a simulation with a general circulation model (GCM) under orbital forcing conditions of 125 kyr BP is compared with a simulation forced with the Late Holocene preindustrial conditions. Consistent with previous GCM simulations for the Eemian, higher northern summer 2-m temperatures are found, which are directly related to the different insolation. Differences in the mean circulation are evident such as, for instance, stronger northern winter westerlies toward Europe, which are associated with warmer temperatures in central and northeastern Europe in the Eemian simulation, while the circulation variability, analyzed by means of a principal component analysis of the sea level pressure (SLP) field, is very similar in both periods. As a consequence of the differences in the mean circulation the simulated Arctic Oscillation (AO) temperature signal in the northern winter, on interannual-to-multidecadal time scales, is weaker during the Eemian than today over large parts of the Northern Hemisphere. Correlations between the AO index and the central European temperature (CET) decrease by about 0.2. The winter and spring SLP anomalies over the North Atlantic/European domain that are most strongly linearly linked to the CET cover a smaller area and are shifted westward over the North Atlantic during the Eemian. However, the strength of the connection between CET and these SLP anomalies is similar in both simulations. The simulated differences in the AO temperature signal and in the SLP anomaly, which is linearly linked to the CET, suggest that during the Eemian the link between the large-scale circulation and temperature-sensitive proxy data from Europe may differ from present-day conditions and that this difference should be taken into account when inferring large-scale climate from temperature-sensitive proxy data.

Description: An atmospheric general circulation model (AGCM) is used to examine the role of Indian Ocean sea surface temperature (SST) anomalies in regional climate variability. In particular, the authors focus on the effect of the basinwide warming that occurs during December through May after the mature phase of El Niño. To elucidate the relative importance of local and remote forcing, model solutions were sought for experiments where SST anomalies are inserted in the (i) tropical Indo-Pacific Oceans, (ii) tropical Pacific Ocean, and (iii) tropical Indian Ocean. A 10-member ensemble simulation is carried out for each of the three forcing scenarios. The model solutions demonstrate that precipitation variations over the southwest Indian Ocean are tied to local SST anomalies and are highly reproducible. Changes in the Indian Ocean–Walker circulation suppress precipitation over the tropical west Pacific–Maritime Continent, contributing to the development of a low-level anticyclone over the Philippine and South China Seas. Our model results indicate that more than 50% of the total precipitation anomalies over the tropical west Pacific–Maritime Continent is forced by remote Indian Ocean SST anomalies, offering an additional mechanism for the Philippine Sea anticyclone apart from Pacific SST. This anticyclone increases precipitation along the East Asian winter monsoon front from December to May. The anomalous subsidence over the Maritime Continent in conjunction with persistent anomalies of SST and precipitation over the Indian Ocean in spring prevent the northwestward migration of the ITCZ and the associated deep moist layer, causing a significant delay in the Indian summer monsoon onset in June by 6–7 days. At time scales of 5 days, however, the reproducibility of the northward progression of the ITCZ during the onset is low. Results indicate that Indian Ocean SST anomalies during December through May that develop in response to both atmospheric and oceanic processes to El Niño need to be considered for a complete understanding of regional climate variability, particularly around the Indian Ocean rim.

Description: Analyses of a 500-yr control integration with the non-flux-adjusted coupled atmosphere–sea ice–ocean model ECHAM5/Max-Planck-Institute Ocean Model (MPI-OM) show pronounced multidecadal fluctuations of the Atlantic overturning circulation and the associated meridional heat transport. The period of the oscillations is about 70–80 yr. The low-frequency variability of the meridional overturning circulation (MOC) contributes substantially to sea surface temperature and sea ice fluctuations in the North Atlantic. The strength of the overturning circulation is related to the convective activity in the deep-water formation regions, most notably the Labrador Sea, and the time-varying control on the freshwater export from the Arctic to the convection sites modulates the overturning circulation. The variability is sustained by an interplay between the storage and release of freshwater from the central Arctic and circulation changes in the Nordic Seas that are caused by variations in the Atlantic heat and salt transport. The relatively high resolution in the deep-water formation region and the Arctic Ocean suggests that a better representation of convective and frontal processes not only leads to an improvement in the mean state but also introduces new mechanisms determining multidecadal variability in large-scale ocean circulation.

Description: The causes of persistent droughts and wet periods, or pluvials, over western North America are examined in model simulations of the period from 1856 to 2000. The simulations used either (i) global sea surface temperature data as a lower boundary condition or (ii) observed data in just the tropical Pacific and computed the surface ocean temperature elsewhere with a simple ocean model. With both arrangements, the model was able to simulate many aspects of the low-frequency (periods greater than 6 yr) variations of precipitation over the Great Plains and in the American Southwest including much of the nineteenth-century variability, the droughts of the 1930s (the “Dust Bowl”) and 1950s, and the very wet period in the 1990s. Results indicate that the persistent droughts and pluvials were ultimately forced by persistent variations of tropical Pacific surface ocean temperatures. It is argued that ocean temperature variations outside of the tropical Pacific, but forced from the tropical Pacific, act to strengthen the droughts and pluvials. The persistent precipitation variations are part of a pattern of global variations that have a strong hemispherically and zonally symmetric component, which is akin to interannual variability, and that can be explained in terms of interactions between tropical ocean temperature variations, the subtropical jets, transient eddies, and the eddy-driven mean meridional circulation. Rossby wave propagation poleward and eastward from the tropical Pacific heating anomalies disrupts the zonal symmetry, intensifying droughts and pluvials over North America. Both mechanisms of tropical driving of extratropical precipitation variations work in summer as well as winter and can explain the year-round nature of the precipitation variations. In addition, land–atmosphere interactions over North America appear important by (i) translating winter precipitation variations into summer evaporation and, hence, precipitation anomalies and (ii) shifting the northward flow of moisture around the North Atlantic subtropical anticyclone eastward from the Plains and Southwest to the eastern seaboard and western Atlantic Ocean.

Description: The Pan-Atlantic sea surface temperature (SST) anomaly pattern that was found in a previous study to have a significant impact on the North Atlantic Oscillation (NAO) in early winter seemed to reflect the nearly uncorrelated influence of a horseshoe SST anomaly in the North Atlantic and an SST anomaly in the eastern equatorial Atlantic. A lagged rotated maximum covariance analysis of a slightly longer dataset shows that the horseshoe SST anomaly influence is robust, but it deemphasizes the center of action southeast of Newfoundland, Canada. On the other hand, it suggests that the link between equatorial SST and the NAO was artificial and due both to ENSO teleconnections and the orthogonality constraint in the maximum covariance analysis.

Description: The response of tropical Pacific SST to increased atmospheric CO2 concentration is reexamined with a new focus on the latitudinal SST gradient. Available evidence, mainly from climate models, suggests that an important tropical SST fingerprint to global warming is an enhanced equatorial warming relative to the subtropics. This enhanced equatorial warming provides a fingerprint of SST response more robust than the traditionally studied El Niño–like response, which is characterized by the zonal SST gradient. Most importantly, the mechanism of the enhanced equatorial warming differs fundamentally from the El Niño–like response; the former is associated with surface latent heat flux, shortwave cloud forcing, and surface ocean mixing, while the latter is associated with equatorial ocean upwelling and wind-upwelling dynamic ocean–atmosphere feedback.

Description: Downward solar irradiance at the sea surface, measured on several voyages of an icebreaker in the Southern Ocean, is used to infer transmittance of solar radiation by clouds. Together with surface albedo estimated from coincident hourly sea ice reports, instantaneous cloud radiative forcing and effective cloud optical depth are obtained. Values of “raw cloud transmittance” (trc), the ratio of downward irradiance under cloud to downward irradiance measured under clear sky, vary from 0.1 to 1.0. Over sea ice, few values of trc were observed between 0.8 and 1.0, possibly due to the threshold nature of the aerosol-to-cloud-droplet transition. This sparsely populated region of transmittances is referred to as the Köhler gap. The instantaneous downward shortwave cloud radiative forcing is computed, as well as the time-averaged net forcing. The net forcing at a solar zenith angle of 60° is typically −250 W m−2 over open ocean, but only half this value over sea ice because of the higher surface albedo and less frequent occurrence of clouds. “Effective” optical depths τ (for a radiatively equivalent horizontally homogeneous cloud) are classified by season and surface type. The frequency distributions of τ are well fitted by decaying exponentials, giving a characteristic optical depth of 15 at 47°S, increasing to 24 in the region of maximum cloud cover at 58°S, and decreasing to 11 at 67°S near the coast of Antarctica.

Description: A new northeast Asian summer monsoon index is introduced to investigate the characteristics of the northeast Asian summer rainfall variation, including Korea, Japan, and northeast China, and its possible connection to the tropical and midlatitude circulations. The summer precipitation over northeast Asia is separated into two components associated with tropical forcing and midlatitude dynamics using this monsoon index. The connection between the northeast Asian summer rainfall and ENSO is clearly identified by separating the Tropics-related component from the northeast Asian summer rainfall. That is, the Tropics-related precipitation over northeast Asia tends to be enhanced after the mature phase of El Niño. On the other hand, it is revealed that the extratropics-related component of summer precipitation is connected to the Eurasian wave pattern with no significant lag correlation. The intensity of the western North Pacific anticyclone modulated by ENSO is a key factor in the variation of the northeast Asian summer precipitation. It is found that the warm SST over the tropical eastern Pacific plays an important role in establishing the western North Pacific anticyclone during the preceding winter of strong northeast Asian summer monsoon years, whereas convective activities over the Bay of Bengal are contributed to the modulation of the anticyclonic circulation in the summer. The warming over the Indian Ocean in the summer of strong monsoon years induces the development of the anticyclone over the western North Pacific and the suppressed convection over the western Pacific tends to enhance the northeast Asian summer rainfall through the Pacific–Japan or East Asia–Pacific teleconnections.

Description: Trends in surface temperature over the last 100, 50, and 30 yr at individual grid boxes in a 5° latitude–longitude grid are compared with model estimates of the natural internal variability of these trends and with the model response to increasing greenhouse gases and sulfate aerosols. Three different climate models are used to provide estimates of the internal variability of trends, one of which appears to overestimate the observed variability of surface temperature at interannual and 5-yr time scales. Significant warming trends are found at a large fraction of the individual grid boxes over the globe, a much larger fraction than can be explained by internal climate variations. The observed warming trends over the last 50 and 30 yr are consistent with the modeled response to increasing greenhouse gases and sulfate aerosols in most of the models. However, in some regions, the observed century-scale trends are significantly larger than the modeled response to increasing greenhouse gases and sulfate aerosols in the atmosphere. Warming trends consistent with the response to anthropogenic forcing are detected at scales on the order of 500 km in many regions of the globe.

Description: Relationships between Gulf of California moisture surges and tropical cyclones (TCs) in the eastern Pacific basin are examined. Standard surface observations are used to identify gulf surge events at Yuma, Arizona, for a multiyear (July–August 1979–2001) period. The surges are related to TCs using National Hurricane Center 6-hourly track data for the eastern Pacific basin. Climate Prediction Center (CPC)- observed daily precipitation analyses and the NCEP Regional Reanalysis are used to examine the relative differences in the precipitation, atmospheric circulation, and moisture fields for several categories of surge events, including those that are directly related to TCs, indirectly related to TCs, and not related to TCs. It is shown that the response to the surge in the southwestern United States and northwestern Mexico is strongly discriminated by the presence or absence of TCs. Surges related to TCs tend to be associated with much stronger and deeper low-level southerly flow, deeper plumes of tropical moisture, and wetter conditions over the core monsoon region than surges that are unrelated to TCs. The response to the surge is also strongly influenced by the proximity of the TC to the Gulf of California (GOC) region. Tropical cyclones that track toward the GOC region exert a stronger, more direct influence on Yuma surges than those that track away from the GOC.

Description: The economic value of seasonal climate forecasting is assessed using a whole-of-chain analysis. The entire system, from sea surface temperature (SST) through pasture growth and animal production to economic and resource outcomes, is examined. A novel statistical forecast method is developed using the partial least squares spatial correlation technique with near-global SST. This method permits forecasts to be tailored for particular regions and industries. The method is used to forecast plant growth days rather than rainfall. Forecast skill is measured by performing a series of retrospective forecasts (hindcasts) over the previous century. The hindcasts are cross-validated to guard against the possibility of artificial skill, so there is no skill at predicting random time series. The hindcast skill is shown to be a good estimator of the true forecast skill obtained when only data from previous years are used in developing the forecast. Forecasts of plant growth, reduced to three categories, are used in several agricultural examples in Australia. For the northeast Queensland grazing industry, the economic value of this forecast is shown to be greater than that of a Southern Oscillation index (SOI) based forecast and to match or exceed the value of a “perfect” category rainfall forecast. Reasons for the latter surprising result are given. Resource degradation, in this case measured by soil loss, is shown to remain insignificant despite increasing production from the land. Two further examples in Queensland, one for the cotton industry and one for wheat, are illustrated in less depth. The value of a forecast is again shown to match or exceed that obtained using the SOI, although further investigation of the decision-making responses to forecasts is needed to extract the maximum benefit for these industries.

Description: From ensembles of 80 AGCM simulations for every December–January–February (DJF) seasonal mean in the 1980–2000 period, interannual variability in atmospheric response to interannual variations in observed sea surface temperature (SST) is analyzed. A unique facet of this study is the use of large ensemble size that allows identification of the atmospheric response to SSTs for each DJF in the analysis period. The motivation of this study was to explore what atmospheric response patterns beyond the canonical response to El Niño–Southern Oscillation (ENSO) SST anomalies exist, and to which SST forcing such patterns may be related. A practical motivation for this study was to seek sources of atmospheric predictability that may lead to improvements in seasonal predictability efforts. This analysis was based on the EOF technique applied to the ensemble mean 200-mb height response. The dominant mode of the atmospheric response was indeed the canonical atmospheric response to ENSO; however, this mode only explained 53% of interannual variability of the ensemble means (often referred to as the external variability). The second mode, explaining 19% of external variability, was related to a general increase (decrease) in the 200-mb heights related to a Tropicwide warming (cooling) in SSTs. The third dominant mode, explaining 12% of external variability, was similar to the mode identified as the “nonlinear” response to ENSO in earlier studies. The realism of different atmospheric response patterns was also assessed from a comparison of anomaly correlations computed between different renditions of AGCM-simulated atmospheric responses and the observed 200-mb height anomalies. For example, the anomaly correlation between the atmospheric response reconstructed from the first mode alone and the observations was compared with the anomaly correlation when the atmospheric response was reconstructed including modes 2 and 3. If the higher-order atmospheric response patterns obtained from the AGCM simulations had observational counterparts, their inclusion in the reconstructed atmospheric response should lead to higher anomaly correlations. Indeed, at some geographical regions, an increase in anomaly correlation with the inclusion of higher modes was found, and it is concluded that the higher-order atmospheric response patterns found in this study may be realistic and may represent additional sources of atmospheric seasonal predictability.

Description: Cloud–climate feedbacks between precipitation, radiation, circulation strength, atmospheric temperature and moisture, and ocean temperature are studied with an idealized model of the Walker circulation in a nonrotating atmosphere coupled to an ocean mixed layer. This study has two main purposes: 1) to formulate a conceptual framework that includes the dominant feedbacks between clouds and a large-scale divergent circulation; and 2) to use this framework to investigate the sensitivity of the climate system to these interactions. Two cloud types—high, convective anvils and low, nonprecipitating stratus—are included and coupled to the large-scale dynamics. The atmosphere is coupled to an ocean mixed layer via a consistent surface energy budget. Analytic approximations with a simplified radiation scheme are derived and used to explain numerical results with a more realistic radiation scheme. The model simplicity allows interactions between different parts of the ocean–atmosphere system to be cleanly elucidated, yet also allows the areal extent of deep convection and the horizontal structure of the Walker circulation to be internally determined by the model. Because of their strong top-of-atmosphere radiative cancellation, high clouds are found to have little overall effect on the circulation strength and convective area fraction. Instead, to leading order, these are set by the horizontally varying ocean heat transport and clear-sky radiative fluxes. Low clouds are found to cool both the ocean and atmosphere, to slightly increase the circulation strength, and to shrink the convective area significantly. The climate is found to be less sensitive to doubled greenhouse gas experiments with low clouds than without.

Description: Projections of future global sea level depend on reliable estimates of changes in the size of polar ice sheets. Calculating this directly from global general circulation models (GCMs) is unreliable because the coarse resolution of 100 km or more is unable to capture narrow ablation zones, and ice dynamics is not usually taken into account in GCMs. To overcome these problems a high-resolution (20 km) dynamic ice sheet model has been coupled to the third Hadley Centre Coupled Ocean–Atmosphere GCM (HadCM3). A novel feature is the use of two-way coupling, so that climate changes in the GCM drive ice mass changes in the ice sheet model that, in turn, can alter the future climate through changes in orography, surface albedo, and freshwater input to the model ocean. At the start of the main experiment the atmospheric carbon dioxide concentration was increased to 4 times the preindustrial level and held constant for 3000 yr. By the end of this period the Greenland ice sheet is almost completely ablated and has made a direct contribution of approximately 7 m to global average sea level, causing a peak rate of sea level rise of 5 mm yr−1 early in the simulation. The effect of ice sheet depletion on global and regional climate has been examined and it was found that apart from the sea level rise, the long-term effect on global climate is small. However, there are some significant regional climate changes that appear to have reduced the rate at which the ice sheet ablates.

Description: The output from an ocean general circulation model driven by observed surface forcing (1958–97) is used to examine the evolution and relative timing of the different branches of the Pacific Subtropical–Tropical Cells (STCs) at both interannual and decadal time scales, with emphasis on the 1976–77 climate shift. The STCs consist of equatorward pycnocline transports in the ocean interior and in the western boundary current, equatorial upwelling, and poleward flow in the surface Ekman layer. The interior pycnocline transports exhibit a decreasing trend after the mid-1970s, in agreement with observational transport estimates, and are largely anticorrelated with both the Ekman transports and the boundary current transports at the same latitudes. The boundary current changes tend to compensate for the interior changes at both interannual and decadal time scales. The meridional transport convergence across 9°S and 9°N as well as the equatorial upwelling are strongly correlated with the changes in sea surface temperature (SST) in the central and eastern equatorial Pacific. However, meridional transport variations do not occur simultaneously at each longitude, so that to understand the phase relationship between transport and SST variations it is important to consider the baroclinic ocean adjustment through westward-propagating Rossby waves. The anticorrelation between boundary current changes and interior transport changes can also be understood in terms of the baroclinic adjustment process. In this simulation, the pycnocline transport variations appear to be primarily confined within the Tropics, with maxima around 10°S and 13°N, and related to the local wind forcing; a somewhat different perspective from previous studies that have emphasized the role of wind variations in the subtropics.

Description: The sea ice motion, area export, and deformation of the Ross Sea ice cover are examined with satellite passive microwave and RADARSAT observations. The record of high-resolution synthetic aperture radar (SAR) data, from 1998 and 2000, allows the estimation of the variability of ice deformation at the small scale (∼10 km) and to assess the quality of the longer record of passive microwave ice motion. Daily and subdaily deformation fields and RADARSAT imagery highlight the variability of motion and deformation in the Ross Sea. With the passive microwave ice motion, the area export at a flux gate positioned between Cape Adare and Land Bay is estimated. Between 1992 and 2003, a positive trend can be seen in the winter (March–November) ice area flux that has a mean of 990 × 103 km2 and ranges from a low of 600 × 103 km2 in 1992 to a peak of 1600 × 103 km2 in 2001. In the mean, the southern Ross Sea produces almost twice its own area of sea ice during the winter. Cross-gate sea level pressure (SLP) gradients explain ∼60% of the variance in the ice area flux. A positive trend in this gradient, from reanalysis products, suggests a “spinup” of the Ross Sea Gyre over the past 12 yr. In both the NCEP–NCAR and ERA-40 surface pressure fields, longer-term trends in this gradient and mean SLP between 1979 and 2002 are explored along with positive anomalies in the monthly cross-gate SLP gradient associated with the positive phase of the Southern Hemisphere annular mode and the extrapolar Southern Oscillation.

Description: The characteristic features of the boreal summer intraseasonal oscillation (BSISO) during its reinitiation period are studied using NCEP–NCAR reanalysis. Based on these observations and with the aid of an anomalous atmospheric general circulation model (AGCM), a possible mechanism responsible for the BSISO reinitiation is elucidated. The western equatorial Indian Ocean along the eastern African coast tends to be a key region for the phase transition of the BSISO from an enhanced to suppressed convective phase, or vise versa. The major precursory feature associated with reinitiation of suppressed convection is found in the divergence and reduced specific humidity in the boundary layer. Numerical experiments indicate that the low-level divergence is caused by the cold horizontal temperature advection and associated adiabatic warming (descending motion) in situ. The summer mean state is found to be important for the cold horizontal temperature advection through the modulation of a Gill-type response to an intraseasonal oscillation (ISO) heating in the eastern equatorial Indian Ocean. The results in this study suggest a self-sustained paradigm in the Indian Ocean for the BSISO; that is, the BSISO could be a basinwide phenomenon instead of a global circumstance system as hypothesized for the boreal winter ISO (i.e., the Madden–Julian oscillation).

Description: A dynamic–thermodynamic sea ice model with the ocean mixed layer forced by atmospheric data is used to investigate spatial and long-term variability of the sea ice cover in the Arctic basin. The model satisfactorily reproduces the averaged main characteristics of the sea ice and its extent in the Arctic Basin, as well as its decrease in the early 1990s. Employment of the average ridge shape for describing the ridging allows the authors to suggest that it occurs in winter and varies from year to year by a factor of 2, depending on an atmospheric circulation pattern. Production and horizontal movement of ridges are the focus in this paper, as they show the importance of interannual variability of the Arctic ice cover. The observed thinning in the 1990s is a result of reduction in ridge formation on the Pacific side during the cyclonic phase of the Arctic Oscillation. The model yields a partial recovery of sea ice cover in the last few years of the twentieth century. In addition to the sea ice cover and average thickness compared with satellite data, the ridge amount is verified with observations taken in the vicinity of the Russian coast. The model results are useful to estimate long-term variability of the probability of ridge-free navigation in different parts of the Arctic Ocean, including the Northern Sea Route area.

Description: The causes for the observed occurrence of anomalous zonally symmetric upper-level pressure ridges in the midlatitude belts of both hemispheres during the year after warm El Niño–Southern Oscillation (ENSO) events have been investigated. Sea surface temperature (SST) anomalies in the Indo–western Pacific (IWP) sector were simulated by allowing an oceanic mixed layer model for that region to interact with local atmospheric changes forced remotely by observed ENSO episodes in the eastern/central tropical Pacific. The spatiotemporal evolution of these SST conditions through a composite ENSO cycle was then inserted as lower boundary conditions within the IWP domain in an ensemble of atmospheric general circulation model (GCM) integrations. This experimental setup is seen to reproduce zonally symmetric geopotential height anomalies with maximum amplitudes being attained over the extratropics in the boreal summer after the peak phase of ENSO. The model evidence hence supports the notion that these global-scale atmospheric changes are primarily responses to SST perturbations in IWP, which are in turn linked to ENSO variability in the equatorial Pacific by the “atmospheric bridge” mechanism. Experimentation with a stationary wave model indicates that the Eastern Hemisphere portion of the aforementioned atmospheric signals are attributable to forcing by tropical heat sources and sinks associated with precipitation anomalies in the IWP region, which are closely related to the underlying SST changes. Diagnosis of the output from the GCM integrations reveals that these circulation changes due to diabatic heating are accompanied by alterations of the propagation path and intensity of the high-frequency eddies at locations farther downstream. The geopotential tendencies associated with the latter disturbances bear some resemblance to the anomalous height pattern in the Western Hemisphere. Such local eddy–mean flow feedbacks hence contribute to the zonal symmetry of the atmospheric response pattern to forcing in the IWP region. Analysis of zonally averaged circulation statistics indicates that the mean meridional circulation induced by divergence of anomalous transient eddy momentum fluxes in ENSO events could also generate zonally symmetric perturbations in midlatitudes. The model-simulated precipitation and surface temperature anomalies in the North American sector in response to SST changes in IWP suggest an increased frequency of droughts and heat waves in that region during the summer season after warm ENSO events.

Description: A coupled global atmosphere–ocean model of intermediate complexity is used to study the influence of glacial boundary conditions on the atmospheric circulation during the Last Glacial Maximum in a systematical manner. A web of atmospheric interactions is disentangled, which involves changes in the meridional temperature gradient and an associated modulation of the atmospheric baroclinicity. This in turn drives anomalous transient eddy momentum fluxes that feed back onto the zonal mean circulation. Moreover, the modified transient activity (weakened in the North Pacific and strengthened in the North Atlantic) leads to a meridional reorganization of the atmospheric heat transport, thereby feeding back onto the meridional temperature structure. Furthermore, positive barotropic conversion and baroclinic production rates over the Laurentide ice sheets and the far eastern North Pacific have the tendency to decelerate the westerlies, thereby feeding back to the stationary wave changes triggered by orographic forcing.

Description: A new physiology-based model of canopy stomatal conductance and photosynthesis is described and included in the latest version of the Goddard Institute for Space Studies (GISS) GCM, ModelE1. The submodel includes responses to atmospheric humidity and CO2 concentration, responses missing from previous GISS GCM land surface schemes. Measurements of moisture, energy, and CO2 fluxes over four vegetation types are used to test and calibrate the submodel. Photosynthetic leaf N is calibrated for each vegetation type from the flux measurements. The new submodel results in surface cooling over many regions previously too warm. Some warm biases of over 2°C are cooled by more than 0.5°C, including over central Eurasia, South America, the western United States, and Australia. In addition, some regions that were previously too cool are warmed, such as northern Eurasia and the Tibetan Plateau. A number of precipitation biases are also reduced, particularly over South America (by up to 1 mm day−1) and the oceanic ITCZs (by over ±1 mm day−1); coastal west Africa becomes significantly wetter. Cloud cover increases over many land areas previously too clear. Higher absolute canopy conductances, and positive feedbacks with atmospheric humidity, are largely responsible for the simulated vegetation influence on the atmosphere. High-latitude climate changes through remote effects of increased tropical latent heating, resulting directly from improved characterization of tropical forest canopy conductance. Realistic representation of the stomatal control on land evaporation is critical for accurate simulation of atmospheric dynamics in the GISS GCM.

Description: The influence of the El Niño–Southern Oscillation (ENSO) on tropical cyclone intensity in the western North Pacific basin is examined. Accumulated cyclone energy (ACE), constructed from the best-track dataset for the region for the period 1950–2002, and other related variables are analyzed. ACE is positively correlated with ENSO indices. This and other statistics of the interannually varying tropical cyclone distribution are used to show that there is a tendency in El Niño years toward tropical cyclones that are both more intense and longer-lived than in La Niña years. ACE leads ENSO indices: during the peak season (northern summer and fall), ACE is correlated approximately as strongly with ENSO indices up to six months later (northern winter), as well as simultaneously. It appears that not all of this lead–lag relationship is easily explained by the autocorrelation of the ENSO indices, though much of it is. Interannual variations in the annual mean lifetime, intensity, and number of tropical cyclones all contribute to the ENSO signal in ACE, though the lifetime effect appears to be the most important of the three.

Description: Observations show the asymmetric nature of El Niño and La Niña sea surface temperature (SST) anomalies. Warm events are often stronger than cold events. This asymmetric behavior is an important feature that can be used to validate coupled models to test their ability to represent the climate system. The asymmetry of El Niño and La Niña SST anomalies has been investigated in a simulation of the Hadley Centre eddy-permitting coupled general circulation model. It is found that the asymmetric behavior is captured by the model with SST anomalies associated with strong El Niño events being greater than those associated with strong La Niña events. The pattern of the SST asymmetry also bears some similar characteristics to those based on observations despite the deficiency that SST anomalies associated with both El Niño and La Niña extend too far westward in the model. Through a heat budget analysis of the ocean mixed layer, it is shown that nonlinear dynamic heating (NDH) is important in generating intense El Niño and the SST asymmetry between El Niño and La Niña events, especially in the eastern tropical Pacific Ocean. This nonlinear dynamic heating enhances the amplitude of El Niño and reduces the amplitude of La Niña, and therefore leads to the asymmetry between El Niño and La Niña events, with El Niño being stronger. However, the skewness and asymmetry in the model are relatively weak, being consistent with a relatively weak nonlinear dynamical heating. It is also shown that the eastward-propagating feature of subsurface anomalies provides a favorable phase relationship between temperature and current anomalies that results in strong nonlinear dynamical heating that tends to produce stronger El Niños. In addition, in the model simulation, the nonlinear nature of zonal wind stress anomalies between El Niño and La Niña events also plays an important role in the central tropical Pacific. These different mechanisms work constructively to determine the asymmetry between El Niño and La Niña events in the model, and they are similar to those proposed in recent studies based on observations. The ability of the model to simulate this asymmetric feature is encouraging and offers hope to the challenge of predicting the amplitude of strong El Niño events.

Description: Diagnostics combining atmospheric reanalysis and station-based temperature data for 1950–2003 indicate that European heat waves can be associated with the occurrence of two specific summertime atmospheric circulation regimes. Evidence is presented that during the record warm summer of 2003, the excitation of these two regimes was significantly favored by the anomalous tropical Atlantic heating related to wetter-than-average conditions in both the Caribbean basin and the Sahel. Given the persistence of tropical Atlantic climate anomalies, their seasonality, and their associated predictability, the suggested tropical–extratropical Atlantic connection is encouraging for the prospects of long-range forecasting of extreme weather in Europe.

Description: The boreal winter Madden–Julian oscillation (MJO) remains very weak and irregular in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model version 2 (CAM2) as in its direct predecessor, the Community Climate Model version 3 (CCM3). The standard version of CAM2 uses the deep convective scheme of Zhang and McFarlane, as in CCM3, with the closure dependent on convective available potential energy (CAPE). Here, sensitivity tests using several versions of the Tiedtke convective scheme are conducted. Typically, the Tiedtke convection scheme gives an improved mean state, intraseasonal variability, space–time power spectra, and eastward propagation compared to the standard version of the model. Coherent eastward propagation of MJO-related precipitation is also much improved, particularly over the Indian–western Pacific Oceans. A composite life cycle of the model MJO indicates that over the Indian Ocean wind-induced surface heat exchange (WISHE) functions, while over the western/central Pacific Ocean aspects of frictional moisture convergence are evident in the maintenance and eastward propagation of the oscillation.

Description: The use of linear statistical methods in building climate prediction models is examined, particularly the use of anomalies. The author’s perspective is that the climate system is a nonlinear interacting system, so the impact of modeling using anomalies rather than observed data directly is considered. With reference to the Lorenz system and a simple model for regime dependence, it is shown that anomalies impair our ability to reconstruct nonlinear dynamics. Some alternative approaches in the literature that offer an attractive way forward are explored, focusing on Bayesian hierarchical methods to construct so-called physical–statistical models. The author’s view is that anomalies should be reserved in most cases as a tool for enhancing graphical representations of climate data. The exceptions are when the implicit assumptions underlying the use of anomalies are met or when an anomaly representation is physically motivated.

Description: A six-member multicoupled model ensemble is created by using six state-of-the-art deep atmospheric convective schemes. The six convective schemes are used inside a single model and make up the ensemble. This six-member ensemble is compared against a multianalysis ensemble, which is created by varying the initial start dates of the atmospheric component of the coupled model. Both ensembles were integrated for seven months (November–May) over a 12-yr period from 1987 to 1998. Examination of the sea surface temperature and precipitation show that while deterministic skill scores are slightly better for the multicoupled model ensemble the probabilistic skill scores favor the multimodel approach. Combining the two ensembles to create a larger ensemble size increases the probabilistic skill score compared to the multimodel. This altering physics approach to create a multimodel ensemble is seen as an easy way for small modeling centers to generate ensembles with better reliability than by only varying the initial conditions.

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

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

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

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

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

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

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

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

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

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

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

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

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