Atlas of susceptibility to pollution in marinas. Application to the Spanish coast

An atlas of susceptibility to pollution of 320 Spanish marinas is provided. Susceptibility is assessed through a simple, fast and low cost empirical method estimating the flushing capacity of marinas. The Complexity Tidal Range Index (CTRI) was selected among eleven empirical methods. The CTRI method was selected by means of statistical analyses because: it contributes to explain the system's variance; it is highly correlated to numerical model results; and, it is sensitive to marinas' location and typology. The process of implementation to the Spanish coast confirmed its usefulness, versatility and adaptability as a tool for the environmental management of marinas worldwide. The atlas of susceptibility, assessed through CTRI values, is an appropriate instrument to prioritize environmental and planning strategies at a regional scale.     

Prevailing-frequency approximation of the coupling ray theory along the SH and SV reference rays in a heterogeneous generally anisotropic medium which is approximately uniaxial

The behaviour of the actual polarization of an electromagnetic wave or elastic S–wave is described by the coupling ray theory, which represents the generalization of both the zero–order isotropic and anisotropic ray theories and provides continuous transition between them. The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. In a generally anisotropic or bianisotropic medium, the actual wave paths may be approximated by the anisotropic–ray–theory rays if these rays behave reasonably. In an approximately uniaxial (approximately transversely isotropic) anisotropic medium, we can define and trace the SH (ordinary) and SV (extraordinary) reference rays, and use them as reference rays for the prevailing–frequency approximation of the coupling ray theory. In both cases, i.e. for the anisotropic–ray–theory rays or the SH and SV reference rays, we have two sets of reference rays. We thus obtain two arrivals along each reference ray of the first set and have to select the correct one. Analogously, we obtain two arrivals along each reference ray of the second set and have to select the correct one. In this paper, we suggest the way of selecting the correct arrivals. We then demonstrate the accuracy of the resulting prevailing–frequency approximation of the coupling ray theory using elastic S waves along the SH and SV reference rays in four different approximately uniaxial (approximately transversely isotropic) velocity models.     

The use of 3D modeling for the prediction of the seismic demands on the gravity dams

Seismic behavior of gravity dams has long been evaluated using a representative two‐dimensional (2D) system. Formulated for the gravity dams built in wide canyons, the assumption is nevertheless utilized extensively for almost all concrete dams due to the established procedures as well as the expected computational costs of a three‐dimensional model. However, a significant number of roller‐compacted concrete dams, characterized as such systems, do not conform to the basic assumptions of these methods by violating the conditions on canyon dimensions and joint‐spacing/details. Based on the premise that the 2D modeling assumption is overstretched for practical purposes in a variety of settings, the purpose of this study is to critically evaluate the use of 2D modeling for the prediction of the seismic demands on these systems. Using a rigorous soil–structure interaction approach, the difference between the two and three‐dimensional response for gravity dams was investigated first in the frequency domain for a range of canyon widths and foundation to dam moduli ratios. Then, the time domain differences between the crest displacements and the maximum principal stress were obtained using 70 different ground motions in order to show the possible bias introduced into the analysis results due to the modeling approach. The results of the study show that even for relatively wide canyons, the 2D analysis can lead to misleading predictions. Copyright © 2017 John Wiley & Sons, Ltd.     

Geophysical evidence of the Eastern Marginal Fault of the Cheb Basin (Czech Republic)

The western part of the Bohemian Massif hosts an intersection of two regional fault zones, the SW-NE trending Oh?e/Eger Graben and the NNW-SSE trending Mariánské Lázně Fault, which has been reactivated several times in the geological history and controlled the formation of the Tertiary Cheb Basin. The broader area of the Cheb Basin is also related to permanent seismic activity of M_L 3+ earthquake swarms. The Eastern Marginal Fault of the Cheb Basin (northern segment of the Mariánské Lázně Fault) separates the basin sediments and underlying granites in the SW from the Kru?né Hory/Erzgebirge Mts. crystalline unit in the NE. We describe a detailed geophysical survey targeted to locating the Eastern Marginal Fault and determining its geometry in the depth. The survey was conducted at the Kopanina site near the Nový Kostel focal zone, which shows the strongest seismic activity of the whole Western Bohemia earthquake swarm region. Complex geophysical survey included gravimetry, electrical resistivity tomography, audiomagnetotellurics and seismic refraction. We found that the rocks within the Eastern Marginal Fault show low resistivity, low seismic velocity and density, which indicates their deep fracturing, weathering and higher water content. The dip of the fault in shallow depths is about 60° towards SW. At greater depths, the slope turns to subvertical with dip angle of about 80°. Results of geoelectrical methods show blocky fabric of the Cheb Basin and deep weathering of the granite bedrock, which is consistent with geologic models based on borehole surveys.     

