Describing long-term trends in precipitation using generalized additive models

With the current concern over climate change, descriptions of how rainfall patterns are changing over time can be useful. Observations of daily rainfall data over the last few decades provide information on these trends. Generalized linear models are typically used to model patterns in the occurrence and intensity of rainfall. These models describe rainfall patterns for an average year but are more limited when describing long-term trends, particularly when these are potentially non-linear. Generalized additive models (GAMS) provide a framework for modelling non-linear relationships by fitting smooth functions to the data. This paper describes how GAMS can extend the flexibility of models to describe seasonal patterns and long-term trends in the occurrence and intensity of daily rainfall using data from Mauritius from 1962 to 2001. Smoothed estimates from the models provide useful graphical descriptions of changing rainfall patterns over the last 40 years at this location. GAMS are particularly helpful when exploring non-linear relationships in the data. Care is needed to ensure the choice of smooth functions is appropriate for the data and modelling objectives. (c) 2008 Elsevier B.V. All rights reserved.

Trends in the start of the wet season over Africa

A quarter of a century of daily rainfall data from the Global Telecommunications System are used to define the temporal and spatial variability of the start of the wet season over Africa and surrounding extreme south of Europe and parts of the Middle East. From 1978 to 2002, the start of the wet season arrived later in the year for the majority of the region, as time progressed. In some parts of the continent, there was an annual increase in the start date of up to 4 days per year. On average, the start of the wet season arrived 9–21 days later from 1978 to 2002, depending on the threshold used to define the start of the rains (varying from 10–30 mm over 2 days, with no dry period in the following 10 days). It is noted that the inter-annual variability of the start of the wet season is high with the range of start dates varying on average from 116 to 142 days dependent on the threshold used to determine the start date. These results may have important implications for agriculturists on all levels (from the individual farmer to those responsible for regional food supply), as knowledge of potential future climate changes starts to play an increasingly important role in the agricultural decision-making process, such as sowing and harvesting times.

Bollène-2002 experiment: radar quantitative precipitation estimation in the Cévennes–Vivarais region, France

The Bollène-2002 Experiment was aimed at developing the use of a radar volume-scanning strategy for conducting radar rainfall estimations in the mountainous regions of France. A developmental radar processing system, called Traitements Régionalisés et Adaptatifs de Données Radar pour l’Hydrologie (Regionalized and Adaptive Radar Data Processing for Hydrological Applications), has been built and several algorithms were specifically produced as part of this project. These algorithms include 1) a clutter identification technique based on the pulse-to-pulse variability of reflectivity Z for noncoherent radar, 2) a coupled procedure for determining a rain partition between convective and widespread rainfall R and the associated normalized vertical profiles of reflectivity, and 3) a method for calculating reflectivity at ground level from reflectivities measured aloft. Several radar processing strategies, including nonadaptive, time-adaptive, and space–time-adaptive variants, have been implemented to assess the performance of these new algorithms. Reference rainfall data were derived from a careful analysis of rain gauge datasets furnished by the Cévennes–Vivarais Mediterranean Hydrometeorological Observatory. The assessment criteria for five intense and long-lasting Mediterranean rain events have proven that good quantitative precipitation estimates can be obtained from radar data alone within 100-km range by using well-sited, well-maintained radar systems and sophisticated, physically based data-processing systems. The basic requirements entail performing accurate electronic calibration and stability verification, determining the radar detection domain, achieving efficient clutter elimination, and capturing the vertical structure(s) of reflectivity for the target event. Radar performance was shown to depend on type of rainfall, with better results obtained with deep convective rain systems (Nash coefficients of roughly 0.90 for point radar–rain gauge comparisons at the event time step), as opposed to shallow convective and frontal rain systems (Nash coefficients in the 0.6–0.8 range). In comparison with time-adaptive strategies, the space–time-adaptive strategy yields a very significant reduction in the radar–rain gauge bias while the level of scatter remains basically unchanged. Because the Z–R relationships have not been optimized in this study, results are attributed to an improved processing of spatial variations in the vertical profile of reflectivity. The two main recommendations for future work consist of adapting the rain separation method for radar network operations and documenting Z–R relationships conditional on rainfall type.

Precipitation diurnal cycle and summer climatology assessment over South America: An evaluation of Regional Climate Model version 3 simulations

