Point process models for fine-resolution rainfall

Abstract

In a recent development in the literature, a new temporal rainfall model, based on the Bartlett-Lewis clustering mechanism and intended for sub-hourly application, was introduced. That model replaced the rectangular rain cells of the original model with finite Poisson processes of instantaneous pulses, allowing greater variability in rainfall intensity over short intervals. In the present paper, the basic instantaneous pulse model is first extended to allow for randomly varying storm types. A systematic comparison of a number of key model variants, fitted to 5-min rainfall data from Germany, then generates further new insights into the models, leading to the development of an additional model extension, which introduces dependence between rainfall intensity and duration in a simple way. The new model retains the original rectangular cells, previously assumed inappropriate for fine-scale data, obviating the need for the computationally more intensive instantaneous pulse model.

Editor D. Koutsoyiannis     

Speleothems in a north Cuban cave register sea‐level changes and Pleistocene uplift rates

A flight of marine terraces along the Cuban coast records Quaternary sea‐level highstands and a general slowly uplifting trend during the Pleistocene. U/Th dating of these limestone terraces is difficult because fossil reef corals have been affected by open system conditions. Terrace ages are thus often based on geological and geomorphological observations. In contrast, the minimum age of the terraces can be constrained by dating speleothems from coastal mixing (flank margin) caves formed during past sea‐level highstands and carving the marine limestones. Speleothems in Santa Catalina Cave have ages >360 ka and show various cycles of subaerial–subaqueous corrosion and speleothem growth. This suggests that the cave was carved during the MIS 11 sea‐level highstand or earlier. Some stalagmites grew during MIS 11 through MIS 8 and were submerged twice, once at the end of MIS 11 and then during MIS 9. Phreatic overgrowths (POS) covering the speleothems suggest anchialine conditions in the cave during MIS 5e. Their altitude at 16 m above present sea level indicates a late Pleistocene uplift rate of <0.1 mm/ka, but modelling also shows uplift to have been insignificant over a long timespan during the middle Pleistocene since the cave was carved. Our study shows that some flank margin caves in the region of Matanzas are older than commonly believed (i.e. MIS 11 rather than MIS 5). These caves not only can be preserved but are good markers of interglacial sea‐level highstands, more reliable than marine abrasion surfaces. Copyright © 2018 John Wiley & Sons, Ltd.     

Factors determining subsidence in urbanized floodplains: evidence from MT‐InSAR in Seville (southern Spain)

Major rivers have traditionally been linked with important human settlements throughout history. The growth of cities over recent river deposits makes necessary the use of multidisciplinary approaches to characterize the evolution of drainage networks in urbanized areas. Since under‐consolidated fluvial sediments are especially sensitive to compaction, their spatial distribution, thickness, and mechanical behavior must be studied. Here, we report on subsidence in the city of Seville (Southern Spain) between 2003 and 2010, through the analysis of the results obtained with the Multi‐Temporal InSAR (MT‐InSAR) technique. In addition, the temporal evolution of the subsidence is correlated with the rainfall, the river water column and the piezometric level. Finally, we characterize the geotechnical parameters of the fluvial sediments and calculate the theoretical settlement in the most representative sectors. Deformation maps clearly indicate that the spatial extent of subsidence is controlled by the distribution of under‐consolidated fine‐grained fluvial sediments at heights comprised in the range of river level variation. This is clearly evident at the western margin of the river and the surroundings of its tributaries, and differs from rainfall results as consequence of the anthropic regulation of the river. On the other hand, this influence is not detected at the eastern margin due to the shallow presence of coarse‐grain consolidated sediments of different terrace levels. The derived results prove valuable for implementing urban planning strategies, and the InSAR technique can therefore be considered as a complementary tool to help unravel the subsidence tendency of cities located over under‐consolidated fluvial deposits. Copyright © 2017 John Wiley & Sons, Ltd.     

Sensitivity analysis of a morphodynamic modelling system applied to a coastal lagoon inlet

This work investigates the recent morphological changes at the inlet of a complex coastal system (Ria de Aveiro lagoon, Portugal). This study was carried out using bathymetric data analysis and numerical simulations obtained with the 2DH morphodynamic modelling system MORSYS2D. The present simulations considered only tidal forcing, and a sensitivity analysis was performed by tuning the formula used to compute the sediment transports. A non-uniform sediment grain size distribution for the Ria de Aveiro inlet is considered in the numerical simulations, based on surveys performed in this area. The model results are analysed to assess if they resemble the observed trends of erosion and deposition, as calculated from bathymetric data. A quantitative analysis of the differences between the bathymetric changes obtained through surveys and the numerical results over a period of 3 years considering different sediment transport formulations shows that the formulations of Ackers and White (1973) and Engelund and Hansen (1967) are the ones that best describe the morphodynamic changes driven by tides in the Ria de Aveiro inlet.     

