Improvement of Nakamura technique by singular spectrum analysis

In this work we investigate the application of the Singular Spectrum Analysis to improve Nakamura Technique. The SSA has a wide and multidisciplinary range of applications; it allows a time series to be decomposed into different components, e.g. the signal itself, as well as various noise components, which can be subsequently removed from the data. Removal of the minor component of the data can lead to significant improvements in the identification of the system. In this paper the use of SSA Technique allows to optimize the signal to noise ratio before computing the classical Nakamura spectral ratios. A number of typical applications are also presented.     

Importance of geochemical transformations in determining submarine groundwater discharge-derived trace metal and nutrient fluxes

Seasonal (Spring and Summer 2002) concentrations of dissolved (<0.22 μm) trace metals (Ag, Al, Co, Cu, Mn, Ni, Pb), inorganic nutrients (NO₃, PO₄, Si), and DOC were determined in groundwater samples from 5 wells aligned along a 30 m shore-normal transect in West Neck Bay, Long Island, NY. Results show that significant, systematic changes in groundwater trace metal and nutrient composition occur along the flowpath from land to sea. While conservative mixing between West Neck Bay water and the groundwaters explains the behavior of Si and DOC, non-conservative inputs for Co and Ni were observed (concentration increases of 10- and 2-fold, respectively) and removal of PO₄ and NO₃ (decreases to about half) along the transport pathway. Groundwater-associated chemical fluxes from the aquifer to the embayment calculated for constituents not exhibiting conservative behavior can vary by orders of magnitude depending on sampling location and season (e.g. Co, 3.4 × 10²– 8.2 × 10³ μmol d⁻¹). Using measured values from different wells as being representative of the true groundwater endmember chemical composition also results in calculation of very different fluxes (e.g., Cu, 6.3 × 10³ μmol d⁻¹ (inland, freshwater well) vs. 2.1 × 10⁵ μmol d⁻¹(seaward well, S = 17 ppt)). This study suggests that seasonal variability and chemical changes occurring within the subterranean estuary must be taken into account when determining the groundwater flux of dissolved trace metals and nutrients to the coastal ocean.     

Petrogenesis of tourmaline rocks associated with Fe-carbonate–graphite metapelite, metabasite and strata-bound polymetallic sulphide mineralisation, Peloritani Mountains, Sicily, Southern Italy

Tourmalinite and tourmaline-rich rocks associated with Fe-carbonate–graphite phyllite, strata-bound polymetallic sulphide deposits, metabasite and marble were studied, for information on the mechanism of tourmaline formation in the pre-Hercynian low-grade metamorphic sequence of the Mandanici Unit in the Peloritani Mountains of Sicily, southern Italy. The major and trace element compositions of the tourmaline rocks suggest the existence of a sedimentary protolith with pre-metamorphic black shale and bedded chert. Boron was interpreted to be accumulated in a restricted sedimentary basin, between platform carbonate formations, with abundant organic matter and Fe–Al–Ti-rich laterite–bauxite soil-derived clastic supply, under a continental volcano-tectonic extensional regime accompanied by a local convective hydrothermal system along faults. Petrographic, crystal–chemical and δ¹¹B isotopic data are compatible with a model of marine sediment dewatering at temperatures below 200 °C, which caused the removal of boron from clay. Metamorphism led to the development of tourmaline in an Al–Ti-rich environment, in equilibrium with other minerals such as ilmenite, albite and muscovite. The upper temperature of metamorphism (almost 375 °C), estimated on the basis of δ¹¹B, fits geothermometric results from Δ¹³C_{carbonate–graphite} on associated rocks. The estimated value of δ¹¹B in the tourmalinite protolith, − 7.5‰ , is also compatible with continental-derived Al-rich sediments.     

Homogeneity of the Geochemical Reference Material BRP-1 (Paraná Basin Basalt) and Assessment of Minimum Mass

Reference materials (RM) are required for quantitative analyses and their successful use is associated with the degree of homogeneity, and the traceability and confidence limits of the values established by characterisation. During the production of a RM, the chemical characterisation can only commence after it has been demonstrated that the material has the required level of homogeneity. Here we describe the preparation of BRP-1, a proposed geochemical reference material, and the results of the tests to evaluate its degree of homogeneity between and within bottles. BRP-1 is the first of two geochemical RM being produced by Brazilian institutions in collaboration with the United States Geological Survey (USGS) and the International Association of Geoanalysts (IAG). Two test portions of twenty bottles of BRP-1 were analysed by wavelength dispersive-XRF spectrometry and major, minor and eighteen trace elements were determined. The results show that for most of the investigated elements, the units of BRP-1 were homogeneous at conditions approximately three times more rigorous than those strived for by the test of "sufficient homogeneity". Furthermore, the within bottle homogeneity of BRP-1 was evaluated using small beam (1 mm²) synchrotron radiation XRF spectrometry and, for comparison, the USGS reference materials BCR-2 and GSP-2 were also evaluated. From our data, it has been possible to assign representative minimum masses for some major constituents (1 mg) and for some trace elements (1-13 mg), except Zr in GSP-2, for which test portions of 74 mg are recommended.     

Sensitivity analysis of seismic hazard for Western Liguria (North Western Italy): A first attempt towards the understanding and quantification of hazard uncertainty

