Witnessing North Atlantic westerlies variability from ships’ logbooks (1685–2008)

A monthly index based on the persistence of the westerly winds over the English Chanel is constructed for 1685–2008 using daily data from ships’ logbooks and comprehensive marine meteorological datasets. The so-called Westerly Index (WI) provides the longest instrumental record of atmospheric circulation currently available. Anomalous WI values are associated with spatially coherent climatic signals in temperature and precipitation over large areas of Europe, which are stronger for precipitation than for temperature and in winter and summer than in transitional seasons. Overall, the WI series accord with the known European climatic history, and reveal that the frequency of the westerlies in the eastern Atlantic during the twentieth century and the Late Maunder Minimum was not exceptional in the context of the last three centuries. It is shown that the WI provides additional and complementary information to the North Atlantic Oscillation (NAO) indices. The analysis of WI series during the industrial era indicates an overall good agreement with the winter and high-summer NAO, with the exception of several multidecadal periods of weakened correlation. These decoupled periods between the frequency and the intensity of the zonal flow are interpreted on the basis of several sources of non-stationarity affecting the centres of the variability of the North Atlantic and their teleconnections. Comparisons with NAO reconstructions and long instrumental indices extending back to the seventeenth century suggest that similar situations have occurred in the past, which call for caution when reconstructing the past atmospheric circulation from climatic proxies. The robustness and extension of its climatic signal, the length of the series and its instrumental nature make the WI an excellent benchmark for proxy calibration in Europe and Greenland.     

Scalar Fluxes Near a Tall Building in an Aligned Array of Rectangular Buildings

Scalar dispersion from ground-level sources in arrays of buildings is investigated using wind-tunnel measurements and large-eddy simulation (LES). An array of uniform-height buildings of equal dimensions and an array with an additional single tall building (wind tunnel) or a periodically repeated tall building (LES) are considered. The buildings in the array are aligned and form long streets. The sensitivity of the dispersion pattern to small changes in wind direction is demonstrated. Vertical scalar fluxes are decomposed into the advective and turbulent parts and the influences of wind direction and of the presence of the tall building on the scalar flux components are evaluated. In the uniform-height array turbulent scalar fluxes are dominant, whereas the tall building produces an increase of the magnitude of advective scalar fluxes that yields the largest component. The presence of the tall building causes either an increase or a decrease to the total vertical scalar flux depending on the position of the source with respect to the tall building. The results of the simulations can be used to develop parametrizations for street-canyon dispersion models and enhance their capabilities in areas with tall buildings.     

North Atlantic atmospheric circulation and surface wind in the Northeast of the Iberian Peninsula: uncertainty and long term downscaled variability

The variability and predictability of the surface wind field at the regional scale is explored over a complex terrain region in the northeastern Iberian Peninsula by means of a downscaling technique based on Canonical Correlation Analysis. More than a decade of observations (1992–2005) allows for calibrating and validating a statistical method that elicits the main associations between the large scale atmospheric circulation over the North Atlantic and Mediterranean areas and the regional wind field. In an initial step the downscaling model is designed by selecting parameter values from practise. To a large extent, the variability of the wind at monthly timescales is found to be governed by the large scale circulation modulated by the particular orographic features of the area. The sensitivity of the downscaling methodology to the selection of the model parameter values is explored, in a second step, by performing a systematic sampling of the parameters space, avoiding a heuristic selection. This provides a metric for the uncertainty associated with the various possible model configurations. The uncertainties associated with the model configuration are considerably dependent on the spatial variability of the wind. While the sampling of the parameters space in the model set up moderately impact estimations during the calibration period, the regional wind variability is very sensitive to the parameters selection at longer timescales. This fact illustrates that downscaling exercises based on a single configuration of parameters should be interpreted with extreme caution. The downscaling model is used to extend the estimations several centuries to the past using long datasets of sea level pressure, thereby illustrating the large temporal variability of the regional wind field from interannual to multicentennial timescales. The analysis does not evidence long term trends throughout the twentieth century, however anomalous episodes of high/low wind speeds are identified.     

