Inter-laboratory calibration of low-field magnetic and anhysteretic susceptibility measurements

Inter-laboratory and absolute calibrations of rock magnetic parameters are fundamental for grounding a rock magnetic database and for semi-quantitative estimates about the magnetic mineral assemblage of a natural sample. Even a dimensionless ratio, such as anhysteretic susceptibility normalized by magnetic susceptibility (K_a/K) may be biased by improper calibration of one or both of the two instruments used to measure K_a and K. In addition, the intensity of the anhysteretic remanent magnetization (ARM) of a given sample depends on the experimental process by which the remanence is imparted. We report an inter-laboratory calibration of these two key parameters, using two sets of artificial reference samples: a paramagnetic rare earth salt, Gd₂O₃ and a commercial “pozzolanico” cement containing oxidized magnetite with grain size of less than 0.1 μm according to hysteresis properties. Using Gd₂O₃ the 10 Kappabridges magnetic susceptibility meters (AGICO KLY-2 or KLY-3 models) tested prove to be cross-calibrated to within 1%. On the other hand, Kappabridges provide a low-field susceptibility value that is ca. 6% lower than the tabulated value for Gd₂O₃, while average high-field susceptibility values measured on a range of instruments are indistinguishable from the tabulated value. Therefore, we suggest that Kappabridge values should be multiplied by 1.06 to achieve absolute calibration. Bartington Instruments magnetic susceptibility meters with MS2B sensors produce values that are 2-13% lower than Kappabridge values, with a strong dependence on sample centering within the sensor. The K_a/K ratio of ca. 11, originally obtained on discrete cement samples with a 2G Enterprises superconducting rock magnetometer and a KLY-2, is consistent with reference parameters for magnetites of grain size <0.1 μm. On the other hand, K_a values from a 2G Enterprises magnetometer and K values from a Bartington Instruments MS2C loop sensor for u-channel and discrete cement samples, will produce average K_a/K values that are unrealistically high if not properly corrected for the nominal volume detected by the sensors for these instruments. Inter-laboratory measurements of K and K_a for standard paleomagnetic plastic cubes filled with cement indicate remarkable differences in the intensity of the newly produced ARMs (with a standard deviation of ca. 21%), that are significantly larger than the differences observed from the calibration of the different magnetometers employed in each laboratory. Differences in the alternating field decay rate are likely the major source of these variations, but cannot account for all the observed variability. With such large variations in experimental conditions, classical interpretation of a “King plot” of K_a versus K would imply significant differences in the determination of grain size of magnetite particles on the same material.     

Rescaling species optima estimated by weighted averaging

The common practice of linear deshrinking in weighted averaging is known to be equivalent to a linear rescaling of the estimated species optima. In published lists of species optima, the use of rescaling is recommended, as it allows values derived from different data sets to be compared and used for new inferences, assuming that taxonomic consistency is assured. Rescaling optima is also shown to influence WA estimates of species tolerances. Non-linear rescaling is also discussed, in the form of cubical rescaling and weighted averaging-partial least squares (WA-PLS). The use of a different deshrinking equation in a small data set did lead to similar prediction errors, probably because of the small size of the data set used.     

Space-time distribution of small earthquakes at Phlegraean Fields,south-central Italy

A preliminary statistical analysis of the space-time distribution of small seismic events in the volcanic area of Phlegraean Fields, south-central Italy, was done on the basis of a catalogue of earthquakes recorded by the local seismic stations in the period January 1, December 31, 1983.The non-random character of the sequence has been tested by matching the observed time-dstribution of seismic events with the theoretical Poisson process.A clustered occurrence of seismic events seems to be the main cause of the departure from a Poisson process as the inter-arrival times distribution clearly shows.Furthermore, the non-random behaviour of the events space-time distribution mainly due to quiescient and clustered recursive seismic phases could be studied by using the method proposed byVon Seggern et al. (1981). The analysis and the space-time diagrams confirm the swarm-type character of the entire seismic sequence.     

Records of environmental and climatic changes during the late Holocene from Svalbard: palaeolimnology of Kongressvatnet

