Enhanced identification of a hydrologic model using streamflow and satellite water storage data: A multicriteria sensitivity analysis and optimization approach

Hydrologic model development and calibration have continued in most cases to focus only on accurately reproducing streamflows. However, complex models, for example, the so‐called physically based models, possess large degrees of freedom that, if not constrained properly, may lead to poor model performance when used for prediction. We argue that constraining a model to represent streamflow, which is an integrated resultant of many factors across the watershed, is necessary but by no means sufficient to develop a high‐fidelity model. To address this problem, we develop a framework to utilize the Gravity Recovery and Climate Experiment's (GRACE) total water storage anomaly data as a supplement to streamflows for model calibration, in a multiobjective setting. The VARS method (Variogram Analysis of Response Surfaces) for global sensitivity analysis is used to understand the model behaviour with respect to streamflow and GRACE data, and the BORG multiobjective optimization method is applied for model calibration. Two subbasins of the Saskatchewan River Basin in Western Canada are used as a case study. Results show that the developed framework is superior to the conventional approach of calibration only to streamflows, even when multiple streamflow‐based error functions are simultaneously minimized. It is shown that a range of (possibly false) system trajectories in state variable space can lead to similar (acceptable) model responses. This observation has significant implications for land‐surface and hydrologic model development and, if not addressed properly, may undermine the credibility of the model in prediction. The framework effectively constrains the model behaviour (by constraining posterior parameter space) and results in more credible representation of hydrology across the watershed.     

Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.     

A method to predict the group fissuring and faulting caused by regional groundwater decline

Ground fissuring is a recurrent problem in many countries where water extraction surpasses the natural recharge of aquifers. Due to differential settlement, the soil layer undergoes deformation and cracks with serious consequences for civil infrastructure. Here, we propose an approximate analysis of the fissuring process that can be used to predict the location of cracks, which increasingly affect some middle- and large-sized cities in the world. For that purpose, the ground loss theory is applied to sediments overlying a sinusoidal-shaped graben. This analysis shows the existence of a tensile zone at the border of the graben with maximal values on its shoulder where tension cracks are more likely to appear. It also shows that soil deformation under differential settlement may evolve into ground faulting if water withdrawal continues. Finally, when a crack has completely developed, the tensile zone shifts towards the center of the graben, creating a new area for potential cracking and faulting.     

Are red tides affecting economically the commercialization of the Galician (NW Spain) mussel farming?

Harmful algal blooms (HABs), more specifically red tides, are among the most critical environmental factors affecting mussel cultivation in Galicia (NW Spain), and they often have been blamed for economic losses for producers. This statement is based on the correlation between days of closure of the production areas and unsold product. The present article shows that such a statement is not always correct, at least in the case of Galician mussel farming, because red tides only cause losses to producers under specific circumstances that arise from the impossibility of placing their product in the market. In addition, this article reveals the importance of finding organizational solutions within the framework of the production sector to counteract the impact of this type of phenomenon.     

The role of bathymetry and directional wave conditions on observed crescentic bar dynamics

Nearshore sandbars are important features in the surf zone of many beaches because they strongly influence the mean circulation and evolving morphology. Due to variations in wave conditions, sandbars can experience cross-shore migration and vary in shape from alongshore uniform (shore-parallel) to alongshore rhythmic (crescentic). Sandbar dynamics have been studied extensively, but existing observational studies usually do not quantify the processes leading to crescentic bar formation and straightening. This study analyses the dynamics of crescentic bar events at the fetch-limited beach of Castelldefels (northwestern Mediterranean Sea, Spain) using 7.5 years of hourly time-exposure video images and detailed wave conditions. The results show that, despite the generally calm wave conditions, the sandbars were very dynamic in the cross-shore and longshore directions. They often migrated rapidly offshore during storms (up to 70 m in one day) and more slowly onshore during post-storm conditions. Crescentic bars were often present at the study site (48% of the time), but only when the sandbar was at least 10 m from the shoreline. They displayed a large variability in wavelengths (100–700 m), alongshore migration speeds (0–50 m/day) and cross-shore amplitudes (5–20 m). Wavelengths increased for larger bar–shoreline distances and the alongshore migration speeds were strongly correlated with the alongshore component of the radiation stresses. Crescentic patterns typically developed during low–medium energetic waves with limited obliquity ( $\theta \lesssim 20$° at 10 m depth), while bar straightening occurred during medium–high energetic waves with strong oblique angles of incidence ( $\theta \gtrsim 15$°). Overall, this study provides further proof for the important role of wave direction in crescentic bar dynamics and highlights the strong dependence of crescentic bar development on the initial bathymetric configuration.     

