Runout analysis of a large rockfall in the Dolomites/Italian Alps using LIDAR derived particle sizes and shapes

LiDAR data were used to quantify and analyse a rockfall event which occurred in 2003 in the Western Dolomites (Italian Alps). In addition to previously existing airborne laserscanning (ALS) data, high resolution terrestrial laserscanning (TLS) data were collected. By using the original point clouds, the volume, axial ratio and runout length of single boulders as well as the surface roughness in the runout zone of the rockfall were derived. The total volume of the rockfall event of approximately 10 000 m³ was estimated by a reconstruction of the pre‐event surface at the detachment zone. The analysis of the laser scanning data of the accumulation zone revealed a power law scaling for boulder volumes larger than 8 m³. The dependence of runout length on boulder volume is complex; it is moderated by particle sphericity. In addition, we used ALS and TLS data to derive the spatial distribution of surface roughness on the talus cone. TLS allow for more accurate roughness mapping than ALS data, but for most applications the point density of ALS data seems to be sufficiently high to derive measures of roughness. Different sampling approaches for plane fitting on the scale of 5 m did not show significant effects besides the computational time. The results of our analyses provide important perspectives for rockfall modelling and process understanding with potential applications in both ‘applied’ (natural hazards) and ‘pure’ geomorphological research. Copyright © 2012 John Wiley & Sons, Ltd.     

Comments on “Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady-state” by O.A. Cirpka and A.J. Valocchi

In this comment we present a re-analysis of the analytical solution presented by Cirpka and Valocchi for steady-state concentrations of dissolved bioreactive compounds and bacterial biomass in porous media. We discuss the validity range of the analytical solution. In particular, the criterion used to determine the sustainability of biomass is revisited. This re-analysis shows that the ω$\omega$ criterion used by Cirpka and Valocchi is only a necessary but not a sufficient criterion to determine the bioreactive zones. As a consequence, the analytical solution does not provide the exact distribution of compounds throughout the domain, but can serve as upper or lower boundaries for species concentrations at a given location. These conclusions are supported by the simulation results obtained from an established reactive transport model.     

Thermal influence of urban groundwater recharge from stormwater infiltration basins

Groundwater warming below cities has become a major environmental issue; but the effect of distinct local anthropogenic sources of heat on urban groundwater temperature distributions is still poorly documented. Our study addressed the local effect of stormwater infiltration on the thermal regime of urban groundwater by examining differences in water temperature beneath stormwater infiltration basins (SIB) and reference sites fed exclusively by direct infiltration of rainwater at the land surface. Stormwater infiltration dramatically increased the thermal amplitude of groundwater at event and season scales. Temperature variation at the scale of rainfall events reached 3 °C and was controlled by the interaction between runoff amount and difference in temperature between stormwater and groundwater. The annual amplitude of groundwater temperature was on average nine times higher below SIB (range: 0·9–8·6 °C) than at reference sites (range: 0–1·2 °C) and increased with catchment area of SIB. Elevated summer temperature of infiltrating stormwater (up to 21 °C) decreased oxygen solubility and stimulated microbial respiration in the soil and vadose zone, thereby lowering dissolved oxygen (DO) concentration in groundwater. The net effect of infiltration on average groundwater temperature depended upon the seasonal distribution of rainfall: groundwater below large SIB warmed up (+0·4 °C) when rainfall occurred predominantly during warm seasons. The thermal effect of stormwater infiltration strongly attenuated with increasing depth below the groundwater table indicating advective heat transport was restricted to the uppermost layers of groundwater. Moreover, excessive groundwater temperature variation at event and season scales can be attenuated by reducing the size of catchment areas drained by SIB and by promoting source control drainage systems. Copyright © 2009 John Wiley & Sons, Ltd.     

It starts at home? Climate policies targeting household consumption and behavioral decisions are key to low-carbon futures

Through their consumption behavior, households are responsible for 72% of global greenhouse gas emissions. Thus, they are key actors in reaching the 1.5 °C goal under the Paris Agreement. However, the possible contribution and position of households in climate policies is neither well understood, nor do households receive sufficiently high priority in current climate policy strategies. This paper investigates how behavioral change can achieve a substantial reduction in greenhouse gas emissions in European high-income countries. It uses theoretical thinking and some core results from the HOPE research project, which investigated household preferences for reducing emissions in four European cities in France, Germany, Norway and Sweden. The paper makes five major points: First, car and plane mobility, meat and dairy consumption, as well as heating are the most dominant components of household footprints. Second, household living situations (demographics, size of home) greatly influence the household potential to reduce their footprint, even more than country or city location. Third, household decisions can be sequential and temporally dynamic, shifting through different phases such as childhood, adulthood, and illness. Fourth, short term voluntary efforts will not be sufficient by themselves to reach the drastic reductions needed to achieve the 1.5 °C goal; instead, households need a regulatory framework supporting their behavioral changes. Fifth, there is a mismatch between the roles and responsibilities conveyed by current climate policies and household perceptions of responsibility. We then conclude with further recommendations for research and policy.     

