Towards adopting open and data-driven science practices in bed form dynamics research,and some steps to this end

In recent years, open and data-driven science has fostered very important scientific breakthroughs. This study describes the challenges and opportunities for the scientific community devoted to bed form dynamics research in adopting such scientific paradigms through, for example, engineered data sharing, formal recognition of scientists who collect the data, and collaborative development of free accessible software. It highlights that once these actions are completed, the potential application of deep learning techniques could substantially improve bed form models and the scientific understanding of bed form dynamics. Likewise, it discusses the potential of Bedforms-ATM, a free available software, to standardize some bed form data analysis techniques. We propose that the technical challenges be tackled by following scholarly accepted/proposed standards (e.g. FAIR Guiding Principles, Geoscience Papers of the Future), using the body of knowledge being built on the matter by some institutions (e.g. Federation of Earth Science Information Partners), and through technical discussions at scientific meetings such as MARID. © 2020 John Wiley & Sons, Ltd.     

The impact of aquifer stratification on saltwater intrusion characteristics. Comprehensive laboratory and numerical study

Laboratory experiments and numerical simulations were utilized in this study to assess the impact of aquifer stratification on saltwater intrusion. Three homogeneous and six layered aquifers were investigated. Image processing algorithms facilitated the precise calculation of saltwater wedge toe length, width of the mixing zone, and angle of intrusion. It was concluded that the length of intrusion in stratified aquifers is predominantly a function of permeability contrast, total aquifer transmissivity and the number of heterogeneous layers, being positively correlated to all three. When a lower permeability layer overlays or underlays more permeable zones its mixing zone widens, while it becomes thinner for the higher permeability strata. The change in the width of the mixing zone (WMZ) is positively correlated to permeability contrast, while it applies to all strata irrespectively of their relative vertical position in the aquifer. Variations in the applied hydraulic head causes the transient widening of WMZ. These peak WMZ values are larger during saltwater retreat and are negatively correlated to the layer's permeability and distance from the aquifer's bottom. Moreover, steeper angles of intrusion are observed in cases where low permeability layers overlay more permeable strata, and milder ones in the inverse aquifer setups. The presence of a low permeability upper layer results in the confinement of the saltwater wedge in the lower part of the stratified aquifer. This occurs until a critical hydraulic head difference is applied to the system. This hydraulic gradient value was found to be a function of layer width and permeability contrast alike.     

Measurements of upper mantle shear wave anisotropy from a permanent network in southern Mexico

Upper mantle shear wave anisotropy under stations in southern Mexico was measured using records of SKS phases. Fast polarization directions where the Cocos plate subducts subhorizontally are oriented in the direction of the relative motion between the Cocos and North American plates, and are trench-perpendicular. This pattern is interpreted as subslab entrained flow, and is similar to that observed at the Cascadia subduction zone. Earlier studies have pointed out that both regions have in common the young age of the subducting lithosphere. Changes in the orientation of the fast axes are observed where the subducting plates change dip and/or are torn, and are thus indicative of 3-D flow around the slab edges. They are consistent with slab rollback, as previously shown by other authors. Some stations located away from the plate boundaries have their fast directions controlled by the absolute motion of the North American plate. The fast axis for station ZAIG, located in the Mesa Central, is oriented WNW-ESE and is different from all the other measurements in this study.     

Dark Matter: A Result of nonadditive gravitational forces

Experimental data in fluids suggest that nonadditive electromagnetic forces between 3 or more molecules account for the existence of critical points, triple states and phase transitions (Robles-Domínguez et al., 2007; Robles-Gutiérrez et al., 2010). Similar nonadditive forces between 3 or more molecules in the gravitational field incorporated into Newton’s universal gravitational law may also explain the existence of dark matter.     

Anthropogenic input of heavy metals in two Audouin’s gull breeding colonies

Anthropogenic activities have become an important source of heavy metals to the marine environments. Biological sentinels like seabirds’ chicks have been widely used to monitorize the levels of some heavy metals. Due to its mainly marine foraging habits, Audouin’s gull fits well for this purpose. Mercury and lead levels were measured in mantle feathers of Audouin’s gull chicks from two colonies in NE Iberian Peninsula: the Ebro Delta and the Llobregat Delta. Both are anthropized areas subject to differential pollutant inputs. Lead levels were significantly higher in the Llobregat Delta probably due to the use of leaded fuel in the nearby Barcelona airport. On the other side, mercury concentrations were higher in the Ebro Delta, in relation with the disposal of the toxic sediments at the Flix site carried down by the Ebro River. These mercury levels in the Ebro chicks reached values that have been described as toxic.     

Factors influencing magnetic susceptibility in the southern Pyrenees

We analyse paramagnetic, ferrimagnetic and frequency-dependent susceptibility in the southern Pyrenees and northern Ebro basin. The analyzed rocks show a wide range of values for the three parameters within the different lithological groups. Sandstones, shales and marls exhibit middle to high values of paramagnetic susceptibility while mudstones usually show low values. However, a general correlation can be established for the different rock types with a wide variability for each group (for example, Eocene marls show average values of paramagnetic susceptibility of 82% but within a range between 12 and 99%, or Flysch materials, with mean values of 47–52% within a range between 15 and 90%). These changes in the paramagnetic contribution to the susceptibility are due to the change in the ferrimagnetic susceptibility. This wide variability is independent of the sedimentological characteristics and can be related to the higher ferrimagnetic contribution to the susceptibility with similar ferrimagnetic mineral phase assemblages. The ferrimagnetic and superparamagnetic contributions, independent of site location with respect to individual macrostructures (thrusts and nappes), show a sharp increase at the occurrence of the macroscopic regional Pyrenean cleavage, probably relating both aspects with changes in pressure-temperature (P-T) conditions.     

