MAST solution of advection problems in irrotational flow fields

A new numerical–analytical Eulerian procedure is proposed for the solution of convection-dominated problems in the case of existing scalar potential of the flow field. The methodology is based on the conservation inside each computational elements of the 0th and 1st order effective spatial moments of the advected variable. This leads to a set of small ODE systems solved sequentially, one element after the other over all the computational domain, according to a MArching in Space and Time technique. The proposed procedure shows the following advantages: (1) it guarantees the local and global mass balance; (2) it is unconditionally stable with respect to the Courant number, (3) the solution in each cell needs information only from the upstream cells and does not require wider and wider stencils as in most of the recently proposed higher-order methods; (4) it provides a monotone solution. Several 1D and 2D numerical test have been performed and results have been compared with analytical solutions, as well as with results provided by other recent numerical methods.     

Water-table and discharge changes associated with the 2016–2017 seismic sequence in central Italy: hydrogeological data and a conceptual model for fractured carbonate aquifers

A seismic sequence in central Italy from August 2016 to January 2017 affected groundwater dynamics in fractured carbonate aquifers. Changes in spring discharge, water-table position, and streamflow were recorded for several months following nine Mw 5.0–6.5 seismic events. Data from 22 measurement sites, located within 100 km of the epicentral zones, were analyzed. The intensity of the induced changes were correlated with seismic magnitude and distance to epicenters. The additional post-seismic discharge from rivers and springs was found to be higher than 9 m³/s, totaling more than 0.1 km³ of groundwater release over 6 months. This huge and unexpected contribution increased streamflow in narrow mountainous valleys to previously unmeasured peak values. Analogously to the L’Aquila 2009 post-earthquake phenomenon, these hydrogeological changes might reflect an increase of bulk hydraulic conductivity at the aquifer scale, which would increase hydraulic heads in the discharge zones and lower them in some recharge areas. The observed changes may also be partly due to other mechanisms, such as shaking and/or squeezing effects related to intense subsidence in the core of the affected area, where effects had maximum extent, or breaching of hydraulic barriers.

     

Patuxent landscape model: 1. Hydrological model development

