Climate change and the incidence of a forest pest in Mediterranean ecosystems: can the North Atlantic Oscillation be used as a predictor?

Many forest pest species strongly depend on temperature in their population dynamics, so that rising temperatures worldwide as a consequence of climatic change are leading to increased frequencies and intensities of insect-pest outbreaks. In the Mediterranean area, the climatic conditions are strongly linked to the effects of the North Atlantic Oscillation (NAO). The aim of this work is to analyze the dynamics of the pine processionary moth (Thaumetopoea pityocampa), a severe pest of Pinus species in the Circunmediterranean, throughout a region of southern Spain, in relation to NAO indices. We related the percentage of forest plots with high defoliation by pine processionary moth each year with NAO values for the present and the three previous winters, using generalized linear models with a binomial error distribution. The time series is 16-year long, and we performed analyses for the whole database and for the five main pine species separately. We found a consistent relationship between the response variable and the NAO index. The relationship is stronger with pine species living at medium-high altitudes, such as Aleppo (P. halepensis), black (P. nigra), and Scots (Pinus sylvestris) pine, which show the higher defoliation intensities up to 3?years after a negative NAO phase. The results highlight, for the first time, the usefulness of using global drivers in order to understand the dynamics of pest outbreaks at a regional scale, and they open the window to the development of NAO-based predictive models as an early-warning signal of severe pest outbreaks.     

Assessment of global warming on the island of Tenerife, Canary Islands (Spain). Trends in minimum, maximum and mean temperatures since 1944

Temperature variation is studied at different altitudes and orientation on the island of Tenerife, according to the trends in the mean, maximum and minimum at 21 meteorological stations. Reference series are obtained by sectors, along with a representative overall series for Tenerife, in which temperature shows a statistically significant growth trend of 0.09?±?0.04°C/decade since 1944. Night-time temperatures have risen most (0.17°C?±?0.04°C/decade), while by day they have been more stable. Consequently, the diurnal temperature range between day and night has narrowed. By regions, warming has been much more intense in the high mountains than the other sectors below the inversion layer between 600 and 1,400?m altitude, and progressively milder towards the coast. The temperature rise on the windward (north-northeast) slopes is greater than on the leeward side and could be related to the increase in cloudiness on the northern side. The general warming of the island is less than in continental areas at between 24 and 44oN, being closer to the sea surface temperature in the same area. This is probably explained largely by the insular conditions. In fact warming is more evident in the high mountains (0.14?±?0.07°C/decade), where the tempering effect of the ocean and the impact of changes in the stratocumulus is weaker, being similar to the mean continental values in the northern hemisphere.     

An upgraded version of the Eta model

Upgrades implemented over a number of years in an open source version of the Eta model, posted at the CPTEC web site http://etamodel.cptec.inpe.br/, are summarized and examples of benefits are shown. The version originates from the NCEP’s Workstation Eta code posted on the NCEP web site http://www.emc.ncep.noaa.gov/mmb/wrkstn_eta, which differs from the NCEP’s latest operational Eta by having the WRF-NMM nonhydrostatic option included. Most of the upgrades made resulted from attention paid to less than satisfactory performance noted in several Eta results, and identification of the reasons for the problem. Others came from simple expectation that including a feature that is physically justified but is missing in the code should help. The most notable of the upgrades are the introduction of the so-called sloping steps, or discretized shaved cells topography; piecewise-linear finite-volume vertical advection of dynamic variables; vapor and hydrometeor loading in the hydrostatic equation, and changes aimed at refining the convection schemes available in the Eta. Several other modifications have to do with the calculation of exchange coefficients, conservation in the vertical diffusion, and diagnostic calculation of 10-m winds. Several examples showing improved performance resulting from the dynamics changes are given. One includes a case of unrealistically low temperatures in several mountain basins generated by a centered vertical advection difference scheme’s unphysical advection from below ground, removed by its replacement with a finite-volume scheme. Another is that of increased katabatic winds in the Terra Nova Bay Antarctica region. Successful forecast of the severe downslope zonda wind case in the lee of the highest peaks of the Andes is also shown, and some of the recent successful verification results of the use of the upgraded model are pointed out. The code is used at numerous places, and along with setup information it is available for outside users at the CPTEC Eta web site given above.     

CO2 dilution in the lower atmosphere from temperature and wind speed profiles

