ALBERT

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

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

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

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

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

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

Description: While it is obvious that the mean diabatic forcing of the atmosphere is crucial for maintaining the mean climate, the importance of diabatic forcing fluctuations is less evident in this regard. Such fluctuations do not appear directly in the equations of the mean climate but affect the mean indirectly through their effects on the time-mean transient-eddy fluxes of heat, momentum, and moisture. How large are these effects? What are the effects of tropical phenomena associated with substantial heating variations such as ENSO and the MJO? To what extent do variations of the extratropical surface heat fluxes and precipitation affect the mean climate? What are the effects of the rapid “stochastic” components of the heating fluctuations? Most current climate models misrepresent ENSO and the MJO and ignore stochastic forcing; they therefore also misrepresent their mean effects. To what extent does this contribute to climate model biases and to projections of climate change? This paper provides an assessment of such impacts by comparing with observations a long simulation of the northern winter climate by a dry adiabatic general circulation model forced only with the observed time-mean diabatic forcing as a constant forcing. Remarkably, despite the total neglect of all forcing variations, the model reproduces most features of the observed circulation variability and the mean climate, with biases similar to those of some state-of-the-art general circulation models. In particular, the spatial structures of the circulation variability are remarkably well reproduced. Their amplitudes, however, are progressively underestimated from the synoptic to the subseasonal to interannual and longer time scales. This underestimation is attributed to the neglect of the variable forcing. The model also excites significant tropical variability from the extratropics on interannual scales, which is overwhelmed in reality by the response to tropical heating variability. It is argued that the results of this study suggest a role for the stochastic, and not only the coherent, components of transient diabatic forcing in the dynamics of climate variability and the mean climate.

Description: The interannual variation in cold-air outbreak activity over the Japan Sea is investigated using Japan Meteorological Agency buoy 21002 and Quick Scatterometer (QuikSCAT) wind data, Japan Oceanographic Data Center sea surface temperature (SST) data, NCEP–NCAR reanalysis surface wind and sea level pressure (SLP) data, and the winter Arctic Oscillation (AO) index of Thompson and Wallace. Cold-air outbreaks occur during the “winter” November–March period, and wind data for this season for the 19-winter period 1981–2000 were analyzed. Wavelet spectra averaged between 5- and 15-day periods were used to evaluate the intensity of cold-air outbreaks quantitatively. The winter mean wavelet spectra exhibited a clear interannual variation and a significant positive correlation with the AO index, indicating that intensive cold-air outbreaks frequently occur during relatively warm winters caused by a quasi-decadal AO. Based on the SST and SLP data, the low atmospheric surface pressure disturbances tend to develop over the warm East China Sea in warm winters in the positive AO phase. As these low SLP disturbances advance toward the northern Japan islands during the positive AO phase, they intensify more, leading to stronger cold-air outbreaks over the Japan Sea and increased sea surface cooling over the northern Japan Sea.

Description: The causes and global context of the North American drought between 1998 and 2004 are examined using atmospheric reanalyses and ensembles of atmosphere model simulations variously forced by global SSTs or tropical Pacific SSTs alone. The drought divides into two distinct time intervals. Between 1998 and 2002 it coincided with a persistent La Niña–like state in the tropical Pacific, a cool tropical troposphere, poleward-shifted jet streams, and, in the zonal mean, eddy-driven descent in midlatitudes. During the winters reduced precipitation over North America in the climate models was sustained by anomalous subsidence and reductions of moisture convergence by the stationary flow and transient eddies. During the summers reductions of evaporation and mean flow moisture convergence drove the precipitation reduction, while transient eddies acted diffusively to oppose this. During these years the North American drought fitted into a global pattern of circulation and hydroclimate anomalies with noticeable zonal and hemispheric symmetry. During the later period of the drought, from 2002 to 2004, weak El Niño conditions prevailed and, while the global climate adjusted accordingly, western North America remained, uniquely among midlatitude regions, in drought. The ensemble mean of the climate model simulations did not simulate the continuation of the drought in these years, suggesting that the termination of the drought was largely unpredictable in terms of global ocean conditions. The global context of the most recent, turn of the century, drought is compared to the five prior persistent North American droughts in the instrumental record from the mid-nineteenth century on. A classic La Niña pattern of ocean temperature in the Pacific is common to all. A cold Indian Ocean, also typical of La Niña, is common to all five prior droughts, but not the most recent one. Except in southern South America the global pattern of precipitation anomalies of the turn of the century drought is similar to that during the five prior droughts. These comparisons suggest that the earlier period of this most recent drought is the latest in a series of multiyear droughts forced by persistent changes in tropical Pacific Ocean temperatures. Warm tropical North Atlantic Ocean temperatures may play a secondary role.

Description: This study provides further evidence of the impacts of tropical Pacific interannual [El Niño–Southern Oscillation (ENSO)] and Northern Pacific decadal–interdecadal [North Pacific index (NPI)] variability on the Pacific–North American (PNA) sector. Both the tropospheric circulation and the North American temperature suggest an enhanced PNA-like climate response and impacts on North America when ENSO and NPI variability are out of phase. In association with this variability, large stationary wave activity fluxes appear in the mid- to high latitudes originating from the North Pacific and flowing downstream toward North America. Atmospheric heating anomalies associated with ENSO variability are confined to the Tropics, and generally have the same sign throughout the troposphere with maximum anomalies at 400 hPa. The heating anomalies that correspond to the NPI variability exhibit a center over the midlatitude North Pacific in which the heating changes sign with height, along with tropical anomalies of comparable magnitudes. Atmospheric heating anomalies of the same sign appear in both the tropical Pacific and the North Pacific with the out-of-phase combination of ENSO and NPI. Both sources of variability provide energy transports toward North America and tend to favor the occurrence of stationary wave anomalies.

Description: Using the coupled climate model of intermediate complexity, CLIMBER-3α, changes in the vertical thermal structure associated with a shutdown of the Atlantic meridional overturning circulation (AMOC) are investigated. When North Atlantic Deep Water formation is inhibited by anomalous freshwater forcing, intermediate depth ventilation can remain active and cool the subsurface water masses (i.e., the “cold case”). However, if intermediate ventilation is completely suppressed, relatively warm water coming from the south penetrates to a high northern latitude beneath the halocline and induces a strong vertical temperature inversion between the surface and intermediate depth (i.e., the “warm case”). Both types of temperature anomalies emerge within the first decade after the beginning of the freshwater perturbation. The sign of subsurface temperature anomaly has a strong implication for the recovery of the AMOC once the anomalous freshwater forcing is removed. While the AMOC recovers from the cold case on centennial time scales, the recovery is much more rapid (decadal time scales) when ventilation is completely suppressed and intermediate depths are anomalously warm. This is explained by a more rapid destabilization of the water column after cessation of the anomalous flux due to a strong vertical temperature inversion. A suite of sensitivity experiments with varying strength and duration of the freshwater perturbation and a larger value of background vertical diffusivity demonstrate robustness of the phenomenon. Implications of the simulated subsurface temperature response to the shutdown of the AMOC for future climate and abrupt climate changes of the past are discussed.

Description: Interannual extremes in New Zealand rainfall and their modulation by modes of Southern Hemisphere climate variability are examined in observations and a coupled climate model. North Island extreme dry (wet) years are characterized by locally increased (reduced) sea level pressure (SLP), cold (warm) sea surface temperature (SST) anomalies in the southern Tasman Sea and to the north of the island, and coinciding reduced (enhanced) evaporation upstream of the mean southwesterly airflow. During extreme dry (wet) years in South Island precipitation, an enhanced (reduced) meridional SLP gradient occurs, with circumpolar strengthened (weakened) subpolar westerlies and an easterly (westerly) anomaly in zonal wind in the subtropics. As a result, via Ekman transport, anomalously cold (warm) SST appears under the subpolar westerlies, while anomalies of the opposite sign occur farther north. The phase and magnitude of the resulting SST and evaporation anomalies cannot account for the rainfall extremes over the South Island, suggesting a purely atmospheric mode of variability as the driving factor, in this case the Southern Annular Mode (SAM). New Zealand rainfall variability is predominantly modulated by two Southern Hemisphere climate modes, namely, the El Niño–Southern Oscillation (ENSO) and the SAM, with a latitudinal gradation in influence of the respective phenomena, and a notable interaction with orographic features. While this heterogeneity is apparent both latitudinally and as a result of orographic effects, climate modes can force local rainfall anomalies with considerable variations across both islands. North Island precipitation is for the most part regulated by both local air–sea heat fluxes and circulation changes associated with the tropical ENSO mode. In contrast, for the South Island the influence of the large-scale general atmospheric circulation dominates, especially via the strength and position of the subpolar westerlies, which are modulated by the extratropical SAM.

Description: This work investigates the dynamic and thermal response of the winter stratosphere to the presence of a weak meridional surface temperature gradient. Previous work suggested that polar stratospheric clouds could have played a decisive role in maintaining high-latitude warmth, especially over continental interiors, during the polar nights of the late Paleocene and early Eocene epochs; both a chemical source of additional water vapor and a dynamical feedback between the surface climate and stratospheric temperatures have been proposed as mechanisms by which such clouds could form. A principal goal of this work is to investigate the latter problem using a general circulation model with stratospheric resolution that is forced with a very weak surface temperature gradient. It is found that temperatures in the lower stratosphere do not deviate significantly from the control run, which results from a robust flux of wave activity into the winter stratosphere. The strength of the stratosphere’s residual circulation increases slightly in the presence of the weak gradient, as wavenumber 3 begins to propagate to stratospheric altitudes. Changes in the zonal wind field that allow for the altered propagation are in balance with a weakened temperature gradient through the full depth of the troposphere. These simulations also suggest that the tropospheric thermal stratification could be maintained by moist convection at all latitudes in warm climate states with a weak temperature gradient.

