Coda-Q and its lapse time dependence analysis in the interaction zone of the Dinarides,the Alps and the Pannonian basin

The single backscattering model was used to estimate total attenuation of coda waves (Q_c) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q₀ and n, parameters of the frequency dependent coda-Q using the relation Q_c = Q₀fⁿ. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Q_c were distributed into classes according to their lapse time, t_L. For t_L = 20–50 s we estimated Q₀ = 45–184 and n = 0.49–0.94, and for t_L = 60–100 s we obtained Q₀ = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q₀, and vice versa. The rates of change of both Q₀ and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.     

UV camera measurements of fumarole field degassing (La Fossa crater,Vulcano Island)

The UV camera is becoming an important new tool in the armory of volcano geochemists to derive high time resolution SO₂ flux measurements. Furthermore, the high camera spatial resolution is particularly useful for exploring multiple-source SO₂ gas emissions, for instance the composite fumarolic systems topping most quiescent volcanoes. Here, we report on the first SO₂ flux measurements from individual fumaroles of the fumarolic field of La Fossa crater (Vulcano Island, Aeolian Island), which we performed using a UV camera in two field campaigns: in November 12, 2009 and February 4, 2010. We derived ~ 0.5 Hz SO₂ flux time-series finding fluxes from individual fumaroles, ranging from 2 to 8.7 t d^?1, with a total emission from the entire system of ~ 20 t d^?1 and ~ 13 t d^?1, in November 2009 and February 2010 respectively. These data were augmented with molar H₂S/SO₂, CO₂/SO₂ and H₂O/SO₂ ratios, measured using a portable MultiGAS analyzer, for the individual fumaroles. Using the SO₂ flux data in tandem with the molar ratios, we calculated the flux of volcanic species from individual fumaroles, and the crater as a whole: CO₂ (684 t d^?1 and 293 t d^?1), H₂S (8 t d^?1 and 7.5 t d^?1) and H₂O (580 t d^?1 and 225 t d^?1).     

Importance of reversible attachment in predicting E. coli transport in saturated aquifers from column experiments

Drinking water wells indiscriminatingly placed adjacent to fecal contaminated surface water represents a significant but difficult to quantify health risk. Here we seek to understand mechanisms that limit the contamination extent by scaling up bacterial transport results from the laboratory to the field in a well constrained setting. Three pulses of Escherichia coli originating during the early monsoon from a freshly excavated pond receiving latrine effluent in Bangladesh were monitored in 6 wells and modeled with a two-dimensional (2-D) flow and transport model conditioned with measured hydraulic heads. The modeling was performed assuming three different modes of interaction of E. coli with aquifer sands: (1) irreversible attachment only (best-fit k_i = 7.6 day⁻¹); (2) reversible attachment only (k_a = 10.5 and k_d = 0.2 day⁻¹); and (3) a combination of reversible and irreversible modes of attachment (k_a = 60, k_d = 7.6, k_i = 5.2 day⁻¹). Only the third approach adequately reproduced the observed temporal and spatial distribution of E. coli, including a 4-log₁₀ lateral removal distance of ∼9 m. In saturated column experiments, carried out using aquifer sand from the field site, a combination of reversible and irreversible attachment was also required to reproduce the observed breakthrough curves and E. coli retention profiles within the laboratory columns. Applying the laboratory-measured kinetic parameters to the 2-D calibrated flow model of the field site underestimates the observed 4-log₁₀ lateral removal distance by less than a factor of two. This is promising for predicting field scale transport from laboratory experiments.     

Evapotranspiration in the Pampean Region using field measurements and satellite data

Evapotranspiration (LE) is an important factor for monitoring crops, water requirements, and water consumption at local and regional scale. In this paper, we applied the semi-empirical model to estimate the daily latent heat flux (LE_d = Rn_d + A − B(Ts − Ta)). LE_d has been estimated using satellite images (Thematic Mapper sensor) and a local dataset (incoming and outgoing short- and long-wave radiation) measured during three years. We first estimated the daily net Radiation (Rn_d) from a linear equation derived from the instantaneous net Radiation (Rn_d = CRn_i + D). Subsequently, coefficients A and B have been estimated for two different cover vegetations (pasture and soybean). For each vegetation cover, an error analysis combining Rn_d, A, B, and surface and air temperatures has been calculated. Results showed that Rn_d had good performance (nonbias and low RMSE). LE_d errors for pasture and soybean were ±28 W m⁻² and ±40 W m⁻² respectively.     

