Geoelectrical constraints on radar probing of shallow water-saturated zones within karstified carbonates in semi-arid environments

Shallow carbonates are of utmost importance as potential sources of groundwater in karstified semi-arid terrains. Ground-Penetrating Radar (GPR) is being increasingly used as a prominent mapping tool in such environments. However, its potential in exploring and identifying shallow water-saturated zones (WSZs) in carbonates is constrained by the geoelectrical properties of carbonate soils as a function of moisture content. We report results of a case study that includes laboratory geoelectrical characterization and their comparison to in situ GPR attenuation measurements performed on Cretaceous Edwards Formation rudist mounds in central Texas, which we hypothesize as analogs for water-bearing formations in semi-arid karstified carbonate terrains. Dielectric measurements on field-collected rock samples carried out in the laboratory under controlled conditions of moisture saturation suggest that real and imaginary parts of dielectric constants of rocks with higher porosity and/or permeability have steeper dependence on pore moisture content; they produce better dielectric contrasts but allow shallower penetration. Our analyses suggest that within carbonates, dielectric contrasts improve with decrease in sounding frequency and/or increase in moisture content; and the relationship between dielectric permittivity and moisture content may be represented by 3rd order polynomial equations. GPR surveys using a wide-band 400 MHz antenna reveal subsurface mound morphologies with heights of ～ 1–2 m and basal diameters of ～ 8–10 m resembling outcrop analogs. Each mound appears to be composed of smaller amalgamated lithounits that seem geoelectrically similar. Amplitudes decays of the backscattered radar signal correlate to moisture distribution. Measuring the differences in signal attenuation allows differentiation between saturated and non-saturated zones. Velocity analyses of GPR profiles enable estimation of moisture distribution in the vicinity of the mounds. Optimal delineation and production of high-resolution GPR data up to a depth of ～ 10 m were observed for a sounding frequency of ～ 250 MHz with moisture content of ～ 5% by weight. Below this moisture level, the dielectric contrast is insufficient to uniquely identify water-saturated zones from the surrounding geoelectrical context, and above it, the radar signal is substantially attenuated leading to a total inefficiency of the method.     

Velocity dispersion in fractured rocks in a wide frequency range

Experimental measurements of fracture-induced seismic waves velocity variations at frequencies ~ 1 kHz, ~ 40 kHz and ~ 1 MHz were performed directly in the field at the rocky outcrop and in the laboratory on specific rock samples collected from the outcrops. The peridotite–lherzolite outcrop appeared macroscopically uniform and contained three systems of visible parallel sub-vertical fractures. This rock has substantial bulk density and higher than average value of seismic wave velocity. The presence of fracture systems gives rise to its velocity anisotropy. The seismic waves passing through the rock fractures are subject to velocity dispersion and frequency dependent attenuation. Our data, obtained from field and laboratory measurements, were compared with theoretical model predictions. In this model we successfully used displacement discontinuity approach. For the velocity dispersion evaluation we used multi-frequency measurements. The a priori observation of orientations and densities of fracture sets allowed evaluation of their stiffness. Our approach revealed that the first arrivals of seismic waves can be used for evaluation of P-wave group velocities, the specific case, in which we expect anomalous velocity dispersion. Our observations contribute to the issue of up-scaling of well-log derived velocities in fractured rock to the scale of standard seismic exploration frequencies.     

Model track studies on fouled ballast using ground penetrating radar and multichannel analysis of surface wave

