A new method to detect transitory signatures and local time/space variability structures in the climate system: the scale-dependent correlation analysis

The study of transitory signals and local variability structures in both/either time and space and their role as sources of climatic memory, is an important but often neglected topic in climate research despite its obvious importance and extensive coverage in the literature. Transitory signals arise either from non-linearities, in the climate system, transitory atmosphere–ocean couplings, and other processes in the climate system evolving after a critical threshold is crossed. These temporary interactions that, though intense, may not last long, can be responsible for a large amount of unexplained variability but are normally considered of limited relevance and often, discarded. With most of the current techniques at hand these typology of signatures are difficult to isolate because the low signal-to-noise ratio in midlatitudes, the limited recurrence of the transitory signals during a customary interval of data considered. Also, there is often a serious problem arising from the smoothing of local or transitory processes if statistical techniques are applied, that consider all the length of data available, rather than taking into account the size of the specific variability structure under investigation. Scale-dependent correlation (SDC) analysis is a new statistical method capable of highlighting the presence of transitory processes, these former being understood as temporary significant lag-dependent autocovariance in a single series, or covariance structures between two series. This approach, therefore, complements other approaches such as those resulting from the families of wavelet analysis, singular-spectrum analysis and recurrence plots. A main feature of SDC is its high-performance for short time series, its ability to characterize phase-relationships and thresholds in the bivariate domain. Ultimately, SDC helps tracking short-lagged relationships among processes that locally or temporarily couple and uncouple. The use of SDC is illustrated in the present paper by means of some synthetic time-series examples of increasing complexity, and it is compared with wavelet analysis in order to provide a well-known reference of its capabilities. A comparison between SDC and companion techniques is also addressed and results are exemplified for the specific case of some relevant El Niño-Southern Oscillation teleconnections.     

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.     

A model that helps explain Sr-isotope disequilibrium between feldspar phenocrysts and melt in large-volume silicic magma systems

Feldspar phenocrysts of silicic volcanic rocks are commonly in Sr-isotopic disequilibrium with groundmass. In some cases the feldspar is more radiogenic, and in others it is less radiogenic. Several explanations have been published previously, but none of these is able to accommodate both senses of disequilibrium. We present a model by which either more- or less-radiogenic feldspar (or even both within a single eruptive unit) can originate. The model requires a magma body open to interaction with biotite- and feldspar-bearing wall rock. Magma is incrementally contaminated as wall rock melts incongruently. Biotite preferentially melts first, followed by feldspar. Such melting behavior, which is supported by both field and experimental studies, first contaminates magma with a relatively radiogenic addition, followed by a less-radiogenic addition. Feldspar phenocrysts lag behind melt (groundmass of volcanic rock) in incorporating the influx of contaminant, thus resulting in Sr-isotopic disequilibrium between the crystals and melt. The sense of disequilibrium recorded in a volcanic rock depends on when eruption quenches the contamination process. This model is testable by isotopic fingerprinting of individual feldspar crystals. For a given set of geologic boundary conditions, specific core-to-rim Sr-isotopic profiles are expectable. Moreover, phenocrysts that nucleate at different times during the contamination process should record different and predictable parts of the history. Initial results of Sr-isotopic fingerprinting of sanidine phenocrysts from the Taylor Creek Rhyolite are consistent with the model. More tests of the model are desirable.     

From subduction to transform motion: a seabeam survey of the Hellenic trench system

Preliminary results of a multi-narrow beam survey of the Hellenic trench system, in the Eastern Mediterranean, are presented. The southwestern Ionian branch is divided in small basins, partly filled with Pleistocene sediments. The morphology suggests that the basins are deformed by a compressional stress acting roughly perpendicularly to the trench along N50°E. This direction is the direction of the regional slip vector of the shallow thrust-type earthquakes. The structure of the southeastern Pliny-Strabo branch is quite different. Narrow en-e´chelon slots, oriented N40°E, have been mapped within the main troughs oriented N60°E. The regional earthquake slip vector is also oriented along N40°E. We conclude that the Hellenic trench system is an active subduction system, dominated by thrust along the Ionian branch and by transform motion along the Pliny-Strabo branch.     

The fit of the continents around the North Atlantic Ocean

A new reassembly of the continents around the North Atlantic Ocean is presented. The first criterion used for this reassembly is the identification of the structural framework related to the opening which consists of marginal fracture zones generated by offsets of the Rift. The Africa—North America, Eurasia—Greenland, Greenland—North America and Eurasia—North America adjustments are successively discussed. It is argued that the adjustments are best made at the 3000-meters isobath between Africa and North America and at the 2000-meters isobath for the younger rifts. The difference is attributed to subsidence and modification of continental margins with time. The importance of the Late Paleozoic tectonic phase in determining the subsequent pattern of Mesozoic rifting is emphasized.     

INTEGRAL/XMM views on the MeV source GRO J1411-64