Saturated area dynamics and streamflow generation from coupled surface–subsurface simulations and field observations

A distributed physically-based model describing coupled surface–subsurface flows is applied to an instrumented catchment to investigate the links between runoff generation processes and the dynamics of saturated areas. The spatial characterization of the system is obtained through geophysical measurements and in situ observations. The model is able to reproduce the dynamics of the system through the calibration of only few parameters with a clear physical interpretation, providing a solid basis for our numerical investigations. Such investigations demonstrate the important control exerted by surface topography on the time evolution of saturated area patterns, mainly mediated by topographic curvature, that dictates both the dominant streamflow generation process at the local scale and the connection-disconnection dynamics of saturated areas. The relation between hillslope water storage and streamflow, Q = f(V), is shown to be highly hysteretical and dependent on the mean saturation of the catchment: higher degrees of saturation tend to yield one-to-one relationships between streamflow and water storage. On the contrary, streamflow-water storage relations are importantly affected by the specific configuration of saturated areas connected to the outlet when the system is far from complete saturation. This observation contradicts common assumptions of a one-to-one relationship Q = f(V) often used to justify widely observed power-law Q vs. dQ/dt recession curves. Furthermore, even when Q = f(V) becomes unique at high degrees of saturation, no power-law form emerged in our runs, speculatively because of the small size of the catchment formed by a single incision and the corresponding hillslope.     

The use of MRS in the determination of hydraulic transmissivity: The case of alluvial aquifers

Magnetic Resonance Sounding (MRS) is nowadays accepted as a new geophysical method that can be used for a reliable determination of the ground water content distribution in the top 150 m. A great effort has also been made in MRS development to deduce the hydraulic transmissivity, based on empiric relationships of the permeability with a factor F which is calculated with NMR parameters measured at laboratory scale. To use this relationship under field conditions a calibration coefficient C_T = T_pt / F has to be previously established, which demands the knowledge of the transmissivity T_pt evaluated in the pumping test. The transmissivity can then be calculated at any other site of the same aquifer using the relation T_mrs = C_TF. The C_T values reported suggest a certain relationship with the lithology, but with a great dispersion and contradictory results. MRS surveys carried out in alluvial aquifers in Spain have shown that the value of C_T evaluated at one site may not be valid at another place of the same aquifer, because of the great heterogeneity of this kind of geological environment. The demand of a pumping test at each site where a MRS is measured invalidates the method actually used for MRS transmissivity evaluation. More than 50 MRS have been used to propose a new methodology. The aquifers visited cover a great range of transmissivities (from 2 × 10^− 6 to 9 × 10^− 3 m²/s). The MRS signal amplitude varies between 20 and 1400 nV, the signal/noise ratio is in the range from 0.6 to 42, and the value of the decay time constant varies from 200 to 800 ms. It has been demonstrated that when the transmissivity increases, the value of F decreases, and C_T increases, except for certain groups of MRS taken at the same aquifer or part of one aquifer, for which F increases with T_pt, keeping C_T constant. A function C_T(F) of the type C_T = mF^− n has been obtained that allows the transmissivity evaluation without the need of T_pt. Considering that both values of transmissivity, T_pt and T_mrs, are subjected to deviations due to the experimental errors as well as due to evaluation errors, the prediction achieved by the proposed equation is rather good. To perform a better evaluation of the values of the coefficients m and n it is necessary to have a greater number of MR soundings of good quality and with a trustworthy inversion at locations where a really comparable and good performed pumping test is available, covering a sufficient range of transmissivities. Though the data we have used do not always fulfil these conditions, the result is promising. Once a trustable function is available, the forecast of the transmissivity using MRS will not need the existence of any pumping test in the area. The general extension of this methodology demands the availability of MRS taken at all kinds of geological and hydrogeological environments, which is impossible without the existence of a universal MRS data base.     

On the effective hydraulic conductivity in mean vertical unsaturated steady flows