Ballast fouling is created by the breakdown of aggregates or outside contamination by coal dust from coal trains, or from soil intrusion beneath rail track. Due to ballast fouling, the conditions of rail track can be deteriorated considerably depending on the type of fouling material and the degree of fouling. So far there is no comprehensive guideline available to identify the critical degree of fouling for different types of fouling materials. This paper presents the identification of degree of fouling and types of fouling using non-destructive testing, namely seismic surface-wave and ground penetrating radar (GPR) survey. To understand this, a model rail track with different degree of fouling has been constructed in Civil engineering laboratory, University of Wollongong, Australia. Shear wave velocity obtained from seismic survey has been employed to identify the degree of fouling and types of fouling material. It is found that shear wave velocity of fouled ballast increases initially, reaches optimum fouling point (OFP), and decreases when the fouling increases. The degree of fouling corresponding after which the shear wave velocity of fouled ballast will be smaller than that of clean ballast is called the critical fouling point (CFP). Ground penetrating radar with four different ground coupled antennas (500 MHz, 800 MHz, 1.6 GHz and 2.3 GHz) was also used to identify the ballast fouling condition. It is found that the 800 MHz ground coupled antenna gives a better signal in assessing the ballast fouling condition. Seismic survey is relatively slow when compared to GPR survey however it gives quantifiable results. In contrast, GPR survey is faster and better in estimating the depth of fouling.     

Combining Electrical Resistivity Tomography and Ground Penetrating Radar to study geological structuring of karst Unsaturated Zone

This paper highlights the efficiency and complementarity of a light package of geophysical techniques to study the structure of karst Unsaturated Zone (UZ) in typical Mediterranean environment where soil cover is thin or absent. Both selected techniques, 2D Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), are widely used in environmental studies and their application is accessible for a lot of scientists/engineers. However, GPR or ERT alone is not able to provide an enhanced characterization of geological features in karst media. In the present study, GPR results supply a near surface high resolution imaging and thus can provide relevant geological information such as stratifications and fractures. Despite the quality of the results GPR's investigation depth remains limited to around 12 m. Apparent and inverted resistivity provided by ERT surveys shows strong lateral and vertical variations. These variations can inform about general geological structuring and feature orientation. ERT is able to prospect down to 40 m but it's a low resolution integrative technique. In the study area the investigated limestone is a commonly electrical resistive formation (more than 2000 Ω.m). However deeper than 5–7 m, the ERT profiles reveal several zones of moderate resistivity (around 900 Ω.m). In these zones a stratification change corresponding to slanted bedding is clearly identified by GPR results. The combination of both GPR and ERT results can allow a well-established geological interpretation. These moderate resistivity zones with slanted beddings can explain the presence of a perennial water flow point 35 m below the surface of the studied site within the underground gallery of the Low-Noise Underground Laboratory (LSBB).     

Scattering coefficients in the mantle revealed from the seismogram envelope analysis based on the multiple isotropic scattering model

At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10^− 4∼1.32 × 10^− 3 km^− 1 and 2.08 × 10^− 4∼6.23 × 10^− 4 km^− 1 at 4 s, and 4.51 × 10^− 4∼7.37 × 10^− 4 km^− 1 and 2.80 × 10^− 5∼2.71 × 10^− 4 km^− 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.     

Application of the Brewster angle to quantify the dielectric properties of ground ice formations

Previous studies of ground ice using moveout type ground-penetrating radar (GPR) surveys indicate that the dielectric permittivity can constrain the type of ground ice present in the subsurface. Due to the high-loss nature of the active layer over permafrost targets, however, the signal strength of GPR signals is often insufficient to resolve the basal boundary required for determining the dielectric permittivity of an underlying unit. We apply a non-conventional antenna orientation and post-processing method to determine the dielectric permittivity of the unit underlying the lowest resolvable boundary. We conduct moveout surveys using a 450 MHz GPR with collinear parallel oriented antennas on two adjacent ground ice formations in the region of Thomas Lee Inlet, Devon Island, Nunavut. We exploit the Brewster angle to calculate the approximate dielectric permittivity of ground ice formations below the active layer. The results agree within 1 dielectric unit with on-ice permittivity measurements made during a complementary study of the site.     

Single-crystal elasticity of stishovite: New experimental data obtained using high-frequency resonant ultrasound spectroscopy and a Gingham check structure model

