The behaviour of incompatible elements during hydrous melting of metasomatized peridotite at 4–6 GPa and 1000 °C–1200 °C

Trace element behavior during hydrous melting of a metasomatized garnet–peridotite was examined at pressures of 4–6 GPa and temperatures of 1000 °C–1200 °C, conditions appropriate for fluid penetrating the mantle wedge atop the subducting slab. Experiments were performed in a rocking multi-anvil apparatus using a diamond-trap setup. The compositions of the fluid and melt phases were measured using the cryogenic LA-ICP-MS technique. The water-saturated solidus of the K-lherzolite composition is located between 900 °C and 1000 °C at 4 GPa and between 1000 °C and 1100 °C at 5 and 6 GPa. The partition coefficients between fluid or melt and clinopyroxene reveal an asymmetric MREE trough with a minimum at Dy. The clinopyroxene in equilibrium with aqueous fluids is characterized by DUfluid–cpx > DThfluid–cpx while DUmelt–cpx tends to be similar to DThmelt–cpx. The partition coefficients between fluid or melt and garnet reveal very strong light to heavy REE fractionation, DLa/DLu from 95 (hydrous melt) to 1600 (aqueous fluid). The LILE are highly incompatible with partition coefficients > 50. The behavior of HFSE are decoupled, with DZr,Hf close to 1 while DNb,Ta > 10. Garnet is characterized by DUmelt/fluid–garnet < DThmelt/fluid–garnet. A comparison of our experimental partitioning results for trivalent cations as well as the results from the literature and the calculations carried out using the lattice strain model adapted to the presence of water in the bulk system indicates that H2O in the fluid or melt phase has a prominent effect on trace element partitioning. Garnet in mantle rocks in equilibrium with an aqueous fluid is characterized by significantly higher Do(3 +) for REE in the X site of the garnet compared with the partitioning values of the optimal cation in garnet in equilibrium with hydrous melts. Our data show for the first time that the change in the nature of the mobile phase (fluid vs. melt) does affect the affinities of trace elements into the garnet crystal at conditions below the second critical endpoint of the system. The same also applies for clinopyroxene, although this is less clear. Consequently, our new data allow for refinements in predictive modeling of element transfer from the slab to the mantle wedge and of possible compositions of metasomatized mantle that sources OIB magmatism.

Density-dependent dynamics of a dominant rain forest tree change with juvenile stage and time of masting

Although negative density dependence (NDD) can facilitate tree species coexistence in forests, the underlying mechanisms can differ, and rarely are the dynamics of seedlings and saplings studied together. Herein we present and discuss a novel mechanism based on our investigation of NDD predictions for the large, grove-forming ectomycorrhizal mast fruiting tree, Microberlinia bisulcata (Caesalpiniaceae), in an 82.5-ha plot at Korup, Cameroon. We tested whether juvenile density, size, growth and survival decreases with increasing conspecific adult basal area for 3245 ‘new’ seedlings and 540 ‘old’ seedlings (< 75-cm tall) during an approximately 4-year study period (2008–2012) and for 234 ‘saplings’ (≥ 75-cm tall) during an approximately 6-year study period (2008–2014). We found that the respective densities of new seedlings, old seedlings and saplings were positively, not and negatively related to increasing BA. Maximum leaf numbers and heights of old seedlings were negatively correlated with increasing basal areas, as were sapling heights and stem diameters. Whereas survivorship of new seedlings decreased by more than one-half with increasing basal area over its range in 2010–2012, that of old seedlings decreased by almost two-thirds, but only in 2008–2010, and was generally unrelated to conspecific seedling density. In 2010–2012 relative growth rates in new seedlings’ heights decreased with increasing basal area, as well as with increasing seedling density, together with increasing leaf numbers, whereas old seedlings’ growth was unrelated to either conspecific density or basal area. Saplings of below-average height had reduced survivorship with increasing basal area (probability decreasing from approx. 0.4 to 0.05 over the basal area range tested), but only sapling growth in terms of leaf numbers decreased with increasing basal area. These static and dynamic results indicate that NDD is operating within this system, possibly stabilizing the M. bisulcata population. However, these NDD patterns are unlikely to be caused by symmetric competition or by consumers. Instead, an alternative mechanism for conspecific adult–juvenile negative feedback is proposed, one which involves the interaction between tree phenology and ectomycorrhizal linkages.

Analysis of recurrent failure times: A time-dependent Yule process approach

