Development of a wind gust model to estimate gust speeds and their return periods

Spatially dense observations of gust speeds are necessary for various applications, but their availability is limited in space and time. This work presents an approach to help to overcome this problem. The main objective is the generation of synthetic wind gust velocities. With this aim, theoretical wind and gust distributions are estimated from 10 yr of hourly observations collected at 123 synoptic weather stations provided by the German Weather Service. As pre-processing, an exposure correction is applied on measurements of the mean wind velocity to reduce the influence of local urban and topographic effects. The wind gust model is built as a transfer function between distribution parameters of wind and gust velocities. The aim of this procedure is to estimate the parameters of gusts at stations where only wind speed data is available. These parameters can be used to generate synthetic gusts, which can improve the accuracy of return periods at test sites with a lack of observations. The second objective is to determine return periods much longer than the nominal length of the original time series by considering extreme value statistics. Estimates for both local maximum return periods and average return periods for single historical events are provided. The comparison of maximum and average return periods shows that even storms with short average return periods may lead to local wind gusts with return periods of several decades. Despite uncertainties caused by the short length of the observational records, the method leads to consistent results, enabling a wide range of possible applications.

Return periods of losses associated with European windstorm series in a changing climate

Possible future changes of clustering and return periods (RPs) of European storm series with high potential losses are quantified. Historical storm series are identified using 40 winters of reanalysis. Time series of top events (1, 2 or 5 year return levels (RLs)) are used to assess RPs of storm series both empirically and theoretically. Additionally, 800 winters of general circulation model simulations for present (1960–2000) and future (2060–2100) climate conditions are investigated. Clustering is identified for most countries, and estimated RPs are similar for reanalysis and present day simulations. Future changes of RPs are estimated for fixed RLs and fixed loss index thresholds. For the former, shorter RPs are found for Western Europe, but changes are small and spatially heterogeneous. For the latter, which combines the effects of clustering and event ranking shifts, shorter RPs are found everywhere except for Mediterranean countries. These changes are generally not statistically significant between recent and future climate. However, the RPs for the fixed loss index approach are mostly beyond the range of pre-industrial natural climate variability. This is not true for fixed RLs. The quantification of losses associated with storm series permits a more adequate windstorm risk assessment in a changing climate.

On the variability of return periods of European winter precipitation extremes over the last three centuries

We investigate the changes of extreme European winter (December-February) precipitation back to 1700 and show for various European regions that return periods of extremely wet and dry winters are subject to significant changes both before and after the onset of anthropogenic influences. Generally, winter precipitation has become more extreme. We also examine the spatial pattern of the changes of the extremes covering the last 300 years where data quality is sufficient. Over central and Eastern Europe dry winters occurred more frequently during the 18th and the second part of the 19th century relative to 1951–2000. Dry winters were less frequent during both the 18th and 19th century over the British Isles and the Mediterranean. Wet winters have been less abundant during the last three centuries compared to 1951–2000 except during the early 18th century in central Europe. Although winter precipitation extremes are affected by climate change, no obvious connection of these changes was found to solar, volcanic or anthropogenic forcing. However, physically meaningful interpretation with atmospheric circulation changes was possible.

A bivariate return period based on copulas for hydrologic dam design: accounting for reservoir routing in risk estimation

A multivariate analysis on flood variables is needed to design some hydraulic structures like dams, as the complexity of the routing process in a reservoir requires a representation of the full hydrograph. In this work, a bivariate copula model was used to obtain the bivariate joint distribution of flood peak and volume, in order to know the probability of occurrence of a given inflow hydrograph. However, the risk of dam overtopping is given by the maximum water elevation reached during the routing process, which depends on the hydrograph variables, the reservoir volume and the spillway crest length. Consequently, an additional bivariate return period, the so-called routed return period, was defined in terms of risk of dam overtopping based on this maximum water elevation obtained after routing the inflow hydrographs. The theoretical return periods, which give the probability of occurrence of a hydrograph prior to accounting for the reservoir routing, were compared with the routed return period, as in both cases hydrographs with the same probability will draw a curve in the peak-volume space. The procedure was applied to the case study of the Santillana reservoir in Spain. Different reservoir volumes and spillway lengths were considered to investigate the influence of the dam and reservoir characteristics on the results. The methodology improves the estimation of the Design Flood Hydrograph and can be applied to assess the risk of dam overtopping

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Estimating present day extreme water level exceedance probabilities around the coastline of Australia : tides, extra-tropical storm surges and mean sea level

