Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

Optimal exponent heat balance and refined integral methods applied to Stefan problems

When using a polynomial approximating function the most contentious aspect of the Heat Balance Integral Method is the choice of power of the highest order term. In this paper we employ a method recently developed for thermal problems, where the exponent is determined during the solution process, to analyse Stefan problems. This is achieved by minimising an error function. The solution requires no knowledge of an exact solution and generally produces significantly better results than all previous HBI models. The method is illustrated by first applying it to standard thermal problems. A Stefan problem with an analytical solution is then discussed and results compared to the approximate solution. An ablation problem is also analysed and results compared against a numerical solution. In both examples the agreement is excellent. A Stefan problem where the boundary temperature increases exponentially is analysed. This highlights the difficulties that can be encountered with a time dependent boundary condition. Finally, melting with a time-dependent flux is briefly analysed without applying analytical or numerical results to assess the accuracy.

Improving the accuracy of heat balance integral methods applied to thermal problems with time dependent boundary conditions

In this paper the two main drawbacks of the heat balance integral methods are examined. Firstly we investigate the choice of approximating function. For a standard polynomial form it is shown that combining the Heat Balance and Refined Integral methods to determine the power of the highest order term will either lead to the same, or more often, greatly improved accuracy on standard methods. Secondly we examine thermal problems with a time-dependent boundary condition. In doing so we develop a logarithmic approximating function. This new function allows us to model moving peaks in the temperature profile, a feature that previous heat balance methods cannot capture. If the boundary temperature varies so that at some time t & 0 it equals the far-field temperature, then standard methods predict that the temperature is everywhere at this constant value. The new method predicts the correct behaviour. It is also shown that this function provides even more accurate results, when coupled with the new CIM, than the polynomial profile. Analysis primarily focuses on a specified constant boundary temperature and is then extended to constant flux, Newton cooling and time dependent boundary conditions.

Adjuvant requirement for successful immunization with recombinant derivatives of Plasmodium vivax merozoite surface protein-1 delivered via the intranasal route

Recently, we generated two bacterial recombinant proteins expressing 89 amino acids of the C-terminal domain of the Plasmodium vivax merozoite surface protein-1 and the hexa-histidine tag (His6MSP1(19)). One of these recombinant proteins contained also the amino acid sequence of the universal pan allelic T-cell epitope (His6MSP1(19)-PADRE). In the present study, we evaluated the immunogenic properties of these antigens when administered via the intra-nasal route in the presence of distinct adjuvant formulations. We found that C57BL/6 mice immunized with either recombinant proteins in the presence of the adjuvants cholera toxin (CT) or the Escherichia coli heat labile toxin (LT) developed high and long lasting titers of specific serum antibodies. The induced immune responses reached maximum levels after three immunizing doses with a prevailing IgG1 subclass response. In contrast, mice immunized by intranasal route with His6MSP1(19)-PADRE in the presence of the synthetic oligonucleotides adjuvant CpG ODN 1826 developed lower antibody titers but when combined to CT, CpG addition resulted in enhanced IgG responses characterized by lower IgG1 levels. Considering the limitations of antigens formulations that can be used in humans, mucosal adjuvants can be a reliable alternative for the development of new strategies of immunization using recombinant proteins of P. vivax.

Caractérisation des microparticules dans les concentrés érythrocytaires : une analyse des antigènes de groupe sanguins par cytométrie de flux