Analysis of recurrent events has been widely discussed in medical, health services, insurance, and engineering areas in recent years. This research proposes to use a nonhomogeneous Yule process with the proportional intensity assumption to model the hazard function on recurrent events data and the associated risk factors. This method assumes that repeated events occur for each individual, with given covariates, according to a nonhomogeneous Yule process with intensity function λx(t) = λ 0(t) · exp( x′β). One of the advantages of using a non-homogeneous Yule process for recurrent events is that it assumes that the recurrent rate is proportional to the number of events that occur up to time t. Maximum likelihood estimation is used to provide estimates of the parameters in the model, and a generalized scoring iterative procedure is applied in numerical computation. ^ Model comparisons between the proposed method and other existing recurrent models are addressed by simulation. One example concerning recurrent myocardial infarction events compared between two distinct populations, Mexican-American and Non-Hispanic Whites in the Corpus Christi Heart Project is examined. ^

Age -period -cohort analysis: An extension of Cox's lifetable regression with time -dependent covariates

It is well known that an identification problem exists in the analysis of age-period-cohort data because of the relationship among the three factors (date of birth + age at death = date of death). There are numerous suggestions about how to analyze the data. No one solution has been satisfactory. The purpose of this study is to provide another analytic method by extending the Cox's lifetable regression model with time-dependent covariates. The new approach contains the following features: (1) It is based on the conditional maximum likelihood procedure using a proportional hazard function described by Cox (1972), treating the age factor as the underlying hazard to estimate the parameters for the cohort and period factors. (2) The model is flexible so that both the cohort and period factors can be treated as dummy or continuous variables, and the parameter estimations can be obtained for numerous combinations of variables as in a regression analysis. (3) The model is applicable even when the time period is unequally spaced.^ Two specific models are considered to illustrate the new approach and applied to the U.S. prostate cancer data. We find that there are significant differences between all cohorts and there is a significant period effect for both whites and nonwhites. The underlying hazard increases exponentially with age indicating that old people have much higher risk than young people. A log transformation of relative risk shows that the prostate cancer risk declined in recent cohorts for both models. However, prostate cancer risk declined 5 cohorts (25 years) earlier for whites than for nonwhites under the period factor model (0 0 0 1 1 1 1). These latter results are similar to the previous study by Holford (1983).^ The new approach offers a general method to analyze the age-period-cohort data without using any arbitrary constraint in the model. ^

The number and timing of time -dependent covariate measurements for the Cox regression model

The problem of analyzing data with updated measurements in the time-dependent proportional hazards model arises frequently in practice. One available option is to reduce the number of intervals (or updated measurements) to be included in the Cox regression model. We empirically investigated the bias of the estimator of the time-dependent covariate while varying the effect of failure rate, sample size, true values of the parameters and the number of intervals. We also evaluated how often a time-dependent covariate needs to be collected and assessed the effect of sample size and failure rate on the power of testing a time-dependent effect.^ A time-dependent proportional hazards model with two binary covariates was considered. The time axis was partitioned into k intervals. The baseline hazard was assumed to be 1 so that the failure times were exponentially distributed in the ith interval. A type II censoring model was adopted to characterize the failure rate. The factors of interest were sample size (500, 1000), type II censoring with failure rates of 0.05, 0.10, and 0.20, and three values for each of the non-time-dependent and time-dependent covariates (1/4,1/2,3/4).^ The mean of the bias of the estimator of the coefficient of the time-dependent covariate decreased as sample size and number of intervals increased whereas the mean of the bias increased as failure rate and true values of the covariates increased. The mean of the bias of the estimator of the coefficient was smallest when all of the updated measurements were used in the model compared with two models that used selected measurements of the time-dependent covariate. For the model that included all the measurements, the coverage rates of the estimator of the coefficient of the time-dependent covariate was in most cases 90% or more except when the failure rate was high (0.20). The power associated with testing a time-dependent effect was highest when all of the measurements of the time-dependent covariate were used. An example from the Systolic Hypertension in the Elderly Program Cooperative Research Group is presented. ^

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Seismic velocities and elastic properties of gabbroic rocks from ODP Hole 118-735B