The occurrence of extreme water levels along low-lying, highly populated and/or developed coastlines can lead to considerable loss of life and billions of dollars of damage to coastal infrastructure. Therefore it is vitally important that the exceedance probabilities of extreme water levels are accurately evaluated to inform risk-based flood management, engineering and future land-use planning. This ensures the risk of catastrophic structural failures due to under-design or expensive wastes due to over-design are minimised. This paper estimates for the first time present day extreme water level exceedence probabilities around the whole coastline of Australia. A high-resolution depth averaged hydrodynamic model has been configured for the Australian continental shelf region and has been forced with tidal levels from a global tidal model and meteorological fields from a global reanalysis to generate a 61-year hindcast of water levels. Output from this model has been successfully validated against measurements from 30 tide gauge sites. At each numeric coastal grid point, extreme value distributions have been fitted to the derived time series of annual maxima and the several largest water levels each year to estimate exceedence probabilities. This provides a reliable estimate of water level probabilities around southern Australia; a region mainly impacted by extra-tropical cyclones. However, as the meteorological forcing used only weakly includes the effects of tropical cyclones, extreme water level probabilities are underestimated around the western, northern and north-eastern Australian coastline. In a companion paper we build on the work presented here and more accurately include tropical cyclone-induced surges in the estimation of extreme water level. The multi-decadal hindcast generated here has been used primarily to estimate extreme water level exceedance probabilities but could be used more widely in the future for a variety of other research and practical applications.

Estimating present day extreme water level exceedance probabilities around the coastline of Australia : tropical cyclone-induced storm surges

The incidence of major storm surges in the last decade have dramatically emphasized the immense destructive capabilities of extreme water level events, particularly when driven by severe tropical cyclones. Given this risk, it is vitally important that the exceedance probabilities of extreme water levels are accurately evaluated to inform risk-based flood and erosion management, engineering and for future land-use planning and to ensure the risk of catastrophic structural failures due to under-design or expensive wastes due to over-design are minimised. Australia has a long history of coastal flooding from tropical cyclones. Using a novel integration of two modeling techniques, this paper provides the first estimates of present day extreme water level exceedance probabilities around the whole coastline of Australia, and the first estimates that combine the influence of astronomical tides, storm surges generated by both extra-tropical and tropical cyclones, and seasonal and inter-annual variations in mean sea level. Initially, an analysis of tide gauge records has been used to assess the characteristics of tropical cyclone-induced surges around Australia. However, given the dearth (temporal and spatial) of information around much of the coastline, and therefore the inability of these gauge records to adequately describe the regional climatology, an observationally based stochastic tropical cyclone model has been developed to synthetically extend the tropical cyclone record to 10,000 years. Wind and pressure fields derived for these synthetically generated events have then been used to drive a hydrodynamic model of the Australian continental shelf region with annual maximum water levels extracted to estimate exceedance probabilities around the coastline. To validate this methodology, selected historic storm surge events have been simulated and resultant storm surges compared with gauge records. Tropical cyclone induced exceedance probabilities have been combined with estimates derived from a 61-year water level hindcast described in a companion paper to give a single estimate of present day extreme water level probabilities around the whole coastline of Australia. Results of this work are freely available to coastal engineers, managers and researchers via a web-based tool (www.sealevelrise.info). The described methodology could be applied to other regions of the world, like the US east coast, that are subject to both extra-tropical and tropical cyclones.

Estimating present day extreme total water level exceedance probabilities around the Australian coastline

