The use of Meteosat satellite data for spatial rainfall estimations and hydrological simulations

Satellite information, in combination with conventional point source measurements, can be a valuable source of information. This thesis is devoted to the spatial estimation of areal rainfall over a region using both the measurements from a dense and sparse network of rain-gauges and images from the meteorological satellites. A primary concern is to study the effects of such satellite assisted rainfall estimates on the performance of rainfall-runoff models. Low-cost image processing systems and peripherals are used to process and manipulate the data. Both secondary as well as primary satellite images were used for analysis. The secondary data was obtained from the in-house satellite receiver and the primary data was obtained from an outside source. Ground truth data was obtained from the local Water Authority. A number of algorithms are presented that combine the satellite and conventional data sources to produce areal rainfall estimates and the results are compared with some of the more traditional methodologies. The results indicate that the satellite cloud information is valuable in the assessment of the spatial distribution of areal rainfall, for both half-hourly as well as daily estimates of rainfall. It is also demonstrated how the performance of the simple multiple regression rainfall-runoff model is improved when satellite cloud information is used as a separate input in addition to rainfall estimates from conventional means. The use of low-cost equipment, from image processing systems to satellite imagery, makes it possible for developing countries to introduce such systems in areas where the benefits are greatest.

Assessing the uncertainty associated with intermittent rainfall fields

[1] In many practical situations where spatial rainfall estimates are needed, rainfall occurs as a spatially intermittent phenomenon. An efficient geostatistical method for rainfall estimation in the case of intermittency has previously been published and comprises the estimation of two independent components: a binary random function for modeling the intermittency and a continuous random function that models the rainfall inside the rainy areas. The final rainfall estimates are obtained as the product of the estimates of these two random functions. However the published approach does not contain a method for estimation of uncertainties. The contribution of this paper is the presentation of the indicator maximum likelihood estimator from which the local conditional distribution of the rainfall value at any location may be derived using an ensemble approach. From the conditional distribution, representations of uncertainty such as the estimation variance and confidence intervals can be obtained. An approximation to the variance can be calculated more simply by assuming rainfall intensity is independent of location within the rainy area. The methodology has been validated using simulated and real rainfall data sets. The results of these case studies show good agreement between predicted uncertainties and measured errors obtained from the validation data.

Statistical modelling of extreme rainfall in Madeira Island

Extreme rainfall events have triggered a significant number of flash floods in Madeira Island along its past and recent history. Madeira is a volcanic island where the spatial rainfall distribution is strongly affected by its rugged topography. In this thesis, annual maximum of daily rainfall data from 25 rain gauge stations located in Madeira Island were modelled by the generalised extreme value distribution. Also, the hypothesis of a Gumbel distribution was tested by two methods and the existence of a linear trend in both distributions parameters was analysed. Estimates for the 50– and 100–year return levels were also obtained. Still in an univariate context, the assumption that a distribution function belongs to the domain of attraction of an extreme value distribution for monthly maximum rainfall data was tested for the rainy season. The available data was then analysed in order to find the most suitable domain of attraction for the sampled distribution. In a different approach, a search for thresholds was also performed for daily rainfall values through a graphical analysis. In a multivariate context, a study was made on the dependence between extreme rainfall values from the considered stations based on Kendall’s τ measure. This study suggests the influence of factors such as altitude, slope orientation, distance between stations and their proximity of the sea on the spatial distribution of extreme rainfall. Groups of three pairwise associated stations were also obtained and an adjustment was made to a family of extreme value copulas involving the Marshall–Olkin family, whose parameters can be written as a function of Kendall’s τ association measures of the obtained pairs.

Spatial intra-annual variability of precipitation based on geostatistics. a case study for the Paraiba do Sul Basin, southeastern Brazil

The emphasis in this research is to evaluate the spatial distribution of the precipitation using a geostatistics approach. Seasonal time scales records considering DJF, MAM, JJA e SON periods performed the analysis. Procedures to evaluate the variogram selection and to produce kriging maps were performed in a GIS environment (ArcGIS®). The results showed that kriging method was very suitable to detect both large changes in the whole area as those local small and subtle changes. Kriging demonstrated be a powerful statistical interpolation method that might be very useful in regions with great complexity in climatology and geomorphology.

A space-time hybrid hourly rainfall model for derived flood frequency analysis