Single-crystal elasticity of stishovite was examined using a new experimental technique and an empirical macroscopic model. Employing high-frequency resonant ultrasound spectroscopy, single-crystal elastic constants of stishovite were determined: C₁₁ = 443(3), C₃₃ = 781(4), C₁₂ = 193(2), C₂₃ = 199(2), C₄₄ = 256(2), and C₃₃ = 316(2) GPa. The frequency range of the resonant ultrasound spectroscopy techniques was 6–20 MHz, which is much lower than the ～10 GHz range of the Brillouin scattering technique. Of the elastic constants, the shear elastic constants C₄₄ and C₆₆ are consistent with the average value of the previously mentioned Brillouin scattering. Conversely, the four elastic constants, C₁₁, C₃₃, C₁₂, and C₂₃, slightly deviate outside the range of previous Brillouin scattering results. The present results, except those for C₁₂, are consistent with recent lattice dynamic analysis of inelastic X-ray scattering data. The adiabatic bulk modulus was calculated as 298 GPa, which is smaller and more consistent with the result of compression experiments than any other Brillouin scattering results (301–312 GPa). The present result shows greater P-wave velocity anisotropy (24.7%) than any preceding work. To understand the unique elastic properties of stishovite, the Gingham check model was proposed and examined. The result shows that the octahedron of 6-coordinated Si in stishovite crystal has stiffness comparable to that of diamond.     

Layer recognition and thickness evaluation of tunnel lining based on ground penetrating radar measurements

Ground penetrating radar (GPR) is a time-dependent geophysical technique for measuring the thickness of second lining, which is one of the most important tunnel quality targets. As the wave transmits from air (or the second lining) to the second lining (or the first lining), the Fresnel reflective coefficient is negative, i.e., the phase of reflective wave is reverse to the incident. In the one dimension time-waveform diagram (A-Scan) of GPR, the lining layers are located on the inflexions which are decided by Fresnel reflection coefficients and the attenuation coefficients of electromagnetic wave in the transmission medium. By towing the antenna over the tunnel surface, two dimension scanning data (B-Scan) is constituted by multiple A-Scan channels, where the grey scale is applied to the amplitude values. In the process of exploration, the lining interfaces are separately plotted by connecting each maximum peak point or second maximum peak point of A-Scan. In the Long Hai Tunnel, the artificial recognition provides 15 sampling values on the exploration line from D9 + 015.195 to D9 + 065.195. However, the automatic recognition can provide more information, such as average thicknesses, standard errors of lining thickness and the qualification rates of lining thickness, etc.     

On the Vp/Vs–Mg# correlation in mantle peridotites: Implications for the identification of thermal and compositional anomalies in the upper mantle

We use thermodynamically self-consistent and hybrid methods to analyze the correlation of important physical parameters (e.g. bulk density, elastic moduli) with bulk Mg# and modal composition in mantle peridotites at upper mantle conditions. Temperature (anharmonic and anelastic), pressure and compositional derivatives for all these parameters are evaluated. The results show that the widely used correlations between Vp/Vs and Mg# in peridotites are strictly valid only for garnet-bearing assemblages at temperatures < 900 °C. The correlation breaks down when: i) spinel is the stable Al-rich phase in the assemblage and ii) when anelastic attenuation of seismic velocities becomes important (T ? 900 °C). This implies that the range of applicability of published Vp/Vs–Mg# correlations for the upper mantle is limited to a depth interval between the spinel–garnet phase transition and the 900 °C isotherm. We use numerical simulations to show that this depth interval is virtually nonexistent in lithospheres thinner than ～ 140 km and can comprise up to ～ 50% of the lithospheric mantle in thick (> 220 km) lithospheric domains. In addition, we show that for most of the upper mantle the expected Δ(Vp/Vs) values associated with compositional variations are smaller than the resolution limit of current seismological methods. All these considerations suggest that the Vp/Vs ratio is not a reliable measure of compositional variations and that for large parts of the upper mantle compositional anomalies cannot be separated from thermal anomalies on the basis of seismological studies only. We further confirm that the only reliable indicator of compositional anomalies in a peridotitic mantle is the ratio of density to shear wave velocities (ρ/Vs). Our results demonstrate that geophysical–petrological models (forward or inverse) that model these two fields (i.e. density and Vs) self-consistently within a robust thermodynamic framework are necessary for characterizing the small-scale thermal and compositional structure of the lithosphere and sublithospheric upper mantle.     

A temperature-dependent multi-relaxation spectroscopic dielectric model for thawed and frozen organic soil at 0.05–15 GHz

