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