Generalized fuzzy clustering for segmentation of multi-spectral magnetic resonance images.

An integrated approach for multi-spectral segmentation of MR images is presented. This method is based on the fuzzy c-means (FCM) and includes bias field correction and contextual constraints over spatial intensity distribution and accounts for the non-spherical cluster's shape in the feature space. The bias field is modeled as a linear combination of smooth polynomial basis functions for fast computation in the clustering iterations. Regularization terms for the neighborhood continuity of intensity are added into the FCM cost functions. To reduce the computational complexity, the contextual regularizations are separated from the clustering iterations. Since the feature space is not isotropic, distance measure adopted in Gustafson-Kessel (G-K) algorithm is used instead of the Euclidean distance, to account for the non-spherical shape of the clusters in the feature space. These algorithms are quantitatively evaluated on MR brain images using the similarity measures.

A meta-analysis of cambium phenology and growth: linear and non-linear patterns in conifers of the northern hemisphere

Background and Aims Ongoing global warming has been implicated in shifting phenological patterns such as the timing and duration of the growing season across a wide variety of ecosystems. Linear models are routinely used to extrapolate these observed shifts in phenology into the future and to estimate changes in associated ecosystem properties such as net primary productivity. Yet, in nature, linear relationships may be special cases. Biological processes frequently follow more complex, non-linear patterns according to limiting factors that generate shifts and discontinuities, or contain thresholds beyond which responses change abruptly. This study investigates to what extent cambium phenology is associated with xylem growth and differentiation across conifer species of the northern hemisphere. Methods Xylem cell production is compared with the periods of cambial activity and cell differentiation assessed on a weekly time scale on histological sections of cambium and wood tissue collected from the stems of nine species in Canada and Europe over 1–9 years per site from 1998 to 2011. Key Results The dynamics of xylogenesis were surprisingly homogeneous among conifer species, although dispersions from the average were obviously observed. Within the range analysed, the relationships between the phenological timings were linear, with several slopes showing values close to or not statistically different from 1. The relationships between the phenological timings and cell production were distinctly non-linear, and involved an exponential pattern. Conclusions The trees adjust their phenological timings according to linear patterns. Thus, shifts of one phenological phase are associated with synchronous and comparable shifts of the successive phases. However, small increases in the duration of xylogenesis could correspond to a substantial increase in cell production. The findings suggest that the length of the growing season and the resulting amount of growth could respond differently to changes in environmental conditions.

A re-examination of petrogenesis and 40Ar/39Ar systematics in the the Chain of Ponds K-feldspar: "diffusion domain" archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, on which thermochronological modeling of Ar diffusion in discrete “domains” was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. 40Ar/39Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest “diffusion domains”. Heterochemical phases also play an important role in producing non-linear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick’s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the "diffusion domain" model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution amongst different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for 40Ar-39Ar hygrochronology as it is for U-Pb geochronology.

Analysis of life expectancy using linear and nonlinear models

Life expectancy has consistently increased over the last 150 years due to improvements in nutrition, medicine, and public health. Several studies found that in many developed countries, life expectancy continued to rise following a nearly linear trend, which was contrary to a common belief that the rate of improvement in life expectancy would decelerate and was fit with an S-shaped curve. Using samples of countries that exhibited a wide range of economic development levels, we explored the change in life expectancy over time by employing both nonlinear and linear models. We then observed if there were any significant differences in estimates between linear models, assuming an auto-correlated error structure. When data did not have a sigmoidal shape, nonlinear growth models sometimes failed to provide meaningful parameter estimates. The existence of an inflection point and asymptotes in the growth models made them inflexible with life expectancy data. In linear models, there was no significant difference in the life expectancy growth rate and future estimates between ordinary least squares (OLS) and generalized least squares (GLS). However, the generalized least squares model was more robust because the data involved time-series variables and residuals were positively correlated. ^

Geochemistry and oxygen and iron isotope ratios of mafic rocks