For derived flood frequency analysis based on hydrological modelling long continuous precipitation time series with high temporal resolution are needed. Often, the observation network with recording rainfall gauges is poor, especially regarding the limited length of the available rainfall time series. Stochastic precipitation synthesis is a good alternative either to extend or to regionalise rainfall series to provide adequate input for long-term rainfall-runoff modelling with subsequent estimation of design floods. Here, a new two step procedure for stochastic synthesis of continuous hourly space-time rainfall is proposed and tested for the extension of short observed precipitation time series. First, a single-site alternating renewal model is presented to simulate independent hourly precipitation time series for several locations. The alternating renewal model describes wet spell durations, dry spell durations and wet spell intensities using univariate frequency distributions separately for two seasons. The dependence between wet spell intensity and duration is accounted for by 2-copulas. For disaggregation of the wet spells into hourly intensities a predefined profile is used. In the second step a multi-site resampling procedure is applied on the synthetic point rainfall event series to reproduce the spatial dependence structure of rainfall. Resampling is carried out successively on all synthetic event series using simulated annealing with an objective function considering three bivariate spatial rainfall characteristics. In a case study synthetic precipitation is generated for some locations with short observation records in two mesoscale catchments of the Bode river basin located in northern Germany. The synthetic rainfall data are then applied for derived flood frequency analysis using the hydrological model HEC-HMS. The results show good performance in reproducing average and extreme rainfall characteristics as well as in reproducing observed flood frequencies. The presented model has the potential to be used for ungauged locations through regionalisation of the model parameters.

Appraising and selecting strategies to combat and mitigate desertification based on stakeholder knowledge and globalbest practices in Cape Verde Archipelago

