Climatic signals in the life histories of insects:the distribution and abundance of heather psyllids (Strophingia spp.) in the UK

1. The population density and age structure of two species of heather psyllid Strophingia ericae and Strophingia cinereae, feeding on Calluna vulgaris and Erica cinerea, respectively, were sampled using standardized methods at locations throughout Britain. Locations were chosen to represent the full latitudinal and altitudinal range of the host plants.

2. The paper explains how spatial variation in thermal environment, insect life-history characteristics and physiology, and plant distribution, interact to provide the mechanisms that determine the range and abundance of Strophingia spp.

3. Strophingia ericae and S. cinereae, despite the similarity in the spatial distribution patterns of their host plants within Britain, display strongly contrasting geographical ranges and corresponding life-history strategies. Strophingia ericae is found on its host plant throughout Britain but S. cinereae is restricted to low elevation sites south of the Mersey-Humber line and occupies only part of the latitudinal and altitudinal range of its host plant. There is no evidence to suggest that S. ericae has reached its potential altitudinal or latitudinal limit in the UK, even though its host plant appears to reach its altitudinal limit.

4. There was little difference in the ability of the two Strophingia spp. to survive shortterm exposure to temperatures as low as - 15 degrees C and low winter temperatures probably do not limit distribution in S. cinereae.

5. Population density of S. ericae was not related to altitude but showed a weak correlation with latitude. The spread of larval instars present at a site, measured as an index of instar homogeneity, was significantly correlated with a range of temperature related variables, of which May mean temperature and length of growing season above 3 degrees C (calculated using the Lennon and Turner climatic model) were the most significant. Factor analysis did not improve the level of correlation significantly above those obtained for single climatic variables. The data confirmed that S. ericae has a I year life cycle at the lowest elevations and a 2 year life cycle at the higher elevations. However, there was no evidence, as previously suggested, for an abrupt change from a one to a 2 year life cycle in S. ericae with increasing altitudes or latitudes.

6. By contrast with S. ericae, S. cinereae had an obligatory 1 year life cycle, its population decreased with altitude and the index of instar homogeneity showed little correlation with single temperature variables. Moreover, it occupied only part of the range of its host plant and its spatial distribution in the UK could be predicted with 96% accuracy using selected variables in discriminant analysis.

7. The life histories of the congeneric heather psyllids reflect adaptations that allow them to exploit host plants with different distributions in climatic and thereby geographical space. Strophingia ericae has the flexible life history that enables it to exploit C. vulgaris throughout its European boreal temperate range. Strophingia cinereae has a less flexible life history and is adapted for living on an oceanic temperate host. While the geographic ranges of the two Strophingia spp. overlap within the UK, the psyllids appear to respond differently to variation in their thermal environment.

A simple regional-scale model for forecasting sorghum yield across North-Eastern Australia

Sorghum is the main dryland summer crop in NE Australia and a number of agricultural businesses would benefit from an ability to forecast production likelihood at regional scale. In this study we sought to develop a simple agro-climatic modelling approach for predicting shire (statistical local area) sorghum yield. Actual shire yield data, available for the period 1983-1997 from the Australian Bureau of Statistics, were used to train the model. Shire yield was related to a water stress index (SI) that was derived from the agro-climatic model. The model involved a simple fallow and crop water balance that was driven by climate data available at recording stations within each shire. Parameters defining the soil water holding capacity, maximum number of sowings (MXNS) in any year, planting rainfall requirement, and critical period for stress during the crop cycle were optimised as part of the model fitting procedure. Cross-validated correlations (CVR) ranged from 0.5 to 0.9 at shire scale. When aggregated to regional and national scales, 78-84% of the annual variation in sorghum yield was explained. The model was used to examine trends in sorghum productivity and the approach to using it in an operational forecasting system was outlined. (c) 2005 Elsevier B.V. All rights reserved.

The potential geographic distribution of the invasive weed Senna obtusifolia in Australia

