Cyclones causing wind storms in the Mediterranean: characteristics, trends and links to large-scale patterns

A climatology of cyclones with a focus on their relation to wind storm tracks in the Mediterranean region (MR) is presented. Trends in the frequency of cyclones and wind storms, as well as variations associated with the North Atlantic Oscillation (NAO), the East Atlantic/West Russian (EAWR) and the Scandinavian variability pattern (SCAND) are discussed. The study is based on the ERA40 reanalysis dataset. Wind storm tracks are identified by tracking clusters of adjacent grid boxes characterised by extremely high local wind speeds. The wind track is assigned to a cyclone track independently identified with an objective scheme. Areas with high wind activity – quantified by extreme wind tracks – are typically located south of the Golf of Genoa, south of Cyprus, southeast of Sicily and west of the Iberian Peninsula. About 69% of the wind storms are caused by cyclones located in the Mediterranean region, while the remaining 31% can be attributed to North Atlantic or Northern European cyclones. The North Atlantic Oscillation, the East Atlantic/West Russian pattern and the Scandinavian pattern all influence the amount and spatial distribution of wind inducing cyclones and wind events in the MR. The strongest signals exist for the NAO and the EAWR pattern, which are both associated with an increase in the number of organised strong wind events in the eastern MR during their positive phase. On the other hand, the storm numbers decrease over the western MR for the positive phase of the NAO and over the central MR during the positive phase of the EAWR pattern. The positive phase of the Scandinavian pattern is associated with a decrease in the number of winter wind storms over most of the MR. A third of the trends in the number of wind storms and wind producing cyclones during the winter season of the ERA40 period may be attributed to the variability of the North Atlantic Oscillation.

Examination of wind storms over Central Europe with respect to circulation weather types and NAO phases

The occurrence of wind storms in Central Europe is investigated with respect to large-scale atmospheric flow and local wind speeds in the investigation area. Two different methods of storm identification are applied for Central Europe as the target region: one based on characteristics of large-scale flow (circulation weather types, CWT) and the other on the occurrence of extreme wind speeds. The identified events are examined with respect to the NAO phases and CWTs under which they occur. Pressure patterns, wind speeds and cyclone tracks are investigated for storms assigned to different CWTs. Investigations are based on ERA40 reanalysis data. It is shown that about 80% of the storm days in Central Europe are connected with westerly flow and that Central European storm events primarily occur during a moderately positive NAO phase, while strongly positive NAO phases (6.4% of all days) account for more than 20% of the storms. A storm occurs over Central Europe during about 10% of the days with a strong positive NAO index. The most frequent pathway of cyclone systems associated with storms over Central Europe leads from the North Atlantic over the British Isles, North Sea and southern Scandinavia into the Baltic Sea. The mean intensity of the systems typically reaches its maximum near the British Isles. Differences between the characteristics for storms identified from the CWT identification procedure (gale days, based on MSLP fields) and those from extreme winds at Central European grid points are small, even though only 70% of the storm days agree. While most storms occur during westerly flow situations, specific characteristics of storms during the other CWTs are also considered. Copyright © 2009 Royal Meteorological Society

Estimation of wind storm impacts over Western Germany under future climate conditions using a statistical-dynamical downscaling approach

A statistical–dynamical regionalization approach is developed to assess possible changes in wind storm impacts. The method is applied to North Rhine-Westphalia (Western Germany) using the FOOT3DK mesoscale model for dynamical downscaling and ECHAM5/OM1 global circulation model climate projections. The method first classifies typical weather developments within the reanalysis period using K-means cluster algorithm. Most historical wind storms are associated with four weather developments (primary storm-clusters). Mesoscale simulations are performed for representative elements for all clusters to derive regional wind climatology. Additionally, 28 historical storms affecting Western Germany are simulated. Empirical functions are estimated to relate wind gust fields and insured losses. Transient ECHAM5/OM1 simulations show an enhanced frequency of primary storm-clusters and storms for 2060–2100 compared to 1960–2000. Accordingly, wind gusts increase over Western Germany, reaching locally +5% for 98th wind gust percentiles (A2-scenario). Consequently, storm losses are expected to increase substantially (+8% for A1B-scenario, +19% for A2-scenario). Regional patterns show larger changes over north-eastern parts of North Rhine-Westphalia than for western parts. For storms with return periods above 20 yr, loss expectations for Germany may increase by a factor of 2. These results document the method's functionality to assess future changes in loss potentials in regional terms.

The Empirical Forcing Function as a tool for the diagnosis of large-scale atmospheric anomalies

