Damage assessment of the 2015 Mw 8.3 Illapel earthquake in the North-Central Chile

Natural Hazards - Destructive megathrust earthquakes, such as the 2015 Mw 8.3 Illapel event, frequently affect Chile. In this study, we assess the damage of the 2015 Illapel Earthquake in the...     

A review of the geodynamic constraints on the development and evolution of geothermal systems in the Central Andean Volcanic Zone (18–28°Lat.S)

ABSTRACT

Geothermal exploration in the Central Andean Volcanic Zone (CAVZ) focuses on heat capacities of known geothermal systems, yet the role of faults, veins, fractures and folds on the evolution and migration of fluids is far from complete. Here, we present a compilation of He and Sr isotope data and newly generated structural maps to examine if particular tectonic configurations are associated with fluids migrating from different crustal levels. Accordingly, we defined three tectono-geothermal environments (T1–T3) depicting specific structural arrangements and spatial relation with geothermal and volcanic manifestations. T1 is dominated by left-lateral strike-slip NW-striking faults, and geothermal and volcanic manifestations occur along the traces of these structures. T2 is dominated by N-striking thrust faults and parallel fault-propagated folds, cut and displaced by NW-striking faults. Here, geothermal manifestations occur at fault intersections and at fold hinges. T3 is defined by left-lateral/normal NW-striking faults, with geothermal and volcanic manifestations lying along fault traces. Each tectono-geothermal environment yields distinctive isotope ratios and geothermal reservoir temperatures. T1 shows high helium and low strontium ratios, and temperatures between 220° and 310°C. T3 shows low helium and high strontium ratios and temperatures between 260° and 320°C. T2 isotope ratios fall between T1 and T3, with lowest (130°-250°C) reservoir temperatures. We argue that these particular isotope signatures are due to a structural control on reservoir location and orientation. The orientation of the fracture mesh genetically associated with each tectono-geothermal environment is a first-order control on the migration pathway of fluids. Our model shows that T1 allows fluxing of deeper fluids, T2 promotes storage and favors longer residence times and T3 enhances subvertical fluid migration. Results here help to explain the evolution of geothermal systems in a wider context, including fault systems and Sr and He isotope variations, thus providing a framework for geothermal exploration in the CAVZ.     

Adaptive POD model reduction for solute transport in heterogeneous porous media

We study the applicability of a model order reduction technique to the solution of transport of passive scalars in homogeneous and heterogeneous porous media. Transport dynamics are modeled through the advection-dispersion equation (ADE) and we employ Proper Orthogonal Decomposition (POD) as a strategy to reduce the computational burden associated with the numerical solution of the ADE. Our application of POD relies on solving the governing ADE for selected times, termed snapshots. The latter are then employed to achieve the desired model order reduction. We introduce a new technique, termed Snapshot Splitting Technique (SST), which allows enriching the dimension of the POD subspace and damping the temporal increase of the modeling error. Coupling SST with a modeling strategy based on alternating over diverse time scales the solution of the full numerical transport model to its reduced counterpart allows extending the benefit of POD over a prolonged temporal window so that the salient features of the process can be captured at a reduced computational cost. The selection of the time scales across which the solution of the full and reduced model are alternated is linked to the Péclet number (P e), representing the interplay between advective and dispersive processes taking place in the system. Thus, the method is adaptive in space and time across the heterogenous structure of the domain through the combined use of POD and SST and by way of alternating the solution of the full and reduced models. We find that the width of the time scale within which the POD-based reduced model solution provides accurate results tends to increase with decreasing P e. This suggests that the effects of local-scale dispersive processes facilitate the POD method to capture the salient features of the system dynamics embedded in the selected snapshots. Since the dimension of the reduced model is much lower than that of the full numerical model, the methodology we propose enables one to accurately simulate transport at a markedly reduced computational cost.     

Human-made disasters and economic impact for a developing economy: evidence from Brazil

Natural Hazards - This article analyzes the short-term economic impacts of the collapse of the ‘Fundão’ mining tailings dam, located in the Brazilian state of Minas Gerais. This...     