Gully-head morphology and implications for gully development on abandoned fields in a semi-arid environment,Sierra de Gata,southeast Spain

In this paper we examine whether gully-head morphology can be used as an indicator for gully development and, hence, for sediment production. A survey was conducted at five hillslopes in the Sierra de Gata where different types of channel heads occur close to each other. The survey included measurements of morphologic and pedologic properties, ground surface, channel and catchment characteristics of every gully head present (n = 59). On the basis of the observed morphologies, the heads were subdivided into four types: gradual, transitional (a short inclined section), abrupt and rilled-abrupt. The analyses showed that it is possible to explain the differences of gully heads and the role of some environmental factors on the basis of their morphologies, at least for the gradual and the abrupt types. The results suggested that steep headcuts (abrupt) were formed from secondary headcuts in the channel, which migrated upstream. The abrupt headcuts were always formed in more than one soil layer of which one was a resistant (stony) layer. However, shear strength measurements (at saturation) showed that the top layer was not always the most resistant one. Width–depth relationships indicated that gradual type headcuts were controlled by fluvial processes and abrupt headcuts by a combination of fluvial and mass-wasting processes. Gradual types occurred more downslope than the abrupt types suggesting that the incisions started by fluvial processes and migrated upwards when knickpoints developed in the channel. The rilled-abrupt types are still actively retreating. Thus, the abrupt types correspond to slower retreat rates. Abrupt gully heads may deteriorate into transitional types when plunge-pool erosion becomes less effective. The conceptual model is supported by data from ephemeral gullies in two other study areas (Sierra de la Torrecilla, Spain, and Alentejo, Portugal). Copyright © 1999 John Wiley & Sons, Ltd.     

Lateral terrestrial water flow contribution to summer precipitation at continental scale – A comparison between Europe and West Africa with WRF-Hydro-tag ensembles

It is well accepted that summer precipitation can be altered by soil moisture condition. Coupled land surface – atmospheric models have been routinely used to quantify soil moisture – precipitation feedback processes. However, most of the land surface models (LSMs) assume a vertical soil water transport and neglect lateral terrestrial water flow at the surface and in the subsurface, which potentially reduces the realism of the simulated soil moisture – precipitation feedback. In this study, the contribution of lateral terrestrial water flow to summer precipitation is assessed in two different climatic regions, Europe and West Africa, for the period June–September 2008. A version of the coupled atmospheric-hydrological model WRF-Hydro with an option to tag and trace land surface evaporation in the modelled atmosphere, named WRF-Hydro-tag, is employed. An ensemble of 30 simulations with terrestrial routing and 30 simulations without terrestrial routing is generated with random realizations of turbulent energy with the stochastic kinetic energy backscatter scheme, for both Europe and West Africa. The ensemble size allows to extract random noise from continental-scale averaged modelled precipitation. It is found that lateral terrestrial water flow increases the relative contribution of land surface evaporation to precipitation by 3.6% in Europe and 5.6% in West Africa, which enhances a positive soil moisture – precipitation feedback and generates more uncertainty in modelled precipitation, as diagnosed by a slight increase in normalized ensemble spread. This study demonstrates the small but non-negligible contribution of lateral terrestrial water flow to precipitation at continental scale.     

Diurnal cycle of surface energy fluxes in high mountain terrain: High-resolution fully coupled atmosphere-hydrology modelling and impact of lateral flow

Water and energy fluxes are inextricably interlinked within the interface of the land surface and the atmosphere. In the regional earth system models, the lower boundary parameterization of land surface neglects lateral hydrological processes, which may inadequately depict the surface water and energy fluxes variations, thus affecting the simulated atmospheric system through land-atmosphere feedbacks. Therefore, the main objective of this study is to evaluate the hydrologically enhanced regional climate modelling in order to represent the diurnal cycle of surface energy fluxes in high spatial and temporal resolution. In this study, the Weather Research and Forecasting model (WRF) and coupled WRF Hydrological modelling system (WRF-Hydro) are applied in a high alpine catchment in Northeastern Tibetan Plateau, the headwater area of the Heihe River. By evaluating and intercomparing model results by both models, the role of lateral flow processes on the surface energy fluxes dynamics is investigated. The model evaluations suggest that both WRF and coupled WRF-Hydro reasonably represent the diurnal variations of the near-surface meteorological fields, surface energy fluxes and hourly partitioning of available energy. By incorporating additional lateral flow processes, the coupled WRF-Hydro simulates higher surface soil moisture over the mountainous area, resulting in increased latent heat flux and decreased sensible heat flux of around 20–50 W/m² in their diurnal peak values during summertime, although the net radiation and ground heat fluxes remain almost unchanged. The simulation results show that the diurnal cycle of surface energy fluxes follows the local terrain and vegetation features. This highlights the importance of consideration of lateral flow processes over areas with heterogeneous terrain and land surfaces.     