Water flow in partially saturated heterogeneous porous formations is modelled by regarding the hydraulic parameters as stationary random space functions (RSFs). As a consequence, the flow variables are also RSFs, and we aim to develop a procedure to derive the effective hydraulic conductivity (EHC). The methodology relies on a perturbation approach which regards the variances of the hydraulic parameters as small quantities. By using the Gardner’s [Gardner WR. Some steady state solutions of unsaturated moisture flow equations with application to evaporation from a water table. Soil Sci 1958;85:228–32] two-parameters (K_s, α) model for the local unsaturated conductivity, we obtain the EHC for any dimensionality d of the flow domain, and arbitrary correlation functions of the input RSFs. Unlike previous studies [e.g. Yeh T-CJ, Gelhar J, Gutjahr A. Stochastic analysis of unsaturated flow in heterogenous soils. 1. Statistically isotropic media. Water Resour Res 1985;21;447–56, Yeh T-CJ, Gelhar J, Gutjahr A. Stochastic analysis of unsaturated flow in heterogenous soils. 2. Statistically anisotropic media with variable α. Water Resour Res 1985:21:457–64], the EHC is represented here as product between the local scale conductivity valid for a domain of mean parameters, and a correction function κ¹ which depends on the medium heterogeneity structure and the mean pressure head.Generally, the correction function κ¹ is expressed by d-fold quadrature. These quadratures are further reduced after adopting specific (i.e. exponential and Gaussian) structure for the (cross) correlation functions involved in the computation of κ¹. We have also focused on some particular formation structures which are relevant for the applications, and permit simplification of the computational aspect, as well.We investigate effects of the heterogeneity formation properties as well as the mean head on the structure of κ¹. Overall, results suggest that, given the formation statistics, the impact of the heterogeneity upon κ¹ is enhanced as the medium becomes drier. This is particularly so when the variability of the fluctuation of Y = ln K_s is small compared with that of ζ = ln α. Conversely, when the heterogeneity of Y is prevalent upon that of ζ, κ¹ is influenced solely by the anisotropic structure of the formation unless the horizontal correlation scales are much greater than the vertical ones.     

Seasonal variability of primary production in a fjord ecosystem of the Chilean Patagonia: Implications for the transfer of carbon within pelagic food webs

We characterized the seasonal cycle of productivity in Reloncaví Fjord (41°30′S), Chilean Patagonia. Seasonal surveys that included measurements of gross primary production, community respiration, bacterioplankton secondary production, and sedimentation rates along the fjord were combined with continuous records of water-column temperature variability and wind forcing, as well as satellite-derived data on regional patterns of wind stress, sea surface temperatures, and surface chlorophyll concentrations. The hydrography and perhaps fjord productivity respond to the timing and intensity of wind forcing over a larger region. Seasonal changes in the direction and intensity of winds, along with a late-winter improvement in light conditions, may determine the timing of phytoplankton blooms and potentially modulate productivity cycles in the region.Depth-integrated gross primary production estimates were higher (0.4–3.8 g C m^?2 d^?1) in the productive season (October, February, and May), and lower (0.1–0.2 g C m^?2 d^?1) in the non-productive season (August). These seasonal changes were also reflected in community respiration and bacterioplankton production rates, which ranged, respectively, from 0.3 to 4.8 g C m^?2 d^?1 and 0.05 to 0.4 g C m^?2 d^?1 during the productive and non-productive seasons and from 0.05 to 0.6 g C m^?2 d^?1 and 0.05 to 0.2 g C m^?2 d^?1 during the same two periods. We found a strong, significant correlation between gross primary production and community respiration (Spearman, r=0.95; p<0.001; n=12), which suggests a high degree of coupling between the synthesis of organic matter and its usage by the planktonic community. Similarly, strong correlations were found between bacterioplankton secondary production and both gross primary production (Spearman, r=0.7, p<0.05, n=9) and community respiration (Spearman, r=0.8, p<0.05, n=9), indicating that bacterioplankton may be processing an important fraction (8–59%) of the organic matter produced by phytoplankton in Reloncaví Fjord. In winter, bacterial carbon utilization as a percentage of gross primary production was >100%, suggesting the use of allochthonous carbon sources by bacterioplankton when the levels of gross primary production are low. Low primary production rates were associated with a greater contribution of small cells to autotrophic biomass, highlighting the importance of small-sized plankton and bacteria for carbon cycling and fluxes during the less productive winter months. Fecal pellet sedimentation was minimal during this period, also suggesting that most of the locally produced organic carbon is recycled within the microbial loop. During the productive season, on the other hand, the area exhibited a great potential to export organic matter, be it to higher trophic levels or vertically towards the bottom.     

Spectral attenuation of solar radiation in Patagonian fjord and coastal waters and implications for algal photobiology