Senna obtusifolia (sicklepod) is an invasive weed of northern Australia, where it significantly impacts agricultural productivity and alters natural ecosystem structure and function. Although currently restricted to northern regions, the potential for S. obtusifolia to spread south is not known. Using the eco-climatic model CLIMEX, this study simulated the potential geographic distribution of S. obtusifolia in Australia under two scenarios. Model parameters for both scenarios were derived from the distribution of S. obtusifolia throughout North and Central America. The first scenario used these base model parameters to predict the distribution of S. obtusifolia in Australia, whilst the second model predicted the distribution of a cold susceptible S. obtusifolia ecotype that is reported to occur in the USA. Both models predicted the potential for an extensive S. obtusifolia distribution, with the first model indicating suitable climatic conditions occurring predominantly in coastal regions from the Northern Territory, to far north Queensland and into northern Victoria. The cold susceptible ecotype displayed a comparatively reduced distribution in the southern parts of Australia, where inappropriate temperatures, a lack of thermal accumulation and cold stress restrict the invasion south to the coastal regions of central New South Wales. The extent of the predicted distribution of both ecotypes of S. obtusifolia reinforces the need for strategic management at a national scale.

Épocas de plantio de milho em função das deficiências hídricas no solo em Cambará-PR

The aim of this study was to try the forecast of corn (Zea mays L.) sowing dates, the understanding of the quantitative effect of water deficits on that crop and crop yield decrease on a basis of a climatic model of water deficit forecast. This model was applied at Cambara (lat. 23 degrees 00'S, long. 50 circle 02'WGr, altitude 450 m), PR, Brazil. The model estimates yield decrease, in relation to potential values, as a function of the sowing dates which determine flowering and grain filling dates, highly critical times in relation to water deficit. The estimates were done from expected values of water deficit, at the 80% probability level and accumulated degrees-days, using several climatological data. Results show that the first ten days of November are the best corn sowing date under dry or irrigated conditions. Under these same conditions, the worst time showed to be August. Estimates of total needs of supplementar irrigation get values of 126 e 226 mm, respectively to the corn sowed at the first ten days of November and August.

Aspectos regionais da variabilidade de precipitação no estado do Pará: estudo observacional e modelagem climática em alta resolução

O objetivo do presente trabalho foi agregar diferentes redes de estações meteorológicas de superfície para a criação de um novo banco de dados integrado, a partir do qual foi gerada uma climatologia recente (1978-2007) para a precipitação do estado do Pará em alta resolução espacial – 30 km, permitindo melhor identificar a variabilidade climática regional, sobretudo influenciada pelos aspectos da fisiografia e em função de mecanismos climáticos de grande escala dos oceanos Pacífico e Atlântico. Buscou-se, ainda, estabelecer uma configuração otimizada do modelo climático RegCM3 utilizando duas diferentes parametrizações de cumulus: RegCM3/Grell e RegCM3/MIT. Foram realizadas 26 simulações (1982/83 a 2007/08) durante a estação chuvosa na Amazônia oriental (dezembro a maio) para cada esquema de parametrização convectiva, utilizando 30 km de resolução espacial. O modelo mostrou-se capaz de capturar os sinais de anomalia na presença de forçantes climáticas extremas, como o El Niño-Oscilação Sul e o dipolo do Atlântico. O RegCM3/MIT obteve ótimo desempenho na região de Altamira/PA e performance razoável nos setores Nordeste (região de Belém), Leste ( região de Marabá), Sudeste (região de Conceição do Araguaia), e Noroeste (região de Tiriós). O RegCM3/Grell destacou-se nas regiões Nordeste, Leste, Sudeste e Noroeste, com desempenho razoável. O setor Norte (região de Macapá) foi o mais problemático, com pouca ou nenhuma sensibilidade apresentada pelo modelo. Embora o RegCM3 tenha obtido resultados razoáveis na maior parte do domínio, foram detectados erros sistemáticos nas simulações, com viés seco para o RegCM3/Grell e viés úmido para o RegCM3/MIT na porção Sul e viés seco na porção Norte. Estas características denotam a necessidade de ajustes às condições regionais dos esquemas de convecção.

Characterization of late Campanian and Maastrichtian planktonic foraminiferal depth habitats and vital activities based on stable isotopes

