Assessment of potential sediment contamination using screening methods (XRF,TGA/MS) taking into account principles of green chemistry,Eastern Slovakia

Mining and milling of metal ores coupled with industries have bequeathed many countries the legacy of wide distribution of metal contaminants in sediments. The aim of this study was to assess potential sediment contamination via useful screening methods (XRF, CHNS, TGA/MS). The sediments were collected from the water reservoir Krompachy Eastern Slovakia in April 2015. Within the frame of evaluation it was found that the concentrations of the study elements (Cu, Zn, As, Pb, Cr, Ni, Cd) exceeded some of the MPC, TV and IV values. Sample c was the most polluted by metals, which evident according to it’s the highest CHNS proportion as well as the highest clay and silt proportion. In the samples studied the best correlation was confirmed between weight losses in the temperature range (400–620 °C) and the following metal concentrations: Cu (r = 0.89), Zn (r = 0.88), As (r = 0.93), Hg (r = 0.83), Pb (r = 0.87). The greatest proportions of m/z 44, m/z 18 were detected at temperatures (400–620 °C) associated with decomposition of minerals such as siderite, barite, and exothermic loss of more refractory aromatic C took also place.     

Possible reasons for the frequency splitting of the harmonics of type II solar radio bursts

AIA/SDO data in the 193 Å channel preceding a coronal mass ejection observed at the solar limb on June 13, 2010 are used to simultaneously identify and examine two different shock fronts. The angular size of each front relative to the CME center was about 20°, and their propagation directions differed by ≈25° (≈4° in position angle). The faster front, called the blast shock, advanced the other front, called the piston shock, by R ≈ (0.02-0.03)R⊙ (R⊙ is the solar radius) and had a maximum initial speed of V_B ≈ 850 km/s (with V_P ≈ 700 km/s for the piston shock). The appearance and motion of these shocks were accompanied by a Type II radio burst observed at the fundamental frequency F and second harmonic H. Each frequency was split into two close frequencies f₁ and f₂ separated by Δf = f₂ - f₁ ? F, H. It is concluded that the observed frequency splitting Δf of the F and H components of the Type II burst could result from the simultaneous propagation of piston and blast shocks moving with different speeds in somewhat different directions displaying different coronal-plasma densities.     

Probabilistic seismic hazard assessment in the Constantine region,Northeast of Algeria

This paper presents a seismic hazard evaluation and develops an earthquake catalogue for the Constantine region over the period from 1357 to 2014. The study contributes to the improvement of seismic risk management by evaluating the seismic hazards in Northeast Algeria. A regional seismicity analysis was conducted based on reliable earthquake data obtained from various agencies (CRAAG, IGN, USGS and ISC). All magnitudes (M _l, m _b) and intensities (I ₀, I _MM, I _MSK and I _EMS) were converted to M _s magnitudes using the appropriate relationships. Earthquake hazard maps were created for the Constantine region. These maps were estimated in terms of spectral acceleration (SA) at periods of 0.1, 0.2, 0.5, 0.7, 0.9, 1.0, 1.5 and 2.0 s. Five seismogenic zones are proposed. This new method differs from the conventional method because it incorporates earthquake magnitude uncertainty and mixed datasets containing large historical events and recent data. The method can be used to estimate the b value of the Gutenberg-Richter relationship, annual activity rate λ(M) of an event and maximum possible magnitude M _max using incomplete and heterogeneous data files. In addition, an earthquake is considered a Poisson with an annual activity rate λ and with a doubly truncated exponential earthquake magnitude distribution. Map of seismic hazard and an earthquake catalogue, graphs and maps were created using geographic information systems (GIS), the Z-map code version 6 and Crisis software 2012.     

Quantifying land use influences on event-based flow frequency,timing, magnitude,and rate of change in an urbanizing watershed of the central USA