The spectral attenuation of solar irradiation was measured during summer in two types of coastal waters in southern Chile, a north Patagonian fjord (Seno Reloncaví) and open coast (Valdivia). In order to relate the light availability with the light requirements of upper subtidal seaweeds, the saturating irradiance for photosynthesis (E_k) from P–I curves was measured. In addition the UV risk was assessed. Based on the z_1% of PAR, the lower limit of the euphotic zone in the studied systems averaged 21 m (K_d 0.24 m^?1) in Seno Reloncaví and 18 m (K_d 0.27 m^?1) in the coast of Valdivia. Photosynthesis of the studied seaweeds was saturated at markedly lower irradiances than found in their natural depths at the time of the study. Solar radiation penetrating into these depths at both locations largely supports the light requirements for the photosynthesis of subtidal species: 50–160 μmol m^?2 s^?1 for seaweeds from Seno Reloncaví (7 m tidal range) and 20–115 μmol m^?2 s^?1 for Valdivia assemblages (2 m tidal range). Optimal light conditions to saturate photosynthesis (E_k) were present at 10–16 m water depth. The attenuation of solar irradiation did not vary significantly between the fjord and coastal sites of this study. However, the underwater light climates to which seaweeds are exposed in these sites vary significantly because of the stronger influence of tidal range affecting the fjord system as compared with the open coastal site. The patterns of UV-B penetration in these water bodies suggest that seaweeds living in upper littoral zones such as the intertidal and shallow subtidal (<3 m) may be at risk.     

Adsorption behaviour of dibutyl phthalate on marine sediments

Laboratory experiments were carried out to investigate the adsorption behaviour of dibutyl phthalate (DBP) on marine sediments collected from five different sites in Victoria Harbour, Hong Kong. DBP adsorption can be well described by the Langmuir isotherm. The maximum DBP adsorption capacity (Q_max) of the marine sediments ranges from 53 to 79 mg g⁻¹, which has a positive correlation with their organic content. Around 90% of the organic can be removed from the sediments with treatment by H₂O₂ oxidation, and the Q_max then decreases to a range between 13 and 22 mg g⁻¹. The black carbon content of the sediments has a much greater DBP adsorption capacity than does the natural organic matter of the sediments. The amount of DBP adsorbed on the sediments increases as the salinity of the marine water increases.     

Estimation of the Free Core Nutation parameters from SG data: Sensitivity study and comparative analysis using linearized least-squares and Bayesian methods

The Free Core Nutation (FCN) is investigated with the help of its resonance effect on the tidal amplitudes in Superconducting Gravimeter (SG) records of the GGP network. The FCN resonance parameters are combined in a resonance equation involving the Earth's interior parameters. The sensitivity of the FCN parameters to the diurnal tidal waves demonstrates that the quality factor of the FCN is strongly dependent on the accuracy of the imaginary part estimates of the gravimetric factors close to the resonance. The weak amplitude of Ψ₁ tidal wave on the Earth, which is the closest in frequency to the FCN, in addition to errors in ocean loading correction, explains the poor determination of the quality factor Q from surface gravimetric data. The inversion of tidal gravimetric factors leads to estimates of the period, Q and resonance strength of the FCN. We show that, by inverting log(Q) instead of Q, the results using the least-squares method optimized using the Levenberg–Marquardt algorithm are in agreement with the Bayesian probabilistic results and agree with the results obtained from VLBI nutation data. Finally, a combined inversion of 7 GGP European SG data is performed giving T = 428 ± 3 days and 7762 < Q < 31,989 (90% C.I.). An experimental estimate of the internal pressure Love number is also proposed.     

Seasonal plankton variability in Chilean Patagonia fjords: Carbon flow through the pelagic food web of Aysen Fjord and plankton dynamics in the Moraleda Channel basin