Depth habitats of 56 late Cretaceous planktonic foraminiferal species from cool and warm climate modes were determined based on stable isotope analyses of deep-sea samples from the equatorial Pacific DSDP Sites 577A and 463, and South Atlantic DSDP Site 525A. The following conclusions can be reached: Planoglobulina multicamerata (De Klasz) and Heterohelix rajagopalani (Govindan) occupied the deepest plankton habitats, followed by Abathomphalus mayaroensis (Bolli), Globotruncanella havanensis (Voorwijk), Gublerina cuvillieri Kikoine, and Laeviheterohelix glabrans (Cushman) also at subthermocline depth. Most keeled globotruncanids, and possibly Globigerinelliodes and Racemiguembelina species, lived at or within the thermocline layer. Heterohelix globulosa (Ehrenberg) and Rugoglobigerina, Pseudotextularia and Planoglobulina occupied the subsurface depth of the mixed layer, and Pseudoguembelina species inhabited the surface mixed layer. However, depth ranking of some species varied depending on warm or cool climate modes, and late Campanian or Maastrichtian age. For example, most keeled globotruncanids occupied similar shallow subsurface habitats as Rugoglobigerina during the warm late Campanian, but occupied the deeper thermocline layer during cool climatic intervals. Two distinct types of "vital effect" mechanisms reflecting photosymbiosis and respiration effects can be recognized by the exceptional delta13C signals of some species. (1) Photosymbiosis is implied by the repetitive pattern of relatively enriched delta13C values of Racemiguembelina (strongest), Planoglobulina, Rosita and Rugoglobigerina species, Pseudoguembelina excolata (weakest). (2) Enriched respiration 12C products are recognized in A. mayaroensis, Gublerina acuta De Klasz, and Heterohelix planata (Cushman). Isotopic trends between samples suggest that photosymbiotic activities varied between localities or during different climate modes, and may have ceased under certain environmental conditions. The appearance of most photosymbiotic species in the late Maastrichtian suggests oligotrophic conditions associated with increased water-mass stratification.

Estudio de la macroflora en yacimientos del Neógeno y Pleistoceno inferior del noroeste peninsular: significación paleoclimática

La presente tesis doctoral se centra en el estudio de dos yacimientos de macrorrestos fosiles vegetales, en especial maderas en diferentes estados de conservacion, y su interpretacion dentro del contexto del Neogeno Iberico. En primer lugar se detallan los objetivos del trabajo y los antecedentes del mismo. A continuacion, se describen detalladamente los dos yacimientos estudiados y se exponen las tecnicas empleadas de acuerdo con los diferentes tipos de materiales. Mas tarde se identifican taxonomicamente todos los restos hallados en los mismos. En el mismo capitulo de resultados se aplica el Metodo de Coexistencia sobre los principales yacimientos neogenos de la peninsula Iberica, para la estimacion de valores climaticos y se elabora un modelo de evolucion climatica para este periodo. Por ultimo, se discuten los resultados obtenidos, estableciendo una hipotesis de composicion y estructura de los bosques que vivieron en el tiempo en el que se depositaron los restos estudiados, se comparan con los datos obtenidos por otros autores y se combinan estos resultados con los datos del modelo de clima para reconstruir una sintesis de la evolucion de la vegetacion de la Peninsula durante el Periodo Neogeno y Gelasiense. Abstract This PhD thesis is based on the study of two different sites of vegetal fossil remains, specially wood in different preservation states, and their meaning in Iberian Neogene context. First of all, the objectives of the work and previous records are explained thoroughly. Secondly, the two deposits are described in detail, as well as the techniques used, according to different types of materials. Later on, all the samples found in each one are identified. In the same chapter of results, the “Coexistence Approach” is applied on the main Neogen sites of the Iberian Peninsula to estimate climatic intervals, and a model of climatic evolution for this period is built. Finally, the achieved results are discussed, setting an hypothesis about the composition and structure of the forests that lived in the time when the studied samples were settled. These results are then compared with the data gathered by different authors, and related to the data of the climatic model, to make a synthesis of the evolution of the vegetation in the Iberian Peninsula, during Neogene and Gelasian.

The use of ZIP and CART to model cryptosporidiosis in relation to climatic variables

This research assesses the potential impact of weekly weather variability on the incidence of cryptosporidiosis disease using time series zero-inflated Poisson (ZIP) and classification and regression tree (CART) models. Data on weather variables, notified cryptosporidiosis cases and population size in Brisbane were supplied by the Australian Bureau of Meteorology, Queensland Department of Health, and Australian Bureau of Statistics, respectively. Both time series ZIP and CART models show a clear association between weather variables (maximum temperature, relative humidity, rainfall and wind speed) and cryptosporidiosis disease. The time series CART models indicated that, when weekly maximum temperature exceeded 31°C and relative humidity was less than 63%, the relative risk of cryptosporidiosis rose by 13.64 (expected morbidity: 39.4; 95% confidence interval: 30.9–47.9). These findings may have applications as a decision support tool in planning disease control and risk management programs for cryptosporidiosis disease.

