Application of K-feldspar-jadeite-quartz barometry to eclogite facies metagranites and metapelites in the Sesia Lanzo Zone (Western Alps, Italy)

The eclogite facies assemblage K-feldspar-jadeite-quartz in metagranites and metapelites from the Sesia-Lanzo Zone (Western Alps, Italy) records the equilibration pressure by dilution of the reaction jadeite + quartz = albite. The metapelites show partial transformation from a pre-Alpine assemblage of garnet (Alm(63)Prp(26)Grs(10))-K-feldspar-plagioclase-biotite +/- sillimanite to the Eo-Alpine high-pressure assemblage garnet (Alm(50)Prp(14)Grs(35))-jadeite (Jd(80-97)Di(0-4)Hd(0-8)Acm(0-7))=zoisite-phengite. Plagioclase is replaced by jadeite-zoisite-kyanite-K-feldspar-quartz and biotite is replaced by garnet-phengite or omphacite-kyanite-phengite. Equilibrium was attained only in local domains in the metapelites and therefore the K-feldspar-jadeite-quartz (KJQ) barometer was applied only to the plagioclase pseudomorphs and K-feldspar domains. The albite content of K-feldspar ranges from 4 to 11 mol% in less equilibrated assemblages from Val Savenca and from 4 to 7 mol% in the partially equilibrated samples from Monte Mucrone and the equilibrated samples from Montestrutto and Tavagnasco. Thermodynamic calculations on the stability of the assemblage K-feldspar-jadeite-quartz using available mixing data for K-feldspar and pyroxene indicate pressures of 15-21 kbar (+/- 1.6-1.9 kbar) at 550 +/- 50 degrees C. This barometer yields direct pressure estimates in high-pressure rocks where pressures are seldom otherwise fixed, although it is sensitive to analytical precision and the choice of thermodynamic mixing model for K-feldspar. Moreover, the KJQ barometer is independent of the ratio P-H2O/P-T. The inferred limiting a(H2O) for the assemblage jadeite-kyanite in the metapelites from Val Savenca is low and varies from 0.2 to 0.6.

A new maximum likelihood method for luminosity calibrations

A new statistical parallax method using the Maximum Likelihood principle is presented, allowing the simultaneous determination of a luminosity calibration, kinematic characteristics and spatial distribution of a given sample. This method has been developed for the exploitation of the Hipparcos data and presents several improvements with respect to the previous ones: the effects of the selection of the sample, the observational errors, the galactic rotation and the interstellar absorption are taken into account as an intrinsic part of the formulation (as opposed to external corrections). Furthermore, the method is able to identify and characterize physically distinct groups in inhomogeneous samples, thus avoiding biases due to unidentified components. Moreover, the implementation used by the authors is based on the extensive use of numerical methods, so avoiding the need for simplification of the equations and thus the bias they could introduce. Several examples of application using simulated samples are presented, to be followed by applications to real samples in forthcoming articles.

An Efficient Segmentation Method for Ultrasound Images based on a Semi-supervised Approach and Patch-based Features

Segmenting ultrasound images is a challenging problemwhere standard unsupervised segmentation methods such asthe well-known Chan-Vese method fail. We propose in thispaper an efficient segmentation method for this class ofimages. Our proposed algorithm is based on asemi-supervised approach (user labels) and the use ofimage patches as data features. We also consider thePearson distance between patches, which has been shown tobe robust w.r.t speckle noise present in ultrasoundimages. Our results on phantom and clinical data show avery high similarity agreement with the ground truthprovided by a medical expert.

Modelling of weather radar echoes from anomalous propagation using a hybrid parabolic equation method and NWP model data

Contamination of weather radar echoes by anomalous propagation (anaprop) mechanisms remains a serious issue in quality control of radar precipitation estimates. Although significant progress has been made identifying clutter due to anaprop there is no unique method that solves the question of data reliability without removing genuine data. The work described here relates to the development of a software application that uses a numerical weather prediction (NWP) model to obtain the temperature, humidity and pressure fields to calculate the three dimensional structure of the atmospheric refractive index structure, from which a physically based prediction of the incidence of clutter can be made. This technique can be used in conjunction with existing methods for clutter removal by modifying parameters of detectors or filters according to the physical evidence for anomalous propagation conditions. The parabolic equation method (PEM) is a well established technique for solving the equations for beam propagation in a non-uniformly stratified atmosphere, but although intrinsically very efficient, is not sufficiently fast to be practicable for near real-time modelling of clutter over the entire area observed by a typical weather radar. We demonstrate a fast hybrid PEM technique that is capable of providing acceptable results in conjunction with a high-resolution terrain elevation model, using a standard desktop personal computer. We discuss the performance of the method and approaches for the improvement of the model profiles in the lowest levels of the troposphere.

Analysis of warm season thunderstorms using an object-oriented tracking method based on radar and total lightning data