Two research cruises (CIMAR 13 Fiordos) were conducted in the N–S oriented macrobasin of the Moraleda Channel (42–47°S), which includes the E–W oriented Puyuhuapi Channel and Aysen Fjord, during two contrasting productive seasons: austral winter (27 July–7 August 2007) and spring (2–12 November 2007). These campaigns set out to assess the spatio-temporal variability, defined by the local topography along Moraleda Channel, in the biological, physical, and chemical oceanographic characteristics of different microbasins and to quantify the carbon budget of the pelagic trophic webs of Aysen Fjord.Seasonal carbon fluxes and fjord-system functioning vary widely in our study area. In terms of spatial topography, two constriction sills (Meninea and Elefantes) define three microbasins along Moraleda Channel, herein the (1) north (Guafo-Meninea), (2) central (Meninea-Elefantes), and (3) south (Elefantes-San Rafael Lagoon) microbasins. In winter, nutrient concentrations were high (i.e. nitrate range: 21–14 μM) and primary production was low (153–310 mgC m^?2 d^?1), suggesting that reduced light radiation depressed the plankton dynamics throughout Moraleda Channel. In spring, primary production followed a conspicuous N–S gradient, which was the highest (5167 mgC m^?2 d^?1) in the north microbasin and the lowest (742 mgC m^?2 d^?1) in the south microbasin. The seasonal pattern of the semi-enclosed Puyuhuapi Channel and Aysen Fjord, however, revealed no significant differences in primary production (～800 mgC m^?2 d^?1), and vertical fluxes of particulate organic carbon were nearly twice as high in spring as in winter (266 vs. 168 mgC m^?2 d^?1).At the time-series station (St. 79), the lithogenic fraction dominated the total sedimented matter (seston). The role of euphausiids in the biological carbon pump of the Patagonian fjords was evident, given the predominance of zooplankton fecal material, mostly euphausiid fecal strings (46% of all fecal material), among the recognizable particles contributing to the particulate organic carbon flux.The topographic constriction sills partially modulated the exchange of oceanic waters (Subantarctic Surface Water) with freshwater river discharges along the Moraleda Channel. This exchange affects salinity and nutrient availability and, thus, the plankton structure. The north microbasin was dominated by a seasonal alternation of the classical (spring) and microbial (winter) food webs. However, in the south microbasin, productivity was low and the system was dominated year-round by large inputs of glacier-derived, silt-rich freshwater carrying predominantly small-sized diatoms (Skeletonema spp) and bacteria. When superimposed upon this scenario, highly variable (seasonal) solar radiation and photoperiods could exacerbate north–south differences along Moraleda Channel.     

Balancing the global oceanic neodymium budget: Evaluating the role of groundwater

The distinctly different, ε_Nd(0) values of the Atlantic, Indian, and Pacific Oceans requires that the residence time of Nd in the ocean (i.e., τ_Nd) be on the order of, or less than, the ocean mixing time of ∼ 500–1500 yr. However, estimates of τ_Nd, based on river influxes, range from 4000 to 15,000 yr, thus exceeding the ocean mixing time. In order to reconcile the oceanic Nd budget and lower the residence time by roughly a factor of 10, an additional, as yet unidentified, and hence “missing Nd flux” to the ocean is necessary. Dissolution of materials deposited on continental margins has previously been proposed as a source of the missing flux. In this contribution, submarine groundwater discharge (SGD) is examined as a possible source of the missing Nd flux. Neodymium concentrations (n = 730) and ε_Nd(0) values (n = 58) for groundwaters were obtained from the literature in order to establish representative groundwater values. Mean groundwater Nd concentrations and ε_Nd(0) values were used along with recent estimates of the terrestrial (freshwater) component of SGD (6% of river discharge on a global basis) to test whether groundwater discharge to the coastal oceans could account for the missing flux. Employing mean Nd concentrations of the compiled data base (i.e., 31.8 nmol/kg for all 730 analyses and 11.3 nmol/kg for 141 groundwater samples from a coastal aquifer), the global, terrestrial-derived SGD flux of Nd is estimated to range between 2.9 × 10⁷ and 8.1 × 10⁷ mol/yr. These estimates are of the same order of magnitude, and within a factor of 2, of the missing Nd flux (i.e., 5.4 × 10⁷ mol/yr). Applying the SGD Nd flux estimates, the global average ε_Nd(0) of SGD is predicted to be − 9.1, which is similar to our estimate for the missing Nd flux (− 9.2), and in agreement with the mean (± S.D.) ε_Nd(0) measured in groundwaters (i.e., ε_Nd(0) = −8.9 ± 4.2). The similarities in the estimated SGD Nd flux and corresponding ε_Nd(0) values to the magnitude and isotope composition of the missing Nd flux are compelling, and suggest that discharge of groundwater to the oceans could account for the missing Nd flux. Future investigations should focus on quantifying the Nd concentrations and isotope compositions of groundwater from coastal aquifers from a variety of coastal settings, as well as the important geochemical reactions that effect Nd concentrations in subterranean estuaries in order to better constrain contributions of SGD to the oceanic Nd budget.     