A dielectric model for thawed and frozen Arctic organic-rich soil (50% organic matter) has been developed. The model is based on soil dielectric measurements that were collected over ranges of gravimetric moisture from 0.03 to 0.55 g/g, dry soil density from 0.72 to 0.87 g/cm³, and temperature from 25 to −30 °C (cooling run) in the frequency range of 0.05–15 GHz. The refractive mixing dielectric model was applied with the Debye multi-relaxation equations to fit the measurements of the soil’s complex dielectric constant as a function of soil moisture and wave frequency. The spectroscopic parameters of the dielectric relaxations for the bound, transient bound, and unbound soil water components were derived and were complimented by the thermodynamic parameters to obtain a complete set of parameters for the proposed temperature-dependent multi-relaxation spectroscopic dielectric model for moist soils. To calculate the complex dielectric constant of the soil, the following input variables must be assigned: (1) density of dry soil, (2) gravimetric moisture, (3) wave frequency, and (4) temperature. The error of the dielectric model was evaluated and yielded RMSE_ε′ values of 0.348 and 0.188 for the soil dielectric constant and the loss factor, respectively. These values are on the order of the dielectric measurement error itself. The proposed dielectric model can be applied in active and passive microwave remote sensing techniques to develop algorithms for retrieving the soil moisture and the freeze/thaw state of organic-rich topsoil in the Arctic regions.     

Shallow fault segmentation of the Alpine fault zone,New Zealand revealed from 2- and 3-D GPR surveying

Where they are preserved, landforms that have been truncated and offset by past fault movements provide potentially valuable quantitative data that can be used to estimate slip rates. At such locations, it is important to investigate the fault zone in sufficient detail to understand how displacements are accommodated on individual fault strands. At a site along a northern section of the Alpine fault zone on the South Island of New Zealand, surface mapping of a series of faulted river terraces and channels has revealed a complicated and poorly understood paleoearthquake history. We have acquired high-resolution 2- and 3-D ground-penetrating radar (GPR) data over a large area (~ 500 × 500 m) of the terraces to map along-strike changes in shallow (<20 m) fault zone morphology. By identifying distinct reflection patterns within the topographically migrated 3-D GPR volumes and extrapolating them to the longer and more widely spaced GPR profiles, we determined the subsurface extent of two main structural/depositional facies that were juxtaposed by three left-stepping en-echelon fault strands. Two regions of warped strata are interpreted to result from transpressive folding between the overlapping strands, where displacement is transferred from one fault to the next. We suggest that diffuse deformation between the overlapping fault tips results in anomalously low estimates for horizontal and vertical fault displacements of some geomorphic features.     

Detection of microwave emission due to rock fracture as a new tool for geophysics: A field test at a volcano in Miyake Island,Japan

This paper describes a field test to verify a newly discovered phenomenon of microwave emission due to rock fracture in a volcano. The field test was carried out on Miyake Island, 150 km south of Tokyo. The main objective of the test was to investigate the applicability of the phenomenon to the study of geophysics, volcanology, and seismology by extending observations of this phenomenological occurrence from the laboratory to the natural field.We installed measuring systems for 300 MHz, 2 GHz, and 18 GHz-bands on the mountain top and mountain foot in order to discriminate local events from regional and global events. The systems include deliberate data subsystems that store slowly sampled data in the long term, and fast sampled data when triggered. We successfully obtained data from January to February 2008. During this period, characteristic microwave pulses were intermittently detected at 300 MHz. Two photographs taken before and after this period revealed that a considerably large-scale collapse occurred on the crater cliff. Moreover, seismograms obtained by nearby observatories strongly suggest that the crater subsidence occurred simultaneously with microwave signals on the same day during the observation period.For confirmation of the microwave emission caused by rock fracture, these microwave signals must be clearly discriminated from noise, interferences, and other disturbances. We carefully discriminated the microwave data taken at the mountaintop and foot, checked the lightning strike data around the island, and consequently concluded that these microwave signals could not be attributed to lightning. Artificial interferences were discriminated by the nature of their waveforms. Thus, we inferred that the signals detected at 300 MHz were due to rock fractures during cliff collapses. This result may provide a useful new tool for geoscientists and for the mitigation of natural hazards.     

