Major element, trace element, and Sr, Nd and Pb isotope studies of Cenozoic basalts from the South China Sea

The whole rock K-Ar ages of basalts from the South China Sea basin vary from 3.8 to 7.9 Ma, which suggest that intra-plate volcanism after the cessation of spreading of the South China Sea (SCS) is comparable to that in adjacent regions around the SCS, i.e., Leiqiong Peninsula, northern margin of the SCS, Indochina block, and so on. Based on detailed petrographic studies, we selected many fresh basaltic rocks and measured their major element, trace element, and Sr-Nd-Pb isotope compositions. Geochemical characteristics of major element and trace element show that these basaltic rocks belong to alkali basalt magma series, and are similar to OIB-type basalt. The extent of partial melting of mantle rock in source region is very low, and magma may experience crystallization differentiation and cumulation during the ascent to or storing in the high-level magma chamber. Sr-Nd-Pb isotopic data of these basaltic rocks imply an inhomogeneous mantle below the South China Sea. The nature of magma origin has a two end-member mixing model, one is EM2 (Enriched Mantle 2) which may be originated from mantle plume, the other is DMM (Depleted MORB Mantle). Pb isotopic characteristics show the Dupal anomaly in the South China Sea, and combined with newly found Dupal anomaly at Gakkel ridge in Arctic Ocean, this implies that Dupal anomaly is not only limited to South Hemisphere. In variation diagrams among Sr, Nd and Pb, the origin nature of mantle below the SCS is similar to those below Leiqiong peninsula, northern margin of the SCS and Indochina peninsula, and is different from those below north and northeast China. This study provides geochemical constraints on Hainan mantle plume.

Genesis of Th-230 excess in basalts from mid-ocean ridges and ocean islands: Constraints from the global U-series isotope database and major and rare earth element geochemistry

Based on Th-230-U-238 disequilibrium and major element data from mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs), this study calculates mantle melting parameters, and thereby investigates the origin of Th-230 excess. (Th-230/U-238) in global MORBs shows a positive correlation with Fe-8, P (o), Na-8, and F-melt (Fe-8 and Na-8 are FeO and Na2O contents respectively after correction for crustal fractionation relative to MgO = 8 wt%, P (o)=pressure of initial melting and F (melt)=degree of melt), while Th-230 excess in OIBs has no obvious correlation with either initial mantle melting depth or the average degree of mantle melting. Furthermore, compared with the MORBs, higher (Th-230/U-238) in OIBs actually corresponds to a lower melting degree. This suggests that the Th-230 excess in MORBs is controlled by mantle melting conditions, while the Th-230 excess in OIBs is more likely related to the deep garnet control. The vast majority of calculated initial melting pressures of MORBs with excess Th-230 are between 1.0 and 2.5 GPa, which is consistent with the conclusion from experiments in recent years that D (U)> D (Th) for Al-clinopyroxene at pressures of > 1.0 GPa. The initial melting pressure of OIBs is 2.2-3.5 GPa (around the spinel-garnet transition zone), with their low excess Ra-226 compared to MORBs also suggesting a deeper mantle source. Accordingly, excess Th-230 in MORBs and OIBs may be formed respectively in the spinel and garnet stability field. In addition, there is no obvious correlation of K2O/TiO2 with (Th-230/U-238) and initial melting pressure (P (o)) of MORBs, so it is proposed that the melting depth producing excess Th-230 does not tap the spinel-garnet transition zone. OIBs and MORBs in both (Th-230/U-238) vs. K2O/TiO2 and (Th-230/U-238) vs. P (o) plots fall in two distinct areas, indicating that the mineral phases which dominate their excess Th-230 are different. Ce/Yb-Ce curves of fast and slow ridge MORBs are similar, while, in comparison, the Ce/Yb-Ce curve for OIBs shows more influence from garnet. The mechanisms generating excess Th-230 in MORBs and OIBs are significantly different, with formation of excess Th-230 in the garnet zone only being suitable for OIBs.

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.

Winkler springs (p-y curves) for liquefied soil from element tests.