The occurrence of extreme water level events along low-lying, highly populated and/or developed coastlines can lead to devastating impacts on coastal infrastructure. Therefore it is very important that the probabilities of extreme water levels are accurately evaluated to inform flood and coastal management and for future planning. The aim of this study was to provide estimates of present day extreme total water level exceedance probabilities around the whole coastline of Australia, arising from combinations of mean sea level, astronomical tide and storm surges generated by both extra-tropical and tropical storms, but exclusive of surface gravity waves. The study has been undertaken in two main stages. In the first stage, a high-resolution (~10 km along the coast) hydrodynamic depth averaged model has been configured for the whole coastline of Australia using the Danish Hydraulics Institute’s Mike21 modelling suite of tools. The model has been forced with astronomical tidal levels, derived from the TPX07.2 global tidal model, and meteorological fields, from the US National Center for Environmental Prediction’s global reanalysis, to generate a 61-year (1949 to 2009) hindcast of water levels. This model output has been validated against measurements from 30 tide gauge sites around Australia with long records. At each of the model grid points located around the coast, time series of annual maxima and the several highest water levels for each year were derived from the multi-decadal water level hindcast and have been fitted to extreme value distributions to estimate exceedance probabilities. Stage 1 provided a reliable estimate of the present day total water level exceedance probabilities around southern Australia, which is mainly impacted by extra-tropical storms. However, as the meteorological fields used to force the hydrodynamic model only weakly include the effects of tropical cyclones the resultant water levels exceedance probabilities were underestimated around western, northern and north-eastern Australia at higher return periods. Even if the resolution of the meteorological forcing was adequate to represent tropical cyclone-induced surges, multi-decadal periods yielded insufficient instances of tropical cyclones to enable the use of traditional extreme value extrapolation techniques. Therefore, in the second stage of the study, a statistical model of tropical cyclone tracks and central pressures was developed using histroic observations. This model was then used to generate synthetic events that represented 10,000 years of cyclone activity for the Australia region, with characteristics based on the observed tropical cyclones over the last ~40 years. Wind and pressure fields, derived from these synthetic events using analytical profile models, were used to drive the hydrodynamic model to predict the associated storm surge response. A random time period was chosen, during the tropical cyclone season, and astronomical tidal forcing for this period was included to account for non-linear interactions between the tidal and surge components. For each model grid point around the coast, annual maximum total levels for these synthetic events were calculated and these were used to estimate exceedance probabilities. The exceedance probabilities from stages 1 and 2 were then combined to provide a single estimate of present day extreme water level probabilities around the whole coastline of Australia.

Idealised simulation of downburst winds

Convectively driven downburst winds pose a threat to structures and communities in many regions of Australia not subject to tropical cyclones. Extreme value analysis shows that for return periods of interest to engineering design these events produce higher gust wind speeds than synoptic scale windstorms. Despite this, comparatively little is known of the near ground wind structure of these potentially hazardous windstorms. With this in mind, a series of idealised three-dimensional numerical simulations were undertaken to investigate convective storm wind fields. A dry, non-hydrostatic, sub-cloud model with parameterisation of the microphysics was used. Simulations were run with a uniform 20 m horizontal grid resolution and a variable vertical resolution increasing from 1 m. A systematic grid resolution study showed further refinement did not alter the morphological structure of the outflow. Simulations were performed for stationary downbursts in a quiescent air field, stationary downbursts embedded within environmental boundary layer winds, and also translating downbursts embedded within environmental boundary layer winds.

Loading of a very tall building in a simulated downburst wind field

Thunderstorm downbursts are important for wind engineers as they have been shown to produce the design wind speeds for mid to high return periods in many regions of Australia [1]. In structural design codes (e.g. AS/NZS1170.02-02) an atmospheric boundary layer (ABL) is assumed, and a vertical profile is interpolated from recorded 10 m wind speeds. The ABL assumption is however inaccurate when considering the complex structure of a thunderstorm outflow, and its effects on engineered structures. Several researchers have shown that the downburst, close to its point of divergence is better represented by an impinging wall jet profile than the traditional ABL. Physical modelling is the generally accepted approach to estimate wind loads on structures and it is therefore important to physically model the thunderstorm downburst so that its effects on engineered structures may be studied. An advancement on the simple impinging jet theory, addressed here is the addition of a pulsing mechanism to the jet which allows not only the divergent characteristics of a downburst to be produced, but also it allows the associated leading ring vortex to be developed. The ring vortex modelling is considered very important for structural design as it is within the horizontal vortex that the largest velocities occur [2]. This paper discusses the flow field produced by a pulsed wall jet, and also discusses the induced pressures that this type of flow has on a scaled tall building.

Delineation of seismic source zones based on seismicity parameters and probabilistic evaluation of seismic hazard using logic tree approach

The delineation of seismic source zones plays an important role in the evaluation of seismic hazard. In most of the studies the seismic source delineation is done based on geological features. In the present study, an attempt has been made to delineate seismic source zones in the study area (south India) based on the seismicity parameters. Seismicity parameters and the maximum probable earthquake for these source zones were evaluated and were used in the hazard evaluation. The probabilistic evaluation of seismic hazard for south India was carried out using a logic tree approach. Two different types of seismic sources, linear and areal, were considered in the present study to model the seismic sources in the region more precisely. In order to properly account for the attenuation characteristics of the region, three different attenuation relations were used with different weightage factors. Seismic hazard evaluation was done for the probability of exceedance (PE) of 10% and 2% in 50 years. The spatial variation of rock level peak horizontal acceleration (PHA) and spectral acceleration (Sa) values corresponding to return periods of 475 and 2500 years for the entire study area are presented in this work. The peak ground acceleration (PGA) values at ground surface level were estimated based on different NEHRP site classes by considering local site effects.