Cette thèse a pour but de caractériser les microparticules isolées à partir des concentrés érythrocytaires et plus précisément de déterminer la présence d'antigènes de groupes sanguins à leur surface. Elle est divisée en trois parties, sous forme d'articles publiés à partir d'un travail de recherche mené au Centre de Transfusion Sanguine d'Epalinges. Dans l'article « Microparticles in stored red blood cells : an approach using flow cytometry and proteomic tools » publié dans Vox Sanguinis en 2008, il est question de la lésion de stockage des globules rouges. Grâce à des techniques alliant la cytométrie de flux et la protéomique, il a été montré que la génération de microparticules augmente au cours du stockage des concentrés érythrocytaires et que leur composition se modifie au cours du temps. L'article de revue « Analysis and clinical relevance of microparticles from red blood cells » publié dans Current Opinion in Hematology en 2010, explique les mécanismes de formation et d'élimination des microparticules de globules rouges. Il fait une revue des implications cliniques liées à la génération de microparticules et discute leur conséquences potentielles dans le domaine de la médecine transfusionnelle. L'article « Red blood cell microparticles and blood group antigens : an analysis by flow cytometry » publié dans Blood Transfusion en 2012, décrit l'étude des antigènes de groupe sanguins à la surface des microparticules générées à partir de concentrés érythrocytaires après ajout de calcium ionophore. Les résultats de cette étude indiquent que les antigènes de groupes sanguins appartenant aux systèmes RH, KEL, JK, FY, MNS, LE et LU sont présents à la surface des microparticules. Ces antigènes pourraient potentiellement être source d'allo-immunisation après transfusion.

States of maximum entropy production in a onedimensional vertical model with convective adjustment

We investigate the hypothesis that the atmosphere is constrained to maximize its entropy production by using a one-dimensional (1-D) vertical model. We prescribe the lapse rate in the convective layer as that of the standard troposphere. The assumption that convection sustains a critical lapse rate was absent in previous studies, which focused on the vertical distribution of climatic variables, since such a convective adjustment reduces the degrees of freedom of the system and may prevent the application of the maximum entropy production (MEP) principle. This is not the case in the radiative–convective model (RCM) developed here, since we accept a discontinuity of temperatures at the surface similar to that adopted in many RCMs. For current conditions, the MEP state gives a difference between the ground temperature and the air temperature at the surface ≈10 K. In comparison, conventional RCMs obtain a discontinuity ≈2 K only. However, the surface boundary layer velocity in the MEP state appears reasonable (≈3 m s-¹). Moreover, although the convective flux at the surface in MEP states is almost uniform in optically thick atmospheres, it reaches a maximum value for an optical thickness similar to current conditions. This additional result may support the maximum convection hypothesis suggested by Paltridge (1978)

Cholesterol and Fatty Acids Regulate Dynamic Caveolin Trafficking through the Golgi Complex and between the Cell Surface and Lipid Bodies

Caveolins are a crucial component of plasma membrane (PM) caveolae but have also been localized to intracellular compartments, including the Golgi complex and lipid bodies. Mutant caveolins associated with human disease show aberrant trafficking to the PM and Golgi accumulation. We now show that the Golgi pool of mainly newly synthesized protein is detergent-soluble and predominantly in a monomeric state, in contrast to the surface pool. Caveolin at the PM is not recognized by specific caveolin antibodies unless PM cholesterol is depleted. Exit from the Golgi complex of wild-type caveolin-1 or -3, but not vesicular stomatitis virus-G protein, is modulated by changing cellular cholesterol levels. In contrast, a muscular dystrophy-associated mutant of caveolin-3, Cav3P104L, showed increased accumulation in the Golgi complex upon cholesterol treatment. In addition, we demonstrate that in response to fatty acid treatment caveolin can follow a previously undescribed pathway from the PM to lipid bodies and can move from lipid bodies to the PM in response to removal of fatty acids. The results suggest that cholesterol is a rate-limiting component for caveolin trafficking. Changes in caveolin flux through the exocytic pathway can therefore be an indicator of cellular cholesterol and fatty acid levels.

Design and construction of precision heat fluxmeters

InAlAs/InGaAs/InP based high electron mobility transistor devices have been structurally and electrically characterized, using transmission electron microscopy and Raman spectroscopy and measuring Hall mobilities. The InGaAs lattice matched channels, with an In molar fraction of 53%, grown at temperatures lower than 530¿°C exhibit alloy decomposition driving an anisotropic InGaAs surface roughness oriented along [1math0]. Conversely, lattice mismatched channels with an In molar fraction of 75% do not present this lateral decomposition but a strain induced roughness, with higher strength as the channel growth temperature increases beyond 490¿°C. In both cases the presence of the roughness implies low and anisotropic Hall mobilities of the two dimensional electron gas.