Magnetic susceptibility and Cs-137 inventory variability as influenced by land use change and slope positions in a hilly,semiarid region of west-central Iran

The objective of this study was to explore the slope position and land use change effects on the variability in magnetic susceptibility and ¹³⁷Cs inventory as the soil redistribution indicators in a hilly semiarid calcareous area in Iran. The selected study area is located in a hilly region with pasture and cultivation land use of Fereydunshahr, Isfahan Province in west-central Iran. In the two mentioned dominant ecosystems, four slope positions including summit, shoulder, backslope and footslope were identified and in each land use and slope position, three cores were selected to collect 72 soil samples from three depths (0–10, 10–20, 20–30 cm) in an area of 15 × 15 cm. Additional 28 soil samples were collected from the reference site for soil loss and deposition calculations by using the Cs-137 measurement. The results of the study with the use of the Cs-137 technique showed that the average soil loss in the pasture land (46.4 t ha^− 1 yr^− 1) was significantly (p < 0.05) lower than the average soil loss in the cultivated land (80.4 t ha^− 1 yr^− 1). The highest soil loss in both land uses was obtained in the shoulder position, 60.1 and 84.4 t ha^− 1 yr^− 1, respectively, for the pasture and cultivated lands. Moreover, the highest rates of soil deposition was observed in a footslope position in both land uses and they were 34 and 32.4 t ha^− 1 yr^− 1 for the pasture and cultivated lands, respectively. Magnetic susceptibility was significantly (p < 0.05) greater in pasture (χ_lf = 41.51 × 10^− 8 m³/kg) than in the cultivated land (χ_lf = 34.90 × 10^− 8 m³/kg). The pasture land with a lower soil loss rate, indicated significantly higher magnetic susceptibility in all landform positions as compared to that in the cultivated land. The results of the correlation analysis showed that among the studied soil physico-chemical properties, χ_lf (r = 0.83, p < 0.01) in the pasture land had the highest correlation with the Cs-137 inventory. Throughout the non-linear regression analysis, χ_lf was introduced for relating soil parameters and the cesium inventory explained 68% and 79% of the total variability of ¹³⁷Cs in the pasture and cultivated lands, respectively. The results implied that the variability in the magnetic susceptibility within the hillslope is consistent with the variation of the Cs-inventory; and the results thus demonstrate the slope and land use effects on soil redistribution.     

Comparing detection and location performance of perpendicular and parallel broadside GPR antenna orientations

GPR (Ground Penetrating Radar) results are shown for perpendicular broadside and parallel broadside antenna orientations. Performance in detection and localization of concrete tubes and steel tanks is compared as a function of acquisition configuration. The comparison is done using 100 MHz and 200 MHz center frequency antennas. All tubes and tanks are buried at the geophysical test site of IAG/USP in São Paulo city, Brazil. The results show that the long steel pipe with a 38-mm diameter was well detected with the perpendicular broadside configuration. The concrete tubes were better detected with the parallel broadside configuration, clearly showing hyperbolic diffraction events from all targets up to 2-m depth. Steel tanks were detected with the two configurations. However, the parallel broadside configuration was generated to a much lesser extent an apparent hyperbolic reflection corresponding to constructive interference of diffraction hyperbolas of adjacent targets placed at the same depth. Vertical concrete tubes and steel tanks were better contained with parallel broadside antennas, where the apexes of the diffraction hyperbolas better corresponded to the horizontal location of the buried target disposition. The two configurations provide details about buried targets emphasizing how GPR multi-component configurations have the potential to improve the subsurface image quality as well as to discriminate different buried targets. It is judged that they hold some applicability in geotechnical and geoscientific studies.     

Polar mesosphere summer echoes (PMSE) studied at Bragg wavelengths of 2.8 m, 67 cm,and 16 cm