Cause of scale inconsistencies in DORIS time series

In this paper we analyze the scale of the DORIS (Doppler Orbitography and Radiopositioning Integrated by Satellite) solutions with respect to DORIS extension of the International Terrestrial Reference Frame (ITRF) for Precise Orbit Determination DPOD2014. The main goal is to explain the scale inconsistencies and to find the optimal solution reaching low-biased and consistent scale time series. Our analysis profits from 4 different strategies based only on the Geodetic Observatory Pecný analysis center solution, using DORIS exchange format data 2.2. A difference in the sequence of the solutions directly corresponds to one of the changes in the solution settings: data elevation dependent weighting, application of data validity indicators and application of phase center - reference point correction. We process multi-satellite and single-satellite solutions for the time period 2011.0–2017.0. Our analysis examines scale inconsistency issues in 2011/2012 and in 2015. The scale increment in 2011/2012 is explained as a result of the concurrence of changes in satellite constellation and change in the provider data validity standards for Cryosat-2 and Jason-2 satellites. The scale increment in 2015 is explained as the effect of change in the standards for phase center - reference center corrections for Saral, Jason-2 and Cryosat-2 satellites. Moreover, comparing the solutions with and without elevation dependent data downweighting using the same elevation cutoff (10°), we found a significant reduction of scale bias and scale variation applying the data downweighting. The data downweighting improved also the station positioning repeatability. We demonstrate that the solution, which is completely free from the additional data associated with observations in DORIS exchange format 2.2 and includes the data downweighting law, eventuates in a consistent scale time series with the lowest offset with respect to DPOD2014 (version 1.0) (12.7 ± 2.3 mm for 2011.0–2017.0). The only remaining scale issue is the part of 2011/2012 increment of around 5 mm, explained by a change in the DORIS satellite constellation.     

Deterioration processes affecting prehistoric rock art engravings in granite in NW Spain

Prehistoric rock art sites are endangered despite conservation efforts. The lack of scientific documentation regarding weathering agents affecting rock art and the absence of specific diagnostic protocols hinder the development of conservation strategies. The aim of this research was to investigate active deterioration processes in a granite petroglyph site located in Mougás (Galicia, NW Spain) by characterizing the granite, conducting a geotechnical study of the outcrop and describing and analysing the main weathering processes. Two main deterioration factors were identified. First, water favours block disjunction at the massif scale and causes pitting and surface erosion at the millimetre scale that affects the readability of the engravings. Second, high temperatures associated with wildfires cause mineral transformations that increase the susceptibility of the rock to weathering. Identifying deterioration factors is a first step in developing appropriate preventive conservation measures, which should aim to reduce rock contact time with water (technically affordable in the short term) and to reduce the probability of wildfire occurrence (technically more complex and possibly with longer‐term results). Copyright © 2018 John Wiley & Sons, Ltd.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

SlugIn 1.0: A Free Tool for Automated Slug Test Analysis

The correct characterization of aquifer parameters is essential for water‐supply and water‐quality investigations. Slug tests are widely used for these purposes. While free software is available to interpret slug tests, some codes are not user‐friendly, or do not include a wide range of methods to interpret the results, or do not include automatic, inverse solutions to the test data. The private sector has also generated several good programs to interpret slug test data, but they are not free of charge. The computer program SlugIn 1.0 is available online for free download, and is demonstrated to aid in the analysis of slug tests to estimate hydraulic parameters. The program provides an easy‐to‐use Graphical User Interface. SlugIn 1.0 incorporates automated parameter estimation and facilitates the visualization of several interpretations of the same test. It incorporates solutions for confined and unconfined aquifers, partially penetrating wells, skin effects, shape factor, anisotropy, high hydraulic conductivity formations and the Mace test for large‐diameter wells. It is available in English and Spanish and can be downloaded from the web site of the Geological Survey of Spain. Two field examples are presented to illustrate how the software operates.     

Extracting water-related features using reflectance data and principal component analysis of Landsat images

This study aimed to map water features using a Landsat image rather than traditional land cover. We involved the original bands, spectral indices and principal components (PCs) of a principal component analysis (PCA) as input data, and performed random forest (RF) and support vector machine (SVM) classification with water, saturated soil and non-water categories. The aim was to compare the efficiency of the results based on various input data. Original bands provided 93% overall accuracy (OA) and bands 4–5–7 were the most informative in this analysis. Except for MNDWI (modified normalized differenced water index, with 98% OA), the performance of all water indices was between 60 and 70% (OA). The PCA-based approach conducted on the original bands resulted in the most accurate identification of all classes (with only 1% error in the case of water bodies). We therefore show that both water bodies and saturated soils can be identified successfully using this approach.     

Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.