The use of logic trees in probabilistic seismic hazard analyses often involves a large number of branches that reflect the uncertainty in the selection of different models and in the selection of the parameter values of each model. The sensitivity analysis, as proposed by Rabinowitz and Steinberg [Rabinowitz, N., Steinberg, D.M., 1991. Seismic hazard sensitivity analysis: a multi-parameter approach. Bull. Seismol. Soc. Am. 81, 796–817], is an efficient tool that allows the construction of logic trees focusing attention on the parameters that have greater impact on the hazard.In this paper the sensitivity analysis is performed in order to identify the parameters that have the largest influence on the Western Liguria (North Western Italy) seismic hazard. The analysis is conducted for six strategic sites following the multi-parameter approach developed by Rabinowitz and Steinberg [Rabinowitz, N., Steinberg, D.M., 1991. Seismic hazard sensitivity analysis: a multi-parameter approach. Bull. Seismol. Soc. Am. 81, 796–817] and accounts for both mean hazard values and hazard values corresponding to different percentiles (e.g., 16%-ile and 84%-ile). The results are assessed in terms of the expected PGA with a 10% probability of exceedance in 50 years for rock conditions and account for both the contribution from specific source zones using the Cornell approach [Cornell, C.A., 1968. Engineering seismic risk analysis. Bull. Seismol. Soc. Am. 58, 1583–1606] and the spatially smoothed seismicity [Frankel, A., 1995. Mapping seismic hazard in the Central and Eastern United States. Seismol. Res. Lett. 66, 8–21]. The influence of different procedures for calculating seismic hazard, seismic catalogues (epicentral parameters), source zone models, frequency–magnitude parameters, maximum earthquake magnitude values and attenuation relationships is considered. As a result, the sensitivity analysis allows us to identify the parameters with higher influence on the hazard. Only these parameters should be subjected to careful discussion or further research in order to reduce the uncertainty in the hazard while those with little or no effect can be excluded from subsequent logic-tree-based seismic hazard analyses.     

Crustal attenuation in the Southern Andean retroarc (38°–39°30′ S) determined from tectonic and gravimetric studies: The Lonco-Luán asthenospheric anomaly

The western retroarc of the Southern Andes between 38° and 40° S is formed by a NNW-elongated ridge not associated with stacked thrust sheets. On the contrary, during the last 4–3 Ma this ridge was affected by extensional deformation, regional uplift and related folding on a very broad scale. Receiver function analysis shows that the drainage divide area and adjacent retroarc lie over an attenuated crust. Expected crustal thickness at these latitudes is around 38 km, whereas in this part of the retroarc the thickness is less than 32 km. The causes for such attenuation have been linked to a moderate steepening of the subducted Nazca plate beneath the South American plate, which is suggested by a westward shift and narrowing of the magmatic arc during the last 4 to 5 Ma. Gravimetric studies show that the upper plate did not react homogeneously to slab steepening, but ancient sutures and lithospheric discontinuities deeply buried under Mesozoic to Cenozoic sequences in the retroarc were locally reactivated. These processes resulted in an asthenospheric anomaly that correlates at the surface with the area of Pliocene to Quaternary doming, widespread extension and three radial troughs. Two of the troughs have accommodated substantial amounts of extension, but the third was probably aborted at an early stage. Moreover, the presence of an anomalous concentration of calderas and large volcanic centers over the proposed asthenospheric anomaly, and their age distribution, may indicate minor migration of the asthenospheric anomaly between 4 and 2 Ma through the western South American plate.     

Mesozoic extension in the Basque–Cantabrian basin (N Spain): Contributions from AMS and brittle mesostructures

In this work we analyse and check the results of anisotropy of magnetic susceptibility (AMS) by means of a comparison with palaeostress orientations obtained from the analysis of brittle mesostructures in the Cabuérniga Cretaceous basin, located in the western end of the Basque–Cantabrian basin, North Spain. The AMS data refer to 23 sites including Triassic red beds, Jurassic and Lower Cretaceous limestones, sandstones and shales. These deposits are weakly deformed, and represent the syn-rift sequence linked to basins formed during the Mesozoic and later inverted during the Pyrenean compression. The observed magnetic fabrics are typical of early stages of deformation, and show oblate, triaxial and prolate magnetic ellipsoids. The magnetic fabric seems to be related to a tectonic overprint of an original, compaction, sedimentary fabric. Most sites display a NE–SW magnetic lineation that is interpreted to represent the stretching direction of the Early Cretaceous extensional stage of the basin, without recording of the Tertiary compressional events, except for sites with compression-related cleavage.Brittle mesostructures include normal faults, calcite and quartz tension gashes and joints, related to the extensional stage. The results obtained from joints and tension gashes show a dominant N–S to NE–SW, and secondary NW–SE, extension direction. Paleostresses obtained from fault analysis (Right Dihedra and stress inversion methods) indicate NW–SE to E–W, and N–S extension direction. The results obtained from brittle mesostructures show a complex pattern resulting from the superposition of several tectonic processes during the Mesozoic, linked to the tectonic activity related to the opening of the Bay of Biscay during the Early Cretaceous. This work shows the potential in using AMS analysis in inverted basins to unravel its previous extensional history when the magnetic fabric is not expected to be modified by subsequent deformational events. Brittle mesostructure analysis seems to be more sensitive to far-field stress conditions and record longer time spans, whereas AMS records deformation on the near distance, during shorter intervals of time.     

On the Intensity Thresholds of the Network and Plage Regions

A period of minimum solar activity in April and May 1996 was used to analyze full-disk CaK spectroheliograms with the aim of deriving the intensity thresholds of the quiet network and the plages by applying the so-called inflexion point method. The average network intensity threshold is found to be influenced mainly by the seeing, whereas the average plage threshold shows an increasing trend from disk center towards the solar limb. Both parameters are compared with the results of other authors.