Classification of seismic signals at Villarrica volcano (Chile) using neural networks and genetic algorithms

Each volcano has its own unique seismic activity. The aim of this work is to construct a system able to classify seismic signals for the Villarrica volcano, one of the most active volcanoes in South America. Since seismic signals are the result of particular processes inside the volcano's structure, they can be used to forecast volcanic activity. This paper describes the different kinds of seismic signals recorded at the Villarrica volcano and their significance. Three kind of signals were considered as most representative of this volcano's activity: the long-period, the tremor, and the energetic tremor signals. A classifier is implemented to read the seismic registers at 30-second intervals, extract the most relevant features of each interval, and classify them into one of the three kinds of signals considered as most representative of this particular volcano. To do so, 1033 different kinds of 30-s signals were extracted and classified by a human expert. A feature extraction process was applied to obtain the main characteristics of each of them. This process was developed using criteria which have been shown by others to effectively classify seismic signals, based on the experience of a human expert. The classifier was implemented with a Multi-Layer Perceptron (MLP) artificial neural network whose architecture and training process were optimized by means of a genetic algorithm. This technique searched for the most adequate MLP configuration to improve the classification performance, optimizing the number of hidden neurons, the transfer functions of the neurons, and the training algorithm. The optimization process also performed a feature selection to reduce the number of signal features, optimizing the number of network inputs. The results show that the optimized classifier reaches more than 93% exactitude. identifying the signals of each kind. The amplitude of the signals is the most important feature for its classification, followed by its frequency content. The described methodology can be used to classify more seismic signals to improve the study of the activity of this volcano or to extend the study to other active volcanoes of the region.     

Cretaceous deformation of the southern Central Andes: synorogenic growth strata in the Neuquén Group (35° 30′–37° S)

The Neuquén Group is an Upper Cretaceous continental sedimentary unit exhumed during the latest Miocene contractional phase occurred in the southern Central Andes, allowing a direct field observation and study of the depositional geometries. The identification of growth strata on these units surrounding the structures of the frontal parts of the Andes, sedimentological analyses and U–Pb dating of detrital components, allowed the definition of a synorogenic unit that coexisted with the uplift of the early Andean orogen since ca. 100 Ma, maximum age obtained in this work, compatible with previous assignments and constrained in the top by the deposition of the Malargüe Group, in the Maastrichtian (ca. 72 Ma). The definition of a wedge top area in this foreland basin system, where growth strata were described, permitted to identify a Late Cretaceous orogenic front and foredeep area, whose location and amplitude contrast with previous hypotheses. This wedge top area was mostly fed from the paleo‐Andes with small populations coming from sources in the cratonic area that are interpreted as a recycling in Jurassic and Lower Cretaceous sections, which contrasts with other analyses performed at the foredeep zone that have mixed sources. In particular, Permian sources are interpreted as coming directly from the cores of the basement structures, where Neopaleozoic sections are exposed, next to the synorogenic sedimentation, implying a strong incision in Late Cretaceous times with an exhumation structural level similar to the present. The maximum recognised advance for this Late Cretaceous deformation in the study area is approximately 500 km east of the Pacific trench, which constitutes an anomaly compared with neighbour segments where Late Cretaceous deformations were found considerably retracted. The geodynamic context of the sedimentation of this unit is interpreted as produced under the westward fast moving of South America, colliding with two consecutive mid‐ocean ridges during a period of important plate reorganisation. The subduction of young, anhydrous, buoyant lithosphere would have produced changes in the subduction geometry, reflected first by an arc waning/gap and subsequently by an arc migration that coexisted with synorogenic sedimentation. These magmatic and deformational processes would be the product of a shallow subduction regime, following previous proposals, which occurred in Late Cretaceous times, synchronous to the sedimentation of the Neuquén Group.     

Digital mapping of coastal boulders – high-resolution data acquisition to infer past and recent transport dynamics