The time-mean quasi-geostrophic potential vorticity equation of the atmospheric flow on isobaric surfaces can explicitly include an atmospheric (internal) forcing term of the stationary-eddy flow. In fact, neglecting some non-linear terms in this equation, this forcing can be mathematically expressed as a single function, called Empirical Forcing Function (EFF), which is equal to the material derivative of the time-mean potential vorticity. Furthermore, the EFF can be decomposed as a sum of seven components, each one representing a forcing mechanism of different nature. These mechanisms include diabatic components associated with the radiative forcing, latent heat release and frictional dissipation, and components related to transient eddy transports of heat and momentum. All these factors quantify the role of the transient eddies in forcing the atmospheric circulation. In order to assess the relevance of the EFF in diagnosing large-scale anomalies in the atmospheric circulation, the relationship between the EFF and the occurrence of strong North Atlantic ridges over the Eastern North Atlantic is analyzed, which are often precursors of severe droughts over Western Iberia. For such events, the EFF pattern depicts a clear dipolar structure over the North Atlantic; cyclonic (anticyclonic) forcing of potential vorticity is found upstream (downstream) of the anomalously strong ridges. Results also show that the most significant components are related to the diabatic processes. Lastly, these results highlight the relevance of the EFF in diagnosing large-scale anomalies, also providing some insight into their interaction with different physical mechanisms.

Assessment of the Wind Gust Estimate Method in mesoscale modelling of storm events over West Germany

A physically based gust parameterisation is added to the atmospheric mesoscale model FOOT3DK to estimate wind gusts associated with storms over West Germany. The gust parameterisation follows the Wind Gust Estimate (WGE) method and its functionality is verified in this study. The method assumes that gusts occurring at the surface are induced by turbulent eddies in the planetary boundary layer, deflecting air parcels from higher levels down to the surface under suitable conditions. Model simulations are performed with horizontal resolutions of 20 km and 5 km. Ten historical storm events of different characteristics and intensities are chosen in order to include a wide range of typical storms affecting Central Europe. All simulated storms occurred between 1990 and 1998. The accuracy of the method is assessed objectively by validating the simulated wind gusts against data from 16 synoptic stations by means of “quality parameters”. Concerning these parameters, the temporal and spatial evolution of the simulated gusts is well reproduced. Simulated values for low altitude stations agree particularly well with the measured gusts. For orographically exposed locations, the gust speeds are partly underestimated. The absolute maximum gusts lie in most cases within the bounding interval given by the WGE method. Focussing on individual storms, the performance of the method is better for intense and large storms than for weaker ones. Particularly for weaker storms, the gusts are typically overestimated. The results for the sample of ten storms document that the method is generally applicable with the mesoscale model FOOT3DK for mid-latitude winter storms, even in areas with complex orography.

Extreme wind storms over Europe in present and future climate: a cluster analysis approach

Boreal winter wind storm situations over Central Europe are investigated by means of an objective cluster analysis. Surface data from the NCEP-Reanalysis and ECHAM4/OPYC3-climate change GHG simulation (IS92a) are considered. To achieve an optimum separation of clusters of extreme storm conditions, 55 clusters of weather patterns are differentiated. To reduce the computational effort, a PCA is initially performed, leading to a data reduction of about 98 %. The clustering itself was computed on 3-day periods constructed with the first six PCs using "k-means" clustering algorithm. The applied method enables an evaluation of the time evolution of the synoptic developments. The climate change signal is constructed by a projection of the GCM simulation on the EOFs attained from the NCEP-Reanalysis. Consequently, the same clusters are obtained and frequency distributions can be compared. For Central Europe, four primary storm clusters are identified. These clusters feature almost 72 % of the historical extreme storms events and add only to 5 % of the total relative frequency. Moreover, they show a statistically significant signature in the associated wind fields over Europe. An increased frequency of Central European storm clusters is detected with enhanced GHG conditions, associated with an enhancement of the pressure gradient over Central Europe. Consequently, more intense wind events over Central Europe are expected. The presented algorithm will be highly valuable for the analysis of huge data amounts as is required for e.g. multi-model ensemble analysis, particularly because of the enormous data reduction.

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

Southern annular mode dynamics in observations and models, Part I: The influence of climatological zonal wind biases in a comprehensive GCM

A common bias among global climate models (GCMs) is that they exhibit tropospheric southern annular mode (SAM) variability that is much too persistent in the Southern Hemisphere (SH) summertime. This is of concern for the ability to accurately predict future SH circulation changes, so it is important that it be understood and alleviated. In this two-part study, specifically targeted experiments with the Canadian Middle Atmosphere Model (CMAM) are used to improve understanding of the enhanced summertime SAM persistence. Given the ubiquity of this bias among comprehensive GCMs, it is likely that the results will be relevant for other climate models. Here, in Part I, the influence of climatological circulation biases on SAM variability is assessed, with a particular focus on two common biases that could enhance summertime SAM persistence: the too-late breakdown of the Antarctic stratospheric vortex and the equatorward bias in the SH tropospheric midlatitude jet. Four simulations are used to investigate the role of each of these biases in CMAM. Nudging and bias correcting procedures are used to systematically remove zonal-mean stratospheric variability and/or remove climatological zonal wind biases. The SAM time-scale bias is not alleviated by improving either the timing of the stratospheric vortex breakdown or the climatological jet structure. Even in the absence of stratospheric variability and with an improved climatological circulation, the model time scales are biased long. This points toward a bias in internal tropospheric dynamics that is not caused by the tropospheric jet structure bias. The underlying cause of this is examined in more detail in Part II of this study.