Regional Climate Model version 3 (RegCM3) simulations of 17 summers (1988-2004) over part of South America south of 5 degrees S were evaluated to identify model systematic errors. Model results were compared to different rainfall data sets (Climate Research Unit (CRU), Climate Prediction Center (CPC), Global Precipitation Climatology Project (GPCP), and National Centers for Environmental Prediction (NCEP) reanalysis), including the five summers mean (1998-2002) precipitation diurnal cycle observed by the Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR). In spite of regional differences, the RegCM3 simulates the main observed aspects of summer climatology associated with the precipitation (northwest-southeast band of South Atlantic Convergence Zone (SACZ)) and air temperature (warmer air in the central part of the continent and colder in eastern Brazil and the Andes Mountains). At a regional scale, the main RegCM3 failures are the underestimation of the precipitation in the northern branch of the SACZ and some unrealistic intense precipitation around the Andes Mountains. However, the RegCM3 seasonal precipitation is closer to the fine-scale analyses (CPC, CRU, and TRMM-PR) than is the NCEP reanalysis, which presents an incorrect north-south orientation of SACZ and an overestimation of its intensity. The precipitation diurnal cycle observed by TRMM-PR shows pronounced contrasts between Tropics and Extratropics and land and ocean, where most of these features are simulated by RegCM3. The major similarities between the simulation and observation, especially the diurnal cycle phase, are found over the continental tropical and subtropical SACZ regions, which present afternoon maximum (1500-1800 UTC) and morning minimum (0900-1200 UTC). More specifically, over the core of SACZ, the phase and amplitude of the simulated precipitation diurnal cycle are very close to the TRMM-PR observations. Although there are amplitude differences, the RegCM3 simulates the observed nighttime rainfall in the eastern Andes Mountains, over the Atlantic Ocean, and also over northern Argentina. The main simulation deficiencies are found in the Atlantic Ocean and near the Andes Mountains. Over the Atlantic Ocean the convective scheme is not triggered; thus the rainfall arises from the grid-scale scheme and therefore differs from the TRMM-PR. Near the Andes, intense (nighttime and daytime) simulated precipitation could be a response of an incorrect circulation and topographic uplift. Finally, it is important to note that unlike most reported bias of global models, RegCM3 does not trigger the moist convection just after sunrise over the southern part of the Amazon.

Origin of Convectively Coupled Kelvin Waves over South America

Convectively coupled Kelvin waves over the South American continent are examined through the use of temporal and spatial filtering of reanalysis, satellite, and gridded rainfall data. They are most prominent from November to April, the season analyzed herein. The following two types of events are isolated: those that result from preexisting Kelvin waves over the eastern Pacific Ocean propagating into the continent, and those that apparently originate over Amazonia, forced by disturbances propagating equatorward from central and southern South America. The events with precursors in the Pacific are mainly upper-level disturbances, with almost no signal at the surface. Those events with precursors over South America, on the other hand, originate as upper-level synoptic wave trains that pass over the continent and resemble the ""cold surges`` documented by Garreaud and Wallace. As the wave train propagates over the Andes, it induces a southerly low-level wind that advects cold air to the north. Precipitation associated with a cold front reaches the equator a few days later and subsequently propagates eastward with the characteristics of a Kelvin wave. The structures of those waves originating over the Pacific are quite similar to those originating over South America as they propagate to eastern South America and into the Atlantic. South America Kelvin waves that originate over neither the Pacific nor the midlatitudes of South America can also be identified. In a composite sense, these form over the eastern slope of the Andes Mountains, close to the equator. There are also cases of cold surges that reach the equator yet do not form Kelvin waves. The interannual variability of the Pacific-originating events is related to sea surface temperatures in the central-eastern Pacific Ocean. When equatorial oceanic conditions are warm, there tends to be an increase in the number of disturbances that reach South America from the Pacific.

Dynamics and spatio-temporal variability of environmental factors in Eastern Australia using functional principal component analysis

This paper introduces a new technique in ecology to analyze spatial and temporal variability in environmental variables. By using simple statistics, we explore the relations between abiotic and biotic variables that influence animal distributions. However, spatial and temporal variability in rainfall, a key variable in ecological studies, can cause difficulties to any basic model including time evolution. The study was of a landscape scale (three million square kilometers in eastern Australia), mainly over the period of 19982004. We simultaneously considered qualitative spatial (soil and habitat types) and quantitative temporal (rainfall) variables in a Geographical Information System environment. In addition to some techniques commonly used in ecology, we applied a new method, Functional Principal Component Analysis, which proved to be very suitable for this case, as it explained more than 97% of the total variance of the rainfall data, providing us with substitute variables that are easier to manage and are even able to explain rainfall patterns. The main variable came from a habitat classification that showed strong correlations with rainfall values and soil types. © 2010 World Scientific Publishing Company.

Stochastic Assessment Of Environmental Flows In Semiarid Environments

