Convectively coupled Kelvin and easterly waves in a regional climate simulation of the tropics

This study evaluates the performance of a regional climate model in simulating two types of synoptic tropical weather disturbances: convectively-coupled Kelvin and easterly waves. Interest in these two wave modes stems from their potential predictability out to several weeks in advance, as well as a strong observed linkage between easterly waves and tropical cyclogenesis. The model is a recent version of the weather research and forecast (WRF) system with 36-km horizontal grid spacing and convection parameterized using a scheme that accounts for key convective triggering and inhibition processes. The domain spans the entire tropical belt between 45°S and 45°N with periodic boundary conditions in the east–west direction, and conditions at the meridional/lower boundaries specified based on observations. The simulation covers 6 years from 2000 to 2005, which is long enough to establish a statistical depiction of the waves through space-time spectral filtering of rainfall data, together with simple lagged-linear regression. Results show that both the horizontal phase speeds and three-dimensional structures of the waves are qualitatively well captured by the model in comparison to observations. However, significant biases in wave activity are seen, with generally overactive easterly waves and underactive Kelvin waves. Evidence is presented to suggest that these biases in wave activity (which are also correlated with biases in time–mean rainfall, as well as biases in the model’s tropical cyclone climatology) stem in part from convection in the model coupling too strongly to rotational circulation anomalies. Nevertheless, the model is seen to do a reasonable job at capturing the genesis of tropical cyclones from easterly waves, with evidence for both wave accumulation and critical layer processes being importantly involved.     

Finite-difference solution of the transport equation: First results

The ray-theoretical transport equation for inhomogeneous isotropic media (2D-SH case) is solved by the method of finite differences on a rectangular grid, both for an incident plane wave (explicit scheme) and a line source (implicit scheme). Results for homogeneous models and for heterogeneous models with structural discontinuities are discussed. First-arrival travel times calculated by various techniques serve as input for the solution of the transport equation and the computation of amplitudes of first arrivals. To obtain correct amplitudes the travel times must be highly accurate and the discontinuities must be smoothed out; the reason is that the spatial second derivatives of the travel time field enter the transport equation. In the simple cases studied, finite differences provide a fast and efficient tool for the computation of first-arrival amplitudes.     

Observation of snowfall with a low-power FM-CW K-band radar (Micro Rain Radar)

Quantifying snowfall intensity especially under arctic conditions is a challenge because wind and snow drift deteriorate estimates obtained from both ground-based gauges and disdrometers. Ground-based remote sensing with active instruments might be a solution because they can measure well above drifting snow and do not suffer from flow distortions by the instrument. Clear disadvantages are, however, the dependency of e.g. radar returns on snow habit which might lead to similar large uncertainties. Moreover, high sensitivity radars are still far too costly to operate in a network and under harsh conditions. In this paper we compare returns from a low-cost, low-power vertically pointing FM-CW radar (Micro Rain Radar, MRR) operating at 24.1?GHz with returns from a 35.5?GHz cloud radar (MIRA36) for dry snowfall during a 6-month observation period at an Alpine station (Environmental Research Station Schneefernerhaus, UFS) at 2,650?m height above sea level. The goal was to quantify the potential and limitations of the MRR in relation to what is achievable by a cloud radar. The operational MRR procedures to derive standard radar variables like effective reflectivity factor (Z _e) or the mean Doppler velocity (W) had to be modified for snowfall since the MRR was originally designed for rain observations. Since the radar returns from snowfall are weaker than from comparable rainfall, the behavior of the MRR close to its detection threshold has been analyzed and a method is proposed to quantify the noise level of the MRR based on clear sky observations. By converting the resulting MRR-Z _e into 35.5?GHz equivalent Z _e values, a remaining difference below 1?dBz with slightly higher values close to the noise threshold could be obtained. Due to the much higher sensitivity of MIRA36, the transition of the MRR from the true signal to noise can be observed, which agrees well with the independent clear sky noise estimate. The mean Doppler velocity differences between both radars are below 0.3?ms^?1. The distribution of Z _e values from MIRA36 are finally used to estimate the uncertainty of retrieved snowfall and snow accumulation with the MRR. At UFS low snowfall rates missed by the MRR are negligible when comparing snow accumulation, which were mainly caused by intensities between 0.1 and 0.8 mm?h^?1. The MRR overestimates the total snow accumulation by about 7%. This error is much smaller than the error caused by uncertain Z _e?Csnowfall rate relations, which would affect the MIRA36 estimated to a similar degree.     

