Statistical Modelling of Extreme Rainfall Indices using Multivariate Extreme Value Distributions

Environmental Modeling & Assessment - Multivariate extreme value models are used to investigate the combined behaviour of several weather variables. To investigate joint dependence of extreme...     

Dynamic extreme values modeling and monitoring by means of sea shores water quality biomarkers and valvometry

Water quality can be evaluated using biomarkers such as tissular enzymatic activities of endemic species. Measurement of molluscs bivalves activity at high frequency (e.g., valvometry) during a long time period is another way to record the animal behavior and to evaluate perturbations of the water quality in real time. As the pollution affects the activity of oysters, we consider the valves opening and closing velocities to monitor the water quality assessment. We propose to model the huge volume of velocity data collected in the framework of valvometry using a new nonparametric extreme values statistical model. The objective is to estimate the tail probabilities and the extreme quantiles of the distribution of valve closing velocity. The tail of the distribution function of valve closing velocity is modeled by a Pareto distribution with parameter ??_t,τ, beyond a threshold τ according to the time t of the experiment. Our modeling approach reveals the dependence between the specific activity of two enzymatic biomarkers (Glutathione-S-transferase and acetylcholinesterase) and the continuous recording of oyster valve velocity, proving the suitability of this tool for water quality assessment. Thus, valvometry allows in real-time in situ analysis of the bivalves behavior and appears as an effective early warning tool in ecological risk assessment and marine environment monitoring.     

One century of hydrological monitoring in two small catchments with different forest coverage

Long-term data on precipitation and runoff are essential to draw firm conclusions about the behavior and trends of hydrological catchments that may be influenced by land use and climate change. Here the longest continuous runoff records from small catchments (<1 km(2)) in Switzerland (and possibly worldwide) are reported. The history of the hydrological monitoring in the Sperbel- and Rappengraben (Emmental) is summarized, and inherent uncertainties in the data arising from the operation of the gauges are described. The runoff stations operated safely for more than 90% of the summer months when most of the major flood events occurred. Nevertheless, the absolute values of peak runoff during the largest flood events are subject to considerable uncertainty. The observed differences in average, base, and peak runoff can only partly be attributed to the substantial differences in forest coverage. This treasure trove of data can be used in various ways, exemplified here with an analysis of the generalized extreme value distributions of the two catchments. These distributions, and hence flood return periods, have varied greatly in the course of one century, influenced by the occurrence of single extreme events. The data will be made publicly available for the further analysis of the mechanisms governing the runoff behavior of small catchments, as well as for testing stochastic and deterministic models.     

downscaling procedures as a tool for integration of multiple air issues

In assessing the risks associated with climate change,downscaling has proven useful in linking surfacechanges, at scales relevant to decision making, tolarge-scale atmospheric circulation derived from GCMoutput. Stochastic downscaling is related to synopticclimatology, weather-typing approaches (classifyingcirculation patterns) such as the Lamb Weather Typesdeveloped for the United Kingdom (UK), the EuropeanGrosswetterlagen (Bardossy and Plate, 1992) and thePerfect Prognosis (Perfect Prog) method from numericalweather prediction. The large-scale atmosphericcirculation is linked with site-specific observationsof atmospheric variables, such as precipitation, windspeed or temperature, within a specified region. Classifying each day by circulation patterns isachieved by clustering algorithms, fuzzy rule bases,neural nets or decision trees. The linkages areextended to GCM output to account for climate change. Stochastic models are developed from the probabilitydistributions for extreme events. Objective analysiscan be used to interpolate values of these models toother locations. The concepts and some applicationsare reviewed to provide a basis for extending thedownscaling approach to assessing the integrated riskof the six air issues: climate change, UV-B radiation,acid rain, transport of hazardous air pollutants, smogand suspended particulates.     

Vegetation NDVI Linked to Temperature and Precipitation in the Upper Catchments of Yellow River

Vegetation in the upper catchment of Yellow River is critical for the ecological stability of the whole watershed. The dominant vegetation cover types in this region are grassland and forest, which can strongly influence the eco-environmental status of the whole watershed. The normalized difference vegetation index (NDVI) for grassland and forest has been calculated and its daily correlation models were deduced by Moderate Resolution Imaging Spectroradiometer products on 12 dates in 2000, 2003, and 2006. The responses of the NDVI values with the inter-annual grassland and forest to three climatic indices (i.e., yearly precipitation and highest and lowest temperature) were analyzed showing that, except for the lowest temperature, the yearly precipitation and highest temperature had close correlations with the NDVI values of the two vegetation communities. The value of correlation coefficients ranged from 0.815 to 0.951 (p?<?0.01). Furthermore, the interactions of NDVI values of vegetation with the climatic indicators at monthly interval were analyzed. The NDVI of vegetation and three climatic indices had strong positive correlations (larger than 0.733, p?<?0.01). The monthly correlations also provided the threshold values for the three climatic indictors, to be used for simulating vegetation growth grassland under different climate features, which is essential for the assessment of the vegetation growth and for regional environmental management.     