Improving the calibration strategy of the physically-based model WaSiM-ETH using critical events

Abstract

The use of a physically-based hydrological model for streamflow forecasting is limited by the complexity in the model structure and the data requirements for model calibration. The calibration of such models is a difficult task, and running a complex model for a single simulation can take up to several days, depending on the simulation period and model complexity. The information contained in a time series is not uniformly distributed. Therefore, if we can find the critical events that are important for identification of model parameters, we can facilitate the calibration process. The aim of this study is to test the applicability of the Identification of Critical Events (ICE) algorithm for physically-based models and to test whether ICE algorithm-based calibration depends on any optimization algorithm. The ICE algorithm, which uses the data depth function, was used herein to identify the critical events from a time series. Low depth in multivariate data is an unusual combination and this concept was used to identify the critical events on which the model was then calibrated. The concept is demonstrated by applying the physically-based hydrological model WaSiM-ETH on the Rems catchment, Germany. The model was calibrated on the whole available data, and on critical events selected by the ICE algorithm. In both calibration cases, three different optimization algorithms, shuffled complex evolution (SCE-UA), parameter estimation (PEST) and robust parameter estimation (ROPE), were used. It was found that, for all the optimization algorithms, calibration using only critical events gave very similar performance to that using the whole time series. Hence, the ICE algorithm-based calibration is suitable for physically-based models; it does not depend much on the kind of optimization algorithm. These findings may be useful for calibrating physically-based models on much fewer data.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Singh, S.K., Liang, J.Y., and Bárdossy, A., 2012. Improving calibration strategy of physically-based model WaSiM-ETH using critical events. Hydrological Sciences Journal, 57 (8), 1487–1505.     

Long-term fluid circulation in extensional faults in the central Catalan Coastal Ranges: P–T constraints from neoformed chlorite and K-white mica

The neoformation of chlorite and K-white mica in fault rocks from two main faults of the central Catalan Coastal Ranges, the Vallès and the Hospital faults, has allowed us to constrain the P–T conditions during fault evolution using thermodynamic modeling. Crystallization of M1 and M2 muscovite and microcline occured as result of deuteric alteration during the exhumation of the pluton (290 °C > T > 370 °C) in the Permian. After that, three tectonic events have been distinguished. The first tectonic event, attributed to the Mesozoic rifting, is characterized by precipitation of M3 and M4 phengite together with chlorite and calcite C1 at temperatures between 190 and 310 °C. The second tectonic event attributed to the Paleogene compression has only been identified in the Hospital fault with precipitation of low-temperature calcite C2. The shortcut produced during inversion of the Vallès fault was probably the responsible for the lack of neoformed minerals within this fault. Finally, the third tectonic event, which is related to the Neogene extension, is characterized in the Vallès fault by a new generation of chlorite, associated with calcite C4 and laumontite, formed at temperatures between 125 and 190 °C in the absence of K-white mica. Differently, the Hospital fault is characterized by the precipitation of calcite C3 during the syn-rift stage at temperatures around 150 °C and by low-temperature fluids precipitating calcites C5, C6 and PC1 during the post-rift stage. During the two extensional events (Mesozoic and Neogene), faults acted as conduits for hot fluids producing anomalous high geothermal gradients (50 °C/km minimum).     

The Son Poc rockfall (Mallorca,Spain) on the 6th of March 2013: 3D simulation

The Son Poc rockfall took place on the 6th of March 2013 in the municipality of Bunyola, on the southern side of the Tramuntana Range (Mallorca) and after a rainy and cold period on the region. A volume of rock of 4.000 m³ was detached from the cliff crowning the peak falling down by toppling. The impact of the boulder caused its fragmentation, and numerous boulders bounced and rolled downslope with volumes from 1 to 35 m³, following two trajectories: southwest (SW) and southeast (SE). The SE trajectory, with a larger runout (376 m), reached an urban area, where some of the boulders hit the roofs and walls of nearby houses, stopping others in their gardening areas. Fortunately, no fatalities occurred despite of the presence of some people at that moment, but the event caused great concern in a region which lives from and for tourism. The Son Poc rockfall has been simulated using RocPro3D software which uses GIS technology to produce 3D rockfall trajectories lines, estimated velocity and energy of falling blocks, as well as bounce heights, impacts, and stopping points. The results are in agreement with field observations and with a very good accuracy between real and modeled outcomes.     

On the application of combined geoelectric weighted inversion in environmental exploration

Geophysical surveying methods are of great importance in environmental exploration. Inversion-based data processing methods are applied for the determination of geometrical and physical parameters of the target model. The use of this geoelectric inversion method is advantageous in environmental research where highly reliable information with large spatial resolution is required. The 2D combined geoelectric inversion (CGI) method performs more accurate parameter estimation than conventional 1D single inversion methods by efficiently decreasing the number of unknowns of the inverse problem (single means that data sets of individual vertical electric sounding stations are inverted separately). The quality improvement in parameter space is demonstrated by comparing the traditional 1D inversion procedure with a 2D series expansion-based inversion technique. The CGI method was further developed by weighting individual direct current geoelectric data sets automatically in order to improve inversion results. The new algorithm was named combined geoelectric weighted inversion, which extracts the solution by a special weighted least squares technique. It is shown that the new inversion methodology is applicable to resolve near-surface structures such as rapidly varying layer boundaries, laterally inhomogeneous formations and pinch-outs.