Future regional projections of extreme temperatures in Europe: a nonstationary seasonal approach

This paper analyzes changes of maximum temperatures in Europe, which are evaluated using two state-of-the-art regional climate models from the EU ENSEMBLES project. Extremes are expressed in terms of return values using a time-dependent generalized extreme value (GEV) model fitted to monthly maxima. Unlike the standard GEV method, this approach allows analyzing return periods at different time scales (monthly, seasonal, annual, etc). The study focuses on the end of the 20th century (1961?C2000), used as a calibration/validation period, and assesses the changes projected for the period 2061?C2100 considering the A1B emission scenario. The performance of the regional models is evaluated for each season of the calibration period against the high-resolution gridded E-OBS dataset, showing a similar South-North gradient with larger values over the Mediterranean basin. The inter-RCM changes in the bias pattern with respect to the E-OBS are larger than the bias resulting from a change in the boundary conditions from ERA-40 to ECHAM5 20c3m. The maximum temperature response to increased green house gases, as projected by the A1B scenario, is consistent for both RCMs. Under that scenario, results indicate that the increments for extremes (e.g. 40-year return values) will be two or three times higher than those for the mean seasonal temperatures, particularly during Spring and Summer in Southern Europe.     

The End-to-End Attribution Problem: From Emissions to Impacts

When a damaging extreme meteorological event occurs, the question often arises as to whether that event was caused by anthropogenic greenhouse gas emissions. The question is more than academic, since people affected by the event will be interested in recurring damages if they find that someone is at fault. However, since this extreme event could have occurred by chance in an unperturbed climate, we are currently unable to properly respond to this question. A solution lies in recognising the similarity with the cause-effect issue in the epidemiological field. The approach there is to consider the changes in the risk of the event occurring as attributable, as against the occurrence of the event itself. Inherent in this approach is a recognition that knowledge of the change in risk as well as the amplitude of the forcing itself are uncertain. Consequently, the fraction of the risk attributable to the external forcing is a probabilistic quantity. Here we develop and demonstrate this methodology in the context of the climate change problem.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

The 2001 Mesoscale Convective Systems over Iberia and the Balearic Islands

Summary This paper characterizes Mesoscale Convective Systems (MCSs) during 2001 over Iberia and the Balearic Islands and their meteorological settings. Enhanced infrared Meteosat imagery has been used to detect their occurrence over the Western Mediterranean region between June and December 2001 according to satellite-defined criteria based on the MCS physical characteristics. Twelve MCSs have been identified. The results show that the occurrence of 2001 MCSs is limited to the August–October period, with September being the most active period. They tend to develop during the late afternoon or early night, with preferred eastern Iberian coast locations and eastward migrations. A cloud shield area of 50.000 km² is rarely exceeded. When our results are compared with previous studies, it is possible to assert that though 2001 MCS activity was moderate, the convective season was substantially less prolonged than usual, with shorter MCS life cycles and higher average speeds. The average MCS precipitation rate was 3.3 mm·h⁻¹ but a wide range of values varying from scarce precipitation to intense events of 130 mm·24 h⁻¹ (6 September) were collected. The results suggest that, during 2001, MCS rainfall was the principal source of precipitation in the Mediterranean region during the convective season, but its impact varied according to the location. Synoptic analysis based on NCEP/NCAR reanalysis show that several common precursors could be identified over the Western Mediterranean Sea when the 2001 MCSs occurred: a low-level tongue of moist air and precipitable water (PW) exceeding 25 mm through the southern portion of the Western Mediterranean area, low-level zonal warm advection over 2 °C·24 h⁻¹ towards eastern Iberia, a modest 1000–850 hPa equivalent potential temperature (θ_e) difference over 20 °C located close to the eastern Iberian coast, a mid level trough (sometimes a cut-off low) over Northern Africa or Southern Spain and high levels geostrophic vorticity advection exceeding 12·10⁻¹⁰ s⁻² over eastern Iberia and Northern Africa. Finally, the results suggest that synoptic, orographic and a warm-air advection were the most relevant forcing mechanisms during 2001.     