Avaliação da vulnerabilidade sísmica de edifícios escolares

A grande motivação para este trabalho resultou na avaliação da vulnerabilidade sísmica das escolas básicas e secundárias que fazem parte integrante do parque escolar de Portugal Continental. Apesar de até ao momento apenas se terem estudado a vulnerabilidade de algumas escolas em algumas zonas particulares do nosso país, para este trabalho de investigação tivemos uma ambição muito maior e obviamente fizemos um esforço muito maior. Estabelecemos assim como meta, a avaliação de todo o parque escolar. Embora todo o parque escolar possua na ordem das três centenas de escolas em todo o território nacional e sendo este projeto de reabilitação, um projeto com a duração de 2007 a 2015. Em 2011, por motivos da crise económica, todo o projeto congelou sendo reabilitadas até à data cerca de apenas um terço das escolas. Esta quantidade foi o número de escolas que avaliamos em todo o país. As escolas, sendo edifícios públicos com uma importância fundamental, tanto pela elevada concentração de jovens, como pela função essencial como centros de aprendizagem para as gerações vindouras, como também pela ameaça que representam na eventualidade dum cenário sísmico pela enorme densidade de utilizadores, e pela vantagem de nesse cenário de catástrofe a importância estrutural ser superior em relação à maior parte dos edifícios correntes, devidamente demonstrado pelos argumentos enumerados, consequentemente as escolas podem servir como instalações de proteção civil perante uma catástrofe sísmica para apoio das populações circundantes afetadas. Portanto para cada uma das escolas cordialmente fornecidas pelo Parque Escolar, E.P.E., foi feito um estudo exaustivo e muito individual de cada uma das escolas, onde cada análise foi desenvolvida por uma metodologia simplificada, sendo cada análise sempre individual, e nunca tendo sida aplicada em série, este fator melhora substancialmente a eficácia da avaliação para a quantificação das vulnerabilidades e da determinação do grau de dano e das frações de perda para os requisitos fundamentais de limitação de danos, de não colapso, e de colapso iminente que correspondem a ações com períodos de retorno de 95, 475 e 975 anos. Este trabalho é fundamental para as entidades competentes terem a consciência da vulnerabilidade das escolas secundárias, para puderem atuar a nível estrutural e diminuir assim a vulnerabilidade sísmica, e mesmo que por impossibilidade económica o poder governamental não intervenha, então pode e principalmente deve, elaborar planos de emergência tanto com engenheiros civis qualificados como com a total colaboração das corporações de bombeiros que fazem parte das forças de operação e socorro da Autoridade Nacional de Proteção Civil (ANPC).

Contribuição do isolamento térmico para o aumento da eficiência energética em edifícios

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Educating for earthquake science and risk in a tectonically slowly deforming region

Over the past decade, scientists have been called to participate more actively in public education and outreach (E&O). This is particularly true in fields of significant societal impact, such as earthquake science. Local earthquake risk culture plays a role in the way that the public engages in educational efforts. In this article, we describe an adapted E&O program for earthquake science and risk. The program is tailored for a region of slow tectonic deformation, where large earthquakes are extreme events that occur with long return periods. The adapted program has two main goals: (1) to increase the awareness and preparedness of the population to earthquake and related risks (tsunami, liquefaction, fires, etc.), and (2) to increase the quality of earthquake science education, so as to attract talented students to geosciences. Our integrated program relies on activities tuned for different population groups who have different interests and abilities, namely young children, teenagers, young adults, and professionals.

Modelación estocástica y trading algorítmico del spread entre acciones mediante procesos de reversión a la media

Las estrategias de inversión pairs trading se basan en desviaciones del precio entre pares de acciones correlacionadas y han sido ampliamente implementadas por fondos de inversión tomando posiciones largas y cortas en las acciones seleccionadas cuando surgen divergencias y obteniendo utilidad cerrando la posición al converger. Se describe un modelo de reversión a la media para analizar la dinámica que sigue el diferencial del precio entre acciones ordinarias y preferenciales de una misma empresa en el mismo mercado. La media de convergencia en el largo plazo es obtenida con un filtro de media móvil, posteriormente, los parámetros del modelo de reversión a la media se estiman mediante un filtro de Kalman bajo una formulación de estado espacio sobre las series históricas. Se realiza un backtesting a la estrategia de pairs trading algorítmico sobre el modelo propuesto indicando potenciales utilidades en mercados financieros que se observan por fuera del equilibrio. Aplicaciones de los resultados podrían mostrar oportunidades para mejorar el rendimiento de portafolios, corregir errores de valoración y sobrellevar mejor periodos de bajos retornos.