BRAMS model performance for significant event involving Brazilian platform ship

Natural Hazards - After a storm displaced the P-70 platform ship located in Guanabara Bay to the coast on January 30, 2020, a numerical investigation was carried out with the BRAMS (Brazilian...     

Synthetic scenario generation of monthly streamflows conditioned to the El Niño–Southern Oscillation: application to operation planning of hydrothermal systems

Stochastic Environmental Research and Risk Assessment - The Brazilian Interconnected Power System is hydro dominated and characterized by large reservoirs presenting multi-year regulation...     

The effectiveness of “clay” liners as basal isolation of landfills: a case study

A sodium bentonite is being used in the basal-isolation of landfills installed on strongly fractured granitic rocks of N Portugal (NW Iberia). To evaluate the performance of this clay as the ultimate impermeable basal barrier, a detailed study was carried out of the granite fracture network; the Na-bentonite was tested to characterize its mechanical and geotechnical behavior; geophysical surveys were made to search for landfill leakage and the waters found around the landfill were geochemically characterized. Fractures in the granite are usually open and devoid of mineral infillings or clay materials and, thus, permeability of the granite is very high. Dispersal of contaminants can be further enhanced by the position of the landfill on a high steep-sided ridge. Geophysical and geochemical data show evidence for seepage and dispersion of pollutants, which means that the clay failed as an ultimate isolation barrier against seepage. This investigation shows that this can be due to fracturing of the clay under load and/or to its non-homogeneous saturation and extreme shrinking character upon drying, which are accompanied by the formation of extensive cracking. Observations and experimental results suggest that the use of synthetic clays in the safe building of landfill bottom liners needs further research, and extreme care should be taken in preventing that clay water content suffers large variations after saturation, as this process considerably degrades the mechanical behavior and sealing properties of the studied Na-bentonite.     

An open source GPS multipath simulator in Matlab/Octave

Multipath is detrimental for both GPS positioning and timing applications. However, the benefits of GPS multipath for reflectometry have become increasingly clear for monitoring soil moisture, snow depth, and vegetation growth. In positioning applications, a simulator can support multipath mitigation efforts in terms of, e.g., site selection, antenna design, receiver performance assessment, and in relating different observations to a common parameterization. For reflectometry, in order to convert observed multipath parameters into useable environmental products, it is important to be able to explicitly link the GPS observables to known characteristics of the GPS receiver/antenna and the reflecting environment. Existing GPS multipath software simulators are generally not readily available for the general scientific community to use and/or modify. Here, a simulator has been implemented in Matlab/Octave and is made available as open source code. It can produce signal-to-noise ratio, carrier phase, and code pseudorange observables, based on L1 and L2 carrier frequencies and C/A, P(Y), and L2C modulations. It couples different surface and antenna types with due consideration for polarization and coherence. In addition to offering predefined material types (water, concrete, soil, etc.), it allows certain dimensional properties to be varied, such as soil moisture and snow density.     

Simulation of Geological Contacts from Interpreted Geological Model Using Multiple-Point Statistics

Applications of multiple-point statistics (mps) algorithms to large non-repetitive geological objects such as those found in mining deposits are difficult because most mps algorithms rely on pattern repetition for simulation. In many cases, an interpreted geological model built from a computer-aided design system is readily available but suffers as a training image due to the lack of patterns repetitiveness. Porphyry copper deposits and iron ore formations are good examples of such mining deposits with non-repetitive patterns. This paper presents an algorithm called contactsim that focuses on reproducing the patterns of the contacts between geological types. The algorithm learns the shapes of the lithotype contacts as interpreted by the geologist, and simulates their patterns at a later stage. Defining a zone of uncertainty around the lithological contact is a critical step in contactsim, because it defines both the zones where the simulation is performed and where the algorithm should focus to learn the transitional patterns between lithotypes. A larger zone of uncertainty results in greater variation between realizations. The definition of the uncertainty zone must take into consideration the geological understanding of the deposit, and the reliability of the contact zones. The contactsim algorithm is demonstrated on an iron ore formation.