Ordinal time series model for forecasting air quality index for ozone in Southern California

Air quality index (AQI) for ozone is currently divided into six states depending on the level of public health concern. Generalized linear type modeling is a convenient and effective way to handle the AQI state, which can be characterized as non-stationary ordinal-valued time series. Various link functions which include cumulative logit, cumulative probit, and complimentary log-log are considered, and the partial maximum likelihood method is used for estimation. For a comparison purpose, the identity link, which yields a multiple regression model on the cumulative probabilities, is also considered. Random time-varying covariates include past AQI states, various meteorological processes, and periodic components. For model selection and comparison, the partial likelihood ratio tests, AIC and SIC are used. The proposed models are applied to 3 years of daily AQI ozone data from a station in San Bernardino County, CA. An independent year-long data from the same station are used to evaluate the performance of day-ahead forecasts of AQI state. The results show that the logit and probit models remove the non-stationarity in residuals, and both models successfully forecast day-ahead AQI states with almost 90 % of the chance.     

Forecasting of cyanobacterial density in Torrão reservoir using artificial neural networks

The ability of general regression neural networks (GRNN) to forecast the density of cyanobacteria in the Torr?o reservoir (Tamega river, Portugal), in a period of 15 days, based on three years of collected physical and chemical data, was assessed. Several models were developed and 176 were selected based on their correlation values for the verification series. A time lag of 11 was used, equivalent to one sample (periods of 15 days in the summer and 30 days in the winter). Several combinations of the series were used. Input and output data collected from three depths of the reservoir were applied (surface, euphotic zone limit and bottom). The model that presented a higher average correlation value presented the correlations 0.991; 0.843; 0.978 for training, verification and test series. This model had the three series independent in time: first test series, then verification series and, finally, training series. Only six input variables were considered significant to the performance of this model: ammonia, phosphates, dissolved oxygen, water temperature, pH and water evaporation, physical and chemical parameters referring to the three depths of the reservoir. These variables are common to the next four best models produced and, although these included other input variables, their performance was not better than the selected best model.     

Development of an Ozone Forecasting Model for Non-Attainment Areas in the State of Ohio

Ground level ozone is responsible for the formation ofsmog, and for a variety of adverse effects on bothhuman and plant life. High concentrations of groundlevel ozone occur during the summer months. This paperdescribes the development of a model to forecast themaximum daily concentration of ozone as a function ofthe maximum surface temperature, for ozonenon-attainment regions in Ohio. The model wasdeveloped by statistical analysis of existing data.Site-specific models were developed initially. Theverification and evaluation of the performancecriteria of the model at each site were explored bycomparing the model with an independent datasetcollected from that site. A generalized statewidemodel was developed from the site-specific models. Theperformance criteria of this model were verified andevaluated by employing the same independent datasetsemployed for the site-specific models. An exceedencemodel to predict the occurrence of ozone exceedencesover 100 ppb has also been presented.     

Adaptation to Climatic Variability and Change: Asian Perspectives on Agriculture and Food Security

The impacts of climate change on potential rice production in Asia are reviewed in the light of the adaptation to climatic variability and change. Collaborative studies carried out by IRRI and US-EPA reported that using process-based crop simulation models increasing temperature may decrease rice potential yield up to 7.4% per degree increment of temperature. When climate scenarios predicted by GCMs were applied it was demonstrated that rice production in Asia may decline by 3.8% under the climates of the next century. Moreover, changes in rainfall pattern and distribution were also found suggesting the possible shift of agricultural lands in the region. The studies however have not taken the impacts of climatic variability into account, which often produce extreme events like that caused by monsoons and El Niño. Shifts in rice-growing areas are likely to be constrained by land-use changes occurring for other developmental reasons, which may force greater cultivation of marginal lands and further deforestation. This should be taken into account and lead to more integrated assessment, especially in developing countries where land-use change is more a top-down policy rather than farmers' decision. A key question is: To what extent will improving the ability of societies to cope with current climatic variability through changing design of agricultural systems and practices help the same societies cope with the likely changes in climate?     

Heat wave hazard classification and risk assessment using artificial intelligence fuzzy logic

