Impact of rock-water interactions and recharge on water resources quality of the Agadir-Essaouira basin,southwestern Morocco

The Agadir-Essaouira area in the occidental High Atlas Mountains of Morocco is characterized by a semi-arid climate. The scarcity and quality of water resources, exacerbated by long drought periods, constitute a major problem for a sustainable development of this region. Groundwater resources of carbonate units within Jurassic and Cretaceous aquifers are requested for drinking and irrigation purposes. In this study, we collected 84 samples from wells, boreholes, springs, and rivers. Hydrochemical and isotopic data were used to examine the mineralization and origin of water, which control groundwater quality. The chemical composition of water seems to be controlled by water-rock interactions, such as dissolution of carbonates (calcite and dolomite), weathering of gypsum, as well as ion exchange processes, which explain the observed variability. Stable isotopes results show that groundwater from the mainly marly Cretaceous aquifer are submitted to an evaporation effect, while samples from the chiefly calcareous Jurassic aquifer indicate a meteoric origin, due to a rapid infiltration of recharge runoff through the karstic outcrops. The low values of δ¹⁸O and δ²H suggest a local recharge from areas with elevations ranging from 400 to 1200 m for the Cretaceous aquifer and from 800 to 1500 m for the Jurassic units.     

Crustal structure of the Southwest Subbasin,South China Sea,from wide-angle seismic tomography and seismic reflection imaging

The Southwest Subbasin (SWSB) is an abyssal subbasin in the South China Sea (SCS), with many debates on its neotectonic process and crustal structure. Using two-dimensional seismic tomography in the SWSB, we derived a detailed P-wave velocity model of the basin area and the northern margin. The entire profile is approximately 311-km-long and consists of twelve oceanic bottom seismometers (OBSs). The average thickness of the crust beneath the basin is 5.3 km, and the Moho interface is relatively flat (10–12 km). No high velocity bodies are observed, and only two thin high-velocity structures (~7.3 km/s) in the layer 3 are identified beneath the northern continent-ocean transition (COT) and the extinct spreading center. By analyzing the P-wave velocity model, we believe that the crust of the basin is a typical oceanic crust. Combined with the high resolution multi-channel seismic profile (MCS), we conclude that the profile shows asymmetric structural characteristics in the basin area. The continental margin also shows asymmetric crust between the north and south sides, which may be related to the large scale detachment fault that has developed in the southern margin. The magma supply decreased as the expansion of the SWSB from the east to the west.     

Temporal Variability of Total Ozone in the Asian Region Inferred from Ground-Based and Satellite Measurement Data

This paper reports investigation data on the temporal variability of total ozone content (TOC) in the Central Asian and Tibet Plateau mountain regions obtained by conventional methods, as well as by spectral, cross-wavelet, and composite analyses. The data of ground-based observation stations located at Huang He, Kunming, and Lake Issyk-Kul, along with the satellite data obtained at SBUV/SBUV2 (SBUV merged total and profile ozone data, Version 8.6) for 1980–2013 and OMI (Ozone Monitoring Instrument) and TOU (Total Ozone Unit) for 2009–2013 have been used. The average relative deviation from the SBUV/SBUV2 data is less than 1% in Kunming and Issyk-Kul for the period of 1980–2013, while the Huang He Station is characterized by an excess of the satellite data over the ground-based information at an average deviation of 2%. According to the Fourier analysis results, the distribution of amplitudes and the periods of TOC oscillations within a range of over 14 months is similar for all series analyzed. Meanwhile, according to the cross-wavelet and composite analyses results, the phase relationships between the series may considerably differ, especially in the periods of 5–7 years. The phase of quasi-decennial oscillations in the Kunming Station is close to the 11-year oscillations of the solar cycle, while in the Huang He and Issyk-Kul stations the TOC variations go ahead of the solar cycle.     

Coseismic Fault Slip of the September 16, 2015 <Emphasis Type="Italic">Mw</Emphasis> 8.3 Illapel,Chile Earthquake Estimated from InSAR Data

The complete surface deformation of 2015 Mw 8.3 Illapel, Chile earthquake is obtained using SAR interferograms obtained for descending and ascending Sentinel-1 orbits. We find that the Illapel event is predominantly thrust, as expected for an earthquake on the interface between the Nazca and South America plates, with a slight right-lateral strike slip component. The maximum thrust-slip and right-lateral strike slip reach 8.3 and 1.5 m, respectively, both located at a depth of 8 km, northwest to the epicenter. The total estimated seismic moment is 3.28 × 10²¹ N.m, corresponding to a moment magnitude Mw 8.27. In our model, the rupture breaks all the way up to the sea-floor at the trench, which is consistent with the destructive tsunami following the earthquake. We also find the slip distribution correlates closely with previous estimates of interseismic locking distribution. We argue that positive coulomb stress changes caused by the Illapel earthquake may favor earthquakes on the extensional faults in this area. Finally, based on our inferred coseismic slip model and coulomb stress calculation, we envision that the subduction interface that last slipped in the 1922 Mw 8.4 Vallenar earthquake might be near the upper end of its seismic quiescence, and the earthquake potential in this region is urgent.     

A Correction on Coastal Heads for Groundwater Flow Models

We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)h_s ? B/2α, in which h_s is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models.     

The Seismotectonics of the Po Plain (Northern Italy): Tectonic Diversity in a Blind Faulting Domain