In practice, piles are most often modelled as "Beams on Non-Linear Winkler Foundation" (also known as “p-y spring” approach) where the soil is idealised as p-y springs. These p-y springs are obtained through semi-empirical approach using element test results of the soil. For liquefied soil, a reduction factor (often termed as p-multiplier approach) is applied on a standard p-y curve for the non-liquefied condition to obtain the p-y curve liquefied soil condition. This paper presents a methodology to obtain p-y curves for liquefied soil based on element testing of liquefied soil considering physically plausible mechanisms. Validation of the proposed p-y curves is carried out through the back analysis of physical model tests.

Applications of finite element computer modelling and thermal infrared remote sensing to the study of geothermal anomalies

Buried heat sources can be investigated by examining thermal infrared images and comparing these with the results of theoretical models which predict the thermal anomaly a given heat source may generate. Key factors influencing surface temperature include the geometry and temperature of the heat source, the surface meteorological environment, and the thermal conductivity and anisotropy of the rock. In general, a geothermal heat flux of greater than 2% of solar insolation is required to produce a detectable thermal anomaly in a thermal infrared image. A heat source of, for example, 2-300K greater than the average surface temperature must be a t depth shallower than 50m for the detection of the anomaly in a thermal infrared image, for typical terrestrial conditions. Atmospheric factors are of critical importance. While the mean atmospheric temperature has little significance, the convection is a dominant factor, and can act to swamp the thermal signature entirely. Given a steady state heat source that produces a detectable thermal anomaly, it is possible to loosely constrain the physical properties of the heat source and surrounding rock, using the surface thermal anomaly as a basis. The success of this technique is highly dependent on the degree to which the physical properties of the host rock are known. Important parameters include the surface thermal properties and thermal conductivity of the rock. Modelling of transient thermal situations was carried out, to assess the effect of time dependant thermal fluxes. One-dimensional finite element models can be readily and accurately applied to the investigation of diurnal heat flow, as with thermal inertia models. Diurnal thermal models of environments on Earth, the Moon and Mars were carried out using finite elements and found to be consistent with published measurements. The heat flow from an injection of hot lava into a near surface lava tube was considered. While this approach was useful for study, and long term monitoring in inhospitable areas, it was found to have little hazard warning utility, as the time taken for the thermal energy to propagate to the surface in dry rock (several months) in very long. The resolution of the thermal infrared imaging system is an important factor. Presently available satellite based systems such as Landsat (resolution of 120m) are inadequate for detailed study of geothermal anomalies. Airborne systems, such as TIMS (variable resolution of 3-6m) are much more useful for discriminating small buried heat sources. Planned improvements in the resolution of satellite based systems will broaden the potential for application of the techniques developed in this thesis. It is important to note, however, that adequate spatial resolution is a necessary but not sufficient condition for successful application of these techniques.

Uniformly convergent finite element and finite difference methods for singularly perturbed ordinary differential equations

This thesis is concerned with uniformly convergent finite element and finite difference methods for numerically solving singularly perturbed two-point boundary value problems. We examine the following four problems: (i) high order problem of reaction-diffusion type; (ii) high order problem of convection-diffusion type; (iii) second order interior turning point problem; (iv) semilinear reaction-diffusion problem. Firstly, we consider high order problems of reaction-diffusion type and convection-diffusion type. Under suitable hypotheses, the coercivity of the associated bilinear forms is proved and representation results for the solutions of such problems are given. It is shown that, on an equidistant mesh, polynomial schemes cannot achieve a high order of convergence which is uniform in the perturbation parameter. Piecewise polynomial Galerkin finite element methods are then constructed on a Shishkin mesh. High order convergence results, which are uniform in the perturbation parameter, are obtained in various norms. Secondly, we investigate linear second order problems with interior turning points. Piecewise linear Galerkin finite element methods are generated on various piecewise equidistant meshes designed for such problems. These methods are shown to be convergent, uniformly in the singular perturbation parameter, in a weighted energy norm and the usual L2 norm. Finally, we deal with a semilinear reaction-diffusion problem. Asymptotic properties of solutions to this problem are discussed and analysed. Two simple finite difference schemes on Shishkin meshes are applied to the problem. They are proved to be uniformly convergent of second order and fourth order respectively. Existence and uniqueness of a solution to both schemes are investigated. Numerical results for the above methods are presented.