We developed a spatially explicit, process-based model of the 2352 km² Patuxent river watershed in Maryland, and its subwatersheds to integrate data and knowledge over several spatial, temporal and complexity scales, and to serve as an aid to regional management. The model was developed using the Library of Hydro-Ecological Modules (LHEM, http://giee.uvm.edu/LHEM), which was designed to create flexible landscape model structures that can be easily modified and extended to suit the requirements of a variety of goals and case studies. The LHEM includes modules that simulate various aspects of ecosystem dynamics. In this paper we consider modules that represent the physical conditions in the environment (climatic factors, geoporphology), and hydrologic processes, both locally and spatially. Where possible the modules are formulated as Stella(R) models, spatial transport processes are presented as C++ code.     

Testing a new rill flow resistance approach using the Water Erosion Prediction Project experimental database

In this paper, a recently theoretically deduced rill flow resistance equation, based on a power‐velocity profile, was tested using the Water Erosion Prediction Project database. This database includes measurements of flow velocity, water depth, cross section area, wetted perimeter, and bed slope that were made in rills shaped on experimental sites distributed across the continental United States. In particular, three different experimental conditions (only rainfall, only flow, and rain with flow) were examined, and for each condition, the theoretically based relationship for estimating the Γ function of the power velocity profile was calibrated. The results established that (a) the Darcy‐Weisbach friction factor can be accurately estimated using the proposed theoretical approach, and (b) the flow resistance increases with the effect of rainfall impact.     

Flow resistance equation for rills

In this paper, a new flow resistance equation for rill flow was deduced applying dimensional analysis and self‐similarity theory. At first, the incomplete self‐similarity hypothesis was used for establishing the flow velocity distribution whose integration gives the theoretical expression of the Darcy–Weisbach friction factor. Then the deduced theoretical resistance equation was tested by some measurements of flow velocity, water depth, cross section area, wetted perimeter, and bed slope carried out in 106 reaches of some rills shaped on an experimental plot. A relationship between the velocity profile, the channel slope, and the flow Froude number was also established. The analysis showed that the Darcy–Weisbach friction factor can be accurately estimated by the proposed theoretical approach based on a power–velocity profile.     

The 5,660?yBP Boquerón explosive eruption, Teide–Pico Viejo complex, Tenerife

Quantitative hazard assessments of active volcanoes require an accurate knowledge of the past eruptive activity in terms of eruption dynamics and the stratified products of eruption. Teide–Pico Viejo (TPV) is one of the largest volcanic complexes in Europe, but the associated eruptive history has only been constrained based on very general stratigraphic and geochronological data. In particular, recent studies have shown that explosive activity has been significantly more frequently common than previously thought. Our study contributes to characterization of explosive activity of TPV by describing for the first time the subplinian eruption of El Boquerón (5,660?yBP), a satellite dome located on the northern slope of the Pico Viejo stratovolcano. Stratigraphic data suggest complex shifting from effusive phases with lava flows to highly explosive phase that generated a relatively thick and widespread pumice fallout deposit. This explosive phase is classified as a subplinian eruption of VEI 3 that lasted for about 9–15?h and produced a plume with a height of up to 9?km above sea level (i.e. 7?km above the vent; MER of 6.9–8.2?×?10⁵?kg/s). The tephra deposit (minimum bulk volume of 4–6?×?10⁷?m³) was dispersed to the NE by up to 10?m/s winds. A similar eruption today would significantly impact the economy of Tenerife (e.g. tourism and aviation), with major consequences mainly for the communities around the Icod Valley, and to a minor extent, the Orotava Valley. This vulnerability shows that a better knowledge of the past explosive history of TPV and an accurate estimate of future potentials to generate violent eruptions is required in order to quantify and mitigate the associated volcanic risk.     

Testing a new sampler for measuring plot soil loss

In order to measure soil loss in equipped plots the estimate of the weight of solid material intercepted at their lower end is required. At the experimental area of Sparacia, Sicily, the runoff produced by an erosive event is collected within storage tanks with a capacity of about 1 m³. In this paper, the use of a new sampler is proposed to measure easily the weight of solid material eroded from an experimental plot and collected into a storage tank. The sampler is a cylinder having a closing valve at the bottom. Two different series of runs were carried out both to test the reliability of the sampler and to establish a sampling procedure, respectively. An analysis of various sampling configurations usable in the field differentiated by the number and location of sampling verticals in the tank cross‐section was finally carried out. The results of the present investigation are that the concentration measurement by the sampler was more accurate than that obtained by other methods involving a collection tank, agitation and sampling of the suspension. This sampler is cheap and usable in combination with a quick field sampling procedure which is particularly advisable when the number of plots equipped at an experimental area is large. The sampler was tested using a clay soil contained within cylinders and a cubic tank, but it appeared also to be usable with coarser sediment than clay and in combination with tanks having a different shape. Copyright © 2015 John Wiley & Sons, Ltd.     

Ejby—A new H5/6 ordinary chondrite fall in Copenhagen,Denmark

On February 6, 2016 at 21:07:19 UT, a very bright fireball was seen over the eastern part of Denmark. The weather was cloudy over eastern Denmark, but many people saw the sky light up—even in the heavily illuminated Copenhagen. Two hundred and thirty three reports of the associated sound and light phenomena were received by the Danish fireball network. We have formed a consortium to describe the meteorite and the circumstances of the fall and the results are presented in this paper. The first fragment of the meteorite was found the day after the fall, and in the following weeks, a total of 11 fragments with a total weight of 8982 g were found. The meteorite is an unbrecciated, weakly shocked (S2), ordinary H chondrite of petrologic type 5/6 (Bouvier et al. 2017). The concentration of the cosmogenic radionuclides suggests that the preatmospheric radius was rather small ~20 cm. The cosmic ray exposure age of Ejby (83 ± 11 Ma) is the highest of an H chondrite and the second highest age for an ordinary chondrite. Using the preatmospheric orbit of the Ejby meteoroid (Spurny et al. 2017) locations of the recovered fragments, and wind data from the date of the fall, we have modeled the dark flight (below 18 km) of the fragments. The recovery location of the largest fragment can only be explained if aerodynamic effects during the dark flight phase are included. The recovery location of all other fragments are consistent with the dark flight modeling.