Temperature and wind speed profiles obtained from 3?years of radio acoustic sounding system sodar measurements at a rural site in the northern Spanish plateau were fitted to polynomial functions. Depending on the extrema of these fits, several groups of profiles were considered. Daily evolution of temperature profiles corresponded to the lower boundary layer evolution. However, wind speed profiles revealed a frequent low-level jet during the whole day. CO₂ surface concentrations were analysed, and surface CO₂ dilution was also considered by selection of thin canopies with variable depth, resulting in dilution rates of 7 and 18?ppm when the layer increased 100?m for the 95th percentile and temperature and wind speed profiles, respectively.     

A seasonal study of the atmospheric dynamics over the Iberian Peninsula based on circulation types

A seasonal analysis of the atmospheric circulation over the Iberian Peninsula (IP) based on circulation types (CTs) obtained from sea level pressure and 500-hPa geopotential height is presented. The study covers the period of 1958–2008, when a high variability and important changes in winter and spring precipitation and temperature have been reported. Frequency, persistence, and the most probable transitions of the circulation types are analyzed. Among the clustering methods available in the literature, two of the most reliable classification methods have been tested, K-means and simulated annealing and diversified randomization. A comparison of both methods over the IP is presented for winter (DJF). The quality of the circulation types obtained through both methods as well as the better stability achieved by K-means suggest this method as more appropriated for our target area. Twelve CTs were obtained for each season and were analyzed. The patterns obtained were regrouped in five general situations: anticyclonic, cyclonic, zonal, summertime, and hybrid-mixed. The analysis of frequencies of these situations offers a similar characterization of the atmospheric circulation that others previously obtained by subjective methods. The analysis of the trends in frequency and persistence for each CT shows few significant trends, mainly in winter and spring with a general decrease of the cyclonic patterns and an increase of the anticyclonic situations. This can be related to the negative precipitation trends reported by other authors. Regarding the persistence, an interesting result is that there is a high interannual variability of the persistence in autumn and spring, when patterns can persist longer than in other seasons. An analysis of the most probable transitions between the CTs has been performed, revealing the existence of cyclic sequences in all seasons. These sequences are related to the high frequency of certain patterns such as the anticyclonic situations in winter. Finally, a clear seasonal dependence of the transitions between cyclonic situations associated with extratropical disturbances was found. This dependence suggests that the transitions of low-pressure systems towards the south of the IP are more likely in spring and autumn than in winter.     

Recent trends in mean maximum and minimum air temperatures over Spain (1961–2006)

This study analyzes the mean maximum and minimum temperature trends on a monthly, seasonal, and annual timescale by applying various statistical tools to data from 476 Spanish weather stations during the period between 1961 and 2006. The magnitude of the trends was derived from the slopes of the regression lines using the least squares method, and the nonparametric Mann–Kendall test was used to determine the statistical significance of the trends. Temperature significantly increased in over 60% of the country in March, June, spring, and summer in the case of maximum temperatures and in March, May, June, August, spring, and summer for minimum temperatures. At the annual resolution, temperatures significantly increased in over 90% of Spain with a rise of around 0.3°C/decade. The maximum temperature increased at a higher rate than the minimum temperature from midsummer to early winter as well as in winter, spring, and summer and also on an annual basis.     

Singular vector decomposition of the internal variability of the Canadian Regional Climate Model

Previous studies have shown that Regional Climate Models (RCM) internal variability (IV) fluctuates in time depending on synoptic events. This study focuses on the physical understanding of episodes with rapid growth of IV. An ensemble of 21 simulations, differing only in their initial conditions, was run over North America using version 5 of the Canadian RCM (CRCM). The IV is quantified in terms of energy of CRCM perturbations with respect to a reference simulation. The working hypothesis is that IV is arising through rapidly growing perturbations developed in dynamically unstable regions. If indeed IV is triggered by the growth of unstable perturbations, a large proportion of the CRCM perturbations must project onto the most unstable singular vectors (SVs). A set of ten SVs was computed to identify the orthogonal set of perturbations that provide the maximum growth with respect to the dry total-energy norm during the course of the CRCM ensemble of simulations. CRCM perturbations were then projected onto the subspace of SVs. The analysis of one episode of rapid growth of IV is presented in detail. It is shown that a large part of the IV growth is explained by initially small-amplitude unstable perturbations represented by the ten leading SVs, the SV subspace accounting for over 70% of the CRCM IV growth in 36?h. The projection on the leading SV at final time is greater than the projection on the remaining SVs and there is a high similarity between the CRCM perturbations and the leading SV after 24–36?h tangent-linear model integration. The vertical structure of perturbations revealed that the baroclinic conversion is the dominant process in IV growth for this particular episode.     