Projection of global wave climate change towards the end of the 21st century

Wind generated waves at the sea surface are of outstanding importance for both their practical relevance in many aspects, such as coastal erosion, protection, or safety of navigation, and for their scientific relevance in modifying fluxes at the air-sea interface. So far long-term changes in ocean wave climate have been studied mostly from a regional perspective with global dynamical studies emerging only recently. Here a global wave climate study is presented, in which a global wave model (WAM) is driven by atmospheric forcing from a global climate model (ECHAM5) for present day and potential future climate conditions represented by the IPCC (Intergovernmental Panel for Climate Change) A1B emission scenario. It is found that changes in mean and extreme wave climate towards the end of the twenty-first century are small to moderate, with the largest signals being a poleward shift in the annual mean and extreme significant wave heights in the mid-latitudes of both hemispheres, more pronounced in the Southern Hemisphere, and most likely associated with a corresponding shift in mid-latitude storm tracks. These changes are broadly consistent with results from the few studies available so far. The projected changes in the mean wave periods, associated with the changes in the wave climate in the mid to high latitudes, are also shown, revealing a moderate increase in the equatorial eastern side of the ocean basins. This study presents a step forward towards a larger ensemble of global wave climate projections required to better assess robustness and uncertainty of potential future wave climate change.

An assessment of a three-beam Doppler lidar wind profiling method for use in urban areas

Currently there are few observations of the urban wind field at heights other than rooftop level. Remote sensing instruments such as Doppler lidars provide wind speed data at many heights, which would be useful in determining wind loadings of tall buildings, and predicting local air quality. Studies comparing remote sensing with traditional anemometers carried out in flat, homogeneous terrain often use scan patterns which take several minutes. In an urban context the flow changes quickly in space and time, so faster scans are required to ensure little change in the flow over the scan period. We compare 3993 h of wind speed data collected using a three-beam Doppler lidar wind profiling method with data from a sonic anemometer (190 m). Both instruments are located in central London, UK; a highly built-up area. Based on wind profile measurements every 2 min, the uncertainty in the hourly mean wind speed due to the sampling frequency is 0.05–0.11 m s−1. The lidar tended to overestimate the wind speed by ≈0.5 m s−1 for wind speeds below 20 m s−1. Accuracy may be improved by increasing the scanning frequency of the lidar. This method is considered suitable for use in urban areas.

High-latitude climate sensitivity to ice-sheet forcing over the last 120kyr

Interpretation of ice-core records is currently limited by paucity of modelling at adequate temporal and spatial resolutions. Several key questions relate to mechanisms of polar amplification and inter-hemispheric coupling on glacial/interglacial timescales. Here, we present the first results from a large set of global ocean–atmosphere climate model ‘snap-shot’ simulations covering the last 120 000 years using the Hadley Centre climate model (HadCM3) at up to 1 kyr temporal resolution. Two sets of simulations were performed in order to examine the roles of orbit and greenhouse gases versus ice-sheet forcing of orbital-scale climate change. A series of idealised Heinrich events were also simulated, but no changes to aerosols or vegetation were prescribed. This paper focuses on high latitudes and inter-hemispheric linkages. The simulations reproduce polar temperature trends well compared to ice-core reconstructions, although the magnitude is underestimated. Polar amplification varies with obliquity, but this variability is dampened by including variations in land ice coverage, while the overall amplification factor increases. The relatively constant amplification of Antarctic temperatures (with ice-sheet forcing included) suggests it is possible to use Antarctic temperature reconstructions to estimate global changes (which are roughly half the magnitude). Atlantic Ocean overturning circulation varies considerably only with the introduction of Northern Hemisphere ice sheets, but only weakens in the North Atlantic in the deep glacial, when ocean–sea-ice feedbacks result in the movement of the region of deep convection to lower latitudes and with the introduction of freshwater to the surface North Atlantic in order to simulate Heinrich events.

Mid-Holocene monsoons: a multi-model analysis of the inter-hemispheric differences in the responses to orbital forcing and ocean feedbacks