Recycling of oceanic crust into the deep mantle via subduction is a widely accepted mechanism for creating compositional heterogeneity in the upper mantle and for explaining the distinct geochemistry of mantle plumes. The oxygen isotope ratios (d18O) of some ocean island basalts (OIB) span values both above and below that of unmetasomatised upper mantle (5.5 ± 0.4 per mil) and provide support for this hypothesis, as it is widely assumed that most variations in d18O are produced by near-surface low-temperature processes. Here we show a significant linear relationship between d18O and stable iron isotope ratios (d57Fe) in a suite of pristine eclogite xenoliths. The d18O values of both bulk samples and garnets range from values within error of normal mantle to significantly lighter values. The observed range and correlation between d18O and d57Fe is unlikely to be inherited from oceanic crust, as d57Fe values determined for samples of hydrothermally altered oceanic crust do not differ significantly from the mantle value and show no correlation with d18O. It is proposed that the correlated d57Fe and d18O variations in this particular eclogite suite are predominantly related to isotopic fractionation by disequilibrium partial melting although modification by melt percolation processes cannot be ruled out. Fractionation of Fe and O isotopes by removal of partial melt enriched in isotopically heavy Fe and O is supported by negative correlations between bulk sample d57Fe and Cr content and bulk sample and garnet d18O and Sc contents, as Cr and Sc are elements that become enriched in garnet- and pyroxene-bearing melt residues. Melt extraction could take place either during subduction, where the eclogites represent the residues of melted oceanic lithosphere, or could take place during long-term residence within the lithospheric mantle, in which case the protoliths of the eclogites could be of either crustal or mantle origin. This modification of both d57Fe and d18O by melting processes and specifically the production of low-d18O signatures in mafic rocks implies that some of the isotopically light d18O values observed in OIB and eclogite xenoliths may not necessarily reflect near-surface processes or components.

Chemical composition and age of metamorphic rocks from the Khavyven Highland, Eastern Kamchatka

Metamorphic rocks of the Khavyven Highland in eastern Kamchatka were determined to comprise two complexes of metavolcanic rocks that have different ages and are associated with subordinate amounts of metasediments. The complex composing the lower part of the visible vertical section of the highland is dominated by leucocratic amphibole-mica (+/-garnet) and epidote-mica (+/-garnet) crystalline schists, whose protoliths were andesites and dacites and their high-K varieties of island-arc calc-alkaline series. The other complex composing the upper part of the vertical section consists of spilitized basaltoids transformed into epidote-amphibole and phengite-epidote-amphibole green schists, which form (together with quartzites, serpentinized peridotites, serpentinites, and gabbroids) a sea-margin ophiolitic association. High LILE concentrations, high K/La, Ba/Th, Th/Ta, and La/Nb ratios, deep Ta-Nb minima, and low (La/Yb)_N and high 87Sr/86Sr ratios of the crystalline schists of the lower unit are demonstrated to testify to their subduction nature and suggest that their protolithic volcanics were produced in the suprasubduction environment of the Ozernoi-Valaginskii (Achaivayam-Valaginskii) island volcanic arc of Campanian-Paleogene age. The green schists of the upper unit show features of depleted MOR tholeiitic melts and subduction melts, which cause the deep Ta-Nb minima, and low K/La and 87Sr/86Sr ratios suggesting that the green schists formed in a marginal basin in front of the Ozernoi-Valaginskaya island arc. Recently obtained K-Ar ages in the Khavyven Highland vary from 32.4 to 39.3 Ma and indicate that metamorphism of the protolithic rocks occurred in Eocene under effect of collision and accretion processes of the arc complexes of the Ozernoi-Valaginskii and Kronotskii island arcs with the Asian continent and the closure of forearc oceanic basins in front of them. The modern position of the collision suture that marks the fossil subduction zone of the Ozernoi-Valaginskii arc and is spatially restricted to the buried Khavyven uplift in the Central Kamchatka Depression characterized by well-pronounced linear gravity anomalies.

Linear sedimentation rates and approximate mass accumulation rates at DSDP Site 61-462