An expressed sequence tag (EST) library for Drosophila serrata, a model system for sexual selection and climatic adaptation studies

The native Australian fly Drosophila serrata belongs to the highly speciose montium subgroup of the melanogaster species group. It has recently emerged as an excellent model system with which to address a number of important questions, including the evolution of traits under sexual selection and traits involved in climatic adaptation along latitudinal gradients. Understanding the molecular genetic basis of such traits has been limited by a lack of genomic resources for this species. Here, we present the first expressed sequence tag (EST) collection for D. serrata that will enable the identification of genes underlying sexually-selected phenotypes and physiological responses to environmental change and may help resolve controversial phylogenetic relationships within the montium subgroup.

Comprehensive evaluation of the WaterGAP3 model across climatic, physiographic, and anthropogenic gradients

In den letzten Jahrzehnten haben sich makroskalige hydrologische Modelle als wichtige Werkzeuge etabliert um den Zustand der globalen erneuerbaren Süßwasserressourcen flächendeckend bewerten können. Sie werden heutzutage eingesetzt um eine große Bandbreite wissenschaftlicher Fragestellungen zu beantworten, insbesondere hinsichtlich der Auswirkungen anthropogener Einflüsse auf das natürliche Abflussregime oder der Auswirkungen des globalen Wandels und Klimawandels auf die Ressource Wasser. Diese Auswirkungen lassen sich durch verschiedenste wasserbezogene Kenngrößen abschätzen, wie z.B. erneuerbare (Grund-)Wasserressourcen, Hochwasserrisiko, Dürren, Wasserstress und Wasserknappheit. Die Weiterentwicklung makroskaliger hydrologischer Modelle wurde insbesondere durch stetig steigende Rechenkapazitäten begünstigt, aber auch durch die zunehmende Verfügbarkeit von Fernerkundungsdaten und abgeleiteten Datenprodukten, die genutzt werden können, um die Modelle anzutreiben und zu verbessern. Wie alle makro- bis globalskaligen Modellierungsansätze unterliegen makroskalige hydrologische Simulationen erheblichen Unsicherheiten, die (i) auf räumliche Eingabedatensätze, wie z.B. meteorologische Größen oder Landoberflächenparameter, und (ii) im Besonderen auf die (oftmals) vereinfachte Abbildung physikalischer Prozesse im Modell zurückzuführen sind. Angesichts dieser Unsicherheiten ist es unabdingbar, die tatsächliche Anwendbarkeit und Prognosefähigkeit der Modelle unter diversen klimatischen und physiographischen Bedingungen zu überprüfen. Bisher wurden die meisten Evaluierungsstudien jedoch lediglich in wenigen, großen Flusseinzugsgebieten durchgeführt oder fokussierten auf kontinentalen Wasserflüssen. Dies steht im Kontrast zu vielen Anwendungsstudien, deren Analysen und Aussagen auf simulierten Zustandsgrößen und Flüssen in deutlich feinerer räumlicher Auflösung (Gridzelle) basieren. Den Kern der Dissertation bildet eine umfangreiche Evaluierung der generellen Anwendbarkeit des globalen hydrologischen Modells WaterGAP3 für die Simulation von monatlichen Abflussregimen und Niedrig- und Hochwasserabflüssen auf Basis von mehr als 2400 Durchflussmessreihen für den Zeitraum 1958-2010. Die betrachteten Flusseinzugsgebiete repräsentieren ein breites Spektrum klimatischer und physiographischer Bedingungen, die Einzugsgebietsgröße reicht von 3000 bis zu mehreren Millionen Quadratkilometern. Die Modellevaluierung hat dabei zwei Zielsetzungen: Erstens soll die erzielte Modellgüte als Bezugswert dienen gegen den jegliche weiteren Modellverbesserungen verglichen werden können. Zweitens soll eine Methode zur diagnostischen Modellevaluierung entwickelt und getestet werden, die eindeutige Ansatzpunkte zur Modellverbesserung aufzeigen soll, falls die Modellgüte unzureichend ist. Hierzu werden komplementäre Modellgütemaße mit neun Gebietsparametern verknüpft, welche die klimatischen und physiographischen Bedingungen sowie den Grad anthropogener Beeinflussung in den einzelnen Einzugsgebieten quantifizieren. WaterGAP3 erzielt eine mittlere bis hohe Modellgüte für die Simulation von sowohl monatlichen Abflussregimen als auch Niedrig- und Hochwasserabflüssen, jedoch sind für alle betrachteten Modellgütemaße deutliche räumliche Muster erkennbar. Von den neun betrachteten Gebietseigenschaften weisen insbesondere der Ariditätsgrad und die mittlere Gebietsneigung einen starken Einfluss auf die Modellgüte auf. Das Modell tendiert zur Überschätzung des jährlichen Abflussvolumens mit steigender Aridität. Dieses Verhalten ist charakteristisch für makroskalige hydrologische Modelle und ist auf die unzureichende Abbildung von Prozessen der Abflussbildung und –konzentration in wasserlimitierten Gebieten zurückzuführen. In steilen Einzugsgebieten wird eine geringe Modellgüte hinsichtlich der Abbildung von monatlicher Abflussvariabilität und zeitlicher Dynamik festgestellt, die sich auch in der Güte der Niedrig- und Hochwassersimulation widerspiegelt. Diese Beobachtung weist auf notwendige Modellverbesserungen in Bezug auf (i) die Aufteilung des Gesamtabflusses in schnelle und verzögerte Abflusskomponente und (ii) die Berechnung der Fließgeschwindigkeit im Gerinne hin. Die im Rahmen der Dissertation entwickelte Methode zur diagnostischen Modellevaluierung durch Verknüpfung von komplementären Modellgütemaßen und Einzugsgebietseigenschaften wurde exemplarisch am Beispiel des WaterGAP3 Modells erprobt. Die Methode hat sich als effizientes Werkzeug erwiesen, um räumliche Muster in der Modellgüte zu erklären und Defizite in der Modellstruktur zu identifizieren. Die entwickelte Methode ist generell für jedes hydrologische Modell anwendbar. Sie ist jedoch insbesondere für makroskalige Modelle und multi-basin Studien relevant, da sie das Fehlen von feldspezifischen Kenntnissen und gezielten Messkampagnen, auf die üblicherweise in der Einzugsgebietsmodellierung zurückgegriffen wird, teilweise ausgleichen kann.