There is a need for research that advances understanding of flow alterations in contemporary watersheds where natural and anthropogenic interactions can confound mitigation efforts. Event-based flow frequency, timing, magnitude, and rate of change were quantified at five-site nested gauging sites in a representative mixed-land-use watershed of the central USA. Statistically independent storms were paired by site (n = 111 × 5 sites) to test for significant differences in event-based rainfall and flow response variables (n = 17) between gauging sites. Increased frequency of small peak flow events (i.e., 64 more events less than 4.0 m³ s^?1) was observed at the rural–urban interface of the watershed. Differences in flow response were apparent during drier periods when small rainfall events resulted in increased flow response at urban sites in the lower reaches. Relationships between rainfall and peak flow were stronger with decreased pasture/crop land use and increased urban land use by approximately 20%. Event-based total rainfall explained 40–68% of the variance in peak flow (p < 0.001). Coefficients of determination (r²) were negatively correlated with pasture/crop land use (r² = 0.92; p = 0.007; n = 5) and positively correlated with urban land use (r² = 0.90; p = 0.008; n = 5). Significant differences in flow metrics were observed between rural and urban sites (p < 0.05; n = 111) that were not explained by differences in rainfall variables and drainage area. An urban influence on flow timing was observed using median time lag to peak centroid and time of maximum precipitation to peak flow. Results highlight the need to establish manageable flow targets in rapidly urbanizing mixed-land-use watersheds.     

First experimental studies of a perturbed zone preceding the front of a coronal mass ejection

The first experimental evidence for a perturbed zone that is likely filled with fast magnetoacoustic oscillations and precedes a coronal mass ejection is presented. When the speed of the coronal mass ejection exceeds the Alfven speed, an outward-moving discontinuity of the plasma density is observed in front of the perturbed zone, on scales comparable to the mean-free path for proton-proton collisions. This suggests that this discontinuity should be interpreted as a collisional shock at distances of R < 30 R _⊙ (where R _⊙ is the solar radius).     

Effect of foliation orientation on the P- and S-wave velocity anisotropies and dynamic elastic constants of the quartz-micaschists metamorphic rocks,Angouran mine,Iran

Laboratory measurements are required to study geophysical properties of the subsurface because of lacking direct observation of Earth’s crust. In this research, compressional (P) and shear (S) wave velocity measurements have been conducted on cylindrical specimens of Quartz-micaschist cored using rock blocks taken from the zinc and lead Angouran mine, Zanjan, northwest of Iran. Cylindrical rock specimens were prepared from the blocks by coring in 0°, 30°, 45°, 60°, and 90° into the foliation direction. P- and S-wave velocities were measured along the cylindrical specimens with different foliation orientations. Percent variations of the P- and S-wave velocities (Thomsen’s anisotropic parameters ε and γ) and constant dynamic modulus of test results have been determined. Percent variations of the P-wave velocity (ε) increase with an increase of the foliation angle with respect to the propagating waves direction by a parabolic function as it shows P-wave velocity differences up to a maximum value of 50 %. Thomsen’s anisotropic parameter of γ has also the same function with the foliation angle. Meanwhile, foliation orientation has a much greater influence on ε than γ for foliation angle from 45° to 90° as \( \frac{\varepsilon }{\gamma } \) ratio increases with an increase of foliation angle. Values of dynamic elastic modulus (E), Poisson’s ratio (ν), shear modulus (μ), bulk modulus (K), and Lamé’s constant (λ) increase with the increase of foliation angle with the parabolic function. The results show that dynamic elastic modulus, Poisson’s ratio, shear modulus, bulk modulus, and Lamé’s constant have anisotropic behavior in relation with the foliation orientation.     

Salinity Tolerance and Competition Drive Distributions of Native and Invasive Submerged Aquatic Vegetation in the Upper San Francisco Estuary