Influence of premetallization surface treatment on the formation of Schottky Au-nGaN contacts

The influence of premetallization surface preparation on the structural, chemical, and electrical properties of Au-nGaN interfaces has been investigated by x-ray photoemission spectroscopy (XPS), current-voltage measurement (I-V) and cross-section transmission electron microscopy (TEM). XPS analysis showed that the three GaN substrate treatments investigated i.e., ex situ hydrofluoric acid etch, in situ anneal in ultrahigh-vacuum (UHV), and in situ Ga reflux cleaning in UHV result in surfaces increasingly free of oxygen contamination. XPS and TEM characterization of Au-nGaN formed after the three premetallization surface treatments show that HF etching and UHV annealing produce abrupt, well-defined interfaces. Conversely, GaN substrate cleaning in a Ga flux results in Au/GaN intermixing. I-V characterization of Au¿nGaN contacts yields a Schottky barrier height of 1.25 eV with a very low-ideality factor and very good contact uniformity for the premetallization UHV anneal, while the Ga reflux cleaning results in a much lower barrier (0.85 eV), with poor ideality and uniformity. I-V and XPS results suggest a high density of acceptor states at the surface, which is further enhanced by UHV annealing. These results are discussed in the context of current models of Schottky barrier formation.

Nitrogen dynamics in soil management systems. I - flux of inorganic nitrogen (NH4+ and NO3-)

In agricultural systems the N-NH4+ and N-NO3- contents is significantly affected by soil management. This study investigated the dynamics of inorganic nitrogen (N; NH4+ and NO3-) in an experimental evaluation of soil management systems (SMSs) adopted in 1988 at the experimental station of the ABC Foundation in Ponta Grossa, in the Central South region of the State of Paraná. The objective of this study was to evaluate the changes in N-NH4+ and N-NO3- flux in the surface layer of a Red Latosol arising from SMSs over a 12-month period. The experiment was arranged in a completely randomized block design in split plots, in three replications. The plots consisted of the following SMSs: 1) conventional tillage (CT); 2) minimum tillage (MT); 3) no-tillage with chisel plow every three years (NT CH); and 4) continuous no-tillage (CNT). To evaluate the dynamics of inorganic N, the subplots represented samplings (11 sampling times, T1 - T11). The ammonium N (N-NH4+) and nitric N (N-NO3-) contents were higher in systems with reduced tillage (MT and NT CH) and without tillage (CNT) than in the CT system. In the period from October 2003 to February 2004, the N-NH4+ was higher than the N-NO3- soil content. Conversely, in the period from May 2004 to July 2004, the N-NO3- was higher than the N-NH4+ content. The greatest fluctuation in the N-NH4+ and N-NO3- contents occurred in the 0-2.5 cm layer, and the highest peak in the N-NH4+ and N-NO3- concentrations occurred after the surface application of N. Both N-NH4+ and N-NO3- were strongly correlated with the soil organic C content, which indicated that these properties vary together in the system.

Soil CO2 flux: a method comparison of closed static chambers in a sugarcane field

A large variety of techniques have been used to measure soil CO2 released from the soil surface, and much of the variability observed between locations must be attributed to the different methods used by the investigators. Therefore, a minimum protocol of measurement procedures should be established. The objectives of this study were (a) to compare different absorption areas, concentrations and volumes of the alkali trapping solution used in closed static chambers (CSC), and (b) to compare both, the optimized alkali trapping solution and the soda-lime trapping using CSC to measure soil respiration in sugarcane areas. Three CO2 absorption areas were evaluated (7; 15 and 20 % of the soil emission area or chamber); two volumes of NaOH (40 and 80 mL) at three concentrations (0.1, 0.25 and 0.5 mol L-1). Three different types of alkaline traps were tested: (a), 80 mL of 0.5 mol L-1 NaOH in glass containers, absorption area 15 % (V0.5); (b) 40 mL of 2 mol L-1 NaOH retained in a sponge, absorption area 80 % (S2) and (c) 40 g soda lime, absorption area 15 % (SL). NaOH concentrations of 0.5 mol L-1 or lower underestimated the soil CO2-C flux or CO2 flux. The lower limit of the alkali trap absorption area should be a minimum of 20 % of the area covered by the chamber. The 2 mol L-1 NaOH solution trap (S2) was the most efficient (highest accuracy and highest CO2 fluxes) in measuring soil respiration.