The average summer temperatures as well as the frequency and intensity of hot days and heat waves are expected to increase due to climate change. Motivated by this consequence, we propose a methodology to evaluate the monthly heat wave hazard and risk and its spatial distribution within large cities. A simple urban climate model with assimilated satellite-derived land surface temperature images was used to generate a historic database of urban air temperature fields. Heat wave hazard was then estimated from the analysis of these hourly air temperatures distributed at a 1-km grid over Athens, Greece, by identifying the areas that are more likely to suffer higher temperatures in the case of a heat wave event. Innovation lies in the artificial intelligence fuzzy logic model that was used to classify the heat waves from mild to extreme by taking into consideration their duration, intensity and time of occurrence. The monthly hazard was subsequently estimated as the cumulative effect from the individual heat waves that occurred at each grid cell during a month. Finally, monthly heat wave risk maps were produced integrating geospatial information on the population vulnerability to heat waves calculated from socio-economic variables.     

A new model for climate system analysis: Outline of the model and application to palaeoclimate simulations

We present a new reduced-form model for climate system analysis. This model, called CLIMBER-2 (for CLIMate and BiosphERe, level 2), fills the current gap between simple, highly parameterized climate models and computationally expensive coupled models of global atmospheric and oceanic circulation. We outline the basic assumptions implicit in CLIMBER-2 and we present examples of climate system analysis including a study of atmosphere–ocean interaction during the last glacial maximum, an analysis of synergism between various components of the climate system during the mid-Holocene around 6000 years ago, and a transient simulation of climate change during the last 8000 years. These studies demonstrate the feasibility of a computationally efficient analysis of climate system dynamics which is a prerequisite for future climate impact research and, more generally, Earth system analysis, i.e., the analysis of feedbacks between our environment and human activities.     

ANN modelling of sediment concentration in the dynamic glacial environment of Gangotri in Himalaya

The present study explores for the first time the possibility of modelling sediment concentration with artificial neural networks (ANNs) at Gangotri, the source of Bhagirathi River in the Himalaya. Discharge, rainfall and temperature have been considered as the main controlling factors of variations in sediment concentration in the dynamic glacial environment of Gangotri. Fourteen feed forward neural networks with error back propagation algorithm have been created, trained and tested for prediction of sediment concentration. Seven models (T1-T7) have been trained and tested in the non-updating mode whereas remaining seven models (T1a-T7a) have been trained in the updating mode. The non-updating mode refers to the scenario where antecedent time (previous time step) values are not used as input to the model. In case of the updating mode, antecedent time values are used as network inputs. The inputs applied in the models are either the variables mentioned above as individual factors (single input networks) or a combination of them (multi-input networks). The suitability of employing antecedent time-step values as network inputs has hence been checked by comparative analysis of model performance in the two modes. The simple feed forward network has been improvised with a series parallel non-linear autoregressive with exogenous input (NARX) architecture wherein true values of sediment concentration have been fed as input during training. In the glacial scenario of Gangotri, maximum sediment movement takes place during the melt period (May–October). Hence, daily data of discharge, rainfall, temperature and sediment concentration for five consecutive melt periods (May–October, 2000–2004) have been used for modelling. High Coefficient of determination values [0.77–0.88] have been obtained between observed and ANN-predicted values of sediment concentration. The study has brought out relationships between variables that are not reflected in normal statistical analysis. A strong rainfall: sediment concentration and temperature: sediment concentration relationship is shown by the models which are not reflected in statistical correlation. It has also been observed that usage of antecedent time-step values as network inputs does not necessarily lead to improvement in model performance.     

Reproducibility and imputation of air toxics data

Ambient air quality datasets include missing data, values below method detection limits and outliers, and the precision and accuracy of the measurements themselves are often unknown. At the same time, many analyses require continuous data sequences and assume that measurements are error-free. While a variety of data imputation and cleaning techniques are available, the evaluation of such techniques remains limited. This study evaluates the performance of these techniques for ambient air toxics measurements, a particularly challenging application, and includes the analysis of intra- and inter-laboratory precision. The analysis uses an unusually complete-dataset, consisting of daily measurements of over 70 species of carbonyls and volatile organic compounds (VOCs) collected over a one year period in Dearborn, Michigan, including 122 pairs of replicates. Analysis was restricted to compounds found above detection limits in > or =20% of the samples. Outliers were detected using the Gumbell extreme value distribution. Error models for inter- and intra-laboratory reproducibility were derived from replicate samples. Imputation variables were selected using a generalized additive model, and the performance of two techniques, multiple imputation and optimal linear estimation, was evaluated for three missingness patterns. Many species were rarely detected or had very poor reproducibility. Error models developed for seven carbonyls showed median intra- and inter-laboratory errors of 22% and 25%, respectively. Better reproducibility was seen for the 16 VOCs meeting detection and reproducibility criteria. Imputation performance depended on the compound and missingness pattern. Data missing at random could be adequately imputed, but imputations for row-wise deletions, the most common type of missingness pattern encountered, were not informative. The analysis shows that air toxics data require significant efforts to identify and mitigate errors, outliers and missing observations, and that these steps are essential and should be performed prior to using these data in receptor, exposure, health and other applications.     