Description: The midsummer drought (MSD) is a diminution in rainfall experienced during the middle of the rainy season in southern Mexico and Central America, as well as in the adjacent Caribbean, Gulf of Mexico, and eastern Pacific seas. The aim of this paper is to describe the regional characteristics of the MSD and to propose some possible forcing mechanisms. Satellite and in situ data are used to form a composite of the evolution of a typical MSD, which highlights its coincidence with a low-level anticyclone centered over the Gulf of Mexico and associated easterly flow across Central America. The diurnal cycle of precipitation over the region is reduced in amplitude during midsummer. The MSD is also coincident with heavy precipitation over the Sierra Madre Occidental (part of the North American monsoon). Reanalysis data are used to show that the divergence of the anomalous low-level flow during the MSD is the main factor governing the variations in precipitation. A linear baroclinic model is used to show that the seasonal progression of the Pacific intertropical convergence zone (ITCZ), which moves northward following warm sea surface temperature (SST) during the early summer, and of the Atlantic subtropical high, which moves westward, are the most important remote factors that contribute toward the low-level easterly flow and divergence during the MSD. The circulation associated with the MSD precipitation deficit helps to maintain the deficit by reinforcing the low-level anticyclonic flow over the Gulf of Mexico. Surface heating over land also plays a role: a large thermal low over the northern United States in early summer is accompanied by enhanced subsidence over the North Atlantic. This thermal low is seen to decrease considerably in midsummer, allowing the high pressure anomalies in the Atlantic and Pacific Oceans to extend into the Gulf of Mexico. These anomalies are maintained until late summer, when an increase in rainfall from the surge in Atlantic tropical depressions induces anomalous surface cyclonic flow with westerlies fluxing moisture from the Pacific ITCZ toward Central America.

Description: Interdecadal climate changes occurring during the latter half of August (LA) in central Japan are described, and the associated changes in rainfall, typhoon tracks, and circulation patterns over East Asia and the western North Pacific (WNP) regions are investigated. Since 1984, rates of sunshine and temperature have increased, while rainfall has decreased significantly during LA in central Japan. In contrast, rainfall over the Southwest Islands, northern part of Taiwan, and in the south of the middle and lower reaches of the Yangtze River in China has increased. Changes in the positions and tracks of typhoons are responsible for these changes. Prior to 1983, many typhoons approached the central-western part of Japan during LA, while after 1984, most typhoons were deflected away from Japan and moved northwestward to Taiwan and China. The North Pacific subtropical high during LA extended more westward after 1984, which affected the interdecadal changes in sunshine, temperature, rainfall, and typhoon tracks, not only in central Japan, but also over East Asia and the WNP regions.

Description: The authors identified an upper-level pressure anomaly pattern corresponding to the interannual variability of the Okhotsk high in midsummer (late July and early August) as a predominant anomaly pattern in the Northern Hemisphere, by using objectively analyzed data. According to the results of empirical orthogonal function (EOF) analyses and composite analyses, a positive pressure anomaly appeared near the tropopause over eastern Siberia in years with strong Okhotsk highs. Examination of the heat budget in the lower troposphere revealed that a negative surface temperature anomaly observed in northern Japan was brought by the advection of the climatological temperature gradient from the anomalous wind associated with the upper-level anticyclonic anomaly. It was also demonstrated that the anomaly field over Siberia does not accompany predominant vorticity forcing or Rossby wave propagation from the west with a specific phase. However, positive kinetic energy conversion from the climatological basic field to the anomaly field is estimated. The energy conversion contributes to maintaining the anomaly pattern. By the numerical experiments using a linear barotropic model, it is suggested that the upper-level anomaly pattern related to the anomalous Okhotsk high appears through the interaction with the climatological basic field, even though the external forcings are homogeneously distributed.

Description: Undocumented changepoints (inhomogeneities) are ubiquitous features of climatic time series. Level shifts in time series caused by changepoints confound many inference problems and are very important data features. Tests for undocumented changepoints from models that have independent and identically distributed errors are by now well understood. However, most climate series exhibit serial autocorrelation. Monthly, daily, or hourly series may also have periodic mean structures. This article develops a test for undocumented changepoints for periodic and autocorrelated time series. Classical changepoint tests based on sums of squared errors are modified to take into account series autocorrelations and periodicities. The methods are applied in the analyses of two climate series.

Description: This study examines the response of the tropical atmospheric and oceanic circulation to increasing greenhouse gases using a coordinated set of twenty-first-century climate model experiments performed for the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). The strength of the atmospheric overturning circulation decreases as the climate warms in all IPCC AR4 models, in a manner consistent with the thermodynamic scaling arguments of Held and Soden. The weakening occurs preferentially in the zonally asymmetric (i.e., Walker) rather than zonal-mean (i.e., Hadley) component of the tropical circulation and is shown to induce substantial changes to the thermal structure and circulation of the tropical oceans. Evidence suggests that the overall circulation weakens by decreasing the frequency of strong updrafts and increasing the frequency of weak updrafts, although the robustness of this behavior across all models cannot be confirmed because of the lack of data. As the climate warms, changes in both the atmospheric and ocean circulation over the tropical Pacific Ocean resemble “El Niño–like” conditions; however, the mechanisms are shown to be distinct from those of El Niño and are reproduced in both mixed layer and full ocean dynamics coupled climate models. The character of the Indian Ocean response to global warming resembles that of Indian Ocean dipole mode events. The consensus of model results presented here is also consistent with recently detected changes in sea level pressure since the mid–nineteenth century.

Description: An observational analysis of the Northern Hemisphere circulation during winter reveals that the leading mode of variability depends on longitude. In particular, the first EOF of the zonal-mean circulation differs over the Atlantic and Pacific sectors. These results provide motivation for a series of model runs where a simple GCM is used to investigate the effects of the zonal jet structure on the leading mode of variability in the Northern Hemisphere. Model results indicate that the leading mode of variability depends on the distance between the eddy-driven and subtropical jets. When the jets are well separated, the leading mode of variability describes latitudinal shifting of the eddy-driven jet. However, when the two jets are nearly collocated, pulsing of the combined jet dominates the variability. This change coincides with a weakening of the positive feedback between the eddies and zonal flow anomalies. These results provide a possible explanation for the reduced amplitude of the Northern Annular Mode in the Pacific sector relative to that in the Atlantic sector during Northern Hemisphere winter.

Description: A recent study of significant events of atmospheric mass depletion from the Northern Hemisphere (NH) during the extended boreal winter indicated that Southeast Asian pressure surges were an important physical mechanism that acted to channel the atmospheric mass equatorward out of the NH on a rapid time scale. This study builds upon this finding and examines both the direct and indirect roles of Southeast Asian pressure surges for a particular event of dry atmospheric mass depletion from the NH. The focus of this study is on the enhanced interhemispheric interactions and associated Southern Hemisphere (SH) tropical and extratropical responses resulting from the pressure surges. First, this study examines the conservation of dry atmospheric mass (i.e., the relationship between the dry meridional winds and the area-integrated dry air surface pressure) in the NCEP reanalysis for the 25 significant events of dry atmospheric mass depletion from the NH. Results indicate that the NCEP dry meridional winds are able to qualitatively capture the dry atmospheric mass evacuation from the NH. In a quantitative sense there is very good agreement between the wind and pressure data in the extratropics of both hemispheres. A distinct negative or southward bias in the NCEP vertically and zonally integrated dry meridional winds is apparent between 5° and 17.5°N. This southward bias was not present in the ECMWF Re-Analysis. The source of the southward bias in NCEP appears to result from a weaker analyzed ITCZ. The particular case of dry atmospheric mass depletion from the NH examined in detail is associated with an intense pressure surge over Southeast Asia. A significant enhancement of convection in the monsoon trough region of northern Australia occurs roughly 4 days after the peak intensity of the Siberian high. A low-level westerly wind burst develops in response to this enhanced zonal pressure gradient caused by the pressure surge as part of the onset of an active phase of the Australian summer monsoon. This study shows that three prominent anticyclonic circulations intensify in the SH extratropics, stretching from the south Indian Ocean to the South Pacific, beneath regions of upper-tropospheric dry atmospheric mass convergence, originating partly from the monsoon convection outflow. These anticyclonic circulations are regional manifestations of the dry atmospheric mass increase in the SH.

Description: Variability in July mean surface air temperatures from 1963 to 2003 accounted for 62% of the variance in the regional annual glacier mass balance signal for the Canadian High Arctic. A regime shift to more negative regional glacier mass balance occurred between 1986 and 1987, and is linked to a coincident shift from lower to higher mean July air temperatures. Both the interannual changes and the regime shifts in regional glacier mass balance and July air temperatures are related to variations in the position and strength of the July circumpolar vortex. In years when the July vortex is “strong” and its center is located in the Western Hemisphere, positive mass balance anomalies prevail. In contrast, highly negative mass balance anomalies occur when the July circumpolar vortex is either weak or strong without elongation over the Canadian High Arctic, and its center is located in the Eastern Hemisphere. The occurrence of westerly positioned July vortices has decreased by 40% since 1987. The associated shift to a dominantly easterly positioned July vortex was associated with an increased frequency of tropospheric ridging over the Canadian High Arctic, higher surface air temperatures, and more negative regional glacier mass balance.