Both abiotic and biotic factors govern distributions of estuarine vegetation, and experiments can reveal effects of these drivers under current and future conditions. In upper San Francisco Estuary (SFE), increased salinity could result from sea level rise, levee failure, or water management. We used mesocosms to test salinity effects on, as well as competition between, the native Stuckenia pectinata (sago pondweed) and invasive Egeria densa (Brazilian waterweed), species with overlapping distributions at the freshwater transition in SFE. Grown alone at a salinity of 5, E. densa decreased fivefold in biomass relative to the freshwater treatment and decomposed within 3 weeks at higher salinities. In contrast, S. pectinata biomass accumulated greatly (~4× initial) at salinities of 0 and 5, doubled at 10, and was unchanged at 15. When grown together in freshwater, S. pectinata produced 75 % less biomass than in monoculture and significantly more nodal roots (suggesting increased nutrient foraging). At a salinity of 5, a decline in E. densa performance coincided with a doubling of S. pectinata shoot density. Additional experiments on E. densa showed elevated temperature (26 and 30 °C) suppressed growth especially at higher salinities (≥5). We conclude that salinity strongly influences distributions of both species and that competition from E. densa may impose limits on S. pectinata abundance in the fresher reaches of SFE. With a salinity increase of 5, S. pectinata is likely to maintain its current distribution while spreading up-estuary at the expense of E. densa, especially if increased temperature also reduces E. densa biomass.     

Impacts of Varying Estuarine Temperature and Light Conditions on <Emphasis Type="Italic">Zostera marina</Emphasis> (Eelgrass) and its Interactions With <Emphasis Type="Italic">Ruppia maritima</Emphasis> (Widgeongrass)

Seagrass populations have been declining globally, with changes attributed to anthropogenic stresses and, more recently, negative effects of global climate change. We examined the distribution of Zostera marina (eelgrass) dominated beds in the York River, Chesapeake Bay, VA over an 8-year time period. Using a temperature-dependent light model, declines in upriver areas were associated with higher light attenuation, resulting in lower light availability relative to compensating light requirements of Z. marina compared with downriver areas. An inverse relationship was observed between SAV growth and temperature with a change between net bed cover increases and decreases for the period of 2004–2011 observed at approximately 23 °C. Z. marina-dominated beds in the lower river have been recovering from a die-off event in 2005 and experienced another near complete decline in 2010, losing an average of 97 % of coverage of Z. marina from June to October. These 2010 declines were attributed to an early summer heat event in which daily mean water temperatures increased from 25 to 30 °C over a 2-week time period, considerably higher than previous years when complete die-offs were not observed. Z. marina recovery from this event was minimal, while Ruppia maritima (widgeongrass) expanded its abundance. Water temperatures are projected to continue to increase in the Chesapeake Bay and elsewhere. These results suggest that short-term exposures to rapidly increasing temperatures by 4–5 °C above normal during summer months can result in widespread diebacks that may lead to Z. marina extirpation from historically vegetated areas, with the potential replacement by other species.     

Potency of constructed wetlands for deportation of pathogens index from rural,urban and industrial wastewater

Pathogen removal is essential for wastewater treatment and its potential reuse in agriculture. Three field-scale wastewater treatment systems consisting of free surface flow were operated around 1.5 years receiving water from urban domestic, rural domestic and industrial sources. The study was conducted to evaluate seasonal performance of constructed wetland systems in removing Escherichia coli, Enterococci and total coliforms under continuous hydraulic flow. Results displayed that all three wetlands gain recognition in removing pathogen load with high removal efficacy till water reaches output ports. Removal efficiencies were even higher, 66–93, 78–92 and 80–94% for E. coli, Enterococci and total coliforms, respectively, within constructed wetlands. Remarkably at shorter temporal scales in CW-A, greater homogeneity of pathogen concentrations was assessed at wetland outlet sites. In outlet ports, results displayed a highly effective removal of E. coli concentration 80–90% (June 2015), 86–92% (October 2015) and 79–92% (February 2016), Enterococci 80–94% (June 2015), 83–94% (October 2015) and 80–94% (February 2016) and total coliforms 85–93% (June 2015), 87–95% (October 2015) and 88–96% (February 2016). Positive correlation was observed between bacterial indicators (E. coli–Enterococci, r = 0.038; p < 0.01 and E. coli–total coliforms, r = 0.142; p < 0.01). Removal of bacterial indicators in constructed wetland was also displayed by PCA in which three-component analysis of variance was 98.39% and showed a clear decrease in measured parameter gradients toward samples from outlet ports. Constructed wetlands provide cost-effective treatment systems for reducing the pathogen load in wastewater in variable agro-climatic conditions and thus improve water quality.     