We present observations of radar volume reflectivities under conditions of polar mesosphere summer echoes (PMSE) at three frequencies, i.e., 53.5, 224, and 930 MHz corresponding to Bragg wavelengths of 2.8, 0.67, and 0.16 m. These measurements were made with the ALWIN radar in Andenes and the EISCAT VHF and UHF radars in Tromsø. Contributions to the signal at 930 MHz by incoherent scatter are used to estimate electron number densities and their gradient at PMSE altitudes, and spectral width measurements of Doppler spectra recorded at 224 MHz are used to estimate the turbulent energy dissipation rate. We further derive a theoretical expression for the radar volume reflectivity for the case of turbulent scatter aided by a large Schmidt number (i.e., the current standard theory of PMSE) and show that our observations quantitatively agree with this theory if Schmidt numbers between 2500 and 5000 are assumed. We then show that these Schmidt numbers correspond to ice particles with radii in the range 20–30 nm which should frequently occur in the polar summer mesopause region. In addition, we show that for the short period when PMSE was observed at UHF frequencies the volume reflectivity is proportional to a factor determined by the turbulent energy dissipation rate, electron number density, and the electron number density gradient in agreement with theory. We consider our findings as strong support that PMSE at all considered frequencies is indeed created by turbulent scatter in the presence of a large Schmidt number. We finally highlight that ultimate proof of this concept will require the direct measurement of ice particle sizes in a PMSE environment probed by radars covering frequencies between 50 MHz and 1 GHz.     

Comparison of OSL ages from young dune sediments with a high-resolution independent age model

In order to test the accuracy of quartz optically stimulated luminescence (OSL) dating for young dune sediments (<100 a), a series of aerial images of a migrating sand dune is used to cross validate OSL ages. The investigated dune is situated on the northern part of the island of Sylt (southern North Sea). Based on aerial images and a map from 1925 to 2009 and the internal architecture of the dune obtained by ground-penetrating radar (GPR), an independent age model has been developed to attribute sedimentary-architectural elements of the dune to time. The annual rate of dune migration is calculated to be around 4.1 m/a. Along a 245 m transect oriented parallel to the direction of dune movement, 13 samples for OSL dating were collected at equidistant locations. Sand-sized quartz (150–250 μm) was used for determining the equivalent dose (D_e) applying a single-aliquot regenerative-dose (SAR) protocol. Results show that the oldest OSL age from the investigated recent dune appeared to be 110 ± 10 a, whereas the modern analogue was dated to 34 ± 3 a. In comparison with the aerial images, the OSL ages show a systematic overestimation of 10–40 a for six out of seven younger samples, which are expected to be younger than ∼60 a. This offset is negligible for older samples, but a substantial error in these younger ages. The overestimation is originated from a combination of small thermal transfer of 4–12 mGy during preheat and incomplete bleaching in medium OSL component causing a residual dose of about 15 mGy. The contribution of the incompletely bleached medium component cannot be removed totally by an early background subtraction approach. Despite the observed offset for youngest samples, this study corroborates the suitability of the OSL technique to date young dune sediments (<100 a).     

Estimation of in-situ density and moisture content in HMA pavements based on GPR trace reflection amplitude using different frequencies

The basic goal of the present research is to investigate the estimation of both the in-situ density and moisture content within the Hot Mix Asphalt (HMA) pavement layer(s) in a non-destructive way using Ground Penetrating Radar (GPR) trace reflection amplitude. For this purpose, an extensive pavement survey was conducted using an air-coupled GPR system, operating at 1 GHz or alternatively with a 2 GHz central frequency. The collected data were analyzed comparatively for the two antennae. The variability of electric permittivity caused by variations in HMA material is discussed, while the effect of the different frequencies is compared on the ability to retrieve permittivity, in-situ density and moisture content of the compacted HMA material using relationships suggested in reviewed international literature. The main finding of the present research is that for the same type of HMA material, the assessment of the material properties appears to be independent from the two central frequencies of investigation. However, there is evidence concerning the variations between the GPR wave data for the two different frequencies. The research highlights that the increased penetration depth of the 1 GHz antenna can provide an increased identification of areas of potential moisture within the body of HMA layer, and suggests that the variations between the permittivity values for the two different frequencies could be used to assess the homogeneity of material density with depth as an indicator of the mixture compaction. Additional findings are included within the paper.     