Description: The authors compare short forecast errors and the balance of terms in the moisture and temperature prediction equations that lead to those errors for the Community Atmosphere Model versions 2 and 3 (CAM2 and CAM3, respectively) at T42 truncation. The comparisons are made for an individual model column from global model forecasts at the Atmospheric Radiation Measurement Program (ARM) Southern Great Plains site for the April 1997 and June–July 1997 intensive observing periods. The goal is to provide insight into parameterization errors in the CAM, which ultimately should lead to improvements in the way processes are modeled. The atmospheric initial conditions are obtained from the 40-yr ECMWF Re-Analysis (ERA-40). The land initial conditions are spun up to be consistent with those analyses. The differences between the model formulations that are responsible for the major differences in the forecast errors and/or parameterization behaviors are identified. A sequence of experiments is performed, accumulating the changes from CAM3 back toward CAM2 to demonstrate the effect of the differences in formulations. In June–July 1997 the CAM3 temperature and moisture forecast errors were larger than those of CAM2. The terms identified as being responsible for the differences are 1) the convective time scale assumed for the Zhang–McFarlane deep convection, 2) the energy associated with the conversion between water and ice of the rain associated with the Zhang–McFarlane convection parameterization, and 3) the dependence of the rainfall evaporation on cloud fraction. In April 1997 the CAM2 and CAM3 temperature and moisture forecast errors are very similar, but different tendencies arising from modifications to one parameterization component are compensated by responding changes in another component to yield the same total moisture tendency. The addition of detrainment of water in CAM3 by the Hack shallow convection to the prognostic cloud water scheme is balanced by a responding difference in the advective tendency. A halving of the time scale assumed for the Hack shallow convection was compensated by a responding change in the prognostic cloud water. Changes to the cloud fraction parameterization affect the radiative heating, which in turn modifies the stability of the atmospheric column and affects the convection. The resulting changes in convection tendency are balanced by responding changes in the prognostic cloud water parameterization tendency.

Description: The earth system model of intermediate complexity ECBilt-CLIO has been used for transient simulations of the last deglaciation and the Holocene. The forcing effects of the ice sheets, greenhouse gas concentrations, and orbital configurations are prescribed as time-varying boundary conditions. In this study two key aspects of the transient simulations are investigated, which are of broader relevance for long-term transient paleoclimate modeling: the effect of using accelerated boundary conditions and of uncertainties in the initial state. Simulations with nonaccelerated boundary conditions and an acceleration factor 10 were integrated. These simulations show that the acceleration can have a significant impact on the local climate history. In the outcropping regions of the high southern latitudes and the convective regions in the North Atlantic, the acceleration leads to damped and delayed temperature response to the boundary conditions. Furthermore, uncertainties in the initial state can strongly bias the climate trajectories in these areas over 500–700 model years. The affected oceanic regions are connected to the large heat capacities of the interior ocean, which cause a strong delay in the response to the forcing. Despite the shown difficulties with the acceleration technique, the accelerated simulations still reproduce the large-scale trend pattern of air temperatures during the Holocene from previous simulations with different models. The accelerated transient model simulation is compared with existing proxy time series at specific sites. The simulation results are in good agreement with those paleoproxies. It is shown that the transient simulations provide valuable insight into whether seasonal or annual signals are recorded in paleoproxies.

Description: The Atlantic warm pool (AWP) is a large body of warm water that comprises the Gulf of Mexico, the Caribbean Sea, and the western tropical North Atlantic. Located to its northeastern side is the North Atlantic subtropical high (NASH), which produces the tropical easterly trade winds. The easterly trade winds carry moisture from the tropical North Atlantic into the Caribbean Sea, where the flow intensifies, forming the Caribbean low-level jet (CLLJ). The CLLJ then splits into two branches: one turning northward and connecting with the Great Plains low-level jet (GPLLJ), and the other continuing westward across Central America into the eastern North Pacific. The easterly CLLJ and its westward moisture transport are maximized in the summer and winter, whereas they are minimized in the fall and spring. This semiannual feature results from the semiannual variation of sea level pressure in the Caribbean region owing to the westward extension and eastward retreat of the NASH. The NCAR Community Atmospheric Model and observational data are used to investigate the impact of the climatological annual mean AWP on the summer climate of the Western Hemisphere. Two groups of the model ensemble runs with and without the AWP are performed and compared. The model results show that the effect of the AWP is to weaken the summertime NASH, especially at its southwestern edge. The AWP also strengthens the summertime continental low over the North American monsoon region. In response to these pressure changes, the CLLJ and its moisture transport are weakened, but its semiannual feature does not disappear. The weakening of the easterly CLLJ increases (decreases) moisture convergence to its upstream (downstream) and thus enhances (suppresses) rainfall in the Caribbean Sea (in the far eastern Pacific west of Central America). Model runs show that the AWP’s effect is to always weaken the southerly GPLLJ. However, the AWP strengthens the GPLLJ’s northward moisture transport in the summer because the AWP-induced increase of specific humidity overcomes the weakening of southerly wind, and vice versa in the fall. Finally, the AWP reduces the tropospheric vertical wind shear in the main development region that favors hurricane formation and development during August–October.

Description: A 124 (1882–2005) summer record of total surface energy content consisting of time series of surface equivalent temperature (TE) and its components T (mean air temperature) and Lq/cp (moist enthalpy, denoted Lq) is developed, quality controlled, and analyzed for Columbus, Ohio, where long records of monthly dewpoint temperature are available. The analysis shows that the highest TE occurs during the summer of 1995 when both T and Lq were very high, associated with a severe midwestern heat wave. That year contrasts with the hot summers of 1930–36, when Lq and TE had relatively low or negative anomalies (low humidity) compared to those of T. Following the 1930–36 summers, T and Lq departures are much more typically the same sign in individual summers, and the two parameters develop a statistically significant high positive correlation into the twenty-first century. Mean T and Lq departures from the long-term normal have opposite signs, however, when summers are stratified either by seasonal total rainfall amounts or by the Palmer drought severity soil moisture index. Normalized trends of T, Lq, and TE are downward from 1940 to 1964 with those of TE exceeding T. Since 1965, however, significant positive T trends slightly exceed TE in magnitude and those of dewpoint temperature and Lq are comparatively lower. A highly significant upward trend in minimum temperatures especially dominates the T variability, creating a significant downward trend in the temperature range that dominates recent summer climate variability more than moisture trends. Regional moisture flux variations are largest away from Columbus, over the upper Midwest and western Atlantic Ocean, during its seasonal extremes in total surface energy.

Description: Combined simultaneous satellite observations are used to evaluate the performance of parameterizations of the microphysical and optical properties of cirrus clouds used for radiative flux computations in climate models. Atmospheric and cirrus properties retrieved from Television and Infrared Observation Satellite (TIROS-N) Operational Vertical Sounder (TOVS) observations are given as input to the radiative transfer model developed for the Met Office climate model to simulate radiative fluxes at the top of the atmosphere (TOA). Simulated cirrus shortwave (SW) albedos are then compared to those retrieved from collocated Scanner for Radiation Budget (ScaRaB) observations. For the retrieval, special care has been given to angular direction models. Three parameterizations of cirrus ice crystal optical properties are represented in the Met Office radiative transfer model. These parameterizations are based on different physical approximations and different hypotheses on crystal habit. One parameterization assumes pristine ice crystals and two ice crystal aggregates. By relating the cirrus ice water path (IWP) retrieved from the effective infrared emissivity to the cirrus SW albedo, differences between the parameterizations are amplified. This study shows that pristine crystals seem to be plausible only for cirrus with IWP less than 30 g m−2. For larger IWP, ice crystal aggregates lead to cirrus SW albedos in better agreement with the observations. The data also indicate that climate models should allow the cirrus effective ice crystal diameter (De) to increase with IWP, especially in the range up to 30 g m−2. For cirrus with IWP less than 20 g m−2, this would lead to SW albedos that are about 0.02 higher than the ones of a constant De of 55 μm.

Description: This study investigated the most recurrent coupled pattern of intraseasonal variability between midlatitude circulation and the Indian summer monsoon (ISM). The leading singular vector decomposition (SVD) pattern reveals a significant, coupled intraseasonal variation between a Rossby wave train across the Eurasian continent and the summer monsoon convection in northwestern India and Pakistan (hereafter referred to as NISM). The wave train associated with an active phase of NISM rainfall displays two high pressure anomalies, one located over central Asia and the other over northeastern Asia. They are accompanied by increased rainfall over the western Siberia plain and northern China and decreased rainfall over the eastern Mediterranean Sea and southern Japan. The circulation of the wave train shows a barotropic structure everywhere except the anomalous central Asian high, located to the northwest of India, where a heat-induced baroclinic circulation structure dominates. The time-lagged SVD analysis shows that the midlatitude wave train originates from the northeastern Atlantic and traverses Europe to central Asia. The wave train enhances the upper-level high pressure and reinforces the convection over the NISM region; meanwhile, it propagates farther toward East Asia along the waveguide provided by the westerly jet. After an outbreak of NISM convection, the anomalous central Asian high retreats westward. Composite analysis suggests a coupling between the central Asian high and the convective fluctuation in the NISM. The significance of the midlatitude–ISM interaction is also revealed by the close resemblance between the individual empirical orthogonal functions and the coupled (SVD) modes of the midlatitude circulation and the ISM. It is hypothesized that the eastward and southward propagation of the wave train originating from the northeastern Atlantic contributes to the intraseasonal variability in the NISM by changing the intensity of the monsoonal easterly vertical shear and its associated moist dynamic instability. On the other hand, the rainfall variations over the NISM reinforce the variations of the central Asian high through the “monsoon–desert” mechanism, thus reenergizing the downstream propagation of the wave train. The coupling between the Eurasian wave train and NISM may be instrumental for understanding their interaction and can provide a way to predict the intraseasonal variations of the Indian summer monsoon and East Asian summer monsoon.

Description: Using station observations of the number of days covered by snow (SCD) and snowfall over the Tibetan Plateau (TP), the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalysis, and precipitation from rain gauge stations in China for the period of 1973–2001, temporal/spatial variations of SCD over the TP and its associations with the hemispheric extratropical atmospheric circulation and East Asian summer monsoon rainfall are investigated. An increase of spring (April–May) SCD over the TP is associated with decreases of local tropospheric temperature and geopotential height in the spring and early summer (June). The anomalies in the tropospheric temperature and geopotential height show a westward propagation from the TP to western Asia as a result of the westward propagation of the anomalous wave energy. These tropospheric anomalies over the TP are connected with changes in the hemispheric extratropical atmospheric circulation along the westerly jet stream that acts as a waveguide. The increase of the spring SCD is also associated with the variation of the East Asian summer monsoon rainfall. Soil moisture in May–June might act as a bridge linking the spring snow anomaly and the subsequent summer monsoon. Corresponding to the increase of SCD, there is a significant decrease of the June 500-mb geopotential height from the TP to the western North Pacific. Meanwhile, the anomalous northeasterlies extend from Japan, through the east coast of China, to central-eastern China, which weaken the East Asian summer monsoon, leading to a decrease of surface air temperature and rainfall in the Yangtze and Hwai Rivers and an increase of rainfall in southeastern China. Additionally, the spring SCD anomaly is likely due to a variation of local synchronous snowfall, rather than previous winter SCD conditions. The spring SCD is not related to previous winter El Niño/La Niña events, but is associated with the equatorial central and eastern Pacific sea surface temperature from the subsequent summer through winter. The climatic implications for this relationship are not clear.