A multi-core, multidisciplinary palaeolimnological study of the partially varved sediment of a deep, meromictic, arctic lake, Kongressvatnet (Svalbard, Western Spitsbergen), provides a record of environmental and climatic changes during last ca. 1800 years. The chronology of sedimentation was established using several dating techniques (¹³⁷Cs, ²¹⁰Pb, varve counts, palaeomagnetic correlation). A multiproxy record of palaeolimnological variability was compiled based on sedimentation rates, magnetic properties, varve thickness, organic matter, geochemistry, pigments from algal and photosynthetic bacteria, mineralogy and biological assemblages (diatoms, Cladocera). The major features recognised in our master core K99-3 include a shift in sediment source and supply (magnetic measurements, geochemistry) probably caused by glaciological changes in the catchment around 38–32 cm core depth (AD 700–820). Additional environmental changes are inferred at 20–18, 8–4.5 and 3–2 cm (AD ca. 1160–1255; 1715–1880; 1940–1963, respectively). During the past ca. 120 years a prominent sedimentological change from brownish-grey, partly laminated silt-clay (varves) to black organic-rich deposits was observed. From AD 1350 to AD1880 the sediment is comprised of a continuous sequence of varves, whereas the earlier sediments are mostly homogeneous with only a few short intercalated laminated sections between AD 860 and 1350. Sedimentation and accumulation rates increased during the last 30 years (modern warming). Pigment concentrations are very low in the lower ca. 32 cm of the core (AD 820) probably because of the high turbidity high energy environment. The high sulphur content in the uppermost 32 cm of sediment has given rise to two horizontally stratified populations of sulphur anaerobic photosynthetic bacteria, as inferred from their specific carotenoids. These bacteria populations are much more abundant during the Little Ice Age (LIA) than during warmer periods (e.g., during the Medieval Warm Period and 20th century). Diatoms are lacking from the core base up to 18 cm (ca. AD 1255); at this level, species indicative of mesotrophic water are present, whereas from 17 cm to the top of the core, oligotrophic taxa such as Staurosira construens/S. pinnata complex dominate, indicating extended ice coverage and more oligotrophic waters during the LIA. The concentration of Cladocera subfossil remains (dominated by Chydorus) are relatively high in the deepest sections (54–32 cm), whereas the upper 32 cm are characterized by a very low concentration of remains, possibly because of the strongly anoxic conditions, and in this upper sediment section rotifer resting eggs become prevalent. We interpret these changes as responses to climate forcing through its impact on glacial melt water, lake ice cover duration and mainly redox conditions in deep water. The observed changes suggest that at least some of our recorded changes may parallel the Greenland Ice core, although our study added more details about the inferred climatic changes. Further aspects are discussed, such as catchment processes, glacial activity, duration of the Medieval Warm Period, the Little Ice Age, local human activity, and limnology.     

The build up of the isotopic signal in skeletons of the stony coral Porites lutea

The build up of the isotopic signal in corals was followed by sampling the newly formed skeleton at a monthly resolution for a period of two years in order to establish the interrelations between the calcification processes and the skeletal isotopic composition. We deployed two underwater sampling schemes, which provide a monitor of the changes in water temperature and δ¹⁸O and in the corresponding newly accreted skeleton of undisturbed Porites lutea colonies under natural conditions and four transplanted colonies, which maintained the genetic identity throughout the experiment. The results indicate that δ¹⁸O of the newly accreted skeleton does not correlate with ambient temperature although the seasonal temperature variability at the site (winter to summer) is in the order of 6 °C and δ¹⁸O of seawater is constant throughout the year. In contrast to the newly formed surface skeleton, the isotopic compositions of the deep and older parts of the skeleton show the predicted annual isotopic pattern with highly significant correlation between δ¹⁸O_s and SST. The transformation between temperature-independent to temperature-dependent isotopic signal occurs several months after the skeleton was formed at the surface. The position of the skeleton in relation to the open sea may generate the difference between δ¹⁸O_s of the surface skeleton and that of the skeleton previously accreted further down the tissue layer. Our data support the general model of a multi-step skeletogenesis process, where the temperature independent skeleton is entails the first step, the production of skeletal scaffold, and the environmental temperature signature is captured by the next two other steps: the thickening and the periodic abrupt uplift occurring at the depth of the tissue layer. However, re-examination and development of the current isotopic models for coral calcification are required in order to explain the observed different temperature dependency during the growth’s sequence.     

Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite‐difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions.     

Depth and morphology of reflectors from the non-linear inversion of arrival-time and waveform semblance data. Part I: method and applications to synthetic data

We propose a two-dimensional, non-linear method for the inversion of reflected/converted traveltimes and waveform semblance designed to obtain the location and morphology of seismic reflectors in a lateral heterogeneous medium and in any source-to-receiver acquisition lay-out. This method uses a scheme of non-linear optimization for the determination of the interface parameters where the calculation of the traveltimes is carried out using a finite-difference solver of the Eikonal equation, assuming an a priori known background velocity model. For the search for the optimal interface model, we used a multiscale approach and the genetic algorithm global optimization technique. During the initial stages of inversion, we used the arrival times of the reflection phase to retrieve the interface model that is defined by a small number of parameters. In the successive steps, the inversion is based on the optimization of the semblance value determined along the calculated traveltime curves. Errors in the final model parameters and the criteria for the choice of the best-fit model are also estimated from the shape of the semblance function in the model parameter space. The method is tested and validated on a synthetic dataset that simulates the acquisition of reflection data in a complex volcanic structure. This study shows that the proposed inversion approach is a valid tool for geophysical investigations in complex geological environments, in order to obtain the morphology and positions of embedded discontinuities.     