Monitoring thunderstorms activity is an essential part of operational weather surveillance given their potential hazards, including lightning, hail, heavy rainfall, strong winds or even tornadoes. This study has two main objectives: firstly, the description of a methodology, based on radar and total lightning data to characterise thunderstorms in real-time; secondly, the application of this methodology to 66 thunderstorms that affected Catalonia (NE Spain) in the summer of 2006. An object-oriented tracking procedure is employed, where different observation data types generate four different types of objects (radar 1-km CAPPI reflectivity composites, radar reflectivity volumetric data, cloud-to-ground lightning data and intra-cloud lightning data). In the framework proposed, these objects are the building blocks of a higher level object, the thunderstorm. The methodology is demonstrated with a dataset of thunderstorms whose main characteristics, along the complete life cycle of the convective structures (development, maturity and dissipation), are described statistically. The development and dissipation stages present similar durations in most cases examined. On the contrary, the duration of the maturity phase is much more variable and related to the thunderstorm intensity, defined here in terms of lightning flash rate. Most of the activity of IC and CG flashes is registered in the maturity stage. In the development stage little CG flashes are observed (2% to 5%), while for the dissipation phase is possible to observe a few more CG flashes (10% to 15%). Additionally, a selection of thunderstorms is used to examine general life cycle patterns, obtained from the analysis of normalized (with respect to thunderstorm total duration and maximum value of variables considered) thunderstorm parameters. Among other findings, the study indicates that the normalized duration of the three stages of thunderstorm life cycle is similar in most thunderstorms, with the longest duration corresponding to the maturity stage (approximately 80% of the total time).

A wavelet-based method for the determination of the relative resolution between remotely sensed images

Spatial resolution is a key parameter of all remote sensing satellites and platforms. The nominal spatial resolution of satellites is a well-known characteristic because it is directly related to the area in ground that represents a pixel in the detector. Nevertheless, in practice, the actual resolution of a specific image obtained from a satellite is difficult to know precisely because it depends on many other factors such as atmospheric conditions. However, if one has two or more images of the same region, it is possible to compare their relative resolutions. In this paper, a wavelet-decomposition-based method for the determination of the relative resolution between two remotely sensed images of the same area is proposed. The method can be applied to panchromatic, multispectral, and mixed (one panchromatic and one multispectral) images. As an example, the method was applied to compute the relative resolution between SPOT-3, Landsat-5, and Landsat-7 panchromatic and multispectral images taken under similar as well as under very different conditions. On the other hand, if the true absolute resolution of one of the images of the pair is known, the resolution of the other can be computed. Thus, in the last part of this paper, a spatial calibrator that is designed and constructed to help compute the absolute resolution of a single remotely sensed image is described, and an example of its use is presented.

New trapping method to survey for presence of the Etruscan shrew Suncus etruscus, the smallest mammal

The presence of the Etruscan shrew Suncus etruscus is hard to prove where its predator, the barn owl Tyto alba, is absent, because most live traps are not triggered by it. I therefore developed a new trapping method involving a feeding period of 1 week followed by one night of trapping using modified Trip Trap traps. I show here in detail how I caught four Etruscan shrews in 2010 with 24 traps in the Valley of Dora Baltea (Piemonte, Italy). In 2011, another 11 Etruscan shrews were caught in Piemonte and Lombardia, Italy, and Ticino, Switzerland. The proposed new method is useful for establishing the presence of the species.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

A new method to measure rolling resistance in treadmill cycling.

The purpose of this study was to evaluate a new method of measuring rolling resistance in treadmill cycling and to establish its sensitivity and reproducibility. One participant was asked to keep a bicycle in equilibrium on a treadmill without pedalling at a constant speed of 5.56 m x s(-1), which was held in place in the front by a dynamometer. For each condition, the method consisted of 11 measurements of the force required to hold the cycle at different treadmill slopes (0-10%, increment 1%). The coefficient of rolling resistance was calculated based on the forces applied to the bicycle in equilibrium. To test the sensitivity of the method, the bicycle was successively equipped with three tyre types (700 x 28, 700 x 23, 700 x 22) and inflation pressure was set at 150, 300, 600, 900, and 1100 kPa. To test the reproducibility of the method, a second experimenter repeated all measurements done with the 700 x 23 tyres. The method was sensitive enough to detect an effect of both tyre type and inflation pressure (P < 0.001: two-way ANOVA). The measurement of the coefficient of rolling resistance by two separate experimenters resulted in a small bias of 0.00029 (95% CI, -0.00011 to 0.00068). In conclusion, the new method is sensitive and reliable, as well as being simple and affordable.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

Maturation of marginal zone and follicular B cells requires B cell activating factor of the tumor necrosis factor family and is independent of B cell maturation antigen.