Description: This study investigated the accuracy and physical representation of air–sea surface heat flux estimates for the Indian Ocean on annual, seasonal, and interannual time scales. Six heat flux products were analyzed, including the newly developed latent and sensible heat fluxes from the Objectively Analyzed Air–Sea Heat Fluxes (OAFlux) project and net shortwave and longwave radiation results from the International Satellite Cloud Climatology Project (ISCCP), the heat flux analysis from the Southampton Oceanography Centre (SOC), the National Centers for Environmental Prediction reanalysis 1 (NCEP1) and reanalysis-2 (NCEP2) datasets, and the European Centre for Medium-Range Weather Forecasts operational (ECMWF-OP) and 40-yr Re-Analysis (ERA-40) products. This paper presents the analysis of the six products in depicting the mean, the seasonal cycle, and the interannual variability of the net heat flux into the ocean. Two time series of in situ flux measurements, one taken from a 1-yr Arabian Sea Experiment field program and the other from a 1-month Joint Air–Sea Monsoon Interaction Experiment (JASMINE) field program in the Bay of Bengal were used to evaluate the statistical properties of the flux products over the measurement periods. The consistency between the six products on seasonal and interannual time scales was investigated using a standard deviation analysis and a physically based correlation analysis. The study has three findings. First of all, large differences exist in the mean value of the six heat flux products. Part of the differences may be attributable to the bias in the numerical weather prediction (NWP) models that underestimates the net heat flux into the Indian Ocean. Along the JASMINE ship tracks, the four NWP modeled mean fluxes all have a sign opposite to the observations, with NCEP1 being underestimated by 53 W m−2 (the least biased) and ECMWF-OP by 108 W m−2 (the most biased). At the Arabian Sea buoy site, the NWP mean fluxes also have an underestimation bias, with the smallest bias of 26 W m−2 (ERA-40) and the largest bias of 69 W m−2 (NCEP1). On the other hand, the OAFlux+ISCCP has the best comparison at both measurement sites. Second, the bias effect changes with the time scale. Despite the fact that the mean is biased significantly, there is no major bias in the seasonal cycle of all the products except for ECMWF-OP. The latter does not have a fixed mean due to the frequent updates of the model platform. Finally, among the four products (OAFlux+ISCCP, ERA-40, NCEP1, and NCEP2) that can be used for studying interannual variability, OAFlux+ISCCP and ERA-40 Qnet have good consistency as judged from both statistical and physical measures. NCEP1 shows broad agreement with the two products, with varying details. By comparison, NCEP2 is the least representative of the Qnet variabilities over the basin scale.

Description: A series of observing system simulation experiments (OSSEs) are performed for the tropical Indian Ocean (±15° from the equator) using a simple analysis system. The analysis system projects an array of observations onto the dominant empirical orthogonal functions (EOFs) derived from an intermediate-resolution (2° × 0.5°) ocean circulation model. This system produces maps of the depth of the 20°C isotherm (D20), representing interannual variability, and the high-pass-filtered mixed layer depth (MLD), representing intraseasonal variability. The OSSEs are designed to assess the suitability of the proposed Indian Ocean surface mooring array for resolving intraseasonal to interannual variability. While the proposed array does a reasonable job of resolving the interannual time scales, it may not adequately resolve the intraseasonal time scales. A procedure is developed to rank the importance of observation locations by determining the observation array that best projects onto the EOFs used in the analysis system. OSSEs using an optimal array clearly outperform the OSSEs using the proposed array. The configuration of the optimal array is sensitive to the number of EOFs considered. The optimal array is also different for D20 and MLD, and depends on whether fixed observations are included that represent an idealized Argo array. Therefore, a relative frequency map of observation locations identified in 24 different OSSEs is compiled and a single, albeit less optimal, array that is referred to as a consolidated array is objectively determined. The consolidated array reflects the general features of the individual optimal arrays derived from all OSSEs. It is found that, in general, observations south of 8°S and off of the Indonesian coast are most important for resolving the interannual variability, while observations a few degrees south of the equator, and west of 75°E, and a few degrees north of the equator, and east of 75°E, are important for resolving the intraseasonal variability. In a series of OSSEs, the consolidated array is shown to outperform the proposed array for all configurations of the analysis system for both D20 and MLD.

Description: Annual mean net heat fluxes from ocean general circulation models (OGCMs) are systematically too low in the tropical Indian Ocean, compared to observations. In the models, only some of the geostrophic inflow replacing southward Ekman outflow is colder than the minimum sea surface temperature (MINSST). Observed heat fluxes imply that much more inflow is colder than MINSST. Since inflow below MINSST can only join the surface Ekman transport after diathermal warming, the OGCMs must underestimate diathermal effects. A crude analog of the annual mean Indian Ocean heat budget was generated, using a rectangular box model with a deep “Indo–Pacific” gap at 7°–10°S in its eastern side. Wind stress was zonal and proportional to the Coriolis parameter, so Ekman transport was spatially constant and equaled Sverdrup transport. For three experiments, zonally integrated Ekman transport was steady and southward at 10 Sv (Sv ≡ 106 m3 s−1). In steady state, a 10 Sv “Indonesian Throughflow” fed a northward western boundary current of 10 Sv, which turned eastward along the northern boundary at 10°N to feed the southward Ekman transport. Most diathermal mixing occurred within an intense eddy in the northwest corner. Some of the geostrophic inflow was at temperatures colder than MINSST (found at the northeast corner of the eddy); it must warm to MINSST via diathermal mixing. Northern boundary upwelling exceeded the 10-Sv Ekman transport. The excess warms as it recirculates around the eddy, apparently supplying the heat to warm inflow below MINSST. In an experiment using the “flux-corrected transport” (FCT) scheme, diathermal mixing occurred in the strongly sheared currents around the eddy. However the Richardson number never became low enough to drive strong diathermal mixing, perhaps because (like that of other published models) the present model’s vertical resolution was too coarse. In three experiments, the dominant mixing was caused by horizontal diffusion, spurious convective overturn, and numerical mixing invoked by the FCT scheme, respectively. All three mixing mechanisms are physically suspect; such model problems (if widespread) must be resolved before the mismatch between observed and modeled heat fluxes can be addressed. However, the fact that the density profile at the western boundary must be hydrostatically stable places a lower limit on the area-integrated heat fluxes. Results from the three main experiments—and from many published OGCMs—are quite close to this lower limit.

Description: In this paper, a series of observing system simulation experiments (OSSEs) are used to study the design of a proposed array of instrumented moorings in the Indian Ocean (IO) outlined by the IO panel of the Climate Variability and Predictability (CLIVAR) Project. Fields of the Ocean Topography Experiment (TOPEX)/Poseidon (T/P) and Jason sea surface height (SSH) and sea surface temperature (SST) are subsampled to simulate dynamic height and SST data from the proposed array. Two different reduced-order versions of the Kalman filter are used to reconstruct the original fields from the simulated observations with the objective of determining the optimal deployment of moored platforms and to address the issue of redundancy and array simplification. The experiments indicate that, in terms of the reconstruction of SSH and SST, the location of the subjectively proposed array compareS favorably with the optimally defined one. The only significant difference between the proposed IO array and the optimal array is the lack of justification for increasing the latitudinal resolution near the equator (i.e., moorings 1.5°S and 1.5°N). An analysis of the redundancy also identifies the equatorial region as the one with the largest amount of redundant information. Thus, in the context of these fields, these results may help define the prioritization of its deployment or redefine the array to extend its latitudinal extent while maintaining the same amount of stations.

Description: The sensitivity of the atmospheric radiation budget to ignoring small ice particles (D ≤ 100 μm) in parameterization of the mean effective size of ice particles was investigated by using the Canadian Centre for Climate Modelling and Analysis (CCCma) third-generation general atmospheric circulation model (AGCM3). The results indicate that small ice particles play two crucial roles in the radiative transfer that influence the simulated climate. First, they inhibit the IR radiation from escaping to space and, second, they enhance the scattering of solar radiation. On average, these two effects tend to partially cancel each other out. However, based on AGCM simulations, the small ice crystals make clouds more opaque to IR radiation. Generally, 5-yr seasonally averaged GCM results suggest that the strongest anomalies in outgoing longwave radiation (OLR) are found in the Tropics, reaching 15 to 25 W m−2 in areas where cold high cirrus anvil clouds are prevalent. The global average change in net cloud radiative forcing was 2.4 W m−2 in June–August (JJA) and 1.7 W m−2 in December–February (DJF). The change in globally averaged 5-yr mean cloud forcing was close to 1.9 W m−2. When the small particles were included, the globally averaged 5-yr mean precipitation decreased by about 8%, but cloudiness increased only slightly (by 2%). The 5-yr averaged global mean surface (screen) temperature also increased slightly (about 0.2°C) when the small ice particles were included.