Estimating uncertainty associated with water stages from a single SAR image

It is the goal of remote sensing to infer information about objects or a natural process from a remote location. This invokes that uncertainty in measurement should be viewed as central to remote sensing. In this study, the uncertainty associated with water stages derived from a single SAR image for the Alzette (G.D. of Luxembourg) 2003 flood is assessed using a stepped GLUE procedure. Main uncertain input factors to the SAR processing chain for estimating water stages include geolocation accuracy, spatial filter window size, image thresholding value, DEM vertical precision and the number of river cross sections at which water stages are estimated. Initial results show that even with plausible parameter values uncertainty in water stages over the entire river reach is 2.8 m on average. Adding spatially distributed field water stages to the GLUE analysis following a one-at-a-time approach helps to considerably reduce SAR water stage uncertainty (0.6 m on average) thereby identifying appropriate value ranges for each uncertain SAR water stage processing factor. For the GLUE analysis a Nash-like efficiency criterion adapted to spatial data is proposed whereby acceptable SAR model simulations are required to outperform a simpler regression model based on the field-surveyed average river bed gradient. Weighted CDFs for all factors based on the proposed efficiency criterion allow the generation of reliable uncertainty quantile ranges and 2D maps that show the uncertainty associated with SAR-derived water stages. The stepped GLUE procedure demonstrated that not all field data collected are necessary to achieve maximum constraining. A possible efficient way to decide on relevant locations at which to sample in the field is proposed. It is also suggested that the resulting uncertainty ranges and flood extent or depth maps may be used to evaluate 1D or 2D flood inundation models in terms of water stages, depths or extents. For this, the extended GLUE approach, which copes with the presence of uncertainty in the observed data, may be adopted.     

Water Table Uncertainties due to Uncertainties in Structure and Properties of an Unconfined Aquifer

We consider two sources of geology‐related uncertainty in making predictions of the steady‐state water table elevation for an unconfined aquifer. That is the uncertainty in the depth to base of the aquifer and in the hydraulic conductivity distribution within the aquifer. Stochastic approaches to hydrological modeling commonly use geostatistical techniques to account for hydraulic conductivity uncertainty within the aquifer. In the absence of well data allowing derivation of a relationship between geophysical and hydrological parameters, the use of geophysical data is often limited to constraining the structural boundaries. If we recover the base of an unconfined aquifer from an analysis of geophysical data, then the associated uncertainties are a consequence of the geophysical inversion process. In this study, we illustrate this by quantifying water table uncertainties for the unconfined aquifer formed by the paleochannel network around the Kintyre Uranium deposit in Western Australia. The focus of the Bayesian parametric bootstrap approach employed for the inversion of the available airborne electromagnetic data is the recovery of the base of the paleochannel network and the associated uncertainties. This allows us to then quantify the associated influences on the water table in a conceptualized groundwater usage scenario and compare the resulting uncertainties with uncertainties due to an uncertain hydraulic conductivity distribution within the aquifer. Our modeling shows that neither uncertainties in the depth to the base of the aquifer nor hydraulic conductivity uncertainties alone can capture the patterns of uncertainty in the water table that emerge when the two are combined.     

Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.     

Photometric survey of binary near-Earth asteroids

Photometric data on 17 binary near-Earth asteroids (15 of them are certain detections, two are probables) were analysed and characteristic properties of the near-Earth asteroid (NEA) binary population were inferred. We have found that binary systems with a secondary-to-primary mean diameter ratio Ds/Dp?0.18 concentrate among NEAs smaller than 2 km in diameter; the abundance of such binaries decreases significantly among larger NEAs. Secondaries show an upper size limit of Ds=0.5-1 km. Systems with Ds/Dp?0.5 are abundant but larger satellites are significantly less common. Primaries have spheroidal shapes and they rotate rapidly, with periods concentrating between 2.2 to 2.8 h and with a tail of the distribution up to ∼4 h. The fast rotators are close to the critical spin for rubble piles with bulk densities about 2 g/cm³. Orbital periods show an apparent cut-off at Porb∼11 h; closer systems with shorter orbital periods have not been discovered, which is consistent with the Roche limit for strengthless bodies. Secondaries are more elongated on average than primaries. Most, but not all, of their rotations appear to be synchronized with the orbital motion; nonsynchronous secondary rotations may occur especially among wider systems with Porb>20 h. The specific total angular momentum of most of the binary systems is similar to within ±20% and close to the angular momentum of a sphere with the same total mass and density, rotating at the disruption limit; this suggests that the binaries were created by mechanism(s) related to rotation near the critical limit and that they neither gained nor lost significant amounts of angular momentum during or since formation. A comparison with six small asynchronous binaries detected in the main belt of asteroids suggests that the population extends beyond the region of terrestrial planets, but with characteristics shifted to larger sizes and longer periods. The estimated mean proportion of binaries with Ds/Dp?0.18 among NEAs larger than 0.3 km is 15±4%. Among fastest rotating NEAs larger than 0.3 km with periods between 2.2 and 2.8 h, the mean proportion of such binaries is (66⁺¹⁰₋₁₂)%.