The nearly continuous recovery of 0.5 km of generally fresh, layer 3 gabbroic rocks at Hole 735B, especially near the bottom of the section, presents scientists an unusual opportunity to study the detailed elastic properties of the lower oceanic crust. Extending compressional-wave and density shipboard measurements at room pressure, Vp and Vs were measured at pressures from 20 to 200 MPa using the pulse transmission method. All of the rocks exhibit significant increases in velocity with increasing pressure up to about 150 MPa, a feature attributed to the closing of microcrack porosity. Measured velocities reflect the mineralogical makeup and microstructures acquired during the tectonic history of Hole 735B. Most of the undeformed and unaltered gabbros are approximately 65:35 plagioclase/clinopyroxene rocks plus olivine or oxide minerals, and the observed densities and velocities are fully consistent with the Voigt-Reuss-Hill (VRH) averages of the component minerals and their proportions. Depending on their olivine content, the predominant olivine gabbros at 200 MPa have average Vp = 7.1 ± 0.2 km/s, Vs = 3.9 ± 0.1 km/s, and grain densities of 2.95 ± 0.5 g/cm3. The less abundant iron-titanium (Fe-Ti) oxide gabbros average Vp = 6.75 ± 0.15 km/s, Vs = 3.70 ± 0.1 km/s, and grain densities of 3.22 ± 0.05 g/cm3, reflecting the higher densities and lower velocities of oxide minerals compared to olivine. About 30% of the core is plastically deformed, and the densities and directionally averaged velocities of these shear-zone tectonites are generally consistent with those of the gabbros, their protoliths. Three sets of observations indicate that the shear-zone metagabbros are elastically anisotropic: (1) directional variations in Vp, both vertical and horizontal and with respect to foliation and lineation; (2) discrepancies among Vp values for the horizontal cores and the VRH averages of the component minerals and their mineral proportions, suggesting preferred crystallographic orientations of anisotropic minerals; and (3) variations of Vs of up to 7%, with polarization directions parallel and perpendicular to foliation. Optical inspection of thin sections of the same samples indicates that plagioclase feldspar, clinopyroxene, and amphibole typically display crystallographic-preferred orientations, and this, plus the elastic anisotropy of these minerals, suggests that preferred orientations are responsible for much of the observed anisotropy, particularly at high pressure. Alteration tends to be localized to brittle faults and brecciated zones, and typical alteration minerals are amphibole and secondary plagioclase, which do not significantly change the velocity-density relationships.

Geochemistry of vein rocks of ODP Hole 118-735B

Rock samples from Hole 735B, Southwest Indian Ridge, were examined to determine the principal vein-related types of alteration that occurred, the nature of fluids that were present, and the temperatures and pressures of these fluids. Samples studied included veined metagabbro, veined mylonitic metagabbro, felsic trondhjemite, and late-stage leucocratic diopside-bearing veins. The methods used were standard petrographic analysis, mineral chemical analysis by electron microprobe, fluid inclusion petrography and analysis by heating/freezing techniques and laser Raman microspectroscopy, and oxygen isotopic analyses of mineral separates. Alteration in lithologic Units I and II (above the level of Core 118-735B-3OR; approximately 140 meters below the seafloor) is dominated by hydration by seawater-derived fluids at high temperature, up to about 700°C, and low water/rock ratio, during and immediately after pervasive ductile deformation. Below Core 118-735B-30R, pervasive deformation is less common, and brittle veining and brecciation are the major alteration styles. Leucocratic centimeter-scale veins, often containing diopside and plagioclase, were produced by interaction of hot (about 500°C) seawater-derived fluid and gabbro. The water/rock ratio was locally high at the veins and breccia zones, but the integrated water/rock ratio for the lower part of the hole is probably low. Accessory hydrous magmatic or deuteric phases formed from magmatic volatiles in some gabbro and in trondhjemite. Most subsequent alteration was affected by fluids that were seawater-derived, based on isotopic and chemical analyses of minerals and analyses of fluid inclusions. Many early-generation fluid inclusions, associated with high-temperature veining, contain appreciable methane as well as saline water. The source of methane is unclear, but it may have formed as seawater was reduced during low water/rock interaction with ultramafic upper mantle or ultramafic and mafic layer 3. Temperatures of alteration were calculated on the basis of coexisting mineral chemistry and isotopic values. Hydrothermal metamorphism commenced at about 720°C and continued to about 550°C. Leucocratic veining took place at about 500°C. Alteration within brecciated horizons was also at about 500° to less than 400°C, and the trondhjemite was altered at about 550° to below 490°C. Pressures calculated from a diopside-bearing vein, based on a combination of fluid inclusion and isotopic analysis, were 90 to 100 MPa. This pressure places the sample, from Core 118-735B-70R in Unit V, at about 2 km below the seafloor.

Grain-size distribution and frictional behaviour of sediment cores from the Nice Slope