Sulfate Reduction and Sulfur Cycles at Two Seagrass Beds Inhabited by Cold Affinity <Emphasis Type="Italic">Zostera marina</Emphasis> and Warm Affinity <Emphasis Type="Italic">Halophila nipponica</Emphasis> in Temperate Coastal Waters

To evaluate the impact of invading seagrass on biogeochemical processes associated with sulfur cycles, we investigated the geochemical properties and sulfate reduction rates (SRRs) in sediments inhabited by invasive warm affinity Halophila nipponica and indigenous cold affinity Zostera marina. A more positive relationship between SRR and below-ground biomass (BGB) was observed at the H. nipponica bed (SRR = 0.6809 × BGB ? 4.3162, r ² = 0.9878, p = 0.0006) than at the Z. marina bed (SRR = 0.3470 × BGB ? 4.0341, r ² = 0.7082, p = 0.0357). These results suggested that SR was more stimulated by the dissolved organic carbon (DOC) exuded from the roots of H. nipponica than by the DOC released from the roots of Z. marina. Despite the enhanced SR in spring-summer, the relatively lower proportion (average, 20%) of acid-volatile sulfur (AVS) in total reduced sulfur and the strong correlation between total oxalate-extractable Fe (Fe_(oxal)) and chromium-reducible sulfur (CRS = 0.2321 × total Fe_(oxal) + 1.8180, r ² = 0.3344, p = 0.0076) in the sediments suggested the rapid re-oxidation of sulfide and precipitation of sulfide with Fe. The turnover rate of the AVS at the H. nipponica bed (0.13 day^?1) was 2.5 times lower than that at the Z. marina bed (0.33 day^?1). Together with lower AVS turnover, the stronger correlation of SRR to BGB in the H. nipponica bed suggests that the extension of H. nipponica resulting from the warming of seawater might provoke more sulfide accumulation in coastal sediments.     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: <Emphasis Type="Italic">b</Emphasis>-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.     

Distribution and abundance of hydrobiid snails in a mixed estuary and a coastal lagoon,Argentina

This study examined the environmental factors influencing the distribution and abundance of hydrobiid snails in two estuaries on the northeastern coast of Argentina in a coastal lagoon (Mar Chiquita, 37°40′S, 57°20′W) and a partially mixed estuary (Quequén Grande, 38°30′S, 58°45′W). Five intertidal study sites in each estuary represented a gradient in environmental conditions. Variations in the main environmental factors and in the abundance of hydrobiids were assessed both spatially and seasonally. The three species wereHeleobia australis, Heleobia conexa, andHeleobia parchappii, and they were primarily distributed across a salinity gradient. This pattern was clearly recognizable in the partially mixed estuary, where the abundance ofH. australis decreased as salinity decreased, and the abundance ofH. conexa gradually increased towards the inner reaches of the estuary.H. parchappii was restricted to areas far away from the influence of the tide. Slight differences in the distribution patterns of these species between Quequén Grande and Mar Chiquita were refated to the different dynamics of environmental factors in each estuary.     

Reconstruction of the metamorphic <Emphasis Type="Italic">P–T</Emphasis> path from the garnet zoning in aluminous schists from the Tsogt Block,Mongolian Altai