The use of Bayesian methods for fitting rating curves,with case studies

Using data from two very large watersheds and five smaller, this paper explores the use of Bayesian methods for fitting rating curves. Posterior distribution of rating-curve parameters were calculated using Markov Chain Monte Carlo (MCMC) methods, and 95% credible intervals were calculated for predicted discharges, given stage. Expected discharge was related to stage using a link function. For the five smaller watersheds, the assumptions were (a) that the distribution of discharge Q, given stage h, is Normal, with variance proportional to h; (b) that a log link function relates μ_Q, the mean of Q∣h, to a function of stage, of the form μ_Q = β(h + α)^λ. For the two large watersheds, however, a better fit was obtained by taking the distribution of Q to be log-Normal, and the link function as ln μ_Q = β₀ + β₁h. For the two large watersheds, priors for all three parameters were taken as uninformative; for the five smaller, the prior for parameter λ was taken as Normally distributed, N(2, 0.5). Acceptable ratings were obtained for all seven sites. It is argued that distributions of derived variables (such as annual maximum discharge) can be derived directly from (a) the posterior distribution of rating-curve parameters, and (b) the stage record, without recourse to additional assumptions. Estimates thus obtained for the T-year event will incorporate rating-curve uncertainty. It is argued that Bayesian methods are appropriate for rating-curve calculation because their inherent flexibility (a) allows the incorporation of prior information about the nature of a rating curve; (b) yields credible intervals for predicted discharges and quantities derived from them; (c) can be extended to allow for uncertainty in stage measurements.     

Inorganic carbon and nutrient fluxes on the Arctic Shelf

Historic data from the Russian-American Hydrochemical Atlas of Arctic Ocean together with data from the TRANSDRIFT II 1994 and TUNDRA 1994 cruises have been used to assess the spatial and inter-annual variability of carbon and nutrient fluxes, as well as air–sea CO₂ exchange in the Laptev and western East Siberian Seas during the summer season. Budget computations using summer data of dissolved inorganic phosphate (DIP), dissolved inorganic nitrogen (DIN) and dissolved inorganic carbon (DIC) gives that the Laptev Sea shelf is a net sink of DIP and DIN of 2.5×10⁶, 23.2×10⁶ mol d⁻¹, respectively, while it is a net source of DIC (excluding air–sea exchange) of 1249×10⁶ mol d⁻¹. In the East Siberian Seas the budget computations give 0.5×10⁶, −11.4×10⁶ and −173×10⁶ mol d⁻¹ (minus being a sink) for DIP, DIN, and DIC, respectively. In summers, the Laptev Sea Shelf is net autotrophic while the East-Siberian Sea Shelf is net heterotrophic, and both systems are weak net denitrifying. The Laptev Sea Shelf takes up 2.1 mmol CO₂ m⁻² d⁻¹ from atmosphere, whereas the western part of the East-Siberian Sea Shelf loose 0.3 mmol CO₂ m⁻² d⁻¹ to the atmosphere. The variability of DIP, DIN and DIC fluxes during summer in the different regions of the Laptev and East Siberian Seas depends on bottom topography, river runoff, exchange with surrounding seas and wind field.     

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.     

Persistent organic pollutants and stable isotopes in pinnipeds from King George Island,Antarctica

In the present work, fat, skin, liver and muscle samples from Leptonychotes weddellii (Weddell seal, n = 2 individuals), Lobodon carcinophagus (crabeater seal, n = 2), Arctocephalus gazella (Antarctic fur seal, n = 3) and Mirounga leonina (southern elephant seal, n = 1) were collected from King George Island, Antarctica, and analysed for POPs (PCBs, organochlorine pesticides and PBDEs) and stable isotopes (δ¹³C and δ¹⁵N in all tissues but fat). PBDEs could be found in only one sample (L. weddellii fat). Generally, PCBs (from 74 to 523 ng g⁻¹ lw), DDTs (from 14 to 168 ng g⁻¹ lw) and chlordanes (from 9 to 78 ng g⁻¹ lw) were the prevailing compounds. Results showed a clear stratification in accordance with ecological data. Nonetheless, stable isotope analyses provide a deeper insight into fluctuations due to migrations and nutritional stress. Correlation between δ¹⁵N and pollutants suggests, to some degree, a considerable ability to metabolize and/or excrete the majority of them.     