The COMPTEL unidentified source GRO J 1411-64 was observed by INTEGRAL and XMM-Newton in 2005. The Circinus Galaxy is the only source detected within the 4σ location error of GRO J1411-64, but in here excluded as the possible counterpart. At soft X-rays, 22 reliable and statistically significant sources (likelihood >10) were extracted and analyzed from XMM-Newton data. Only one of these sources, XMMU J141255.6-635932, is spectrally compatible with GRO J1411-64 although the fact the soft X-ray observations do not cover the full extent of the COMPTEL source position uncertainty make an association hard to quantify and thus risky. At the best location of the source, detections at hard X-rays show only upper limits, which, together with MeV results obtained by COMPTEL suggest the existence of a peak in power output located somewhere between 300–700 keV for the so-called low state. Such a spectrum resembles those in blazars or microquasars, and might suggest at work by the models accordingly. However, an analysis using a microquasar model consisting on a magnetized conical jet filled with relativistic electrons, shows that it is hard to comply with all observational constrains. This fact and the non-detection at hard X-rays introduce an a-posteriori question mark upon the physical reality of this source, what is discussed here.     

Integrative taxonomic description of Plakina kanaky,a new polychromatic sponge species from New Caledonia (Porifera: Homoscleromorpha)

Four similar sponges of different colors, all unknown to science, were collected in submarine caves of New Caledonia. We aimed at determining whether the four chromotypes represented different species or phenotypic variations of a unique new species. We used an integrative taxonomic approach combining morphologic, molecular and metabolomic analyses. The main traits that define these specimens are a skeleton made of monolophose, trilophose and tetralophose calthrops only, high chemical diversity and a high abundance and diversity of prokaryotic symbionts. The symbiotic community includes two unique prokaryote morphotypes, which are described for the first time in Homoscleromorpha, and appeared to be vertically transmitted. Although several features slightly differ among chromotypes, the most parsimonious conclusion was to propose a single new species Plakina kanaky sp. nov. Our phylogenetic analysis indicated the paraphyly of the Plakina genus, with P. kanaky sp. nov. belonging to a clade that includes Plakina jani and Plakina trilopha. The present work demonstrates that integrative taxonomy should be used in order to revise the entire Plakinidae family and especially the non‐monophyletic genus Plakina.     

Tagging-the-triangle algorithm for partitioning features with inconsistent boundaries

In spatial data sets, gaps or overlaps among features are frequently found in spatial tessellations due to the non-abutting edges with adjacent features. These non-abutting edges in loose tessellations are also called inconsistent boundaries or slivers; polygons containing at least one inconsistent boundary are called inconsistent polygons or sliver polygons. The existing algorithms to solve topological inconsistencies in sliver polygons suffer from one or more of three major issues, namely determination of tolerances, excessive CPU processing time for large data sets and loss of vertex history. In this article, we introduce a new algorithm that mitigates these three issues. Our algorithm efficiently searches the features with inconsistent polygons in a given spatial data set and logically partitions them among adjacent features. The proposed algorithm employs the constrained Delaunay triangulation technique to generate labelled triangles from which inconsistent polygons with gaps and overlaps are identified using label counts. These inconsistent polygons are then partitioned using the straight skeleton method. Moreover, each of these partitioned gaps or overlaps is distributed among the adjacent features to improve the topological consistency of the spatial data sets. We experimentally verified our algorithm using the real land cadastre data set. The comparison results show that the proposed algorithm is four times faster than the existing algorithm for data sets with 200,000 edges.     

Geophysical exploration on the subsurface geology of La Garrotxa monogenetic volcanic field (NE Iberian Peninsula)

We applied self-potential (SP) and electrical resistivity tomography (ERT) to the exploration of the uppermost part of the substrate geology and shallow structure of La Garrotxa monogenetic volcanic field, part of the European Neogene–Quaternary volcanic province. The aim of the study was to improve knowledge of the shallowest part of the feeding system of these monogenetic volcanoes and of its relationship with the subsurface geology. This study complements previous geophysical studies carried out at a less detailed scale and aimed at identifying deeper structures, and together will constitute the basis to establish volcanic susceptibility in La Garrotxa. SP study complemented previous smaller-scale studies and targeted key areas where ERT could be conducted. The main new results include the generation of resistivity models identifying dykes and faults associated with several monogenetic cones. The combined results confirm that shallow tectonics controlling the distribution of the foci of eruptive activity in this volcanic zone mainly correspond to NNW–SSE and accessorily by NNE–SSW Neogene extensional fissures and faults and concretely show the associated magmatic intrusions. These structures coincide with the deeper ones identified in previous studies, and show that previous Alpine tectonic structures played no apparent role in controlling the loci of this volcanism. Moreover, the results obtained show that the changes in eruption dynamics occurring at different vents located at relatively short distances in this volcanic area are controlled by shallow stratigraphical, structural and hydrogeological differences underneath these monogenetic volcanoes.     

Can a Training Image Be a Substitute for a Random Field Model?

In most multiple-point simulation algorithms, all statistical features are provided by one or several training images (TI) that serve as a substitute for a random field model. However, because in practice the TI is always of finite size, the stochastic nature of multiple-point simulation is questionable. This issue is addressed by considering the case of a sequential simulation algorithm applied to a binary TI that is a genuine realization of an underlying random field. At each step, the algorithm uses templates containing the current target point as well as all previously simulated points. The simulation is validated by checking that all statistical features of the random field (supported by the simulation domain) are retrieved as an average over a large number of outcomes. The results are as follows. It is demonstrated that multiple-point simulation performs well whenever the TI is a complete (infinitely large) realization of a stationary, ergodic random field. As soon as the TI is restricted to a limited domain, the statistical features cannot be obtained exactly, but integral range techniques make it possible to predict how much the TI should be extended to approximate them up to a prespecified precision. Moreover, one can take advantage of extending the TI to reduce the number of disruptions in the execution of the algorithm, which arise when no conditioning template can be found in the TI.