Cape Verde is considered part of Sahelian Africa, where drought and desertification are common occurrences. The main activity of the rural population is rain-fed agriculture, which over time has been increasingly challenged by high temporal and spatial rainfall variability, lack of inputs, limited land area, fragmentation of land, steep slopes, pests, lack of mechanization and loss of top soil by water erosion. Human activities, largely through poor farming practices and deforestation (Gomez, 1989) have accelerated natural erosion processes, shifting the balance between soil erosion and soil formation (Norton, 1987). According to previous studies, vegetation cover is one of the most important factors in controlling soil loss (Cyr et al., 1995; Hupy, 2004; Zhang et al., 2004; Zhou et al., 2006). For this reason, reforestation is a touchstone of the Cape Verdean policy to combat desertification. After Independence in 1975, the Cape Verde government had pressing and closely entangled environmental and socio-economic issues to address, as long-term desertification had resulted in a lack of soil cover, severe soil erosion and a scarcity of water resources and fuel wood. Across the archipelago, desertification was resulting from a variety of processes including poor farming practices, soil erosion by water and wind, soil and water salinity in coastal areas due to over pumping and seawater intrusion, drought and unplanned urbanization (DGA-MAAP, 2004). All these issues directly affected socio-economic vulnerability in rural areas, where about 70% of people depended directly or indirectly on agriculture in 1975. By becoming part of the Inter- State Committee for the Fight against Drought in the Sahel in 1975, the government of Cape Verde gained structured support to address these issues more efficiently. Presentday policies and strategies were defined on the basis of rational use of resources and human efforts and were incorporated into three subsequent national plans: the National Action Plan for Development (NDP) (1982–1986), the NDP (1986–1990) and the NDP (1991–1995) (Carvalho

Characteristics of 2D convective structures in Catalonia (NE Spain): an analysis using radar data and GIS

Flood simulation studies use spatial-temporal rainfall data input into distributed hydrological models. A correct description of rainfall in space and in time contributes to improvements on hydrological modelling and design. This work is focused on the analysis of 2-D convective structures (rain cells), whose contribution is especially significant in most flood events. The objective of this paper is to provide statistical descriptors and distribution functions for convective structure characteristics of precipitation systems producing floods in Catalonia (NE Spain). To achieve this purpose heavy rainfall events recorded between 1996 and 2000 have been analysed. By means of weather radar, and applying 2-D radar algorithms a distinction between convective and stratiform precipitation is made. These data are introduced and analyzed with a GIS. In a first step different groups of connected pixels with convective precipitation are identified. Only convective structures with an area greater than 32 km2 are selected. Then, geometric characteristics (area, perimeter, orientation and dimensions of the ellipse), and rainfall statistics (maximum, mean, minimum, range, standard deviation, and sum) of these structures are obtained and stored in a database. Finally, descriptive statistics for selected characteristics are calculated and statistical distributions are fitted to the observed frequency distributions. Statistical analyses reveal that the Generalized Pareto distribution for the area and the Generalized Extreme Value distribution for the perimeter, dimensions, orientation and mean areal precipitation are the statistical distributions that best fit the observed ones of these parameters. The statistical descriptors and the probability distribution functions obtained are of direct use as an input in spatial rainfall generators.

Seasonality of westerly moisture transport in the East Asian summer monsoon and its implications for interpreting precipitation δ18O

East Asian summer monsoon (EASM) rainfall impacts the world's most populous regions. Accurate EASM rainfall prediction necessitates robust paleoclimate reconstructions from proxy data and quantitative linkage to modern climatic conditions. Many precisely dated oxygen isotope records from Chinese stalagmites have been interpreted as directly reflecting past EASM rainfall amount variability, but recent research suggests that such records instead integrate multiple hydroclimatic processes. Using a Lagrangian precipitation moisture source diagnostic, we demonstrate that EASM rainfall is primarily derived from the Indian Ocean. Conversely, Pacific Ocean moisture export peaks during winter, and the moisture uptake area does not differ significantly between summer and winter and is thus a minor contributor to monsoonal precipitation. Our results are substantiated by an accurate reproduction of summer and winter spatial rainfall distributions across China. We also correlate modern EASM rainfall oxygen isotope ratios with instrumental rainfall amount and our moisture source data. This analysis reveals that the strength of the source effect is geographically variable, and differences in atmospheric moisture transport may significantly impact the isotopic signature of EASM rainfall at the Hulu, Dongge, and Wanxiang Cave sites. These results improve our ability to isolate the rainfall amount signal in paleomonsoon reconstructions and indicate that precipitation across central and eastern China will directly respond to variability in Indian Ocean moisture supply.

Precipitação sazonal sobre a Amazônia oriental no período chuvoso: observações e simulações regionais com o RegCM3

O presente trabalho apresenta uma contribuição aos estudos de modelagem climática com ênfase na variabilidade pluviométrica sazonal da Amazônia oriental, durante as estações de verão e outono (DJF e MAM). Baseado nos resultados das simulações regionais do RegCM3 para um período de 26 anos (1982/83 a 2007/08) e usando domínio em alta resolução espacial (30 Km) e dois diferentes esquemas de convecção (Grell e MIT), foi investigado o desempenho do modelo em simular a distribuição regional de precipitação sazonal na Amazônia oriental, com referência a um novo conjunto de dados observacional compilado com informações de uma ampla rede integrada de estações pluviométricas. As análises quantitativas evidenciaram que o RegCM3 apresenta erros sistemáticos, sobretudo aqueles relacionados com viés seco no Amapá e norte/nordeste do Pará usando ambos os esquemas Grell e MIT, os quais apontam que o modelo não reproduz as características da ZCIT sobre o Atlântico equatorial. As simulações usando MIT, também apresentaram viés úmido no sudoeste/sul/sudeste do Pará e norte do Tocantins. Além disso, através da técnica de composições, também foi investigado o desempenho do RegCM3 em reproduzir os padrões espaciais anômalos de precipitação sazonal em associação aos episódios ENOS, e as fases do gradiente térmico sobre o Atlântico intertropical. Os resultados demonstraram que o modelo conseguiu representar realisticamente bem o padrão espacial das anomalias pluviométricas acima (abaixo) do normal em grande parte da Amazônia oriental, durante os conhecidos cenários favoráveis, i.e., condições de La Niña e gradiente de aTSM para o Atlântico sul (desfavoráveis, i.e., El Niño e gradiente de aTSM para o Atlântico norte.

Data assimilation for precipitation nowcasting using Bayesian inference

This work introduces a new variational Bayes data assimilation method for the stochastic estimation of precipitation dynamics using radar observations for short term probabilistic forecasting (nowcasting). A previously developed spatial rainfall model based on the decomposition of the observed precipitation field using a basis function expansion captures the precipitation intensity from radar images as a set of ‘rain cells’. The prior distributions for the basis function parameters are carefully chosen to have a conjugate structure for the precipitation field model to allow a novel variational Bayes method to be applied to estimate the posterior distributions in closed form, based on solving an optimisation problem, in a spirit similar to 3D VAR analysis, but seeking approximations to the posterior distribution rather than simply the most probable state. A hierarchical Kalman filter is used to estimate the advection field based on the assimilated precipitation fields at two times. The model is applied to tracking precipitation dynamics in a realistic setting, using UK Met Office radar data from both a summer convective event and a winter frontal event. The performance of the model is assessed both traditionally and using probabilistic measures of fit based on ROC curves. The model is shown to provide very good assimilation characteristics, and promising forecast skill. Improvements to the forecasting scheme are discussed

Estimation of Convective Extreme Events over Germany based on Large Scale Parameters

Strong convective events can produce extreme precipitation, hail, lightning or gusts, potentially inducing severe socio-economic impacts. These events have a relatively small spatial extension and, in most cases, a short lifetime. In this study, a model is developed for estimating convective extreme events based on large scale conditions. It is shown that strong convective events can be characterized by a Weibull distribution of radar-based rainfall with a low shape and high scale parameter value. A radius of 90km around a station reporting a convective situation turned out to be suitable. A methodology is developed to estimate the Weibull parameters and thus the occurrence probability of convective events from large scale atmospheric instability and enhanced near-surface humidity, which are usually found on a larger scale than the convective event itself. Here, the probability for the occurrence of extreme convective events is estimated from the KO-index indicating the stability, and relative humidity at 1000hPa. Both variables are computed from ERA-Interim reanalysis. In a first version of the methodology, these two variables are applied to estimate the spatial rainfall distribution and to estimate the occurrence of a convective event. The developed method shows significant skill in estimating the occurrence of convective events as observed at synoptic stations, lightning measurements, and severe weather reports. In order to take frontal influences into account, a scheme for the detection of atmospheric fronts is implemented. While generally higher instability is found in the vicinity of fronts, the skill of this approach is largely unchanged. Additional improvements were achieved by a bias-correction and the use of ERA-Interim precipitation. The resulting estimation method is applied to the ERA-Interim period (1979-2014) to establish a ranking of estimated convective extreme events. Two strong estimated events that reveal a frontal influence are analysed in detail. As a second application, the method is applied to GCM-based decadal predictions in the period 1979-2014, which were initialized every year. It is shown that decadal predictive skill for convective event frequencies over Germany is found for the first 3-4 years after the initialization.

Maize yield and rainfall on different spatial and temporal scales in Southern Brazil

This study aimed to establish relationships between maize yield and rainfall on different temporal and spatial scales, in order to provide a basis for crop monitoring and modelling. A 16-year series of maize yield and daily rainfall from 11 municipalities and micro-regions of Rio Grande do Sul State was used. Correlation and regression analyses were used to determine associations between crop yield and rainfall for the entire crop cycle, from tasseling to 30 days after, and from 5 days before tasseling to 40 days after. Close relationships between maize yield and rainfall were found, particularly during the reproductive period (45-day period comprising the flowering and grain filling). Relationships were closer on a regional scale than at smaller scales. Implications of the crop-rainfall relationships for crop modelling are discussed.

Spatial variability of the rainfall erosive potential in the State of Mato Grosso do Sul, Brazil

Information about rainfall erosivity is important during soil and water conservation planning. Thus, the spatial variability of rainfall erosivity of the state Mato Grosso do Sul was analyzed using ordinary kriging interpolation. For this, three pluviograph stations were used to obtain the regression equations between the erosivity index and the rainfall coefficient EI30. The equations obtained were applied to 109 pluviometric stations, resulting in EI30 values. These values were analyzed from geostatistical technique, which can be divided into: descriptive statistics, adjust to semivariogram, cross-validation process and implementation of ordinary kriging to generate the erosivity map.Highest erosivity values were found in central and northeast regions of the State, while the lowest values were observed in the southern region. In addition, high annual precipitation values not necessarily produce higher erosivity values.

Maize yield and rainfall on different spatial and temporal scales in Southern Brazil

This study aimed to establish relationships between maize yield and rainfall on different temporal and spatial scales, in order to provide a basis for crop monitoring and modelling. A 16-year series of maize yield and daily rainfall from 11 municipalities and micro-regions of Rio Grande do Sul State was used. Correlation and regression analyses were used to determine associations between crop yield and rainfall for the entire crop cycle, from tasseling to 30 days after, and from 5 days before tasseling to 40 days after. Close relationships between maize yield and rainfall were found, particularly during the reproductive period (45-day period comprising the flowering and grain filling). Relationships were closer on a regional scale than at smaller scales. Implications of the crop-rainfall relationships for crop modelling are discussed.

Spatial disaggregation and intensity correction of TRMM-based rainfall time series for hydrological applications in dryland catchments

A novel approach is presented for combining spatial and temporal detail from newly available TRMM-based data sets to derive hourly rainfall intensities at 1-km spatial resolution for hydrological modelling applications. Time series of rainfall intensities derived from 3-hourly 0.25° TRMM 3B42 data are merged with a 1-km gridded rainfall climatology based on TRMM 2B31 data to account for the sub-grid spatial distribution of rainfall intensities within coarse-scale 0.25° grid cells. The method is implemented for two dryland catchments in Tunisia and Senegal, and validated against gauge data. The outcomes of the validation show that the spatially disaggregated and intensity corrected TRMM time series more closely approximate ground-based measurements than non-corrected data. The method introduced here enables the generation of rainfall intensity time series with realistic temporal and spatial detail for dynamic modelling of runoff and infiltration processes that are especially important to water resource management in arid regions.