A 560-meter-thick sequence of Cenomanian through Pleistocene sediments cored at DSDP Site 462 in the Nauru Basin overlies a 500-meter-thick complex unit of altered basalt flows, diabase sills, and thin intercalated volcaniclastic sediments. The Upper Cretaceous and Cenozoic sediments contain a high proportion of calcareous fossils, although the site has apparently been below the calcite compensation depth (CCD) from the late Mesozoic to the Pleistocene. This fact and the contemporaneous fluctuations of the calcite and opal accumulation rates suggest an irregular influx of displaced pelagic sediments from the shallow margins of the basin to its center, resulting in unusually high overall sedimentation rates for such a deep (5190 m) site. Shallow-water benthic fossils and planktonic foraminifers both occur as reworked materials, but usually are not found in the same intervals of the sediment section. We interpret this as recording separate erosional interludes in the shallow-water and intermediate-water regimes. Lower and upper Cenozoic hiatuses also are believed to have resulted from mid-water events. High accumulation rates of volcanogenic material during Santonian time suggest a corresponding significant volcanic episode. The coincidence of increased carbonate accumulation rates during the Campanian and displacement of shallow-water fossils during the late Campanian-early Maestrichtian with the volcanic event implies that this early event resulted in formation of the island chains around the Nauru Basin, which then served as platforms for initial carbonate deposition.

Geochemistry and mineralogy of ODP Leg 107 basalts

During ODP Leg 107, the basement of the Tyrrhenian Sea was drilled at Site 650, located in the Marsili basin, and at Sites 651 and 655, both located in the Vavilov basin. In addition, a lava flow was drilled at Site 654 on the Sardinia rifted margin. Mineral and whole rock major and trace element chemistry, including rare earth element (REE) and Sr and Nd isotopic ratios, were determined in samples of these rocks. Site 654 lava was sampled within uppermost Pliocene postrift sediments. This lava is a basaltic andesite of intraplate affinity, and is analogous to some Plio-Pleistocene tholeiitic lavas from Sardinia. Site 650 basalts, drilled beneath 1.7-1.9-Ma-old basal sediment, are strongly altered and vesicular suggesting a rapid subsidence of the Marsili basin. Based on incompatible trace elements, these basalts show calc-alkaline affinity like some products of the Marsili Seamount and the Eolian arc. The basement of the two sites drilled within Vavilov basin shows contrasting petrologies. Site 655, located along the Gortani ridge in the western part of the basin, drilled a 116-m-thick sequence of basalt flows beneath 3.4-3.6-Ma-old basal sediments. These basalts are chemically relatively homogeneous and show affinity to transitional MORB. Four units consisting of slightly differentiated basaltic lavas, have been identified. Site 655 basalts are geochemically similar to the high Ti lavas from DSDP Leg 42, Site 373 (Vavilov Basin). The basement at Site 651, overlain by 40 m of metalliferous dolostone covered by fossiliferous sediments with an age of 2 Ma, consists of two basalt units separated by a dolerite-albitite intrusive body; serpentinized harzburgites were drilled for 30 m at the base of the hole. The two basalt units of Site 651 are distinct petrochemically, though both show incompatible elements affinity with high-K calc-alkaline/calc-alkaline magmas from Eolian arc. The cpx chemistry and high K/Na ratio of the lower unit lavas suggest a weak alkaline tendency of potassic lineage. Leg 107 basement rock data, together with data from DSDP Site 373 and from dredged samples, indicate that the deepest basins of the central Tyrrhenian Sea are underlain by a complex back-arc basin crust produced by magmas with incompatible element affinities to transitional MORB (Site 655 and DSDP Site 373), and to calc-alkaline and high-K calc-alkaline converging plate margin basalts (Sites 650 and 651). This petrogenetic complexity is in accordance with the back-arc setting of the Vavilov and Marsili basins. Other back-arc basin basalts, particularly those from ensialic basins such as the Bransfield Strait (Antarctica), show a comparable petrogenetic complexity (cf., Sounders and Tarney, 1984).

Deep water formation in the North Pacific and deglacial CO2 rise