Impact of anthropogenic CO 2 on the next glacial cycle

The model of Paillard and Parrenin (Earth Planet Sci Lett 227(3–4):263–271, 2004) has been recently optimized for the last eight glacial cycles, leading to two different relaxation models with model-data correlations between 0.8 and 0.9 (García-Olivares and Herrero (Clim Dyn 1–25, 2012b)). These two models are here used to predict the effect of an anthropogenic CO ₂ pulse on the evolution of atmospheric CO ₂, global ice volume and Antarctic ice cover during the next 300 kyr. The initial atmospheric CO ₂ condition is obtained after a critical data analysis that sets 1300 Gt as the most realistic carbon Ultimate Recoverable Resources (URR), with the help of a global compartmental model to determine the carbon transfer function to the atmosphere. The next 20 kyr will have an abnormally high greenhouse effect which, according to the CO ₂ values, will lengthen the present interglacial by some 25 to 33 kyr. This is because the perturbation of the current interglacial will lead to a delay in the future advance of the ice sheet on the Antarctic shelf, causing that the relative maximum of boreal insolation found 65 kyr after present (AP) will not affect the developing glaciation. Instead, it will be the following insolation peak, about 110 kyr AP, which will find an appropriate climatic state to trigger the next deglaciation.     

Vertical behaviour and meteorological properties of air masses in the southwest of the Iberian Peninsula (1997–2007)

Air masses are characterized by physical (temperature, humidity) and chemical (transported gases and aerosols) properties, being associated their arrival to different meteorological scenarios. The knowledge of the air masses over a region is fundamental as complementary information in several atmospheric studies, being the calculation of back-trajectory the most widely used tool whenever air masses are analyzed. A study of air masses has been carried out in southwestern Iberian Peninsula using 5-day kinematic back trajectories computed by the HYSPLIT model at three heights (500, 1,500 and 3,000 m) from 1997 to 2007. The main aims have been to characterize their vertical behaviour and their thermal and humidity properties. Thirteen trajectory clusters have been defined, showing the northerly and westerly clusters a high coupling degree. Seasonal daily variation of potential temperature and specific humidity has been analyzed, obtaining higher differences among clusters in the cold season.     

A densification mechanism to model the mechanical effect of methane hydrates in sandy sediments

Recent pore-scale observations and geomechanical investigations suggest the lack of true cohesion in methane hydrate-bearing sediments (MHBSs) and propose that their mechanical behavior is governed by kinematic constrictions at pore-scale. This paper presents a constitutive model for MHBS, which does not rely on physical bonding between hydrate crystals and sediment grains but on the densification effect that pore invasion with hydrate has on the sediment mechanical properties. The Hydrate-CASM extends the critical state model Clay and Sand Model (CASM) by implementing the subloading surface model and introducing the densification mechanism. The model suggests that the decrease of the sediment available void volume during hydrate formation stiffens its structure and has a similar mechanical effect as the increase of sediment density. In particular, the model attributes stress-strain changes observed in MHBS to the variations in sediment available void volume with hydrate saturation and its consequent effect on isotropic yield stress and swelling line slope. The model performance is examined against published experimental data from drained triaxial tests performed at different confining stress and with distinct hydrate saturation and morphology. Overall, the simulations capture the influence of hydrate saturation in both the magnitude and trend of the stiffness, shear strength, and volumetric response of synthetic MHBS. The results are validated against those obtained from previous mechanical models for MHBS that examine the same experimental data. The Hydrate-CASM performs similarly to previous models, but its formulation only requires one hydrate-related empirical parameter to express changes in the sediment elastic stiffness with hydrate saturation.