The regimen of environmental flows (EF) must be included as terms of environmental demand in the management of water resources. Even though there are numerous methods for the computation of EF, the criteria applied at different steps in the calculation process are quite subjective whereas the results are fixed values that must be meet by water planners. This study presents a friendly-user tool for the assessment of the probability of compliance of a certain EF scenario with the natural regimen in a semiarid area in southern Spain. 250 replications of a 25-yr period of different hydrological variables (rainfall, minimum and maximum flows, ...) were obtained at the study site from the combination of Monte Carlo technique and local hydrological relationships. Several assumptions are made such as the independence of annual rainfall from year to year and the variability of occurrence of the meteorological agents, mainly precipitation as the main source of uncertainty. Inputs to the tool are easily selected from a first menu and comprise measured rainfall data, EF values and the hydrological relationships for at least a 20-yr period. The outputs are the probabilities of compliance of the different components of the EF for the study period. From this, local optimization can be applied to establish EF components with a certain level of compliance in the study period. Different options for graphic output and analysis of results are included in terms of graphs and tables in several formats. This methodology turned out to be a useful tool for the implementation of an uncertainty analysis within the scope of environmental flows in water management and allowed the simulation of the impacts of several water resource development scenarios in the study site.

Hydrodynamic Simulations Of Vacuum-Driven Storage Tanks

As a result of urbanization, stormwater runoff flow rates and volumes are significantly increased due to increasing impervious land cover and the decreased availability of depression storage. Storage tanks are the basic devices to efficiently control the flow rate in drainage systems during wet weather. Presented in the paper conception of vacuum-driven detention tanks allows to increase the storage capacity by usage of space above the free surface water elevation at the inlet channel. Partial vacuum storage makes possible to gain cost savings by reduction of both the horizontal area of the detention tank and necessary depth of foundations. Simulation model of vacuum-driven storage tank has been developed to estimate potential profits of its application in urban drainage system. Although SWMM5 has no direct options for vacuum tanks an existing functions (i.e. control rules) have been used to reflect its operation phases. Rainfall data used in simulations were recorded at raingage in Czestochowa during years 2010÷2012 with time interval of 10minutes. Simulation results gives overview to practical operation and maintenance cost (energy demand) of vacuum driven storage tanks depending of the ratio: vacuum-driven volume to total storage capacity. The following conclusion can be drawn from this investigations: vacuum-driven storage tanks are characterized by uncomplicated construction and control systems, thus can be applied in newly developed as well as in the existing urban drainage systems. the application of vacuum in underground detention facilities makes possible to increase of the storage capacity of existing reservoirs by usage the space above the maximum depth. Possible increase of storage capacity can achieve even a few dozen percent at relatively low investment costs. vacuum driven storage tanks can be included in existing simulation software (i.e. SWMM) using options intended for pumping stations (including control and action rules ).

Estudo estatístico sobre eventos de precipitação intensa no nordeste do Brasil

The Northeast of Brazil (NEB) shows high climate variability, ranging from semiarid regions to a rainy regions. According to the latest report of the Intergovernmental Panel on Climate Change, the NEB is highly susceptible to climate change, and also heavy rainfall events (HRE). However, few climatology studies about these episodes were performed, thus the objective main research is to compute the climatology and trend of the episodes number and the daily rainfall rate associated with HRE in the NEB and its climatologically homogeneous sub regions; relate them to the weak rainfall events and normal rainfall events. The daily rainfall data of the hydrometeorological network managed by the Agência Nacional de Águas, from 1972 to 2002. For selection of rainfall events used the technique of quantiles and the trend was identified using the Mann-Kendall test. The sub regions were obtained by cluster analysis, using as similarity measure the Euclidean distance and Ward agglomerative hierarchical method. The results show that the seasonality of the NEB is being intensified, i.e., the dry season is becoming drier and wet season getting wet. The El Niño and La Niña influence more on the amount of events regarding the intensity, but the sub-regions this influence is less noticeable. Using daily data reanalysis ERAInterim fields of anomalies of the composites of meteorological variables were calculated for the coast of the NEB, to characterize the synoptic environment. The Upper-level cyclonic vortex and the South atlantic convergene zone were identified as the main weather systems responsible for training of EPI on the coastland

Potencial de erosão da bacia do Rio Uberaba