Description: In this paper the authors explore the spatial and temporal patterns of lightning strikes in northern Australia for the first time. In particular, the possible relationships between lightning strikes and elevation, vegetation type, and fire scars (burned areas) are examined. Lightning data provided by the Bureau of Meteorology were analyzed for a 6-yr period (1998–2003) over the northern, southern, and coastal regions of the Northern Territory (NT) through the use of Geographical Information Systems (GIS) to determine the spatial and temporal characteristics of lightning strikes. It was determined that the highest densities of lightning strikes occurred during the monsoon transitional period (dry to wet) and during the active monsoon periods, when atmospheric moisture is highest. For the period of this study, lightning was far more prevalent over the northern region (1.21 strikes per km2 yr−1) than over the southern (0.58 strikes per km2 yr−1) and coastal regions (0.71 strikes per km2 yr−1). Differences in vegetation cover were suggested to influence the lightning distribution over the northern region of the NT, but no relationship was found in the southern region. Lightning strikes in the southern region showed a positive relationship with elevations above 800 m, but no relationship was found in the northern region, which could be due to the low-lying topography of the area. A comparison of lightning densities between burned and unburned areas showed high variability; however, the authors suggest that, under ideal atmospheric conditions, large-scale fire scars (〉500 m) could produce lightning strikes triggered by either enhanced free convection or mesoscale circulations.

Description: Present-day OGCMs give low values of annual mean net heat flux (AMNHF) in the tropical Indian Ocean, compared to climatologies. AMNHF generation is examined in an open-boundary model of this region with realistic coastlines. In the first two of three experiments only annual mean wind stresses were applied so that a modified form of the “minimum depth criterion” of the previous paper would be applicable. Area-integrated AMNHF was well below observed values, despite the fact that western boundary inflow was substantially deeper and colder than was expected from the modified minimum depth estimate. The model showed large “spikes” in the gradient of “depth-integrated steric height” (DISH) along the western boundary, coinciding with coastline steps (which were absent in the previous paper). Most diapycnal entrainment occurred next to the coast, near these steps. In a third experiment a seasonal cycle of wind stress was added to the same annual mean. Annual mean diapycnal mixing and entrainment increased and the western boundary inflow deepened, resulting in substantially greater AMNHF for the same annual mean Ekman transports. However, area-integrated AMNHF was still well below the mean of directly observed surface fluxes. The recirculation around the “Great Whirl” doubled, permitting more cold water crossing the equator in one year to mix with recirculated water generated in a previous year. Entrainment up to the surface thus went by stages, over more than one year. The increased Great Whirl was related to stronger annual mean curls of nonlinear terms in the momentum equation, while the deeper entrainment was caused by stronger annual mean diapycnal mixing. In all experiments, diapycnal mixing was primarily due to the “flux corrected transport” (FCT) advective scheme, which in effect replaces spurious convective overturn by numerical diffusion. More research is needed to solve such problems, but sensitivity of AMNHF in OGCMs to time-varying forcing—due to seasonal, intraseasonal, or baroclinic instability—may offer a new source of climate predictability.

Description: The mean characteristics and trends of Southern Hemisphere (SH) winter extratropical cyclones occurring at six levels of the troposphere over the period 1979–2001 have been investigated using the 40-yr ECMWF Re-Analysis (ERA-40) data. Cyclonic systems were identified with the Melbourne University cyclone finding and tracking scheme. This study shows that mean sea level pressure (MSLP) cyclones are more numerous, more intense, smaller, deeper, and slower moving than higher-level cyclones. The novel vertical tracing scheme devised for this research revealed that about 52% of SH winter MSLP cyclones have a vertically well organized structure, extending through to the 500-hPa level. About 80% of these vertically coherent SH cyclones keep their westward tilt until the surface cyclones reach their maximum depths, and the mean distance is 300 km between the surface and the 500-hPa level cyclone centers when the surface cyclones obtain their maturity. According to the authors’ definition of vertical organization, explosively developing cyclones are vertically very well organized systems, whose surface development is antecedent to their 500-hPa level counterpart. Over 1979–2001 cyclones have increased in their system density, intensity, and translational velocity but decreased in their scale at almost all levels. However, some of the trends are not statistically significant. The proportion of vertically well organized systems in the entire population of SH winter extratropical cyclones has considerably increased over the last 23 yr, and the mean distance between the surface and the 500-hPa- level cyclone centers has decreased. Such changes in vertical organization of extratropical cyclones are statistically significant at the 95% confidence level.

Description: The automatic tracking technique used by Thorncroft and Hodges has been used to identify coherent vorticity structures at 850 hPa over West Africa and the tropical Atlantic in the 40-yr ECMWF Re-Analysis. The presence of two dominant source regions, north and south of 15°N over West Africa, for storm tracks over the Atlantic was confirmed. Results show that the southern storm track provides most of the storms that reach the main development region where most tropical cyclones develop. There exists marked seasonal variability in location and intensity of the storms leaving the West African coast, which may influence the likelihood of downstream intensification and longevity. There exists considerable year-to-year variability in the number of West African storm tracks, both in numbers over the land and continuing out over the tropical Atlantic Ocean. While the low-frequency variability is well correlated with Atlantic tropical cyclone activity, West African rainfall, and SSTs, the interannual variability is found to be uncorrelated with these. In contrast, variance of the 2–6-day-filtered meridional wind, which provides a synoptic-scale measure of African easterly wave activity, shows a significant, positive correlation with tropical cyclone activity at interannual time scales.

Description: Two atmospheric general circulation models (AGCMs), differing in numerics and physical parameterizations, are employed to test the hypothesis that El Niño–induced sea surface temperature (SST) anomalies in the tropical Indian Ocean impact considerably the Northern Hemisphere extratropical circulation anomalies during boreal winter [January–March +1 (JFM +1)] of El Niño years. The hypothesis grew out of recent findings that ocean dynamics influence SST variations over the southwest Indian Ocean (SWIO), and these in turn impact local precipitation. A set of ensemble simulations with the AGCMs was carried out to assess the combined and individual effects of tropical Pacific and Indian Ocean SST anomalies on the extratropical circulation. To elucidate the dynamics responsible for the teleconnection, solutions were sought from a linear version of one of the AGCMs. Both AGCMs demonstrate that the observed precipitation anomalies over the SWIO are determined by local SST anomalies. Analysis of the circulation response shows that over the Pacific–North American (PNA) region, the 500-hPa height anomalies, forced by Indian Ocean SST anomalies, oppose and destructively interfere with those forced by tropical Pacific SST anomalies. The model results validated with reanalysis data show that compared to the runs where only the tropical Pacific SST anomalies are specified, the root-mean-square error of the height anomalies over the PNA region is significantly reduced in runs in which the SST anomalies in the Indian Ocean are prescribed in addition to those in the tropical Pacific. Among the ensemble members, both precipitation anomalies over the SWIO and the 500-hPa height over the PNA region show high potential predictability. The solutions from the linear model indicate that the Rossby wave packets involved in setting up the teleconnection between the SWIO and the PNA region have a propagation path that is quite different from the classical El Niño–PNA linkage. The results of idealized experiments indicate that the Northern Hemisphere extratropical response to Indian Ocean SST anomalies is significant and the effect of this response needs to be considered in understanding the PNA pattern during El Niño years. The results presented herein suggest that the tropical Indian Ocean plays an active role in climate variability and that accurate observation of SST there is of urgent need.

Description: Subtropical shallow cumulus convection is shown to play an important role in tropical climate dynamics, in which convective mixing between the atmospheric boundary layer and the free troposphere initiates a chain of large-scale feedbacks. It is found that the presence of shallow convection in the subtropics helps set the width and intensity of oceanic ITCZs, a mechanism here termed the shallow cumulus humidity throttle because of the control exerted on the moisture supply to the deep convection zones. These conclusions are reached after investigations based on a tropical climate model of intermediate complexity, with sufficient vertical degrees of freedom to capture (i) the effects of shallow convection on the boundary layer moisture budget and (ii) the dependency of deep convection on the free-tropospheric humidity. An explicit shallow cumulus mixing time scale in this simple parameterization is varied to assess sensitivity, with moist static energy budget analysis aiding to identify how the local effect of shallow convection is balanced globally. A reduction in the mixing efficiency of shallow convection leads to a more humid atmospheric mixed layer, and less surface evaporation, with a drier free troposphere outside of the convecting zones. Advection of drier free-tropospheric air from the subtropics by transients such as dry intrusions, as well as by mean inflow, causes a substantial narrowing of the convection zones by inhibition of deep convection at their margins. In the tropical mean, the reduction of convection by this narrowing more than compensates for the reduction in surface evaporation. Balance is established via a substantial decrease in tropospheric temperatures throughout the Tropics, associated with the reduction in convective heating. The temperature response—and associated radiative contribution to the net flux into the column—have broad spatial scales, while the reduction of surface evaporation is concentrated outside of the convecting zones. This results in differential net flux across the convecting zone, in a sense that acts to destabilize those areas that do convect. This results in stronger large-scale convergence and more intense convection within a narrower area. Finally, mixed layer ocean experiments show that in a coupled ocean–atmosphere system this indirect feedback mechanism can lead to SST differences up to +2 K between cases with different shallow cumulus mixing time, tending to counteract the direct radiative impact of low subtropical clouds on SST.

Description: Shortly after the advent of the first imaging passive microwave sensor on board a research satellite an anomalous climate feature was observed within the Weddell Sea. During the years 1974–1976, a 250 × 103 km2 area within the seasonal sea ice cover was virtually free of winter sea ice. This feature, the Weddell Polynya, was created as sea ice formation was inhibited by ocean convection that injected relatively warm deep water into the surface layer. Though smaller, less persistent polynyas associated with topographically induced upwelling at Maud Rise frequently form in the area, there has not been a reoccurrence of the Weddell Polynya since 1976. Archived observations of the surface layer salinity within the Weddell gyre suggest that the Weddell Polynya may have been induced by a prolonged period of negative Southern Annular Mode (SAM). During negative SAM the Weddell Sea experiences colder and drier atmospheric conditions, making for a saltier surface layer with reduced pycnocline stability. This condition enables Maud Rise upwelling to trigger sustained deep-reaching convection associated with the polynya. Since the late 1970s SAM has been close to neutral or in a positive state, resulting in warmer, wetter conditions over the Weddell Sea, forestalling repeat of the Weddell Polynya. A contributing factor to the Weddell Polynya initiation may have been a La Niña condition, which is associated with increased winter sea ice formation in the polynya area. If the surface layer is made sufficiently salty due to a prolonged negative SAM period, perhaps aided by La Niña, then Maud Rise upwelling meets with positive feedback, triggering convection, and a winter persistent Weddell Polynya.

