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.

Deriving global flood hazard maps of fluvial floods through a physical model cascade

Global flood hazard maps can be used in the assessment of flood risk in a number of different applications, including (re)insurance and large scale flood preparedness. Such global hazard maps can be generated using large scale physically based models of rainfall-runoff and river routing, when used in conjunction with a number of post-processing methods. In this study, the European Centre for Medium Range Weather Forecasts (ECMWF) land surface model is coupled to ERA-Interim reanalysis meteorological forcing data, and resultant runoff is passed to a river routing algorithm which simulates floodplains and flood flow across the global land area. The global hazard map is based on a 30 yr (1979–2010) simulation period. A Gumbel distribution is fitted to the annual maxima flows to derive a number of flood return periods. The return periods are calculated initially for a 25×25 km grid, which is then reprojected onto a 1×1 km grid to derive maps of higher resolution and estimate flooded fractional area for the individual 25×25 km cells. Several global and regional maps of flood return periods ranging from 2 to 500 yr are presented. The results compare reasonably to a benchmark data set of global flood hazard. The developed methodology can be applied to other datasets on a global or regional scale.

Identification and ranking of extraordinary rainfall events over northwest Italy: the role of Atlantic moisture

An objective identification and ranking of extraordinary rainfall events for Northwest Italy is established using time series of annual precipitation maxima for 1938–2002 at over 200 stations. Rainfall annual maxima are considered for five reference durations (1, 3, 6, 12, and 24 h). In a first step, a day is classified as an extraordinary rainfall day when a regional threshold calculated on the basis of a two-components extreme value distribution is exceeded for at least one of the stations. Second, a clustering procedure taking into account the different rainfall durations is applied to the identified 163 events. Third, a division into six clusters is chosen using Ward's distance criteria. It is found that two of these clusters include the seven strongest events as quantified from a newly developed measure of intensity which combines rainfall intensities and spatial extension. Two other clusters include the weakest 72% historical events. The obtained clusters are analyzed in terms of typical synoptic characteristics. The two top clusters are characterized by strong and persistent upper air troughs inducing not only moisture advection from the North Atlantic into the Western Mediterranean but also strong northward flow towards the southern Alpine ranges. Humidity transports from the North Atlantic are less important for the weaker clusters. We conclude that moisture advection from the North Atlantic plays a relevant role in the magnitude of the extraordinary events over Northwest Italy.

Estudo de Fractalidade e Evolução Dinâmica de Sistemas Complexos

In this work, the study of some complex systems is done with use of two distinct procedures. In the first part, we have studied the usage of Wavelet transform on analysis and characterization of (multi)fractal time series. We have test the reliability of Wavelet Transform Modulus Maxima method (WTMM) in respect to the multifractal formalism, trough the calculation of the singularity spectrum of time series whose fractality is well known a priori. Next, we have use the Wavelet Transform Modulus Maxima method to study the fractality of lungs crackles sounds, a biological time series. Since the crackles sounds are due to the opening of a pulmonary airway bronchi, bronchioles and alveoli which was initially closed, we can get information on the phenomenon of the airway opening cascade of the whole lung. Once this phenomenon is associated with the pulmonar tree architecture, which displays fractal geometry, the analysis and fractal characterization of this noise may provide us with important parameters for comparison between healthy lungs and those affected by disorders that affect the geometry of the tree lung, such as the obstructive and parenchymal degenerative diseases, which occurs, for example, in pulmonary emphysema. In the second part, we study a site percolation model for square lattices, where the percolating cluster grows governed by a control rule, corresponding to a method of automatic search. In this model of percolation, which have characteristics of self-organized criticality, the method does not use the automated search on Leaths algorithm. It uses the following control rule: pt+1 = pt + k(Rc − Rt), where p is the probability of percolation, k is a kinetic parameter where 0 < k < 1 and R is the fraction of percolating finite square lattices with side L, LxL. This rule provides a time series corresponding to the dynamical evolution of the system, in particular the likelihood of percolation p. We proceed an analysis of scaling of the signal obtained in this way. The model used here enables the study of the automatic search method used for site percolation in square lattices, evaluating the dynamics of their parameters when the system goes to the critical point. It shows that the scaling of , the time elapsed until the system reaches the critical point, and tcor, the time required for the system loses its correlations, are both inversely proportional to k, the kinetic parameter of the control rule. We verify yet that the system has two different time scales after: one in which the system shows noise of type 1 f , indicating to be strongly correlated. Another in which it shows white noise, indicating that the correlation is lost. For large intervals of time the dynamics of the system shows ergodicity

A method to assess annual average renewable groundwater reserves for large regions in Spain