Deep water formation in the North Atlantic and Southern Ocean is widely thought to influence deglacial CO2 rise and climate change; here we suggest that deep water formation in the North Pacific may also play an important role. We present paired radiocarbon and boron isotope data from foraminifera from sediment core MD02-2489 at 3640 m in the North East Pacific. These show a pronounced excursion during Heinrich Stadial 1, with benthic-planktic radiocarbon offsets dropping to ~350 years, accompanied by a decrease in benthic d11B. We suggest this is driven by the onset of deep convection in the North Pacific, which mixes young shallow waters to depth, old deep waters to the surface, and low-pH water from intermediate depths into the deep ocean. This deep water formation event was likely driven by an increase in surface salinity, due to subdued atmospheric/monsoonal freshwater flux during Heinrich Stadial 1. The ability of North Pacific Deep Water (NPDW) formation to explain the excursions seen in our data is demonstrated in a series of experiments with an intermediate complexity Earth system model. These experiments also show that breakdown of stratification in the North Pacific leads to a rapid ~30 ppm increase in atmospheric CO2, along with decreases in atmospheric d13C and D14C, consistent with observations of the early deglaciation. Our inference of deep water formation is based mainly on results from a single sediment core, and our boron isotope data are unavoidably sparse in the key HS1 interval, so this hypothesis merits further testing. However we note that there is independent support for breakdown of stratification in shallower waters during this period, including a minimum in d15N, younging in intermediate water 14C, and regional warming. We also re-evaluate deglacial changes in North Pacific productivity and carbonate preservation in light of our new data, and suggest that the regional pulse of export production observed during the Bølling-Allerød is promoted by relatively stratified conditions, with increased light availability and a shallow, potent nutricline. Overall, our work highlights the potential of NPDW formation to play a significant and hitherto unrealized role in deglacial climate change and CO2 rise.

40Ar/39Ar single-crystal dating of detrital muscovite and K-feldspar of ODP Holes 116-717A and 116-718C

Detrital K-feldspars and muscovites from Ocean Drilling Program Leg 116 cores that have depositional ages from 0 to 18 Ma have been dated by the 40Ar/39Ar technique. Four to thirteen individual K-feldspars have been dated from seven stratigraphic levels, each of which have a very large range, up to 1660 Ma. At each level investigated, at least one K-feldspar yielded an age minimum which is, within uncertainty, identical to the age of deposition. One to twelve single muscovite crystals from each of six levels have also been studied. The range of muscovite ages is less than that of the K-feldspars and, with one exception, reveal only a 20-Ma spread in ages. As with the K-feldspars, each level investigated contains muscovites with mineral ages essentially identical to depositional ages. These results indicate that a significant portion of the material in the Bengal Fan is first-cycle detritus derived from the Himalayas. Therefore, the significant proportion of sediment deposited in the distal fan in the early to mid Miocene can be ascribed to a significant pulse of uplift and erosion in the collision zone. Moreover, these data indicate that during the entire Neogene, some portion of the Himalayan orogen was experiencing rapid erosion (<= uplift). The lack of granulite facies rocks in the eastern Himalayas and Tibetan Plateau suggests that very rapid uplift must have been distributed in brief pulses in different places in the mountain belt. We suggest that the great majority of the crystals with young apparent ages have been derived from the southern slope of the Himalayas, predominantly from near the main central thrust zone. These data provide further evidence against tectonic models in which the Himalayas and Tibetan plateaus are uplifted either uniformly during the past 40 m.y. or mostly within the last 2 to 5 m.y.

Strontium isotope ratios of planktonic foraminifera from sediment samples acroos the K/T boundary

Seawater 87Sr/86Sr values increase abruptly by 28 * 10**-6 across the Cretaceous/Tertiary boundary (KTB). This small, but rapid shift is superimposed on the larger scale structure of the seawater Sr isotope curve. The time scale of radiogenic Sr addition appears to be too rapid to reconcile with sources associated with volcanism, and we show that the amount of Sr required to produce even this small increase is too large to be derived from: (1) a KT bolide of the size constrained by the Ir anomaly, (2) continental crust ejecta from the impact of such a bolide, (3) soot from global wildfires initiated by an impact, or (4) any combination of these sources. The probable source of the radiogenic Sr is enhanced continental weathering, but the high rate of increase appears to rule out processes such as sea level regression, glaciation or tectonism. A plausible mechanism for rapid addition of radiogenic Sr to the oceans is enhanced weathering associated with globally distributed acid rain (pH c. 1) which is a proposed by-product of a bolide impact (Prinn and Fegley, 1987, doi:10.1016/0012-821X(87)90046-X).