Optical properties and light penetration in a deep,naturally acidic,iron rich lake: Lago Caviahue (Patagonia,Argentina)

The optical properties and light climate in the deep and extremely acid Lake Caviahue have been studied in order to better understand its characteristics and the possible influence upon the phytoplankton community. The absorption coefficients for the dissolved fraction were maximal in the ultraviolet (UV) region and the water absorption spectra showed a shoulder around 300 nm, which was attributed to the concentration of Fe(III). No radiation was detected in the water column below 360 nm. The depth of the 1% incident radiation was dependent of wavelength, showing its maximum of 13.3 m at 565 nm, compared to 1.7 m and 4.8 m at 400 nm and 700 nm, respectively. Phytoplankton biomass was low and showed an almost constant profile with depth despite the relative darkness of the water column. Optical climate of Lake Caviahue is not typical of high elevation lakes but is more similar to low elevation shallow lakes of the Andean region. The chemical composition of the water, mainly Fe oxidation state and concentration, is the responsible for the high attenuation of the UV radiation (UVR). Living organisms are protected of UVR because Lake Caviahue waters are a shield against UV-B.     

Experimental and numerical evaluation of the effectiveness of a stiff wave barrier in the soil

This paper discusses the design, the installation, and the experimental and numerical evaluation of the effectiveness of a stiff wave barrier in the soil as a mitigation measure for railway induced vibrations. A full scale in situ experiment has been conducted at a site in El Realengo (Spain), where a barrier consisting of overlapping jet grout columns has been installed along a railway track. This barrier is stiff compared to the soil and has a depth of 7.5 m, a width of 1 m, and a length of 55 m. Geophysical tests have been performed prior to the installation of the barrier for the determination of the dynamic soil characteristics. Extensive measurements have been carried out before and after installation of the barrier, including free field vibrations during train passages, transfer functions between the track and the free field, and the track receptance. Measurements have also been performed at a reference section adjacent to the test section in order to verify the effect of changing train, track, and soil conditions over time. The in situ measurements show that the barrier is very effective: during train passages, a reduction of vibration levels by 5 dB is already obtained from 8 Hz upwards, while a peak reduction of about 12 dB is observed near 30 Hz immediately behind the barrier. The performance decreases further away from the jet grouting wall, but remains significant. The experimental results are also compared to numerical simulations based on a coupled finite element–boundary element methodology. A reasonable agreement between experiments and predictions is found, largely confirming the initially predicted reduction. This in situ test hence serves as a ‘proof of concept׳, demonstrating that stiff wave barriers are capable of significantly reducing vibration levels, provided that they are properly designed.     

Clay content evaluation in soils through GPR signal processing

The mechanical behavior of soils is partly affected by their clay content, which arises some important issues in many fields of employment, such as civil and environmental engineering, geology, and agriculture. This work focuses on pavement engineering, although the method applies to other fields of interest. Clay content in bearing courses of road pavement frequently causes damages and defects (e.g., cracks, deformations, and ruts). Therefore, the road safety and operability decreases, directly affecting the increase of expected accidents. In this study, different ground-penetrating radar (GPR) methods and techniques were used to non-destructively investigate the clay content in sub-asphalt compacted soils. Experimental layout provided the use of typical road materials, employed for road bearing courses construction. Three types of soils classified by the American Association of State Highway and Transportation Officials (AASHTO) as A1, A2, and A3 were used and adequately compacted in electrically and hydraulically isolated test boxes. Percentages of bentonite clay were gradually added, ranging from 2% to 25% by weight. Analyses were carried out for each clay content using two different GPR instruments. A pulse radar with ground-coupled antennae at 500 MHz centre frequency and a vector network analyzer spanning the 1–3 GHz frequency range were used. Signals were processed in both time and frequency domains, and the consistency of results was validated by the Rayleigh scattering method, the full-waveform inversion, and the signal picking techniques. Promising results were obtained for the detection of clay content affecting the bearing capacity of sub-asphalt layers.