We present a systematic and updated overview of a seismotectonic model for the Po Plain (northern Italy). This flat and apparently quiet tectonic domain is, in fact, rather active as it comprises the shortened foreland and foredeep of both the Southern Alps and the Northern Apennines. Assessing its seismic hazard is crucial due to the concentration of population, industrial activities, and critical infrastructures, but it is also complicated because (a) the region is geologically very diverse, and (b) nearly all potential seismogenic faults are buried beneath a thick blanket of Pliocene–Pleistocene sediments, and thus can be investigated only indirectly. Identifying and parameterizing the potential seismogenic faults of the Po Plain requires proper consideration of their depth, geometry, kinematics, earthquake potential and location with respect to the two confronting orogens. To this end, we subdivided them into four main, homogeneous groups. Over the past 15 years we developed new strategies for coping with this diversity, resorting to different data and modeling approaches as required by each individual fault group. The most significant faults occur beneath the thrust fronts of the Ferrara-Romagna and Emilia arcs, which correspond to the most advanced and buried portions of the Northern Apennines and were the locus of the destructive May 2012 earthquake sequence. The largest known Po Plain earthquake, however, occurred on an elusive reactivated fault cutting the Alpine foreland south of Verona. Significant earthquakes are expected to be generated also by a set of transverse structures segmenting the thrust system, and by the deeper ramps of the Apennines thrusts. The new dataset is intended to be included in the next version of the Database of Individual Seismogenic Sources (DISS; http://diss.rm.ingv.it/diss/, version 3.2.0, developed and maintained by INGV) to improve completeness of potential sources for seismic hazard assessment.     

Large Earthquakes at the Ibero-Maghrebian Region: Basis for an EEWS

Large earthquakes (M_w > 6, I_max > VIII) occur at the Ibero-Maghrebian region, extending from a point (12ºW) southwest of Cape St. Vincent to Tunisia, with different characteristics depending on their location, which cause considerable damage and casualties. Seismic activity at this region is associated with the boundary between the lithospheric plates of Eurasia and Africa, which extends from the Azores Islands to Tunisia. The boundary at Cape St. Vincent, which has a clear oceanic nature in the westernmost part, experiences a transition from an oceanic to a continental boundary, with the interaction of the southern border of the Iberian Peninsula, the northern border of Africa, and the Alboran basin between them, corresponding to a wide area of deformation. Further to the east, the plate boundary recovers its oceanic nature following the northern coast of Algeria and Tunisia. The region has been divided into four zones with different seismic characteristics. From west to east, large earthquake occurrence, focal depth, total seismic moment tensor, and average seismic slip velocities for each zone along the region show the differences in seismic release of deformation. This must be taken into account in developing an EEWS for the region.     

Spatial Periodic Boundary Condition for MODFLOW

Small‐scale hyporheic zone (HZ) models often use a spatial periodic boundary (SPB) pair to simulate an infinite repetition of bedforms. SPB's are common features of commercially available multiphysics modeling packages. MODFLOW's lack of this boundary type has precluded it from being effectively utilized in this area of HZ research. We present a method to implement the SPB in MODFLOW by development of the appropriate block‐centered finite‐difference expressions. The implementation is analogous to MODFLOW's general head boundary package. The difference is that the terms on the right hand side of the solution equations must be updated with each iteration. Consequently, models that implement the SPB converge best with solvers that perform both inner and outer iterations. The correct functioning of the SPB condition in MODFLOW is verified by two examples. This boundary condition allows users to build HZ‐bedform models in MODFLOW, facilitating further research using related codes such as MT3DMS and PHT3D.     

An Overview of the 28 October 2012 M<Subscript>w</Subscript> 7.7 Earthquake in Haida Gwaii,Canada: A Tsunamigenic Thrust Event Along a Predominantly Strike-Slip Margin

The boundary between the Pacific and North America plates along Canada’s west coast is one of the most seismically active regions of Canada, and is where Canada’s two largest instrumentally recorded earthquakes have occurred. Although this is a predominantly strike-slip transform fault boundary, there is a component of oblique convergence between the Pacific and North America plates off Haida Gwaii. The 2012 M_w 7.7 Haida Gwaii earthquake was a thrust event that generated a tsunami with significant run up of over 7 m in several inlets on the west coast of Moresby Island (several over 6 m, with a maximum of 13 m). Damage from this earthquake and tsunami was minor due to the lack of population and vulnerable structures on this coast.     

Impact of Three-Parameter Weibull Models in Probabilistic Assessment of Earthquake Hazards

This paper investigates the suitability of a three-parameter (scale, shape, and location) Weibull distribution in probabilistic assessment of earthquake hazards. The performance is also compared with two other popular models from same Weibull family, namely the two-parameter Weibull model and the inverse Weibull model. A complete and homogeneous earthquake catalog (Yadav et al. in Pure Appl Geophys 167:1331–1342, 2010) of 20 events (M ≥ 7.0), spanning the period 1846 to 1995 from north–east India and its surrounding region (20°–32°N and 87°–100°E), is used to perform this study. The model parameters are initially estimated from graphical plots and later confirmed from statistical estimations such as maximum likelihood estimation (MLE) and method of moments (MoM). The asymptotic variance–covariance matrix for the MLE estimated parameters is further calculated on the basis of the Fisher information matrix (FIM). The model suitability is appraised using different statistical goodness-of-fit tests. For the study area, the estimated conditional probability for an earthquake within a decade comes out to be very high (≥0.90) for an elapsed time of 18 years (i.e., 2013). The study also reveals that the use of location parameter provides more flexibility to the three-parameter Weibull model in comparison to the two-parameter Weibull model. Therefore, it is suggested that three-parameter Weibull model has high importance in empirical modeling of earthquake recurrence and seismic hazard assessment.