Urban Thermodynamic Island in a Coastal City Analysed from an Optimized Surface Network

Within the framework of ESCOMPTE, a French experiment performed in June and July 2001 in the south-east of France to study the photo-oxidant pollution at the regional scale, the urban boundary layer (UBL) program focused on the study of the urban atmosphere over the coastal city of Marseille. A methodology developed to optimize a network of 20 stations measuring air temperature and moisture over the city is presented. It is based on the analysis of a numerical simulation, performed with the non-hydrostatic, mesoscale Meso-NH model, run with four nested-grids down to a horizontal resolution of 250 m over the city and including a specific parametrization for the urban surface energy balance. A three-day period was modelled and evaluated against data collected during the preparatory phase for the project in summer 2000. The simulated thermodynamic surface fields were analysed using an empirical orthogonal function (EOF) decomposition in order to determine the optimal network configuration designed to capture the dominant characteristics of the fields. It is the first attempt of application of this kind of methodology to the field of urban meteorology. The network, of 20 temperature and moisture sensors, was implemented during the UBL-ESCOMPTE experiment and continuously recorded data from 12 June to 14 July 2001. The measurements were analysed in order to assess the urban thermodynamic island spatio-temporal structure, also using EOF decomposition. During nighttime, the influence of urbanization on temperature is clear the field is characterized by concentric thermo-pleths around the old core of the city, which is the warmest area of the domain. The moisture field is more influenced by proximity to the sea and airflow patterns. During the day, the sea breeze often moves from west or south-west and consequently the spatial pattern for both parameters is characterized by a gradient perpendicular to the shoreline. Finally, in order to assess the methodology adopted, the spatial structures extracted from the simulation of the 2000 preparatory campaign and observations gathered in 2001 have been compared. They are highly correlated, which is a relevant validation of the methodology proposed. The relations between these spatial structures and geographical characteristics of the site have also been studied. High correlations between temperature spatial structure during nighttime and urban cover fraction or street aspect ratio are observed and simulated. For temperature during daytime or moisture during both daytime and nighttime these geographical factors are not correlated with thermodynamic fields spatial structures.     

A method for climate and vegetation reconstruction through the inversion of a dynamic vegetation model

Climate reconstructions from data sensitive to past climates provide estimates of what these climates were like. Comparing these reconstructions with simulations from climate models allows to validate the models used for future climate prediction. It has been shown that for fossil pollen data, gaining estimates by inverting a vegetation model allows inclusion of past changes in carbon dioxide values. As a new generation of dynamic vegetation model is available we have developed an inversion method for one model, LPJ-GUESS. When this novel method is used with high-resolution sediment it allows us to bypass the classic assumptions of (1) climate and pollen independence between samples and (2) equilibrium between the vegetation, represented as pollen, and climate. Our dynamic inversion method is based on a statistical model to describe the links among climate, simulated vegetation and pollen samples. The inversion is realised thanks to a particle filter algorithm. We perform a validation on 30 modern European sites and then apply the method to the sediment core of Meerfelder Maar (Germany), which covers the Holocene at a temporal resolution of approximately one sample per 30 years. We demonstrate that reconstructed temperatures are constrained. The reconstructed precipitation is less well constrained, due to the dimension considered (one precipitation by season), and the low sensitivity of LPJ-GUESS to precipitation changes.     

Sensitivity study and validation of a land surface parameterization using the HAPEX-MOBILHY data set

A simple parameterization of land surface processes, amenable to the structure of a two-layer soil model, including a representation of the vegetation, has been designed for use in meteorological models. Prior to implementation in a mesoscale model, it is necessary to check the components and to verify the good working order of the parameterization as a whole. The aims of this paper then are: (i) evaluation and a sensitivity study of the various components of the model, specifying the needed accuracy for the parameters; (ii) micrometeorological validation of the model against the HAPEX-MOBILHY data set.First, we present the basic scheme. The focus is on the parameterization of surface resistance, and especially on its relationship with soil moisture.A sensitivity study is then performed through a set of one-dimensional simulations which allow a full interaction between the ground and the atmosphere. Above bare ground, it is shown that both soil texture and initial moisture greatly influence the outcome of the simulation. Latent heat flux ranges from that associated with potential evaporation through a switch-like behavior to that of dry soil. Next, the effects of transpiring vegetation canopies on the physical processes involved and the surface energy balance are examined. The sensitivity of the latent heat flux to changes in the soil and canopy parameters is emphazised; the major influence of the initial mean soil moisture and of the vegetation cover is pointed out. Finally, the evolution of the boundary layer in response to various surface conditions is studied.A validation of the land surface scheme is conducted through daily cycles during cloudless days. Simulated turbulent fluxes are successfully compared to micrometeorological measurements over a maize field at different growth stages. Over a pine forest, the correct simulation of the turbulent fluxes is obtained with an adequate parameterization of the surface resistance accounting for the atmospheric moisture deficit.