Coastal boulder fields provide clues to long-term frequency-magnitude patterns of coastal flooding events and have the potential to play an important role in coastal hazard assessment. Mapping boulders in the field is time and labour-intensive, and work on intertidal reef platforms, as in the present study, is physically challenging. By addressing coastal scientists who are not specialists in remote sensing, this contribution reports on the possibilities and limitations of digital applications in boulder mapping in Eastern Samar, Philippines, where recent supertyphoons Haiyan and Hagupit induced high waves, coastal flooding and boulder transport. It is demonstrated how satellite imagery of sub-metre resolution (from Pléiades and WorldView-3 imagery) enables efficient analysis of transport vectors and distances of larger boulders, reflecting variation in latitudes of both typhoon tracks and approaching angles of typhoon-generated waves. During the investigated events, boulders with a-axes of up to 8 m were clearly identified to have been shifted for up to 32 m, mostly along the seaward margin of the boulder field. It is, however, hard to keep track of smaller boulders, and the length of a-axes and b-axes including their orientation is often impossible to map with sufficient accuracy. Orthophotographs and digital surface models created through the application of an unmanned aerial vehicle and the ‘Structure from Motion’ technique provide ultra-high-resolution data, and have the potential to not only improve the results of satellite image analysis, but also those from field mapping and may significantly reduce overall time in the field. Orthophotographs permit unequivocal mapping of a-axes and b-axes including their orientation, while precise values for c-axes can be derived from the respective digital surface models. Volume of boulders is best inferred from boulder-specific Structure from Motion-based three-dimensional models. Battery power, flight speed and altitude determine the limits of the area covered, while patches shielded by the boulders are difficult to resolve. For some tasks, field mapping remains mandatory and cannot be replaced by currently available remote sensing tools: for example, sampling for rock type, density and age dating, recording of lithological separation of boulders from the underlying geological unit and of geomorphic features on a millimetre to decimetre-scale, or documentation of fine-grained sediment transport in between the boulders in supratidal settings. In terms of future events, the digital products presented here will provide a valuable reference to track boulder transport on a centimetre to decimetre-scale and to better understand the hydrodynamics of extreme-wave events on a fringing reef coastline.     

The Ushnus in the astronomy of the Inca culture

Inferences are made about the relationship that existed between the Ushnus, pyramid-shaped, terraced structures used by the Incas in the most important ceremonies of the Tawantinsuyo, and Inka Astronomy. We draw attention to Ayni, Kawsaypacha, Duality, and Tinkuy principles, multidimensional codes of conduct and wisdom that are at the root of the Andean cosmovision and on their perception of the world and the Cosmos. These principles, examined as postulates, allow to elaborate axiomatic propositions to identify the Ushnus with ancient Astronomy practices. In a complementary statement, starting from a bi-conditional proposition, we may infer through reciprocal corollaries that the Inka earliest roots to a holistic learning and educational ambient in the Tawantinsuyo was not elitist, instead it was based on a epistemological construct that differs from the corresponding Western educational ambients. An epistemological and cognitive approach allows to identify an ancient elaborate process of knowledge construction, based on the four fundamental principles, corresponding to different levels of assimilation and comprehension. As a complementary aspect, we identify some of the most preserved Ushnus of the Inka “Empire.” Then we complement this contribution with a broader interpretation for the Ushnus.     

The Loncopué Trough: A Cenozoic basin produced by extension in the southern Central Andes

The Loncopué Trough is located in the hinterland Andean zone between 36°30′ and 39°S. It constitutes a topographic low bounded by normal faults and filled by lavas and sediments less than 5 Ma old. Reprocessed seismic lines show wedge-like depocenters up to 1700 m deep associated with high-angle faults, correlated with the 27–17 Ma Cura Mallín basin deposits, and buried beneath Pliocene to Quaternary successions and Late Miocene foreland sequences. The southern Central Andes seem to have been under extension in the hinterland zone some 27 Ma ago and again at approximately 5 Ma ago. This last extensional period could have been the product of slab steepening after a shallow subduction cycle in the area, although other alternatives are discussed. Orogenic wedge topography, altered by the first extensional stage in the area, was recovered through Late Miocene inversion, and was associated with foreland sequences. However, since the last extension (<5 Ma) the Andes have not recovered their characteristic contractional behavior that controlled past orogenic growth.