The role of temperature on treeline migration for an eastern African mountain during the Last Glacial Maximum

Paleo-data suggest that East African mountain treelines underwent an altitudinal shift during the Last Glacial Maximum (LGM). Understanding the ecological and physiological processes underlying treeline response to such past climate change will help to improve forecasts of treeline change under future global warming. In spite of significant improvements in paleoclimatic reconstruction, the climatic conditions explaining this migration are still debated and important factors such as atmospheric CO₂ concentration, the impact of lapse rate decreasing temperature along altitudinal gradients and rainfall modifications due to elevation have often been neglected or simplified. Here, we assess the effects of these different factors and estimate the influence of the most dominant factors controlling changes in past treeline position using a multi-proxy approach based on simulations from BIOME4, a coupled biogeography and biogeochemistry model, modified to account for the effect of elevation on vegetation, compared with pollen, and isotopic data. The results indicate a shift in mountain vegetation at the LGM was controlled by low pCO₂ and low temperatures promoting species morphologically and physiologically better adapted to LGM conditions than many trees composing the forest belt limit. Our estimate that the LGM climate was cooler than today’s by ?4.5 °C (range: ?4.3 to ?4.6 °C) at the upper limit of the treeline, whereas at 831 m it was cooler by ?1.4 °C (range: ?2.6 to ?0.6 °C), suggests that a possible lapse rate modification strongly constrained the upper limit of treeline, which may limit its potential extension under future global warming.     

Trends in the frequency of synoptic types in central-southern Chile in the period 1961–2012 using the Jenkinson and Collison synoptic classification

Atmospheric circulation patterns in southern Chile (42° 30′ S) were studied in order to determine and analyse the most characteristic synoptic types and their recent trends, as well as to gain an understanding of how they are associated with low-frequency variability patterns. According to the Jenkinson and Collison (J&C) classification method, a 16-point grid of sea-level pressure data was employed. The findings reveal that some synoptic types show statistically significant trends with a 95% confidence level, positively for anticyclonic westerly hybrids (AW) and advective types for third and fourth quadrant wind flows (W, NW, and N) and negatively for SW and cyclonic hybrids (CS and CSW). A model has been constructed of the linear regression of some weather types with teleconnections that most affect Chile: the undetermined types (U), AW were associated with El Niño or the warm phase of the Pacific Decadal Oscillation (PDO), whereas the cyclonic northerly and cyclonic northeasterly types (CN and CNE) were associated with La Niña or cool phase of the PDO. The weather types associated with Antarctic Oscillation (AAO) in its positive phase are anticyclonic northerly and northeasterly and northerly advection types, while in its negative phase are cyclonic southwesterly and advection types.     

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.     

Impacts of Indian Ocean SST biases on the Indian Monsoon: as simulated in a global coupled model

In this study, the impact of the ocean–atmosphere coupling on the atmospheric mean state over the Indian Ocean and the Indian Summer Monsoon (ISM) is examined in the framework of the SINTEX-F2 coupled model through forced and coupled control simulations and several sensitivity coupled experiments. During boreal winter and spring, most of the Indian Ocean biases are common in forced and coupled simulations, suggesting that the errors originate from the atmospheric model, especially a dry islands bias in the Maritime Continent. During boreal summer, the air-sea coupling decreases the ISM rainfall over South India and the monsoon strength to realistic amplitude, but at the expense of important degradations of the rainfall and Sea Surface Temperature (SST) mean states in the Indian Ocean. Strong SST biases of opposite sign are observed over the western (WIO) and eastern (EIO) tropical Indian Ocean. Rainfall amounts over the ocean (land) are systematically higher (lower) in the northern hemisphere and the south equatorial Indian Ocean rainfall band is missing in the control coupled simulation. During boreal fall, positive dipole-like errors emerge in the mean state of the coupled model, with warm and wet (cold and dry) biases in the WIO (EIO), suggesting again a significant impact of the SST errors. The exact contributions and the distinct roles of these SST errors in the seasonal mean atmospheric state of the coupled model have been further assessed with two sensitivity coupled experiments, in which the SST biases are replaced by observed climatology either in the WIO (warm bias) or EIO (cold bias). The correction of the WIO warm bias leads to a global decrease of rainfall in the monsoon region, which confirms that the WIO is an important source of moisture for the ISM. On the other hand, the correction of the EIO cold bias leads to a global improvement of precipitation and circulation mean state during summer and fall. Nevertheless, all these improvements due to SST corrections seem drastically limited by the atmosphere intrinsic biases, including prominently the unimodal oceanic position of the ITCZ (Inter Tropical Convergence Zone) during summer and the enhanced westward wind stress along the equator during fall.     

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.     

Atmospheric boundary-layer structure observed during a haze event due to forest-fire smoke

During a haze event in Baltimore, U.S.A. from July 6 to 8, 2002, smoke from forest fires in the Québec region (Canada), degraded air quality and impacted upon local climate, decreasing solar radiation and air temperature. The smoke particles in and above the atmospheric boundary layer (ABL) served as a tracer and provided a unique opportunity to investigate the ABL structure, especially entrainment. Elastic backscatter lidar measurements taken during the haze event distinctly reveal the downward sweeps (or wisps) of smoke-laden air from the free atmosphere into the ABL. Visualisations of mechanisms such as dry convection, the entrainment process, detrainment, coherent entrainment structures, and mixing inside the ABL, are presented. Thermals overshooting at the ABL top are shown to create disturbances in the form of gravity waves in the free atmosphere aloft, as evidenced by a corresponding ripple structure at the bottom of the smoke layer. Lidar data, aerosol ground-based measurements and supporting meteorological data are used to link free atmosphere, mixed-layer and ground-level aerosols. During the peak period of the haze event (July 7, 2002), the correlation between time series of elastic backscatter lidar data within the mixed layer and the scattering coefficient from a nephelometer at ground level was found to be high (R=0.96 for z =324 m, and R=0.89 for z=504 m). Ground-level aerosol concentration was at a maximum about 2 h after the smoke layer intersected with the growing ABL, confirming that the wisps do not initially reach the ground.     

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.     

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).     

Climate Suitability for Stable Malaria Transmission in Zimbabwe Under Different Climate Change Scenarios

Climate is one factor that determines the potential range of malaria. As such, climate change may work with or against efforts to bring malaria under control. We developed a model of future climate suitability for stable Plasmodium falciparum malaria transmission in Zimbabwe. Current climate suitability for stable malaria transmission was based on the MARA/ARMA model of climatic constraints on the survival and development of the Anopheles vector and the Plasmodium falciparum malaria parasite. We explored potential future geographic distributions of malaria using 16 projections of climate in 2100. The results suggest that, assuming no future human-imposed constraints on malaria transmission, changes in temperature and precipitation could alter the geographic distribution of malaria in Zimbabwe, with previously unsuitable areas of dense human population becoming suitable for transmission. Among all scenarios, the highlands become more suitable for transmission, while the lowveld and areas with low precipitation show varying degrees of change, depending on climate sensitivity and greenhouse gas emission stabilization scenarios, and depending on the general circulation model used. The methods employed can be used within or across other African countries.