Downscaling of South America present climate driven by 4-member HadCM3 runs

The objective of this work is to evaluate climate simulations over South America using the regional Eta Model driven by four members of an ensemble of the UK Met Office Hadley Centre HadCM3 global model. The Eta Model has been modified with the purpose of performing long-term decadal integrations and has shown to reproduce “present climate”—the period 1961–1990—reasonably well when forced by HadCM3. The global model lateral conditions with a resolution of 2.5° latitude?×?3.75° longitude were provided at a frequency of 6?h. Each member of the global model ensemble has a different climate sensitivity, and the four members were selected to span the range of uncertainty encompassed by the ensemble. The Eta Model nested in the HadCM3 global model was configured with 40-km horizontal resolution and 38 layers in the vertical. No large-scale internal nudging was applied. Results are shown for austral summer and winter at present climate defined as 1961–90. The upper and low-level circulation patterns produced by the Eta-CPTEC/HadCM3 experiment set-up show good agreement with reanalysis data and the mean precipitation and temperature with CRU observation data. The spread in the downscaled mean precipitation and temperature is small when compared against model errors. On the other hand, the benefits in using an ensemble is clear in the improved representation of the seasonal cycle by the ensemble mean over any one realization. El Ni?o and La Ni?a years were identified in the HadCM3 member runs based on the NOAA Climate Prediction Center criterion of sea surface temperature anomalies in the Ni?o 3.4 area. The frequency of the El Ni?o and La Ni?a events in the studied period is underestimated by HadCM3. The precipitation and temperature anomalies typical of these events are reproduced by most of the Eta-CPTEC/HadCM3 ensemble, although small displacements of the positions of the anomalies occur. This experiment configuration is the first step on the implementation of Eta-CPTEC/HadCM3 upcoming experiments on climate change studies that are discussed in a companion paper.     

Interannual to decadal variations of precipitation and daily maximum and daily minimum temperatures in southern Brazil

The interannual to decadal variability of precipitation and daily maximum and daily minimum temperature (TMAX and TMIN) in southern Brazil for the 1913–2006 period is investigated using indices for these variables. Relations of these indices and the sea surface temperature (SST) variability in the eastern equatorial Pacific (Niño 3.4 index) and southwestern subtropical Atlantic (SSA index) are also investigated. Analyses are based on the Morlet wavelet transform. The interannual precipitation variability during the 1913–2006 period is mostly due to the high variances that occurred in specific sub-periods. The TMAX and TMIN indices also show significant interannual variability. The coherency and phase difference analyses of the precipitation index and the SST indices show higher coherency with the Niño 3.4 than with the SSA index. On the other hand, examining the coherencies and phase difference between the temperature indices and the SST indices, weaker coherencies with the Niño 3.4 index than with the SSA index are noted. These differences are discussed in the context of previous results. The new results here, with possible application for climate monitoring tasks, refer to the spectral differences between TMIN and TMAX.     

Future climate in the Pacific Northwest

Climate models used in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) on the whole reproduce the observed seasonal cycle and twentieth century warming trend of 0.8°C (1.5°F) in the Pacific Northwest, and point to much greater warming for the next century. These models project increases in annual temperature of, on average, 1.1°C (2.0°F) by the 2020s, 1.8°C (3.2°F) by the 2040s, and 3.0°C (5.3°F) by the 2080s, compared with the average from 1970 to 1999, averaged across all climate models. Rates of warming range from 0.1°C to 0.6°C (0.2°F to 1.0°F) per decade. Projected changes in annual precipitation, averaged over all models, are small (+1% to +2%), but some models project an enhanced seasonal cycle with changes toward wetter autumns and winters and drier summers. Changes in nearshore sea surface temperatures, though smaller than on land, are likely to substantially exceed interannual variability, but coastal upwelling changes little. Rates of twenty-first century sea level rise will depend on poorly known factors like ice sheet instability in Greenland and Antarctica, and could be as low as twentieth century values (20 cm, 8^″) or as large as 1.3 m (50^″).