Turbidity limits gas exchange in a large macrotidal estuary

In estuaries, the gas transfer velocity (k) is driven by a combination of two major physical drivers, wind and water current. The k values for CO₂ in the macrotidal Gironde Estuary were obtained from 159 simultaneous pCO₂ and floating chamber flux measurements. Values of k increased with wind speed and were significantly greater when water currents and wind were in opposing directions. At low wind speeds (<1 m s⁻¹), k increased with water current velocities (0–1.5 m s⁻¹) following an exponential trend. The latter was a good proxy for the Y-intercept in a generic equation for k versus wind speed in estuaries. We also found that, in this turbid estuary, k was significantly lower at high turbidity. The presence of suspended material in great concentrations (TSS > 0.2 g L⁻¹) had a significant role in attenuating turbulence and therefore gas exchange. This result has important consequences for modeling water oxygenation in estuarine turbidity maxima. For seven low turbidity estuaries previously described in the literature, the slope of the linear regression between k and wind speed correlates very well with the estuary surface area due to a fetch effect. In the Gironde Estuary, this slope follows the same trend at low turbidity (TSS < 0.2 g L⁻¹), but is on average significantly lower than in other large estuaries and decreases linearly with the TSS concentration. A new generic equation for estuaries is proposed that gives k as a function of water current velocity, wind speed, estuarine surface area and TSS concentration.     

A Critical Analysis of Transverse Dispersivity Field Data

Transverse dispersion, or tracer spreading orthogonal to the mean flow direction, which is relevant e.g, for quantifying bio-degradation of contaminant plumes or mixing of reactive solutes, has been studied in the literature less than the longitudinal one. Inferring transverse dispersion coefficients from field experiments is a difficult and error-prone task, requiring a spatial resolution of solute plumes which is not easily achievable in applications. In absence of field data, it is a questionable common practice to set transverse dispersivities as a fraction of the longitudinal one, with the ratio 1/10 being the most prevalent. We collected estimates of field-scale transverse dispersivities from existing publications and explored possible scale relationships as guidance criteria for applications. Our investigation showed that a large number of estimates available in the literature are of low reliability and should be discarded from further analysis. The remaining reliable estimates are formation-specific, span three orders of magnitude and do not show any clear scale-dependence on the plume traveled distance. The ratios with the longitudinal dispersivity are also site specific and vary widely. The reliability of transverse dispersivities depends significantly on the type of field experiment and method of data analysis. In applications where transverse dispersion plays a significant role, inference of transverse dispersivities should be part of site characterization with the transverse dispersivity estimated as an independent parameter rather than related heuristically to longitudinal dispersivity.     

Decoupled kinematics of two neighbouring permafrost creeping landforms in the Eastern Italian Alps

An overall acceleration of rock glacier displacement rates in the Alps has been observed in recent decades, with several cases of destabilization leading to potential geomorphological hazards. This behaviour has been attributed to the rising permafrost temperature, induced by atmospheric warming and regulated by thermo-hydrological processes. Landforms derived from the interaction of glacier remnants and permafrost are widespread in mountain areas, but are less studied and monitored than talus rock glaciers. This work presents a comparative study of a talus rock glacier and a glacial-permafrost composite landform (GPCL) in the Eastern Italian Alps. The two landforms are only 10 km apart, but have rather different elevation ranges and main slope aspects. The kinematics and ground thermal conditions were monitored from 2001 to 2015 along with geomorphological surveys, analyses of historical maps and remote sensing data. The dynamic behaviour of the rock glacier was similar to the majority of monitored rock glaciers in the Alps, with an acceleration after 2008 and a velocity peak in 2015. In contrast, the GPCL had a nearly unchanged displacement rate during the observation period. Statistical analyses of kinematic vs. nivo-meteorological variables revealed a dynamic decoupling of the two landforms after 2008 that corresponds with increased winter snow accumulation. Although the kinematics of both landforms respond to ground surface temperature variations, the collected evidence suggests a different reaction of ground surface temperature to variations in the precipitation regime. This different reaction is likely due to local topo-climatic conditions that affect snow redistribution by wind. The different reactions of the two systems to the same climatic forcing is likely a legacy of their different origins. GPCL dynamics result from interaction of permafrost and residual glacial dynamics that are associated with possible peculiarities in the internal/basal meltwater circulation, whose future response is uncertain and requires improved understanding. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.