The paper presents original authors’ data on aluminous schists in the Tsogt tectonic plate in the Southern Altai Metamorphic Belt. The nappe includes a medium-temperature/medium-pressure zonal metamorphic complex, whose metamorphic grade varies from the greenschist to epidote-amphibolite facies. The garnet and garnet–staurolite schists contain three garnet generations of different composition and morphology. The P–T metamorphic parameters estimated by mineralogical geothermometers and geobarometers and by numerical modeling with the PERPLEX 668 software provide evidence of two successive metamorphic episodes: high-gradient (of the andalusite–sillimanite type, geothermal gradient approximately 40–50°/km) and low-gradient (kyanite–sillimanite type, geothermal gradient approximately 27°/km). The P-T parameters of the older episode are T = 545–575°C and P = 3.1–3.7 kbar. Metamorphism during the younger episode was zonal, and its peak parameters were T = 560–565°C, P = 6.4–7.2 kbar for the garnet zone and T = 585–615°C, P = 7.1–7.8 kbar for the staurolite zone. The metamorphism evolved according to a clockwise P–T path: the pressure increased during the first episode at a practically constant temperature, and then during the second episode, the temperature increased at a nearly constant pressure. Such trends are typical of metamorphism related to collisional tectonic settings and may be explained by crustal thickening due to overthrusting. The regional crustal thickening reached at least 15–18 km.     

Magnetic field models for HD 116458 and HD 126515

We have modeled the magnetic fields of the slowly rotating stars HD 116458 and HD 126515 using the “magnetic charge” technique. HD 116458 has a small angle between its rotation axis and dipole axis (β = 12°), whereas this angle is large for HD 126515 (β = 86°). Both stars can be described with a decentered-dipole model, with the respective displacements being r = 0.07 and r = 0.24 in units of the stellar radius. The decentered-dipole model is able to satisfactorily explain the phase relations for the effective field, B_e(P), and the mean surface field, B_s(P), for both stars, along with the fact that the B_e(P) phase relation for HD 126515 is anharmonic. We discuss the role of systematic measurement errors possibly resulting from instrumental or methodical effects in one or both of the phase relations. The displacement of the dipole probably reflects real asymmetry of the stellar field structure, and is not due to measurement errors. Using both phase relations, B_e(P) and B_s(P), in the modeling considerably reduces the influence of the nonuniform distribution of chemical elements on the stellar surface.     

The Effect of High Temperature on Tensile Strength and Damage Characteristics of Limestone

In this paper, the limestone specimens are heated from room temperature 25 to 800 °C in a high temperature furnace and then are subjected to Brazilian test with the AG-I250 electronic precision material testing machine. The physical properties, mechanical properties, disc failure pattern, energy absorbed per unit area and damage characteristics of disc are comprehensively investigated. The results show that: with the increase of temperature, the changing trends of tensile strength, peak strain, tensile modulus and accumulated energy absorbed per unit area of disc are similar, they are first increases, then decrease, the energy consumption index is consistent with the macroscopic damage characteristics; the value of ε_s increase first and then reduce, reaches the maximum at 600 °C. The value of n is increasing and fluctuating, but the change trend of D_c is opposite, which is decreasing and fluctuating. The slope of the damage variable-strain curves decreases first and then increases, the minimum value at 600 °C. This study is of significance to the prediction and evaluation of the stability and safety of rock mass post-high temperature.     

Variations in the orbital periods of the Algol-type eclipsing binaries RZ Cas and Z Dra

A detailed study of variations of the orbital periods of the Algol-type eclipsing binary systems RZ Cas and Z Dra is presented. The fairly complex variations of the periods of both systems can be represented as a superposition of a secular increase of the period, slow periodic fluctuations, and quasiperiodic oscillations with a small amplitude occurring on timescales of decades. The secular increase of the period can be explained by the steady mass transfer from the less massive to the more massive component with conservation of the total angular momentum. The mass-transfer rate is 5.7 × 10^?9M_⊙/yr for RZ Cas and 3.0×10^?8M_⊙/yr for Z Dra. To explain the long-period cyclic variations of the orbital periods of RZCas and Z Dra, it must be assumed that the eclipsing binaries move in long-period orbits. RZ Cas moves with a period of 133 yr around a third body withmass M₃ > 0.55M_⊙, while Z Dra moves with a period of 60 yr around a third body with mass M₃ > 0.7M_⊙. The residual fluctuations of the periods may be due to a superposition of variations due to magnetic cycles and non-stationary ejections of matter.     

Detailed stratigraphy of the Jurassic-Cretaceous boundary beds of the Lena River lower reached based on ammonites and Buchiids