Glacier shrinkage and climatic change in the Russian Altai from the mid-20th century: An assessment using remote sensing and PRECIS regional climate model

This paper examines changes in the surface area of glaciers in the North and South Chuya Ridges, Altai Mountains in 1952-2004 and their links with regional climatic variations. The glacier surface areas for 2004 were derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. Data from the World Glacier Inventory (WGI)dating to 1952 and aerial photographs from 1952 were used to estimate the changes. 256 glaciers with a combined area of 253±5.1 km2 have been identified in the region in 2004. Estimation of changes in extent of 126 glaciers with the individual areas not less than 0.5 km2 in 1952 revealed a 19.7±5.8% reduction. The observed glacier retreat is primarily driven by an increase in summer temperatures since the 1980s when air temperatures were increasing at a rate of 0.10 - 0.13oC a-1 at the glacier tongue elevation. The regional climate projections for A2 and B2 CO2 emission scenarios developed using PRECIS regional climate model indicate that summer temperatures will increase in the Altai in 2071-2100 by 6-7oC and 3-5oC respectively in comparison with 1961-1990 while annual precipitation will increase by 15% and 5%. The length of the ablation season will extend from June-August to the late April – early October. The projected increases in precipitation will not compensate for the projected warming and glaciers will continue to retreat in the 21st century under both B2 and A2 scenarios.

Climatic impacts of fresh water hosing under Last Glacial Maximum conditions: a multi-model study