Climate drift in an ocean model coupled to a simple,perfectly matched atmosphere

A very simple, diffusive energy balance atmosphere is coupled to the GFDL ocean circulation model. This provides a useful tool for analyzing climate drift in the ocean model after coupling, and may be used to assess various schemes for minimizing such drift. In the experiment reported here, the atmosphere is constructed in such a way that it provides the ocean model at the moment of coupling with the same fluxes as during spinup. The experiment is therefore equivalent to coupling a perfectly flux-corrected atmosphere model, and is used to investigate the response of the ocean model under these conditions. In spite of the steady, passive, flux-corrected atmosphere, the ocean model drifts to a new equilibrium state after coupling. The transition takes about 2000 years; the new state is characterized by different sites of deep convection and resulting changes in high-latitude SST and global deep temperatures. The mechanism for the transition is an instability of the oceanic convection patterns under the new feedback, felt after coupling. A similar state transition of the ocean model may be triggered by the coupling shock in fully coupled GCMs. If this is so, the transition would contaminate the results of climate scenario experiments, and it would explain part of the residual drift observed in coupled models in spite of the use of flux corrections.     

Tropical versus high latitude freshwater influence on the Atlantic circulation

We investigate the model sensitivity of the Atlantic meridional overturning circulation (AMOC) to anomalous freshwater flux in the tropical and northern Atlantic. Forcing in both locations leads to the same qualitative response: a positive freshwater anomaly induces a weakening of the AMOC and a negative freshwater anomaly strengthens the AMOC. Strong differences arise in the temporal characteristics and amplitude of the response. The advection of the tropical anomaly up to the deep water formation area leads to a time delayed response compared to a northern forcing. Thus, in its transient response, the AMOC is less sensitive to a constant anomalous freshwater flux in the tropics than in the north. This difference decreases with time and practically vanishes in equilibrium with constant freshwater forcing. The equilibrium response of the AMOC shows a non-linear dependence on freshwater forcing in both locations, with a stronger sensitivity to positive freshwater forcing. As a consequence, competitive forcing in both regions is balanced when the negative forcing is about 1.5 times larger than the positive forcing. The relaxation time of the AMOC after termination of a freshwater perturbation depends significantly on the AMOC strength itself. A strong overturning exhibits a faster relaxation to its unperturbed state. By means of a set of complementary experiments (pulse-perturbations, constant and stochastic forcing) we quantify these effects and discuss the corresponding time scales and physical processes.     

Circumglobal wave train and the summer monsoon over northwestern India and Pakistan: the explicit role of the surface heat low

This study examines the influence of the mid-latitude circulation on the surface heat low (HL) and associated monsoon rainfall over northwestern India and Pakistan using the ERA40 data and high resolution (T106L31) climate model ECHAM5 simulation. Special emphasis is given to the surface HL which forms over Pakistan and adjoining areas of India, Iran and Afghanistan during the summer season. A heat low index (HLI) is defined to depict the surface HL. The HLI displays significant correlations with the upper level mid-latitude circulation over western central Asia and low level monsoon circulation over Arabian Sea and acts as a bridge connecting the mid-latitude wave train to the Indian summer monsoon. A time-lagged singular value decomposition analysis reveals that the eastward propagation of the mid-latitude circumglobal wave train (CGT) influences the surface pressure anomalies over the Indian domain. The largest low (negative) pressure anomalies over the western parts of the HL region (i.e., Iran and Afghanistan) occur in conjunction with the upper level anomalous high that develops over western-central Asia during the positive phase of the CGT. The composite analysis also reveals a significant increase in the low pressure anomalies over Iran and Afghanistan during the positive phase of CGT. The westward increasing low pressure anomalies with its north?Csouth orientation provokes enormous north?Csouth pressure gradient (lower pressure over land than over sea). This in turn enables the moist southerly flow from the Arabian Sea to penetrate farther northward over northwestern India and Pakistan. A monsoon trough like conditions develops over northwestern India and Pakistan where the moist southwesterly flow from the Arabian Sea and the Persian Gulf converge. The convergence in association with the orographic uplifting expedites convection and associated precipitation over northwestern India and Pakistan. The high resolution climate model ECHAM5 simulation also underlines the proposed findings and mechanism.     