The response of monsoon circulation in the northern and southern hemisphere to 6 ka orbital forcing has been examined in 17 atmospheric general circulation models and 11 coupled ocean–atmosphere general circulation models. The atmospheric response to increased summer insolation at 6 ka in the northern subtropics strengthens the northern-hemisphere summer monsoons and leads to increased monsoonal precipitation in western North America, northern Africa and China; ocean feedbacks amplify this response and lead to further increase in monsoon precipitation in these three regions. The atmospheric response to reduced summer insolation at 6 ka in the southern subtropics weakens the southern-hemisphere summer monsoons and leads to decreased monsoonal precipitation in northern South America, southern Africa and northern Australia; ocean feedbacks weaken this response so that the decrease in rainfall is smaller than might otherwise be expected. The role of the ocean in monsoonal circulation in other regions is more complex. There is no discernable impact of orbital forcing in the monsoon region of North America in the atmosphere-only simulations but a strong increase in precipitation in the ocean–atmosphere simulations. In contrast, there is a strong atmospheric response to orbital forcing over northern India but ocean feedback reduces the strength of the change in the monsoon although it still remains stronger than today. Although there are differences in magnitude and exact location of regional precipitation changes from model to model, the same basic mechanisms are involved in the oceanic modulation of the response to orbital forcing and this gives rise to a robust ensemble response for each of the monsoon systems. Comparison of simulated and reconstructed changes in regional climate suggest that the coupled ocean–atmosphere simulations produce more realistic changes in the northern-hemisphere monsoons than atmosphere-only simulations, though they underestimate the observed changes in precipitation in all regions. Evaluation of the southern-hemisphere monsoons is limited by lack of quantitative reconstructions, but suggest that model skill in simulating these monsoons is limited.

Observations of wind speed profiles over Greater London, UK, using a Doppler lidar

To calculate the potential wind loading on a tall building in an urban area, an accurate representation of the wind speed profile is required. However, due to a lack of observations, wind engineers typically estimate the characteristics of the urban boundary layer by translating the measurements from a nearby reference rural site. This study presents wind speed profile data obtained from a Doppler lidar in central London, UK, during an 8 month observation period. Used in conjunction with wind speed data measured at a nearby airport, the data have been used to assess the accuracy of the predictions made by the wind engineering tools currently available. When applied to multiple changes in surface roughness identified from morphological parameters, the non-equilibrium wind speed profile model developed by Deaves (1981) provides a good representation of the urban wind speed profile. For heights below 500 m, the predicted wind speed remains within the 95% confidence interval of the measured data. However, when the surface roughness is estimated using land use as a proxy, the model tends to overestimate the wind speed, particularly for very high wind speed periods. These results highlight the importance of a detailed assessment of the nature of the surface when estimating the wind speed above an urban surface.

Observations of diurnal and spatial variability of radiative forcing by equatorial deep convective clouds

From geostationary satellite observations of equatorial Africa and the equatorial east Atlantic during May and June 2000 we explore the radiative forcing by deep convective cloud systems in these regions. Deep convective clouds (DCCs) are associated with a mean radiative forcing relative to non–deep convective areas of −39 W m−2 over the Atlantic Ocean and of +13 W m−2 over equatorial Africa (±10 W m−2 in both cases). We show that over land the timing of the daily cycle of convection relative to the daily cycle in solar illumination and surface temperature significantly affects the mean radiative forcing by DCCs. Displacement of the daily cycle of DCC coverage by 2 hours changes their overall radiative effect by ∼10 W m−2, with implications for the simulation of the radiative balance in this region. The timing of the minimum DCC cover over land, close to noon local time, means that the mean radiative forcing is nearly maximized.

20th century intraseasonal Asian monsoon dynamics viewed from Isomap

The Asian summer monsoon is a high dimensional and highly nonlinear phenomenon involving considerable moisture transport towards land from the ocean, and is critical for the whole region. We have used daily ECMWF reanalysis (ERA-40) sea-level pressure (SLP) anomalies to the seasonal cycle, over the region 50-145°E, 20°S-35°N to study the nonlinearity of the Asian monsoon using Isomap. We have focused on the two-dimensional embedding of the SLP anomalies for ease of interpretation. Unlike the unimodality obtained from tests performed in empirical orthogonal function space, the probability density function, within the two-dimensional Isomap space, turns out to be bimodal. But a clustering procedure applied to the SLP data reveals support for three clusters, which are identified using a three-component bivariate Gaussian mixture model. The modes are found to appear similar to active and break phases of the monsoon over South Asia in addition to a third phase, which shows active conditions over the Western North Pacific. Using the low-level wind field anomalies the active phase over South Asia is found to be characterised by a strengthening and an eastward extension of the Somali jet whereas during the break phase the Somali jet is weakened near southern India, while the monsoon trough in northern India also weakens. Interpretation is aided using the APHRODITE gridded land precipitation product for monsoon Asia. The effect of large-scale seasonal mean monsoon and lower boundary forcing, in the form of ENSO, is also investigated and discussed. The outcome here is that ENSO is shown to perturb the intraseasonal regimes, in agreement with conceptual ideas.