Fresh water hosing simulations, in which a fresh water flux is imposed in the North Atlantic to force fluctuations of the Atlantic Meridional Overturning Circulation, have been routinely performed, first to study the climatic signature of different states of this circulation, then, under present or future conditions, to investigate the potential impact of a partial melting of the Greenland ice sheet. The most compelling examples of climatic changes potentially related to AMOC abrupt variations, however, are found in high resolution palaeo-records from around the globe for the last glacial period. To study those more specifically, more and more fresh water hosing experiments have been performed under glacial conditions in the recent years. Here we compare an ensemble constituted by 11 such simulations run with 6 different climate models. All simulations follow a slightly different design, but are sufficiently close in their design to be compared. They all study the impact of a fresh water hosing imposed in the extra-tropical North Atlantic. Common features in the model responses to hosing are the cooling over the North Atlantic, extending along the sub-tropical gyre in the tropical North Atlantic, the southward shift of the Atlantic ITCZ and the weakening of the African and Indian monsoons. On the other hand, the expression of the bipolar see-saw, i.e., warming in the Southern Hemisphere, differs from model to model, with some restricting it to the South Atlantic and specific regions of the southern ocean while others simulate a widespread southern ocean warming. The relationships between the features common to most models, i.e., climate changes over the north and tropical Atlantic, African and Asian monsoon regions, are further quantified. These suggest a tight correlation between the temperature and precipitation changes over the extra-tropical North Atlantic, but different pathways for the teleconnections between the AMOC/North Atlantic region and the African and Indian monsoon regions.

Climatic impacts of land-use change due to crop yield increases and a universal carbon tax from a scenario model*

Future land cover will have a significant impact on climate and is strongly influenced by the extent of agricultural land use. Differing assumptions of crop yield increase and carbon pricing mitigation strategies affect projected expansion of agricultural land in future scenarios. In the representative concentration pathway 4.5 (RCP4.5) from phase 5 of the Coupled Model Intercomparison Project (CMIP5), the carbon effects of these land cover changes are included, although the biogeophysical effects are not. The afforestation in RCP4.5 has important biogeophysical impacts on climate, in addition to the land carbon changes, which are directly related to the assumption of crop yield increase and the universal carbon tax. To investigate the biogeophysical climatic impact of combinations of agricultural crop yield increases and carbon pricing mitigation, five scenarios of land-use change based on RCP4.5 are used as inputs to an earth system model [Hadley Centre Global Environment Model, version 2-Earth System (HadGEM2-ES)]. In the scenario with the greatest increase in agricultural land (as a result of no increase in crop yield and no climate mitigation) there is a significant -0.49 K worldwide cooling by 2100 compared to a control scenario with no land-use change. Regional cooling is up to -2.2 K annually in northeastern Asia. Including carbon feedbacks from the land-use change gives a small global cooling of -0.067 K. This work shows that there are significant impacts from biogeophysical land-use changes caused by assumptions of crop yield and carbon mitigation, which mean that land carbon is not the whole story. It also elucidates the potential conflict between cooling from biogeophysical climate effects of land-use change and wider environmental aims.

Radio-wave propagation model for UHF band in different climatic conditions with dyadic green’s function

This article proposes a deterministic radio propagation model using dyadic Green's function to predict the value of the electric field. Dyadic is offered as an efficient mathematical tool which has symbolic simplicity and robustness, as well as taking account of the anisotropy of the medium. The proposed model is an important contribution for the UHF band because it considers climatic conditions by changing the constants of the medium. Most models and recommendations that include an approach for climatic conditions, are designed for satellite links, mainly Ku and Ka bands. The results obtained by simulation are compared and validated with data from a Digital Television Station measurement campaigns conducted in the Belém city in Amazon region during two seasons. The proposed model was able to provide satisfactory results by differentiating between the curves for dry and wet soil and these corroborate the measured data, (the RMS errors are between 2-5 dB in the case under study).

A population dynamical model of Operophtera brumata, L. extended by climatic factors

Setting out from the database of Operophtera brumata, L. in between 1973 and 2000 due to the Light Trap Network in Hungary, we introduce a simple theta-logistic population dynamical model based on endogenous and exogenous factors, only. We create an indicator set from which we can choose some elements with which we can improve the fitting results the most effectively. Than we extend the basic simple model with additive climatic factors. The parameter optimization is based on the minimized root mean square error. The best model is chosen according to the Akaike Information Criterion. Finally we run the calibrated extended model with daily outputs of the regional climate model RegCM3.1, regarding 1961-1990 as reference period and 2021-2050 with 2071-2100 as future predictions. The results of the three time intervals are fitted with Beta distributions and compared statistically. The expected changes are discussed.