Sustainability of winter tourism in a changing climate over Kashmir Himalaya

Mountain areas are sensitive to climate change. Implications of climate change can be seen in less snow, receding glaciers, increasing temperatures, and decreasing precipitation. Climate change is also a severe threat to snow-related winter sports such as skiing, snowboarding, and cross-country skiing. The change in climate will put further pressure on the sensitive environment of high mountains. Therefore, in this study, an attempt has been made to know the impact of climate change on the snow precipitation, water resources, and winter tourism in the two famous tourist resorts of the Kashmir Valley. Our findings show that winters are getting prolonged with little snow falls on account of climate change. The average minimum and maximum temperatures are showing statistically significant increasing trends for winter months. The precipitation is showing decreasing trends in both the regions. A considerable area in these regions remains under the snow and glacier cover throughout the year especially during the winter and spring seasons. However, time series analysis of LandSat MODIS images using Normalized Difference Snow Index shows a decreasing trend in snow cover in both the regions from past few years. Similarly, the stream discharge, comprising predominantly of snow- and glacier-melt, is showing a statistically significant declining trend despite the melting of these glaciers. The predicted futuristic trends of temperature from Predicting Regional Climates for Impact Studies regional climate model are showing an increase which may enhance snow-melting in the near future posing a serious threat to the sustainability of winter tourism in the region. Hence, it becomes essential to monitor the changes in temperature and snow cover depletion in these basins in order to evaluate their effect on the winter tourism and water resources in the region.     

Prediction of extreme ozone levels in Barcelona, Spain

Barcelona is one of the most polluted cities in Western Europe, although our levels of air pollution are within the World Health Organisation air quality guidelines. However, high concentrations of air pollution have not been studied yet. Ground ozone levels is a topic of considerable environmental concern, since excessive level of ozone are taken as indicative of high pollution. In terms of the air quality guidelines ozone levels higher than 100 µg m^–3 can start to be health-hazards for human health. Our objective is to report a detailed analysis of ozone data exceeding the thresholds established by the air quality guidelines. Data analysed were collected in two measurement stations in Barcelona, for the reference period 1991–1996. Applying statistical techniques commonly used in the analysis of extreme values, mainly the Peak Over Threshold method was used for in this study. The analysis reveal that the ozone threshold values for the protection of human health has exceeded many times in both stations. The estimated return values for 3, 10, and 40 yr exceed the threshold value for information to the public of almost once in both stations, also it seems to be unlikely that the threshold value for warning to the public will be exceeded in 40 yr.     

Biophysical and socioeconomic characterization of a water-stressed area and simulating agri-production estimates and land use planning under normal and extreme climatic events: a case study

This study aims in linking the biophysical and socioeconomic data base layers with the technical coefficients or simulation models for agri-production estimates and land use planning under normal and extreme climatic events, and exploring the resource and inputs management options in village Shikohpur, Gurgaon district located in the northwest part of India. The socioeconomic profile of Shikohpur is highly skewed with mostly small and marginal farmers. Though the areas under wheat in Shikohpur are increasing, the productivity is declining or remaining stagnant over the years. Most of the area during kharif season (June-September) remains fallow. Pearl millet based cropping systems (pearl millet-mustard and pearl millet-wheat) are predominant. Soils are mostly loamy sand to sandy loam with average of 70-80% sand content. Organic C content in soil is less than 0.3%, due to high prevailing temperature with little rainfall and also intensive agriculture followed in this region. Though the annual average seasonal rainfall in Gurgaon did not have much variation over the years, occurrence of extreme climate events has increased in the last two decades. The crop intensity is low and the water table is declining. Water and nitrogen production functions were developed for the important crops of the region, for their subsequent use in scheduling of the inputs. InfoCrop, WTGROWS and technical coefficients were used for crop planning and resource management under climate change and its variability, extreme events, limited resource availability and crop intensification. These will help in disseminating necessary agro-advisories to the farmers so that they will be able to manipulate the inputs and agronomic management practices for sustained agricultural production under normal as well as extreme climatic conditions.     

Comparing impacts across climate models

In this paper we combine a climate-forecasting model, COSMIC, with a global impact model, GIM, to compare the market impacts of climate change projected by 14 general circulation models. Given a specific date (2100), carbon dioxide concentration (612 ppmv), and global temperature sensitivity (2.5°C), predicted impacts to economies are calculated using climate-response functions from Experimental and Cross-sectional evidence. The Cross-sectional impact model predicts small global benefits across all climate models, whereas the Experimental impact model predicts a range from small benefits to small damages. High-latitude countries are less sensitive to temperature increases than low-latitude countries because they are currently cool. Uniform global temperature changes overestimate global damages because they underestimate the benefits in polar regions and overestimate the damages in tropical regions compared to the GCM predictions.