Magnetization versus heat treatment in rapidly solidified NdFeB alloys

NdFeB melt-spun amorphous or partially amorphous alloys of four compositions were prepared. Their crystallization kinetics induced by thermal treatment was studied by differential scanning calorimetry and scanning and transmission electron microscopy. Scanning electron microscopy demonstrated that heterogeneous nucleation occurs preferentially at the ribbon surface which was in contact with the wheel. The explicit form of the kinetic equation that best describes the first stage of crystallization under high undercooling conditions was obtained for each alloy. From the crystallization results, the lower part of the experimental time-temperature-transformation curves was deduced for each alloy and extrapolated up to the high-temperature limit of their validity. Microstructural observations showed a typical size of the microcrystals obtained by heat treatment of ~100 nm. From the magnetic properties measured with a vibrating sample magnetometer, the same magnetic behavior of partially crystallized alloys is observed regardless of the temperature of annealing provided the same crystallization fraction, x, is achieved, at least for small values of x (typically ~10%).

Modification of Ca-Montmorillonite by Low-Temperature Heat Treatment, HR-106, 1964

Objectives of this investigation were to measure the effects of moderate heat treatments (below the dehydroxylation temperature) on physical and chemical properties of a calcium-montmorillonite clay. Previous workers have noted the reduction in cation exchange capacity and swelling property after heating in the range 200 to 400°C, and have suggested several possible explanations, such as hysteresis effect, increased inter-layer attractions due to removal of inter-layer water, or changes in the disposition of inter-layer or layer surface ions. The liquid limits of Ca-montmorillonite were steadily decreased with increased temperature of treatment, levelling at about 450°C. The plastic limit decreased slightly up to 350°C, above which samples could no longer be rolled into threads. The gradual change is in contrast with sudden major changes noted for weight loss (maximum rates of change at l00°C and 500°C), glycol retention surface area (520°C), and d001 diffraction peak intensity (17.7 A spacing) and breadth after glycolation (530°C). Other properties showing more gradual reductions with heat treatment were amount of exchangeable calcium (without water soaking), cation exchange capacity by NH4AC method, and d001 intensity (21 A spacing) after storing at 100% r.h. one month and re-wetting with water. Previous water soaking allowed much greater release of fixed Ca++ up to 450°C. Similar results were obtained with cation exchange capacities when samples were treated with N CaCl2 solution. The 21.0 A peak intensity curve showed close similarity to the liquid limit and plastic index curves in the low temperature range, and an explanation is suggested.

Learning form Nature to improve the heat generation of iron-oxide nanoparticles for magnetic hyperthermia applications.

The performance of magnetic nanoparticles is intimately entwined with their structure, mean size and magnetic anisotropy. Besides, ensembles offer a unique way of engineering the magnetic response by modifying the strength of the dipolar interactions between particles. Here we report on an experimental and theoretical analysis of magnetic hyperthermia, a rapidly developing technique in medical research and oncology. Experimentally, we demonstrate that single-domain cubic iron oxide particles resembling bacterial magnetosomes have superior magnetic heating efficiency compared to spherical particles of similar sizes. Monte Carlo simulations at the atomic level corroborate the larger anisotropy of the cubic particles in comparison with the spherical ones, thus evidencing the beneficial role of surface anisotropy in the improved heating power. Moreover we establish a quantitative link between the particle assembling, the interactions and the heating properties. This knowledge opens new perspectives for improved hyperthermia, an alternative to conventional cancer therapies.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.