The upper shelf of the landslide-prone Ligurian Margin (Western Mediterranean Sea) off Nice well-known for the 1979 Airport Landslide is a natural laboratory to study preconditioning factors and trigger mechanisms for submarine landslides. For this study low-stress ring shear experiments have been carried out on a variety of sediments from >50 gravity cores to characterise the velocity-dependent frictional behaviour. Mean values of the peak coefficient of friction vary from 0.46 for clay-dominated samples (53 % clay, 46 % silt, 1 %) sand up to 0.76 for coarse-grained sediments (26 % clay, 57 % silt, 17 % sand). The majority of the sediments tested show velocity strengthening regardless of the grain size distribution. For clayey sediments the peak and residual cohesive strength increases with increasing normal stress, with values from 1.3 to 10.6 kPa and up to 25 % of all strength supported by cohesive forces in the shallowmost samples. A pseudo-static slope stability analysis reveals that the different lithologies (even clay-rich material with clay content >=50 %) tested are stable up to slope angles <26° under quasi-drained conditions.

Temperature and time dependent threshold voltage characterization of AlGaN/GaN high electron mobility transistors

Relacionado con línea de investigación del GDS del ISOM ver http://www.isom.upm.es/dsemiconductores.php

A method for revealing phase space structures of general time dependent dynamical systems

In this paper we develop new techniques for revealing geometrical structures in phase space that are valid for aperiodically time dependent dynamical systems, which we refer to as Lagrangian descriptors. These quantities are based on the integration, for a finite time, along trajectories of an intrinsic bounded, positive geometrical and/or physical property of the trajectory itself. We discuss a general methodology for constructing Lagrangian descriptors, and we discuss a “heuristic argument” that explains why this method is successful for revealing geometrical structures in the phase space of a dynamical system. We support this argument by explicit calculations on a benchmark problem having a hyperbolic fixed point with stable and unstable manifolds that are known analytically. Several other benchmark examples are considered that allow us the assess the performance of Lagrangian descriptors in revealing invariant tori and regions of shear. Throughout the paper “side-by-side” comparisons of the performance of Lagrangian descriptors with both finite time Lyapunov exponents (FTLEs) and finite time averages of certain components of the vector field (“time averages”) are carried out and discussed. In all cases Lagrangian descriptors are shown to be both more accurate and computationally efficient than these methods. We also perform computations for an explicitly three dimensional, aperiodically time-dependent vector field and an aperiodically time dependent vector field defined as a data set. Comparisons with FTLEs and time averages for these examples are also carried out, with similar conclusions as for the benchmark examples.

Lagrangian descriptors: a method for revealing phase space structures of general time dependent dynamical systems

Lagrangian descriptors are a recent technique which reveals geometrical structures in phase space and which are valid for aperiodically time dependent dynamical systems. We discuss a general methodology for constructing them and we discuss a "heuristic argument" that explains why this method is successful. We support this argument by explicit calculations on a benchmark problem. Several other benchmark examples are considered that allow us to assess the performance of Lagrangian descriptors with both finite time Lyapunov exponents (FTLEs) and finite time averages of certain components of the vector field ("time averages"). In all cases Lagrangian descriptors are shown to be both more accurate and computationally efficient than these methods.

Nonlinear seismic behaviour of concrete arch bridges