Young eclogite from the Greater Himalayan Sequence,Arun Valley,eastern Nepal: P–T–t path and tectonic implications

Garnet geochronology was used to provide the first direct measurement of the timing of eclogitization in the central Himalaya. Lu–Hf dates from garnet separates in one relict eclogite from the Arun River Valley in eastern Nepal indicate an age of 20.7 ± 0.4 Ma, significantly younger than ultra-high pressure eclogites from the western Himalaya, reflecting either different origins or substantial time lags in tectonics along strike. Four proximal garnet amphibolites from structurally lower horizons are 14–15 Ma, similar to post-eclogitization ages published for rocks along strike in southern Tibet. P–T calculations indicate three metamorphic episodes for the eclogite: i) eclogite-facies metamorphism at ～ 670 °C and ≥ 15 kbar at 23–16 Ma; ii) a peak-T granulite event at ～ 780 °C and 12 kbar; and iii) late-stage amphibolite-facies metamorphism at ～ 675 °C and 6 kbar at ～ 14 Ma. The garnet amphibolites were metamorphosed at ～ 660 °C. Three models are considered to explain the observed P–T–t evolution. The first assumes that the Main Himalayan Thrust (basal thrust of the Himalayan thrust system) cuts deeper at Arun than elsewhere. While conceptually the simplest, this model has difficulty explaining both the granulite-facies overprint and the pulse of exhumation between 25 and 14 Ma. A second model assumes that (aborted) subduction, slab breakoff, and ascent of India's leading edge occurred diachronously: ～ 50 Ma in the western Himalaya, ～ 25 Ma in the central Himalaya of Nepal, and presumably later in the eastern Himalaya. This model explains the P–T–t path, particularly heating during initial exhumation, but implies significant along-strike diachroneity, which is generally lacking in other features of the Himalaya. A third model assumes repeated loss of mantle lithosphere, first by slab breakoff at ～ 50 Ma, and again by delamination at ～ 25 Ma; this model explains the P–T–t path, but requires geographically restricted tectonic behavior at Arun. The P–T–t history of the Arun eclogites may imply a change in the physical state of the Himalayan metamorphic wedge at 16–25 Ma, ultimately giving rise to the Main Central Thrust by 15–16 Ma.     

Luminescence dating of archaeosediments: A comparison of K-feldspar and plagioclase IRSL ages

IRSL single-aliquot dating of three feldspar fractions from a Howieson's Poort (SA) industry site has been investigated. In the case of the plagioclase fraction (ρ > 2.62 g cm^?3) equivalent dose is similar to that of the other fractions (i.e. ρ < 2.58 g cm^?3, high-K-feldspar; and 2.58 < ρ < 2.62 g cm^?3; low-K-feldspar). However, the g values measured for this fraction is such that the fading-corrected age is overestimated. The K content of each fraction was estimated by SEM. The DRC-corrected ages for the two lighter feldspar fractions are consistent with the archaeological interpretation. This study shows that reliable single-aliquot luminescence ages can be obtained using K-feldspar and Na-plagioclase.     

Monitoring of Guadalentín valley (southern Spain) through a fast SAR Interferometry method

Ground subsidence is a critical problem in populated areas, which requires attention and routine monitoring. Specifically, subsidence related to the exploitation of aquifer systems is of particular interest in terms of risk for the possible damage to buildings and availability of water resources. We propose to realise a survey of the aquifer systems near the city of Lorca, southern Spain, using Synthetic Aperture Radar Interferometry (InSAR) with a light method necessitating a small number of interferograms rather than analysing hundreds of radar images. We computed nine interferograms with good coherence over the years 2004 and 2005. We analysed them by a two-step filtering process in order to characterise the atmospheric contribution and the ground deformation. We are then able to identify two deforming areas and quantify subsidence rates at 0.20 ± 0.02 mm/day south of Lorca and southeast to the city of Alhama de Murcia with the same efficiency as the studies using sophisticated InSAR methods. Using continuous records of water levels in boreholes, we deduced a mean storage coefficient (S_k) at 4.1 × 10^− 3 which is very close to the previous published values, and coincides in this case to an inelastic deformation mode of all the system. With this value, we determine the groundwater level variations from the measurements of the ground subsidence.