Applying probabilistic projections of climate change with impact models: a case study for sub-arctic palsa mires in Fennoscandia

A comparison of two approaches for determining probabilistic climate change impacts is presented. In the first approach, ensemble climate projections are applied directly as inputs to an impact model and the risk of impact is computed from the resulting ensemble of outcomes. As this can involve large numbers of projections, the approach may prove to be impractical when applied to complex impact models with demanding input requirements. The second approach is to construct an impact response surface based on a sensitivity analysis of the impact model with respect to changes in key climatic variables, and then to superimpose probabilistic projections of future climate onto the response surface to assess the risk of impact. To illustrate this comparison, an impact model describing the spatial distribution of palsas in Fennoscandia was applied to estimate the risk of palsa disappearance. Palsas are northern mire complexes with permanently frozen peat hummocks, located at the outer limit of the permafrost zone and susceptible to rapid decline due to regional warming. Probabilities of climate changes were derived from an ensemble of coupled atmosphere–ocean general circulation model (AOGCM) projections using a re-sampling method. Results indicated that the response surface approach, though introducing additional uncertainty, gave risk estimates of area decline for palsa suitability that were comparable to those obtained using multiple simulations with the original palsa model. It was estimated as very likely (>90% probability) that a decline of area suitable for palsas to less than half of the baseline distribution will occur by the 2030s and likely (>66%) that all suitable areas will disappear by the end of the twenty-first century under scenarios of medium (A1B) and moderately high (A2) emissions. For a low emissions (B1) scenario, it was more likely than not (>50%) that conditions over a small fraction of the current palsa distribution would remain suitable until the end of the twenty-first century.     

Stratospheric aerosol sulfuric acid: First direct in situ measurements using a novel balloon-based mass spectrometer apparatus

The first direct in situ measurements of the sulfuric acid contained in stratospheric aerosol particles were made using a novel balloon-based Ion Molecule Reaction Mass Spectrometer instrument (IMRMS) equipped with an aerosol vaporizer. The IMRMS method employed offers direct aerosol sulfuric acid measurements with high spatial resolution. The balloon flight took place on 23 October 1995 in middle latitudes (44°N) and reached a maximum altitude of 24 km. Measured molecular number densities of aerosol sulfuric acid decreased with increasing altitude from about 2.9 × 109 cm-3 at 15km altitude to about 2.4 × 108 cm-3 at 21 km. Corresponding mass mixing ratios are 2.5 and 0.6 ppbm, respectively. Calculated sulfuric acid mixing ratios from aerosol volumes inferred from aerosol size distribution measurements on the same balloon agree well with the IMRMS data using standard assumption aerosol composition.     

Documentary evidence of climate variability during cold seasons in Lesotho, southern Africa, 1833–1900

This study presents the first 19th century cold season climate chronology for the Kingdom of Lesotho in southern Africa. The chronology is constructed using a variety of documentary sources including letters, diaries, reports, monographs and newspaper articles obtained from southern African and British archives. Information relating to cold season weather phenomena during the austral autumn, winter and early spring months were recorded verbatim. Each of the cold seasons from 1833 to 1900 was then classified as “very severe”, “severe” or “normal/mild”, with a confidence rating ranging from low (1) to high (3) awarded against each annual classification. The accuracy of the document-derived chronology was verified against temperature data for Maseru for the period 1893–1900. Excellent correspondence of the document-derived chronology with the Maseru instrumental data and also with other global proxy temperature records for the 19th century is achieved. The results indicate 12 (18% of the total) very severe, 16 (23%) severe and 40 (59%) normal/mild cold seasons between 1833 and 1900. The overall trend is for more severe and snow-rich cold seasons during the early part of the study period (1833–1854) compared with the latter half of the 19th century (with the exception of the 1880s). A reduction in the duration of the frost season by over 20 days during the 19th century is also tentatively identified. Several severe to very severe cold seasons in Lesotho follow after major tropical and SH volcanic eruptions; such years are usually characterized by early frosts, and frequent and heavy snowfalls. The blocking of solar radiation and the enhanced northward displacement of polar fronts that are directly or indirectly associated with volcanic events, may account for many of the most severe Lesotho winters during the 19th century.