Arch bridge structural solution has been known for centuries, in fact the simple nature of arch that require low tension and shear strength was an advantage as the simple materials like stone and brick were the only option back in ancient centuries. By the pass of time especially after industrial revolution, the new materials were adopted in construction of arch bridges to reach longer spans. Nowadays one long span arch bridge is made of steel, concrete or combination of these two as "CFST", as the result of using these high strength materials, very long spans can be achieved. The current record for longest arch belongs to Chaotianmen bridge over Yangtze river in China with 552 meters span made of steel and the longest reinforced concrete type is Wanxian bridge which also cross the Yangtze river through a 420 meters span. Today the designer is no longer limited by span length as long as arch bridge is the most applicable solution among other approaches, i.e. cable stayed and suspended bridges are more reasonable if very long span is desired. Like any super structure, the economical and architectural aspects in construction of a bridge is extremely important, in other words, as a narrower bridge has better appearance, it also require smaller volume of material which make the design more economical. Design of such bridge, beside the high strength materials, requires precise structural analysis approaches capable of integrating the combination of material behaviour and complex geometry of structure and various types of loads which may be applied to bridge during its service life. Depend on the design strategy, analysis may only evaluates the linear elastic behaviour of structure or consider the nonlinear properties as well. Although most of structures in the past were designed to act in their elastic range, the rapid increase in computational capacity allow us to consider different sources of nonlinearities in order to achieve a more realistic evaluations where the dynamic behaviour of bridge is important especially in seismic zones where large movements may occur or structure experience P - _ effect during the earthquake. The above mentioned type of analysis is computationally expensive and very time consuming. In recent years, several methods were proposed in order to resolve this problem. Discussion of recent developments on these methods and their application on long span concrete arch bridges is the main goal of this research. Accordingly available long span concrete arch bridges have been studied to gather the critical information about their geometrical aspects and properties of their materials. Based on concluded information, several concrete arch bridges were designed for further studies. The main span of these bridges range from 100 to 400 meters. The Structural analysis methods implemented in in this study are as following: Elastic Analysis: Direct Response History Analysis (DRHA): This method solves the direct equation of motion over time history of applied acceleration or imposed load in linear elastic range. Modal Response History Analysis (MRHA): Similar to DRHA, this method is also based on time history, but the equation of motion is simplified to single degree of freedom system and calculates the response of each mode independently. Performing this analysis require less time than DRHA. Modal Response Spectrum Analysis (MRSA): As it is obvious from its name, this method calculates the peak response of structure for each mode and combine them using modal combination rules based on the introduced spectra of ground motion. This method is expected to be fastest among Elastic analysis. Inelastic Analysis: Nonlinear Response History Analysis (NL-RHA): The most accurate strategy to address significant nonlinearities in structural dynamics is undoubtedly the nonlinear response history analysis which is similar to DRHA but extended to inelastic range by updating the stiffness matrix for every iteration. This onerous task, clearly increase the computational cost especially for unsymmetrical buildings that requires to be analyzed in a full 3D model for taking the torsional effects in to consideration. Modal Pushover Analysis (MPA): The Modal Pushover Analysis is basically the MRHA but extended to inelastic stage. After all, the MRHA cannot solve the system of dynamics because the resisting force fs(u; u_ ) is unknown for inelastic stage. The solution of MPA for this obstacle is using the previously recorded fs to evaluate system of dynamics. Extended Modal Pushover Analysis (EMPA): Expanded Modal pushover is a one of very recent proposed methods which evaluates response of structure under multi-directional excitation using the modal pushover analysis strategy. In one specific mode,the original pushover neglect the contribution of the directions different than characteristic one, this is reasonable in regular symmetric building but a structure with complex shape like long span arch bridges may go through strong modal coupling. This method intend to consider modal coupling while it take same time of computation as MPA. Coupled Nonlinear Static Pushover Analysis (CNSP): The EMPA includes the contribution of non-characteristic direction to the formal MPA procedure. However the static pushovers in EMPA are performed individually for every mode, accordingly the resulted values from different modes can be combined but this is only valid in elastic phase; as soon as any element in structure starts yielding the neutral axis of that section is no longer fixed for both response during the earthquake, meaning the longitudinal deflection unavoidably affect the transverse one or vice versa. To overcome this drawback, the CNSP suggests executing pushover analysis for governing modes of each direction at the same time. This strategy is estimated to be more accurate than MPA and EMPA, moreover the calculation time is reduced because only one pushover analysis is required. Regardless of the strategy, the accuracy of structural analysis is highly dependent on modelling and numerical integration approaches used in evaluation of each method. Therefore the widely used Finite Element Method is implemented in process of all analysis performed in this research. In order to address the study, chapter 2, starts with gathered information about constructed long span arch bridges, this chapter continuous with geometrical and material definition of new models. Chapter 3 provides the detailed information about structural analysis strategies; furthermore the step by step description of procedure of all methods is available in Appendix A. The document ends with the description of results and conclusion of chapter 4.

Control of time-dependent biological processes by temporally patterned input

Temporal patterning of biological variables, in the form of oscillations and rhythms on many time scales, is ubiquitous. Altering the temporal pattern of an input variable greatly affects the output of many biological processes. We develop here a conceptual framework for a quantitative understanding of such pattern dependence, focusing particularly on nonlinear, saturable, time-dependent processes that abound in biophysics, biochemistry, and physiology. We show theoretically that pattern dependence is governed by the nonlinearity of the input–output transformation as well as its time constant. As a result, only patterns on certain time scales permit the expression of pattern dependence, and processes with different time constants can respond preferentially to different patterns. This has implications for temporal coding and decoding, and allows differential control of processes through pattern. We show how pattern dependence can be quantitatively predicted using only information from steady, unpatterned input. To apply our ideas, we analyze, in an experimental example, how muscle contraction depends on the pattern of motorneuron firing.