B cells undergo a complex series of maturation and selection steps in the bone marrow and spleen during differentiation into mature immune effector cells. The tumor necrosis factor (TNF) family member B cell activating factor of the TNF family (BAFF) (BLyS/TALL-1) plays an important role in B cell homeostasis. BAFF and its close homologue a proliferation-inducing ligand (APRIL) have both been shown to interact with at least two receptors, B cell maturation antigen (BCMA) and transmembrane activator and cyclophilin ligand interactor (TACI), however their relative contribution in transducing BAFF signals in vivo remains unclear. To functionally inactivate both BAFF and APRIL, mice transgenic for a soluble form of TACI were generated. They display a developmental block of B cell maturation in the periphery, leading to a severe depletion of marginal zone and follicular B2 B cells, but not of peritoneal B1 B cells. In contrast, mice transgenic for a soluble form of BCMA, which binds APRIL, have no detectable B cell phenotype. This demonstrates a crucial role for BAFF in B cell maturation and strongly suggests that it signals via a BCMA-independent pathway and in an APRIL-dispensable way.

Synorogenic extension: Quantitative constraints on the age and displacement of the Zanskar shear zone (northwest Himalaya)

subsequent extension-induced exhumation. Geochronological dating of various Structural, thermobarometric, and geochronological data place limits on the age and tectonic displacement along the Zanskar shear zone, a major north-dipping synorogenic extensional structure separating the high-grade metamorphic sequence of the High Himalayan Crystalline Sequence from the overlying low-grade sedimentary rocks of the Tethyan Himalaya, A complete Barrovian metamorphic succession, from kyanite to biotite zone mineral assemblages, occurs within the I-km-thick Zanskar shear zone. Thermobarometric data indicate a difference In equilibration depths of 12 +/- 3 km between the lower kyanite zone and the garnet zone, which is Interpreted as a minimum estimate for the finite vertical displacement accommodated by the Zanskar shear zone. For the present-day dip of the structure (20 degrees), a simple geometrical model shows that a net slip of 35 +/- 9 km is required to regroup these samples to the same structural level. Because the kyanite to garnet zone rocks represent only part of the Zanskar shear zone, and because its original dip may have been less than the present-day dip, these estimates fur the finite displacement represent minimum values. Field relations and petrographic data suggest that migmatization and associated leucogranite intrusion in the footwall of the Zanskar shear zone occurred as a continuous profess starting at the Barrovian metamorphic peak and lasting throughout the subsequent extension-induced exhumation. Geochronological dataing of various leucogranitic plutons and dikes in the Zanskar shear zone footwall indicates that the main ductile shearing along the structure ended by 19.8 Ma and that extension most likely initiated shortly before 22.2 Ma.

The central Sesia-Lanzo zone (Western Italian Alps) - New field observations and lithostratigraphic subdivisions

An extensive study of the central part of the Sesia Lanzo Zone has been undertaken to identify pre-Alpine protoliths and to reconstruct the lithologic and tectonic setting of this part of the Western Alps. Three main complexes have been defined: 1) the Polymetamorphic Basement Complex, corresponding to the lower unit of the Sesia Lanzo Zone after COMPAGNONI et al. (1977), is further subdivided into the three following units: a) an Internal Unit characterized by eo-Alpine high pressure (HP) assemblages (DAL PIAZ et al., 1972) (Eclogitic Micaschists); b) an Intermediate Unit where HP parageneses are partially re-equilibrated under greenschist conditions and c) an External Unit where the main foliation is defined by a greenschist paragenesis (Gneiss Minuti auct.). 2) the Monometamorphic Cover Complex, subdivided into the followings: a) the Bonze Unit, composed of sheared metagabbros, eclogitized metabasalts with MORB geochemical affinity and related metasediments (micaschists, quartzites and Mn-cherts) and b) the Scalaro Unit, containing predominantly metasediments of supposed Permo-Triassic age (yellow dolomitic marbles, calcschists and conglomeratic limestones, micaschists and quartzites with thin levels of basic rocks with within plate basalts [WPB] geochimical affinity). Multiple lithostratigraphic sequences for the Monometamorphic Cover Complex are proposed. The contact between the Bonze and Scalaro Units is defined by repetitions of dolomitic marbles and metabasalts; the ages of the metasediments have been assigned solely by analogy with other sediments of the Western Alps, due to the absence of fossils. The Monometamorphic Cover Complex can be considered as the autochthonous cover of the Sesia Lanzo Zone because of the primary contacts with the basement and because of the presence of preAlpine HT basement blocks in the cover sequences. 3) The pre-Alpine high temperature (HT) Basement Complex (or `'Seconda Zona Diorito-Kinzigitica''), comprises HT Hercynian rocks like kinzigites, amphibolites, granulites and calcite marbles; this Complex is always located between the Internal and the External Units and can be followed continuously for several kilometers south of the Gressoney Valley to the Orco Valley. A schematic evolution for the Sesia Lanzo Zone is proposed; based on available data together with new geochronological data, this study shows that the internal and external parts of the polymetamorphic basement of the Sesia Zone experienced different cooling histories .