Description: The North American monsoon (NAM) is a prominent summertime feature over northwestern Mexico and the southwestern United States. It is characterized by a distinct shift in midlevel winds from westerly to easterly as well as a sharp, marked increase in rainfall. This maximum in rainfall accounts for 60%–80% of the annual precipitation in northwestern Mexico and nearly 40% of the yearly rainfall over the southwestern United States. Gulf surges, or coastally trapped disturbances that occur over the Gulf of California, are important mechanisms in supplying the necessary moisture for the monsoon and are hypothesized in previous studies to be initiated by the passage of a tropical easterly wave (TEW). Since the actual number of TEWs varies from year to year, it is possible that TEWs are responsible for producing some of the interannual variability in the moisture flux and rainfall seen in the NAM. To explore the impact of TEWs on the NAM, four 1-month periods are chosen for study that represent a reasonable variability in TEW activity. Two continuous month-long simulations are produced for each of the selected months using the Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model. One simulation is a control run that uses the complete boundary condition data, whereas a harmonic analysis is used to remove TEWs with periods of approximately 3.5 to 7.5 days from the model boundary conditions in the second simulation. These simulations with and without TEWs in the boundary conditions are compared to determine the impact of the waves on the NAM. Fields such as meridional moisture flux, rainfall totals, and surge occurrences are examined to define similarities and differences between the model runs. Results suggest that the removal of TEWs not only reduces the strength of gulf surges, but also rearranges rainfall over the monsoon region. Results further suggest that TEWs influence rainfall over the Southern Plains of the United States, with TEWs leading to less rainfall in this region. While these results are only suggestive, since rainfall is the most difficult model forecast parameter, it may be that TEWs alone can explain part of the inverse relationship between NAM and Southern Plains rainfall.

Description: Idealized experiments are conducted using a GCM coupled to a 20-m slab ocean model to examine the short-term response to an initial localized positive equatorial SST anomaly, or “hot spot.” A hot spot is imposed upon an aquaplanet with globally uniform 28°C SST, insolation, and trace gas concentrations designed to mimic tropical warm pool conditions. No boundary condition or external parameter other than the Coriolis parameter varies with latitude. A 15-member ensemble is initiated using random atmospheric initial conditions. A 2°C equatorial warm anomaly is switched on, along with ocean coupling (day 0). Enhanced deep convection rapidly develops near the hot spot, forcing an anomalous large-scale circulation that resembles the linear response of a dry atmosphere to a localized heating, as in the Gill model. Enhanced convection, the anomalous large-scale circulation, and enhanced wind speed peak in amplitude at about day 15. Enhanced latent heat fluxes driven primarily by an increase in vector mean wind damp the anomalous heat content of the ocean near the hot spot before day 20. Between day 20 and day 50, suppressed latent heat fluxes due to suppressed synoptic eddy variance cause a warming of the remote Tropics in regions of anomalous low-level easterly flow. This wind-driven evaporative atmosphere–ocean exchange results in a 60–70-day oscillation in tropical mean oceanic heat content, accompanied by a compensating out-of-phase oscillation in vertically integrated atmospheric moist static energy. Beyond day 70 of the simulation, positive SST anomalies are found across much of the tropical belt. These slowly decay toward the 28°C background state.

Description: The authors demonstrate that variability in vegetation cover can potentially influence oceanic variability through the atmospheric bridge. Experiments aimed at isolating the impact of variability in forest cover along the poleward side of the Asian boreal forest on North Pacific SSTs are performed using the fully coupled model, Fast Ocean Atmosphere Model–Lund Potsdam Jena (FOAM-LPJ), with dynamic atmosphere, ocean, and vegetation. The northern edge of the simulated Asian boreal forest is characterized by substantial variability in annual forest cover, with an east–west dipole pattern marking its first EOF mode. Simulations in which vegetation cover is allowed to vary over north/central Russia exhibit statistically significant greater SST variance over the Kuroshio Extension. Anomalously high forest cover over North Asia supports a lower surface albedo with higher temperatures and lower sea level pressure, leading to a reduction in cold advection into northern China and in turn a decrease in cold air transport into the Kuroshio Extension region. Variability in the large-scale circulation pattern is indirectly impacted by the aforementioned vegetation feedback, including the enhancement in upper-level jet wind variability along the north–south flanks of the East Asian jet stream.

Description: The extent to which soil water storage can support an average dry season evapotranspiration (ET) is investigated using observations from the Rebio Jarú site for the period of 2000 to 2002. During the dry season, when total rainfall is less than 100 mm, the soil moisture storage available to root uptake in the top 3-m layer is sufficient to maintain the ET rate, which is equal to or higher than that in the wet season. With a normal or less-than-normal dry season rainfall, more than 75% of the ET is supplied by soil water below 1 m, whereas during a rainier dry season, about 50% of ET is provided by soil water from below 1 m. Soil moisture below 1-m depth is recharged by rainfall during the previous wet season: dry season rainfall rarely infiltrates to this depth. These results suggest that, even near the southern edge of the Amazon forest, seasonal and moderate interannual rainfall deficits can be mitigated by an increase in root uptake from deeper soil. How dry season ET varies geographically within the Amazon and what might control its geographic distribution are examined by comparing in situ observations from 10 sites from different areas of Amazonia reported during the last two decades. Results show that the average dry season ET varies less than 1 mm day−1 or 30% from the driest to nearly the wettest parts of Amazonia and is largely correlated with the change of surface net radiation of 25% and 30%. Thus the geographic variation of the average dry season ET appears to be mainly determined by the surface radiation.

Description: The 227-yr daily precipitation record gathered for Seoul, South Korea, represents one of the longest instrumental measurements, which provides an exceptional opportunity for detecting climate singularity (a property of phase locking to annual cycle) of extreme weather events and multidecadal–centennial variability of the rainy season structure. From late June to early September, the occurrence of heavy rain events shows a climatological quasi-biweekly oscillation. The rainy season characteristics, including the dates of onset, retreat, summit, and the duration, all show significant centennial variations. The rainy season summit shows a tendency toward delayed occurrence, which changed from the 37th pentad (P37; 30 June–4 July) during the 1778–1807 period to P44 (4–8 August) during the 1975–2004 period. The amplitude of the interannual (2–6 yr) variation of summer precipitation shows a prominent fluctuation with a 50-yr rhythm. A notable climatological break (around 9–13 August) divides the rainy season into a changma (Korean for continuous rain period) and a post-changma period. The major modes of subseasonal variability of the rainy season are characterized by an advanced changma and an enhanced post-changma, respectively. The former is dominated by biennial variation, whereas the latter has a major 5-yr spectral peak, suggesting that the processes leading to their variability are different. The occurrence of severe drought events exhibits a 4-yr spectral peak along with large power on a centennial time scale, while the severe flood events have a spectral peak at 3 and 19 yr, respectively. The remarkable climate variability in Seoul rainfall suggests that trends detected by using a 50-yr-or-shorter precipitation record likely reflect natural variability.

Description: Daily variations in Australian rainfall and surface temperature associated with the Southern Hemisphere annular mode (SAM) are documented using observations for the period 1979–2005. The high index polarity of the SAM is characterized by a poleward contraction of the midlatitude westerlies. During winter, the high index polarity of the SAM is associated with decreased daily rainfall over southeast and southwest Australia, but during summer it is associated with increased daily rainfall on the southern east coast of Australia and decreased rainfall in western Tasmania. Variations in the SAM explain up to ∼15% of the weekly rainfall variance in these regions, which is comparable to the variance accounted for by the El Niño–Southern Oscillation, especially during winter. The most widespread temperature anomalies associated with the SAM occur during the spring and summer seasons, when the high index polarity of the SAM is associated with anomalously low maximum temperature over most of central/eastern subtropical Australia. The regions of decreased maximum temperature are also associated with increased rainfall. Implications for recent trends in Australian rainfall and temperature are discussed.

Description: The physical processes associated with the ∼70-yr period climate mode, known as the Atlantic multidecadal oscillation (AMO), are examined. Based on analyses of observational data, a deterministic mechanism relying on atmosphere–ocean–sea ice interactions is proposed for the AMO. Variations in the thermohaline circulation are reflected as uniform sea surface temperature anomalies in the North Atlantic. These anomalies are associated with a hemispheric wavenumber-1 sea level pressure (SLP) structure in the atmosphere that is amplified through atmosphere–ocean interactions in the North Pacific. The SLP pattern and its associated wind field affect the sea ice export through Fram Strait, the freshwater balance in the northern North Atlantic, and consequently the strength of the large-scale ocean circulation. It generates sea surface temperature anomalies with opposite signs in the North Atlantic and completes a negative feedback. The authors find that the time scale of the cycle is associated with the thermohaline circulation adjustment to freshwater forcing, the SST response to it, the oceanic adjustment in the North Pacific, and the sea ice response to the wind forcing. Finally, it is argued that the Great Salinity Anomaly in the late 1960s and 1970s is part of AMO.

Description: A recently developed variance decomposition approach is applied to predict seasonal mean 500-hPa geopotential height anomalies in the Southern Hemisphere. In terms of predictability of both the winter and summer height fields, the Southern Oscillation and the Southern Annular Mode are identified as the first and second most important factors affecting the variability. Based on this study, a statistical prediction scheme has been developed. The linear trend in the leading empirical orthogonal function of the height field, the November Southern Annular Mode index, the austral spring Niño-3 index, and the November Coral Sea index are identified as the main predictors for the summer height field, while the March–May Southern Annular Mode index, the May Niño-4 index, and the austral autumn central Indian Ocean index are the main predictors for the winter height field. The predictive skill in forecasts of National Centers for Environmental Prediction–National Center for Atmospheric Research and European Centre for Medium-Range Weather Forecasts reanalysis 500-hPa geopotential height anomaly fields, in terms of a spatiotemporal anomaly correlation, is considerably higher than a single prediction achieved by a coupled general circulation seasonal forecast model.

