Anthropogenic climate change over the Mediterranean region simulated by a global variable resolution model

Two 30-year simulations corresponding to 1960-1989 and 2070-2099 have been performed with a variable resolution atmospheric model. The model has a maximum horizontal resolution of 0.5° over the Mediterranean Sea. Simulations are driven by IPCC-B2 scenario radiative forcing. Sea surface temperatures (SSTs) are prescribed from monthly observations for the present climate simulation, and from a blend of observations and coupled simulations for the scenario. Another pair of forced atmospheric simulations has been run under these forcings with the same uniform low resolution as the coupled model. Comparisons with observations show that the variable resolution model realistically reproduces the main climate characteristics of the Mediterranean region. At a global scale, changes in latitudinal temperature profiles are similar for the forced and coupled models, justifying the time-slice approach. The 2 m temperature and precipitation responses predict a warming and drying of the Mediterranean region. A comparison with the coupled simulation and forced low-resolution simulation shows that this pattern is robust. The decrease in mean precipitation is associated with a significant decrease in soil wetness, and could involve considerable impact on water resources around the Mediterranean basin.     

A Simple And General Subgrid Model Suitable Both For Surface Layer And Free-Stream Turbulence

A new and general approach is presented to allow standard subgrid schemes to besuitable both for surface layer and free-stream turbulence. Simple modificationsto subgrid schemes are proposed and derived for any vertical stabilityconditions. They are simple to implement in models and are suitable for morecomplicated simulations such as large-eddy simulation with inhomogeneoussurface conditions or complex topography. They are also applicable to mesoscaleand large-scale models. These modifications are physically justified by recentmeasurements of spectra close to the ground. The spectral analysis presentedshows how the energy deficit of blocked turbulence for a given dissipation(`anomalous dissipation') dramatically affects the coefficients to be used insubgrid schemes. As shown for neutral and convective cases with wind shear,these changes allow us to substantially improve the prediction of profiles for themean quantities in the surface layer. Agreement with similarity laws in the unstablecase is found up to about 0.2z_i, for simulated shear, stabilityprofiles and dissipation rates of turbulent kinetic energy.     

Mean flow-transient perturbation interaction in the Southern Hemisphere: a simulation using a variable-resolution GCM

The ability of an atmospheric general circulation model to reproduce fundamental features of the wintertime extratropical Southern Hemisphere (SH) circulation is evaluated with emphasis on the daily variability of the SH mean flow and the mean flow-transient perturbations interaction. Two 10-year simulations using a new version of the LMDZ GCM with a stretched grid scheme centered at 45 °S and forced by climatological SST are performed: a high (144Ꮡ) and low (64Ꭹ) horizontal resolution runs. The performance of both simulations was determined by comparing several simulated fields (zonal wind, temperature, kinetic energy, transient eddy momentum and heat fluxes, Eliassen-Palm fluxes, Eady growth rate and baroclinic conversion term) against the European Centre for Medium Range Weather Forecast reanalyses (ERA). High and low-resolution simulations are similar in many respects; in particular, both experiments reproduce the main patterns of the southern extratropical large-scale circulation satisfactorily. Increasing resolution does not improve universally some spurious aspects of the low resolution simulation (e.g. the cold bias in the high polar troposphere, the debilitated subtropical jet, the low baroclinic conversion rate). Those aspects present little sensitivity to the model resolution. The interaction between transient eddies and zonal mean flow are examined. The low-resolution experiment is able to qualitatively represent the acceleration/deceleration of the mean flow by transient perturbations, south/north of 30 °S with an accuracy similar to that of the high-resolution experiment. Although both experiments represent the baroclinic structure of the mean flow satisfactorily, the model underestimates some transient properties due to the underestimation of the baroclinic conversion term in middle latitudes. Such misrepresentation does not improve with increasing resolution and is related to the relatively weak meridional temperature gradient and the inadequate geographical distribution of the eddy heat fluxes. In particular, the eddy kinetic energy is always underestimated. Eddy kinetic energy does not improve convincingly with increasing resolution, suggesting that the adequate representation of the storm tracks is highly influenced by the physical parametrizations.     

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.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

Synoptic and mesoscale study of a severe convective outbreak with the nonhydrostatic Global Environmental Multiscale (GEM) model

Summary ?A nonhydrostatic 4-km version of the Global Environmental Multiscale (GEM) model, with detailed microphysics included, was used to forecast the initiation, development, and structure of a tornado-producing supercell storm that occurred near Pine Lake (Alberta, Canada) on 15 July 2000. Examination of observations and comparison with conceptual models indicate that this storm is a good example of supercell storms that regularly produce summertime severe weather over Alberta. It was found that the high-resolution model was able to reproduce the early initiation of convective activity along the Rocky Mountains foothills, as well as the rapid northeastward propagation towards the Pine Lake area and the subsequent intensification into a supercell storm. The general structures of the forecasted convective system correspond well with conceptual representations of such events. Overall, this high-resolution forecast of the Pine Lake supercell storm was a significant improvement over the current operational version of the GEM model (24 km), which was not able to intensify the foothills’ convection into a supercell storm. Finally, it was found that the nonhydrostatic version of the model produces better trajectory and propagation speed of the convective system, as compared with the hydrostatic one. Received March 20, 2001; revised August 24, 2001     

Projections of climate change impacts on central America tropical rainforest

Tropical rainforest plays an important role in the global carbon cycle, accounting for a large part of global net primary productivity and contributing to CO₂ sequestration. The objective of this work is to simulate potential changes in the rainforest biome in Central America subject to anthropogenic climate change under two emissions scenarios, RCP4.5 and RCP8.5. The use of a dynamic vegetation model and climate change scenarios is an approach to investigate, assess or anticipate how biomes respond to climate change. In this work, the Inland dynamic vegetation model was driven by the Eta regional climate model simulations. These simulations accept boundary conditions from HadGEM2-ES runs in the two emissions scenarios. The possible consequences of regional climate change on vegetation properties, such as biomass, net primary production and changes in forest extent and distribution, were investigated. The Inland model projections show reductions in tropical forest cover in both scenarios. The reduction of tropical forest cover is greater in RCP8.5. The Inland model projects biomass increases where tropical forest remains due to the CO₂ fertilization effect. The future distribution of predominant vegetation shows that some areas of tropical rainforest in Central America are replaced by savannah and grassland in RCP4.5. Inland projections under both RCP4.5 and RCP8.5 show a net primary productivity reduction trend due to significant tropical forest reduction, temperature increase, precipitation reduction and dry spell increments, despite the biomass increases in some areas of Costa Rica and Panama. This study may provide guidance to adaptation studies of climate change impacts on the tropical rainforests in Central America.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

Strong ground motion record selection for the reliable prediction of the mean seismic collapse capacity of a structure group

How to select a limited number of strong ground motion records (SGMRs) is an important challenge for the seismic collapse capacity assessment of structures. The collapse capacity is considered as the ground motion intensity measure corresponding to the drift‐related dynamic instability in the structural system. The goal of this paper is to select, from a general set of SGMRs, a small number of subsets such that each can be used for the reliable prediction of the mean collapse capacity of a particular group of structures, i.e. of single degree‐of‐freedom systems with a typical behaviour range. In order to achieve this goal, multivariate statistical analysis is first applied, to determine what degree of similarity exists between each selected small subset and the general set of SGMRs. Principal Component analysis is applied to identify the best way to group structures, resulting in a minimum number of SGMRs in a proposed subset. The structures were classified into six groups, and for each group a subset of eight SGMRs has been proposed. The methodology has been validated by analysing a first‐mode‐dominated three‐storey‐reinforced concrete structure by means of the proposed subsets, as well as the general set of SGMRs. The results of this analysis show that the mean seismic collapse capacity can be predicted by the proposed subsets with less dispersion than by the recently developed improved approach, which is based on scaling the response spectra of the records to match the conditional mean spectrum. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluation of methods of spatial interpolation for monthly rainfall data over the state of Rio de Janeiro,Brazil

Five deterministic methods of spatial interpolation of monthly rainfall were compared over the state of Rio de Janeiro, southeast Brazil. The methods were the inverse distance weight (IDW), nearest neighbor (NRN), triangulation with linear interpolation (TLI), natural neighbor (NN), and spline tension (SPT). A set of 110 weather stations was used to test the methods. The selection of stations had two criteria: time series longer than 20 years and period of data from 1960 to 2009. The methods were evaluated using cross-validation, linear regression between values observed and interpolated, root mean square error (RMSE), coefficient of determination (r ²), coefficient of variation (CV, %), and the Willmott index of agreement (d). The results from different methods are influenced by the meteorological systems and their seasonality, as well as by the interaction with the topography. The methods presented higher precision (r ²) and accuracy (d, RMSE) during the summer and transition to autumn, in comparison with the winter or spring months. The SPT had the highest precision and accuracy in relation to other methods, in addition to having a good representation of the spatial patterns expected for rainfall over the complex terrain of the state and its high spatial variability.