The Jurassic-Cretaceous (Volgian-Lower Valanginian) boundary deposits in the lower reaches of the Lena River (near the village of Chekurovka and on the Cape Chucha, North Yakutia) are described bed-by-bed. The taxonomic composition of ammonites and bivalves (mainly Buchia) is determined. The species assemblages are biostratigraphically analyzed and biostratigraphic units of the bed and zone ranks are recognized based on by Buchia and ammonites. Analogues of the Buchia Zones (buchiazones) Unschensis, Okensis, Tolmatschewi, Inflata and ammonite zones Exoticus, Okensis, Sibiricus, Analogus, Mesezhnikowi and ? Klimovskiensis are established in this region for the first time. A correlation of the distinguished buchiazones with the buchiazones of the Nordvik Peninsula, New Siberian Islands (Stolbovoy Island), the basin of the Anyuy River, and Northern California is proposed. The ammonite assemblages constantly contain phylloceratids; some stratigraphic intervals also contain lytoceratids; this is probably connected with a Paleo-Pacific influence. Unlike the ammonite assemblages of the same age of the Khatanga depression, the subzones cannot be distinguished in the Sibiricus and Kochi Zones of the Lower Lena due to differences in the stratigraphic ranges of some ammonite species. Hectoroceras in the Kochi zone is very rare and Praetollia predominates here.     

Correlation between geological units and mineralized zones using fractal modeling in Zarshuran gold deposit (NW Iran)

In this study, the relation between ore grade and geological characteristic has been studied as a principle and also important conceptual in Zarshuran gold deposit in NW Iran. The main subject in this study was identifying a correlation among the ore grade populations and rock types which could be used in other steps of local estimation in the deposit concentration–number (C–N) fractal model and logratio matrix. The C–N log–log plot reveals six geochemical zones defined by Au?<?0.02 ppm as non-mineralized zone and Au?>?0.02 ppm as mineralized zones. According to geological logging and field geology inspection, black gauge, jasperoid, fault gauge and breccia, and carbonaceous rocks are considered as main rock types which contain major Au mineralized zones. The correlation between geological and fractal modeling by logratio matrix shows that there is a good correlation between geological assumed host rocks and C–N fractal modeling. Black gauge rock type with 93.48 % of overall accuracy shows a significant correlation with supergene zone of fractal model, and jasperoid with 92.5 % and carbonaceous rock type with 52.90 % have a decent correlation with highly and lowly mineralized zone of fractal model relatively. Black gauge, jasperoid, and fault gauge and breccia have an approximately near cooperation in this zone for mineralization.     

Modelling the effects of porous media deformation on the propagation of water-table waves in a sandy unconfined aquifer

This paper examines the influence of porous media deformation on water-table wave dispersion in an unconfined aquifer using a numerical model which couples Richards’ equation to the poro-elastic model. The study was motivated by the findings of Shoushtari et al. (J Hydrol 533:412–440, 2016) who were unable to reproduce the observed wave dispersion in their sand flume data with either numerical Richards’ equation models (assuming rigid porous media) or existing analytic solutions. The water-table wave dispersion is quantified via the complex wave number extracted from the predicted amplitude and phase profiles. A sensitivity analysis was performed to establish the influence of the main parameters in the poro-elastic model, namely Young’s modulus (E) and Poisson’s ratio (ν). For a short oscillation period (T?=?16.4 s), the phase lag increase rate (k _i) is sensitive to the chosen values of E and ν, demonstrating an inverse relationship with both parameters. Changes in the amplitude decay rate (k _r), however, were negligible. For a longer oscillation period (T?=?908.6 s), variations in the values of E and ν resulted in only small changes in both k _r and k _i. In both the short and long period cases, the poro-elastic model is unable to reproduce the observed wave dispersion in the existing laboratory data. Hence porous media deformation cannot explain the additional energy dissipation in the laboratory data. Shoushtari SMH, Cartwright N, Perrochet P, Nielsen P (2016) The effects of oscillation period on groundwater wave dispersion in a sandy unconfined aquifer: sand flume experiments and modelling. J Hydrol 533:412–440.