Description: The adjustment of the tropical climate outside the Pacific (the “remote Tropics”) to the abrupt onset of El Niño conditions is examined in a tropical atmosphere model that assumes simplified vertical structure and quasi-equilibrium (QE) convective closure. The El Niño signal is rapidly (∼1 week) communicated to the remote Tropics via an eastward-propagating Kelvin-like wave that induces both anomalous subsidence and tropospheric warming. Widespread reductions in convective precipitation occur in conjunction with the spreading of the temperature and subsidence anomalies. The remote rainfall suppression peaks roughly 5–15 days after the initiation of El Niño conditions, after which the anomalous remote rainfall field recovers to a state characterized by a smaller remote areal mean rainfall deficit and the appearance of localized positive rainfall anomalies. The initial remote precipitation reduction after El Niño onset is tied to both tropospheric warming (i.e., stabilization of the troposphere to deep convection) and the suppression of remote humidity levels; recovery of the initial deficits occurs as feedbacks modulate the subsequent evolution of humidity anomalies in the tropospheric column. Apart from the short-term response, there is a longer-term adjustment of the remote climate related to the thermal inertia of the underlying surface: surface thermal disequilibrium, which is related to the depth of the ocean mixed layer, maintains larger precipitation deficits than would be expected for equilibrated conditions. This result supports a previous prediction by one of the authors for a significant disequilibrium mechanism in the precipitation teleconnection to El Niño resulting from the local vertical coupling of the troposphere to the surface through moist convection.

Description: The correlation between parameters characterizing observed westerly wind bursts (WWBs) in the equatorial Pacific and the large-scale SST is analyzed using singular value decomposition. The WWB parameters include the amplitude, location, scale, and probability of occurrence for a given SST distribution rather than the wind stress itself. This approach therefore allows for a nonlinear relationship between the SST and the wind signal of the WWBs. It is found that about half of the variance of the WWB parameters is explained by only two large-scale SST modes. The first mode represents a developed El Niño event, while the second mode represents the seasonal cycle. More specifically, the central longitude of WWBs, their longitudinal extent, and their probability seem to be determined to a significant degree by the ENSO-driven signal. The amplitude of the WWBs is found to be strongly influenced by the phase of the seasonal cycle. It is concluded that the WWBs, while partially stochastic, seem an inherent part of the large-scale deterministic ENSO dynamics. Implications for ENSO predictability and prediction are discussed.

Description: This paper presents a framework based on Bayesian regression and constrained least squares methods for incorporating prior beliefs in a linear regression problem. Prior beliefs are essential in regression theory when the number of predictors is not a small fraction of the sample size, a situation that leads to overfitting—that is, to fitting variability due to sampling errors. Under suitable assumptions, both the Bayesian estimate and the constrained least squares solution reduce to standard ridge regression. New generalizations of ridge regression based on priors relevant to multimodel combinations also are presented. In all cases, the strength of the prior is measured by a parameter called the ridge parameter. A “two-deep” cross-validation procedure is used to select the optimal ridge parameter and estimate the prediction error. The proposed regression estimates are tested on the Development of a European Multimodel Ensemble System for Seasonal to Interannual Prediction (DEMETER) hindcasts of seasonal mean 2-m temperature over land. Surprisingly, none of the regression models proposed here can consistently beat the skill of a simple multimodel mean, despite the fact that one of the regression models recovers the multimodel mean in a suitable limit. This discrepancy arises from the fact that methods employed to select the ridge parameter are themselves sensitive to sampling errors. It is plausible that incorporating the prior belief that regression parameters are “large scale” can reduce overfitting and result in improved performance relative to the multimodel mean. Despite this, results from the multimodel mean demonstrate that seasonal mean 2-m temperature is predictable for at least three months in several regions.

Description: Rainfall in South America as simulated by a 24-ensemble member of the ECHAM 4.5 atmospheric general circulation model is compared and contrasted with observations (in areas in which data are available) for the period 1976–2001. Emphasis is placed on determining the onset and end of the rainy season, from which its length and rain rate are determined. It is shown that over large parts of the domain the onset and ending dates are well simulated by the model, with biases of less than 10 days. There is a tendency for model onset to occur early and ending to occur late, resulting in a simulated rainy season that is on average too long in many areas. The model wet season rain rate also tends to be larger than observed. To estimate the relative importance of errors in wet season length and rain rate in determining biases in the annual total, adjusted totals are computed by substituting both the observed climatological wet season length and rate for those of the model. Problems in the rain rate generally are more important than problems in the length. The wet season length and rain rate also contribute substantially to interannual variations in the annual total. These quantities are almost independent, and it is argued that they are each associated with different mechanisms. The observed onset dates almost always lie within the range of onset of the ensemble members, even in the areas with a large model onset bias. In some areas, though, the model does not perform well. In southern Brazil the model ensemble average onset always occurs in summer, whereas the observations show that winter is often the wettest period. Individual members, however, do occasionally show a winter rainfall peak. In southern Northeast Brazil the model has a more distinct rainy season than is observed. In the northwest Amazon the model annual cycle is shifted relative to that observed, resulting in a model bias. No interannual relationship between model and observed onset dates is expected unless onset in the model and observations has a mutual relationship with SST anomalies. In part of the near-equatorial Amazon, there does exist an interannual relationship between onset dates. Previous studies have shown that in this area there is a relationship between SST anomalies and variations in seasonal total rainfall.

Description: The skill of state-of-the-art operational dynamical models in predicting the two most important modes of variability in the Northern Hemisphere extratropical atmosphere, the North Atlantic Oscillation (NAO) and Pacific–North American (PNA) teleconnection patterns, is investigated at time scales ranging from daily to seasonal. Two uncoupled atmospheric models used for deterministic forecasting in the short to medium range as well as eight fully coupled atmosphere–land–ocean forecast models used for monthly and seasonal forecasting are examined and compared. For the short to medium range, the level of forecast skill for the two indices is higher than that for the entire Northern Hemisphere extratropical flow. The forecast skill of the PNA is higher than that of the NAO. The forecast skill increases with the magnitude of the NAO and PNA indices, but the relationship is not pronounced. The probability density function (PDF) of the NAO and PNA indices is negatively skewed, in agreement with the distribution of skewness of the geopotential field. The models maintain approximately the observed PDF, including the negative skewness, for the first week. Extreme negative NAO/PNA events have larger absolute values than positive extremes in agreement with the negative skewness of the two indices. Recent large extreme events are generally well forecasted by the models. On the intraseasonal time scale it is found that both NAO and PNA have lingering forecast skill, in contrast to the Northern Hemisphere extratropical flow as a whole. This fact offers some hope for extended range forecasting, even though the skill is quite low. No conclusive positive benefit is seen from using higher horizontal resolution or coupling to the oceans. On the monthly and seasonal time scales, the level of forecast skill for the two indices is generally quite low, with the exception of winter predictions at short lead times.

Description: Sea ice variability in the North Pacific and its associations with the east Asia–North Pacific winter climate were investigated using observational data. Two dominant modes of sea ice variability in the North Pacific were identified. The first mode features a dipole pattern between the Sea of Okhotsk and the Bering Sea. The second mode is characterized by more uniform ice changes throughout the North Pacific. Using the principal components of the two dominant modes as the indices (PC1 and PC2), analyses show that the positive phases of PC1 feature a local warming (cooling) in the Sea of Okhotsk (the Bering Sea), which is associated with the formation of the anomalous anticyclone extending from the northern Pacific to Siberia, accompanied by a weakening of the east Asian jet stream and trough. The associated anomalous southeasterlies/easterlies reduce the climatological northwesterlies/westerlies, leading to warm and wet conditions in northeast China and central Siberia. The positive phases of PC2 are characterized by a strong local warming in the northern Pacific that coincides with the anomalous cyclone occupying the entire North Pacific, accompanied by a strengthening of the east Asia jet stream and trough. The associated anomalous northerlies intensify the east Asian winter monsoon (EAWM), leading to cold and dry conditions in the east coast of Asia. The intensified EAWM also strengthens the local Hadley cell, which in turn strengthens the east Asian jet stream and leads to a precipitation deficit over subtropical east Asia. The linkages between PC1 and PC2 and large-scale modes of climate variability were also discussed. It is found that PC1 is a better indicator than the Arctic Oscillation of the recent Siberian warming, whereas PC2 may be a valuable predictor of EAWM.

Description: This study assesses the impact of two different remote sensing–derived leaf area index (RSLAI) datasets retrieved from the same source (i.e., Advanced Very High Resolution Radiometer measurements) on a general circulation model’s (GCM) seasonal climate simulations as well as the mechanisms that lead to the improvement in simulations over several regions. Based on the analysis of these two RSLAI datasets for 17 yr from 1982 to 1998, their spatial distribution patterns and characteristics are discussed. Despite some disagreements in the RSLAI magnitudes and the temporal variability between these two datasets over some areas, their effects on the simulation of near-surface climate and the regions with significant impact are generally similar to each other. Major disagreements in the simulated climate appear in a few limited regions. The GCM experiment using the RSLAI and other satellite-derived land surface products showed substantial improvements in the near-surface climate in the East Asian and West African summer monsoon areas and boreal forests of North America compared to the control experiment that used LAI extrapolated from limited ground surveys. For the East Asia and northwest U.S. regions, the major role of RSLAI changes is in partitioning the net radiative energy into latent and sensible heat fluxes, which results in discernable warming and decrease of precipitation due to the smaller RSLAI values compared to the control. Meanwhile, for the West African semiarid regions, where the LAI difference between RSLAI and control experiments is negligible, the decrease in surface albedo caused by the high vegetation cover fraction in the satellite-derived dataset plays an important role in altering local circulation that produces a positive feedback in land/atmosphere interaction.