This paper proposes a method for assessing the groundwater renewable reserves of large regions for an average year, based on the integration of the recession curves for their basins springs or the natural base flow of their rivers. In this method, the hydrodynamic volume (or renewable reserves), were estimated from the baseflow equation. It was assumed that the flow was the same as the natural recharge, and that the recession coefficients were derived by the hydrogeological parameters and geometrical characteristics of aquifers, and adjusted to fit the recession curves at gauging stations. The method was applied to all the aquifers of Spain, which have a total groundwater renewable reserve of 86,895 hm3 four times the mean annual recharge. However, the distribution of these reserves is very variable; 18.6% of the country aquifers contain 94.7% of the entire reserve.

Background effects in the induced polarization method of geophysical exploration : annual progress report for 1956-57 /

"June 28, 1957."

The apostolical method of preaching the gospel. A sermon, preached at the annual meeting of the Northern Education Society, August 31st, 1814.

Mode of access: Internet.

Estimation of average annual daily traffic (AADT) using geographically weighted regression (GWR) method and geographic information system (GIS)

Annual average daily traffic (AADT) is important information for many transportation planning, design, operation, and maintenance activities, as well as for the allocation of highway funds. Many studies have attempted AADT estimation using factor approach, regression analysis, time series, and artificial neural networks. However, these methods are unable to account for spatially variable influence of independent variables on the dependent variable even though it is well known that to many transportation problems, including AADT estimation, spatial context is important. ^ In this study, applications of geographically weighted regression (GWR) methods to estimating AADT were investigated. The GWR based methods considered the influence of correlations among the variables over space and the spatially non-stationarity of the variables. A GWR model allows different relationships between the dependent and independent variables to exist at different points in space. In other words, model parameters vary from location to location and the locally linear regression parameters at a point are affected more by observations near that point than observations further away. ^ The study area was Broward County, Florida. Broward County lies on the Atlantic coast between Palm Beach and Miami-Dade counties. In this study, a total of 67 variables were considered as potential AADT predictors, and six variables (lanes, speed, regional accessibility, direct access, density of roadway length, and density of seasonal household) were selected to develop the models. ^ To investigate the predictive powers of various AADT predictors over the space, the statistics including local r-square, local parameter estimates, and local errors were examined and mapped. The local variations in relationships among parameters were investigated, measured, and mapped to assess the usefulness of GWR methods. ^ The results indicated that the GWR models were able to better explain the variation in the data and to predict AADT with smaller errors than the ordinary linear regression models for the same dataset. Additionally, GWR was able to model the spatial non-stationarity in the data, i.e., the spatially varying relationship between AADT and predictors, which cannot be modeled in ordinary linear regression. ^

Development of a Rational Characterization Method for Iowa Fly Ash, HR-286, Annual Progress Report 1, 1986

The following report summarizes research activities on the project for the period December 1, 1985 through November 31, 1986. Research efforts for the first year have proceeded basically as outlined in the project proposal.

Development of a Rational Characterization Method for Iowa Fly Ash, HR-286, Annual Progress Report 2, 1987

The following report summarizes research activities on the project for the period December 1, 1986 to November 30, 1987. Research efforts for the second year deviated slightly from those described in the project proposal. By the end of the second year of testing, it was possible to begin evaluating how power plant operating conditions influenced the chemical and physical properties of fly ash obtained from one of the monitored power plants (Ottumwa Generating Station, OGS). Hence, several of the tasks initially assigned to the third year of the project (specifically tasks D, E, and F) were initiated during the second year of the project. Manpower constraints were balanced by delaying full scale implementation of the quantitative X-ray diffraction and differential thermal analysis tasks until the beginning of the third year of the project. Such changes should have little bearing on the outcome of the overall project.

Hybrid Method to Assess Sensitive Process Interruption Costs Due to Faults in Electric Power Distribution Networks

This paper shows a new hybrid method for risk assessment regarding interruptions in sensitive processes due to faults in electric power distribution systems. This method determines indices related to long duration interruptions and short duration voltage variations (SDVV), such as voltage sags and swells in each customer supplied by the distribution network. Frequency of such occurrences and their impact on customer processes are determined for each bus and classified according to their corresponding magnitude and duration. The method is based on information regarding network configuration, system parameters and protective devices. It randomly generates a number of fault scenarios in order to assess risk areas regarding long duration interruptions and voltage sags and swells in an especially inventive way, including frequency of events according to their magnitude and duration. Based on sensitivity curves, the method determines frequency indices regarding disruption in customer processes that represent equipment malfunction and possible process interruptions due to voltage sags and swells. Such approach allows for the assessment of the annual costs associated with each one of the evaluated power quality indices.