Este trabalho teve por objetivo identificar qualitativamente as áreas suscetíveis à erosão laminar na bacia do Rio Uberaba, localizada em Uberaba -MG, apoiado no modelo matemático da Equação Universal de Perda de Solo (EUPS). Foram utilizadas cartas de: solos, uso e ocupação das terras, redes de drenagem, declividade e dados pluviográficos, utilizando-se de um Sistema de Informação Geográfica (SIG -IDRISI). A espacialização do potencial de erosão só foi possível a partir da estimativa da tolerância às perdas laminares para cada tipo de solo da bacia, e da profundidade dos solos, por entender que as perdas são mais significativas em solos mais rasos do que em solos muito profundos. Na análise dos resultados, verificou-se que 37% da área total da bacia do Rio Uberaba (905,24 km²) sofrem perdas de solos acima do limite de tolerância, sendo 12% em solos profundos e 25% em muito profundos, e a espacialização deste evento favorece a adoção de ações efetivas quanto à conservação dos solos da bacia.

Análise de uso e cobertura do solo no município de Parelhas/RN

The human interference in the semiarid region of Seridó Potiguar has promoted the increase of degraded areas. The economic dynamic that was established in the Seridó territory, especially after the fall of the trinomial cattle-cotton-mining in the 70s and 80s of the 20th century as pillars of the regional economy, resulted in an accelerated process of erosion of natural resources. The municipalities of the Seridó region have been spatially reordered by this new economic dynamic, marked by the growth of existing enterprises, and the development of new agricultural practices. One of the municipalities in the region that restructured its territorial space with the emergence of new agro-industrial activities was the town of Parelhas. With the demise of the trinomial cattle-cotton-mining in the 1980s, other productive activities were intensified from the 1990s, amongst them, pottery, responsible for the vegetal extraction for use as energy source. This recent economic and spatial restructuring in the region, reflected in the Parelhense municipal territory, required new productive ingredients responsible for the modification of past production relations that were based on cattle, cotton and mining. By that a process of exploring the environment was unleashed, especially the native vegetation, in an uncontrolled manner. In this context, the objective of this study was to survey and detect deforestation in the areas of Caatinga vegetation, used indiscriminately as energy supply for new agricultural practices, using remote sensing techniques based on the quantification of the Normalized Difference Vegetation Index / NDVI, Soil-Adjusted Vegetation Index / SAVI, surface temperature and rainfall data in the years 1990 and 2010. The results indicated that SAVI values above 0.2 in 1990 and 2010 represent the areas with the highest density of vegetation that occur exclusively along the major drainages in the town and areas of higher elevations. The areas between the ranges of values from 0.5 to 0.15 SAVI are areas with poor vegetation. On the other hand the highest values of temperature are distributed in the western and southeastern parts of the township, usually in places where the soil is exposed or there is sparse vegetation. The areas of bare soil decreased in extension in 2010 at 11, 6% when related to 1990, this was caused by a higher rainfall intensity in the first half of 2010, but no regeneration of vegetation occurred in some places in the western and southeastern areas of the municipality today, due to the extraction of firewood to fuel the furnaces of industries in town

Probabilidade de ocorrencia de periodos secos e chuvosos em Rio Claro - SP.

Using the daily rainfall data from 1961 through 1980 the frequency of dry and wet periods was determined. The results of the frequency distribution of dry and wet periods indicated that observed data fit very closely an equation of the type Y = aebn. -after English summary

Proposição e validações de equações para estimativa da erosividade de dois municípios do estado de São Paulo

The equations and extrapolation use to localities whose characteristics of soil and climate, even if partial, distinguish the town to which they were generated, still permeate in studies to estimate the rainfall erosivity (EI 30). This work has objective to propose and validate mathematical equations to estimate the rainfall erosivity of two cities of Sao Paulo State's. The adjusted to estimate obtaining and validate data of equations of erosivity (EI 30) according to values of coefficient of rain (Rc) were obtained from pluviographic and pluviometric rainfall data, respectively, using of distinct historical rainfall series. Mutiple comparisions test and confidence intervals were performed to compare absolute average of EI 30, pluviometric data (Pp), and Rc. The correlation between EI 30 and Rc was verified by of Pearson correlation coefficient. Test of the hypothesis of equality between population variance was used to compare the equations. Pluviometrics data of historical series rainfall data different than those that the models were generated were used to validate and to assess the performance of the equations, proposed of this study and compare them with another equation already consolidated in literature. The results show that for the conditions under which the study was conducted, the simple linear equations, shown to be the most appropriate to estimate the rainfall erosivity these two cities. According to the test of the hypothesis of equality variances between populations, the equations adjusted for each city differ statistically so that the rainfall erosivity of each city must be estimated by their respective equation.

As chuvas no estado do Tocantins: distribuição geográfica e gênese das variações rítmicas

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

As chuvas na bacia do Paraná: aspectos temporais, espaciais e rítmicos. -

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)