Description: Climate records over the last millennium place the twentieth-century warming in a longer historical context. Reconstructions of millennial temperatures show a wide range of variability, raising questions about the reliability of currently available reconstruction techniques and the uniqueness of late-twentieth-century warming. A calibration method is suggested that avoids the loss of low-frequency variance. A new reconstruction using this method shows substantial variability over the last 1500 yr. This record is consistent with independent temperature change estimates from borehole geothermal records, compared over the same spatial and temporal domain. The record is also broadly consistent with other recent reconstructions that attempt to fully recover low-frequency climate variability in their central estimate. High variability in reconstructions does not hamper the detection of greenhouse gas–induced climate change, since a substantial fraction of the variance in these reconstructions from the beginning of the analysis in the late thirteenth century to the end of the records can be attributed to external forcing. Results from a detection and attribution analysis show that greenhouse warming is detectable in all analyzed high-variance reconstructions (with the possible exception of one ending in 1925), and that about a third of the warming in the first half of the twentieth century can be attributed to anthropogenic greenhouse gas emissions. The estimated magnitude of the anthropogenic signal is consistent with most of the warming in the second half of the twentieth century being anthropogenic.

Description: How unstable is the tropical ocean–atmosphere system? Are two successive El Niño events independent, or are they part of a continual (perhaps weakly damped) cycle sustained by random atmospheric disturbances? How important is energy dissipation for ENSO dynamics? These closely related questions are frequently raised in connection with several climate problems ranging from El Niño predictability to the impact of atmospheric “noise” on ENSO. One of the factors influencing the system’s stability and other relevant properties is the damping (decay) time scale for the thermocline anomalies associated with the large-scale oceanic motion. Here this time scale is estimated by considering energy balance and net energy dissipation in the tropical ocean and it is shown that there are two distinct dissipative regimes: in the interannual frequency band the damping rate is approximately (2.3 yr)−1; however, in a near-annual frequency range the damping appears to be much stronger, roughly (8 months)−1.

Description: A methodology is proposed in which a few prognostic variables of a regional climate model (RCM) are strongly constrained at certain wavelengths to what is prescribed from the bias-corrected atmospheric general circulation model (AGCM; driver model) integrations. The goal of this strategy is to reduce the systematic errors in a RCM that mainly arise from two sources: the lateral boundary conditions and the RCM errors. Bias correction (which essentially corrects the climatology) of the forcing from the driving model addresses the former source while constraining the solution of the RCM beyond certain relatively large wavelengths in the regional domain [also termed as scale-selective bias correction (SSBC)] addresses the latter source of systematic errors in RCM. This methodology is applied to experiments over the South American monsoon region. It is found that the combination of bias correction and SSBC on the nested variables of divergence, vorticity, and the log of surface pressure of an RCM yields a major improvement in the simulation of the regional climate variability over South America from interannual to intraseasonal time scales. The basis for such a strategy is derived from a systematic empirical approach that involved over 100 regional seasonal climate integrations.

Description: The isotopic composition of present-day Antarctic snow is simulated for the period September 1980–August 2002 using a Rayleigh-type isotope distillation model in combination with backward trajectory calculations with 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) data as meteorological input. Observed spatial isotopic gradients are correctly reproduced, especially in West Antarctica and in the coastal areas. However, isotopic depletion of snow on the East Antarctic plateau is underestimated, a problem that is also observed in general circulation models equipped with isotope tracers. The spatial isotope–temperature relation varies strongly, which indicates that this widely used relation is not applicable to all sites and temporal scales. Spatial differences in the seasonal amplitude are identified, with maximum values in the Antarctic interior and hardly any seasonal isotope signature in Marie Byrd Land, West Antarctica. The modeled signature of deuterium excess remains largely preserved during the last phase of transport, though the simulated relation of deuterium excess with δ18O suggests that parameterizations of kinetic isotopic fractionation can be improved.

Description: From a dataset of weather observations from land stations worldwide, about 5400 stations were selected as having long periods of record with cloud-type information; they cover all continents and many islands. About 185 million synoptic reports were analyzed for total cloud cover and the amounts of nine different cloud types, for the 26-yr period 1971–96. Monthly and seasonal averages were formed for day and night separately. Time series of total-cloud-cover anomalies for individual continents show a large decrease for South America, small decreases for Eurasia and Africa, and no trend for North America. The largest interannual variations (2.7%) are found for Australia, which is strongly influenced by ENSO. The zonal average trends of total cloud cover are positive in the Arctic winter and spring, 60°–80°N, but negative in all seasons at most other latitudes. The global average trend of total cloud cover over land is small, −0.7% decade−1, offsetting the small positive trend that had been found for the ocean, and resulting in no significant trend for the land–ocean average. Significant regional trends are found for many cloud types. The night trends agree with day trends for total cloud cover and for all cloud types except cumulus. Cirrus trends are generally negative over all continents. A previously reported decline in total cloud cover over China and its neighbors appears to be largely attributable to high and middle clouds. Global trends of the cloud types exhibit trade-offs, with convective cloud types increasing at the expense of stratiform clouds, in both the low and middle levels. Interannual variations over Europe, particularly of nimbostratus, are well correlated with the North Atlantic Oscillation; significant correlations are also found across northern Asia. Interannual variations in many parts of the Tropics are well correlated with an ENSO index. Little correlation was found with an index of smoke aerosol, in seven regions of seasonal biomass burning. In the middle latitudes of both hemispheres, seasonal anomalies of cloud cover are positively correlated with surface temperature in winter and negatively correlated in summer, as expected if the direction of causality is from clouds to temperature.

Description: Observations of the development of recent El Niño events suggest a pivotal role for the Madden–Julian oscillation (MJO). Previous attempts to uncover a systematic relationship between MJO activity and the El Niño–Southern Oscillation (ENSO), however, have yielded conflicting results. In this study the MJO–ENSO relationship is stratified by season, and the focus is on MJO activity in the equatorial western Pacific. The results demonstrate that MJO activity in late boreal spring leads El Niño in the subsequent autumn–winter for the period 1979–2005. Spring is the season when MJO activity is least asymmetric with respect to the equator and displays the most sensitivity to SST variations at the eastern edge of the warm pool. Enhanced MJO activity in the western Pacific in spring is associated with an eastward-expanded warm pool and low-frequency westerly surface zonal wind anomalies. These sustained westerly anomalies in the western Pacific are hypothesized to project favorably onto a developing El Niño in spring.

Description: A lagged maximum covariance analysis (MCA) of monthly anomaly data from the NCEP–NCAR reanalysis shows significant relations between the large-scale atmospheric circulation in two seasons and prior North Pacific sea surface temperature (SST) anomalies, independent from the teleconnections associated with the ENSO phenomenon. Regression analysis based on the SST anomaly centers of action confirms these findings. In late summer, a hemispheric atmospheric signal that is primarily equivalent barotropic, except over the western subtropical Pacific, is significantly correlated with an SST anomaly mode up to at least 5 months earlier. Although the relation is most significant in the upper troposphere, significant temperature anomalies are found in the lower troposphere over North America, the North Atlantic, Europe, and Asia. The SST anomaly is largest in the Kuroshio Extension region and along the subtropical frontal zone, resembling the main mode of North Pacific SST anomaly variability in late winter and spring, and it is itself driven by the atmosphere. The predictability of the atmospheric signal, as estimated from cross-validated correlation, is highest when SST leads by 4 months because the SST anomaly pattern is more dominant in the spring than in the summer. In late fall and early winter, a signal resembling the Pacific–North American (PNA) pattern is found to be correlated with a quadripolar SST anomaly during summer, up to 4 months earlier, with comparable statistical significance throughout the troposphere. The SST anomaly changes shape and propagates eastward, and by early winter it resembles the SST anomaly that is generated by the PNA pattern. It is argued that this results via heat flux forcing and meridional Ekman advection from an active coupling between the SST and the PNA pattern that takes place throughout the fall. Correspondingly, the predictability of the PNA-like signal is highest when SST leads by 2 months. In late summer, the maximum atmospheric perturbation at 250 mb reaches 35 m K−1 in the MCA and 20 m K−1 in the regressions. In early winter, the maximum atmospheric perturbation at 250 mb ranges between 70 m K−1 in the MCA and about 35 m K−1 in the regressions. This suggests that North Pacific SST anomalies have a substantial impact on the Northern Hemisphere climate. The back interaction of the atmospheric response onto the ocean is also discussed.

Description: This study examines the capabilities and limitations of the fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model (MM5) in predicting the precipitation and circulation features that accompanied the 2004 North American monsoon (NAM). When the model is reinitialized every 5 days to restrain the growth of modeling errors, its results for precipitation checked at subseasonal time scales (not for individual rainfall events) become comparable with ground- and satellite-based observations as well as with the NAM’s diagnostic characteristics. The modeled monthly precipitation illustrates the evolution patterns of monsoon rainfall, although it underestimates the rainfall amount and coverage area in comparison with observations. The modeled daily precipitation shows the transition from dry to wet episodes on the monsoon onset day over the Arizona–New Mexico region, and the multiday heavy rainfall (〉1 mm day−1) and dry periods after the onset. All these modeling predictions agree with observed variations. The model also accurately simulated the onset and ending dates of four major moisture surges over the Gulf of California during the 2004 monsoon season. The model reproduced the strong diurnal variability of the NAM precipitation, but did not predict the observed diurnal feature of the precipitation peak’s shift from the mountains to the coast during local afternoon to late night. In general, the model is able to reproduce the major, critical patterns and dynamic variations of the NAM rainfall at intraseasonal time scales, but still includes errors in precipitation quantity, pattern, and timing. The numerical study suggests that these errors are due largely to deficiencies in the model’s cumulus convective parameterization scheme, which is responsible for the model’s precipitation generation.