On the strength of electric currents and zonal magnetic fields at the top of the Earth''s core: Methodology and preliminary estimates

A new method is described for estimating: (a) the meridional electric current density, j_θ, (b) the vertical growth rate of the zonal magnetic field, ∂B_φ/∂r, or its scale-height, B_φ/∂B_φ/∂r) and (c) the vertical growth rate of the vertical current density, ∂j_r/∂r, at a few isolated points on the top surface of the Earth's core from observations of the internal geomagnetic field at the Earth's surface. The theoretical technique rests on combining unaccelerated, gravitationally-driven Boussinesq fluid dynamics of the core with frozen-flux electromagnetism, the mantle being treated as a spherically symmetric insulator.

Insertion into this theory of main field models for epochs 1965, 1975 leads to preliminary values for these quantities of magnitude: (a) j_θ 1 A/m², (b) ∂B_φ/∂r 10⁻⁶ T/m or B_φ/(∂B_φ/∂r) 10 m, (c) ∂j_r/∂r 10⁻⁶ A/m³. Some geophysical implications of these estimates are discussed.     

Response spectra for differential motion of columns paper II: Out-of-plane response

It is shown that the common response spectrum method for synchronous ground motion can be extended to make it applicable for earthquake response analyses of extended structures experiencing differential out-of-plane ground motion. A relative displacement spectrum for design of first-story columns SDC (T, T_T, ζ, ζ_T, τ, δ) is defined. In addition to the natural period of the out-of-plane response, T, and the corresponding fraction of critical damping, ζ, this spectrum also depends on the fundamental period of torsional vibrations, T_T, and the corresponding fraction of critical damping, ζ_T, on the “travel time,” τ (of the waves in the soil over a distance of about one-half the length of the structure), and on a dimensionless factor δ, describing the relative response of the first floor. The new spectrum, SDC, can be estimated by using the empirical scaling equations for relative displacement spectra, SD, and for peak ground velocity, v_max. For recorded strong-motion acceleration, and for symmetric buildings, the new spectrum can be computed from Duhamel's integrals of two uncoupled equations for dynamics equilibrium describing translation and rotation of a two-degree-of-freedom system. This representation is accurate when the energy of the strong-motion is carried by waves in the ground the wavelengths of which are one order of magnitude or more longer than the characteristic length of the structure.     

Imposex in Morula Granulata as Bioindicator of Tributyltin (TBT) Contamination in the Dampier Archipelago, Western Australia

This study was undertaken to examine the use of imposex in the tropical neogastropod whelk Morula granulata as bioindicator of tributyltin (TBT) contamination. Samples were collected from sites throughout the Dampier Archipelago; both impact sites close to intensive shipping activity and remote control sites. TBT in this area originates mainly from commercial vessels using the Port of Dampier that are >25 m in length and exempt from legislation controlling the use of TBT in antifouling paints. A field survey was undertaken in July 1997; 100 M. granulata were taken at each of 18 sites with varying vessel activity, and examined for imposex. The level of TBT contamination in the Archipelago was low; ranging from <0.3 to 25 ng Sn l⁻¹ in water,<0.3 to 80 ng Sn g−¹ in oyster tissue (Saccostrea cuccullata – prey of M. granulata), and <0.3 to 33 ng Sn g⁻¹ in whelks. Percentage of imposex in M. granulata ranged from 0% to 57%, and was correlated with distance to vessel activity (r=0.646, p<0.05). The only significant relationship of imposex with butyltin concentrations was with TBT in oysters (r=0.826, p<0.05), and dibutyltin (DBT) in whelks (r=0.783, p<0.05). It is postulated that this is because of interspecies variation in ability to metabolize TBT to DBT and MBT (monobutyltin), as well as induction of imposex in earlier life stages, so that imposex in adults could be a reflection of prior contamination. M. granulata is a suitable bioindicator of TBT contamination in the Dampier Archipelago, which could be used throughout the Indo-Pacific given its abundance in the region.     

Steady state analysis of unsaturated flow above a shallow water-table aquifer drained by ditches

The one-dimensional pressure head profile above a fixed water-table was studied for different steady infiltration rates. As shown in previous studies, when the infiltration rate (q_in) is less than the soil saturated hydraulic conductivity (Ks), this profile can be divided into two parts: (1) from the water-table surface (z₀) to an elevation z_γ, the pressure head varies from 0 to a value h_γ such as K(h_γ)=−q_in; (2) above the elevation z_γ the pressure head is constant and equal to h_γ. Above the water-table the zone where the pressure head is variable has been called ‘transition zone’. Its height is shown to be highly variable and to depend on soil properties as well as on the infiltration rate. This transition zone is not the ‘capillary fringe’ as defined by Gillham (Gillham R.W., 1984. The capillary fringe and its effect on water-table response. J. Hydrol. 67, 307–324). Numerical experiments performed with HYDRUS-2D^® for the case of a water-table drained by parallel ditches have shown that the height of the transition zone is similar in the one-dimensional profile and in the two-dimensional system as long as the local slope of the water-table is small. This result is important since in a two-dimensional system, the transition zone is the place where all the horizontal component of the unsaturated flow occurs. The ratio of the horizontal component of the unsaturated flow vs the total horizontal flow in both the unsaturated and saturated zones has been computed. For a given soil, this ratio decreases as the infiltration rate increases; for a given infiltration rate, the soil with the thinnest transition zone transfers the largest amount of water above the water-table.     

Histopathological effects of petroleum hydrocarbons and heavy metals on the American oyster (Crassostrea virginica) from Tabasco, Mexico

Oyster tissues may be affected by the concentration of ions in the water (i.e. salinity) and by contaminants such as petroleum hydrocarbons and heavy metals. Oyster populations from three coastal lagoons (Mecoacan, Camen and Machona) in the Mexican state of Tabasco, in the southern Gulf of Mexico, were sampled for pollution studies during June, September and November 1992 and May 1993. No statistically significant relationships were found between the concentration of contaminants in the body tissues of organisms and their shell length, soft tissue weight and particulate matter; however, a significant correlation was found between condition index and salinity (r = −0.72). Generally, the lowest weights, shell lengths and indices of condition were found in Mecoacan. The concentrations of cadmium and zinc were inversely related to salinity (r = −0.52 and r = −0.32, respectively), so a riverine input is suspected. On average, 63% of the individual oysters showed histopathological lesions, which can be related to salinity and to the concentrations of cadmium and the unresolved fraction of hydrocarbons (UCM). The response of each of the tissues analysed was different. The percentage of individuals with damage in the digestive diverticulum increased linearly with UCM (r = 0.71), but in a saturation-response fashion (Y = B₀ − B₀/(1 + X^B1) r = 0.66) with cadmium. The percentage of individuals with damage in the gills increased linearly with weight (r = 0.68), cadmium (r = 0.60) and UCM (r = 0.60). The lesions in the connective tissue decreased linearly (r = −0.82) with salinity, but increased in a saturation-response way (r = 0.83) with cadmium. Finally, the percentage of individuals with lesions in the digestive tube decreased linearly with salinity (r = −0.59). Only the damage to the gills and digestive diverticula were dependent on gonadal maturity, while damage to the connective tissue was dependent on the sex of the individual.     

The ground-state conductivity profile of the core-mantle system of the Earth

The ground state of the core-mantle conductivity system is defined as a step function transition from the constant non-zero conductivity of the core to zero conductivity in the mantle. This ground state is reached by letting the thickness of a transition region Δ → 0, where Δ appears as a parameter in the model conductivity of the system. By a transformation of the equation governing the behavior of the electrostatic potential φ, it is shown that the function V(x) = (1/σ^1/2)(d²σ^1/2/dx²) acts as a potential barrier in the quantum mechanical sense, and that for certain conductivity profiles σ(x), where x is the usual Cartesian coordinate, the electrostatic potential is screened out in regions where σ(x) → 0 as Δ → 0. Consequently, E = −φ also vanishes in these regions. The results generalize to the time-dependent case. Conditions that the conductivity function σ(x) must satisfy to qualify as the ground-state conductivity are defined and an example is provided.     

Assessment of a Mediterranean shallow lentic ecosystem (Lake Pamvotis, Greece) using benthic community diversity: Response to environmental parameters

Macroinvertebrates play a key role in freshwater lentic and lotic ecosystems. The macroinvertebrate benthic community of a shallow Mediterranean lake (Lake Pamvotis, NW Greece) was studied. The benthic assemblage was sampled monthly at five sites during a period of 1 year (Apr. 1998–Mar. 1999). In addition hypolimnetic water quality variables were monitored over the same period at each site.

The aim of the study was (a) to describe the intra-annual and spatial variability in benthic communities, (b) to relate possible community changes to environmental conditions and (c) to evaluate the responses of the lake's ecological status on community indices.

The benthic fauna of Lake Pamvotis was found to be very limited with a total of 10 species belonging to five taxonomic groups. The oligochaete community comprised 80% of the total benthic fauna with Potamothrix bavaricus as a new record for the Lake Pamvotis and Potamothrix hammoniensis, being the dominant benthic species represented more than 61% of the total benthic fauna. Chironomus plumosus was the most abundant chironomid species contributing with about 6% of the total benthic fauna, and Chaoborus flavicans with 19% was the important dipteran. Almost all benthic species showed the same intra-annual seasonal pattern, with peak population densities during spring and early summer except P. hammoniensis which predominated during the whole sampling period. Dissolved oxygen and temperature seemed to be the main environmental factors affecting community indices.

Benthic communities are affected by human disturbances in Lake Pamvotis shifting their composition to more tolerant taxa, reflecting also the eutrophic to hypertophic character of the lake.     

Mantle and crustal contribution in the genesis of Recent basalts from off-rift zones in Iceland: Constraints from Th, Sr and O isotopes

Along the two volcanic off-rift zones in Iceland, the Snfellsnes volcanic zone (SNVZ) and the South Iceland volcanic zone (SIVZ), geochemical parameters vary regularly along the strike towards the centre of the island. Recent basalts from the SNVZ change from alkali basalts to tholeiites where the volcanic zone reaches the active rift axis, and their⁸⁷Sr/⁸⁶Sr andTh/U ratios decrease in the same direction. These variations are interpreted as the result of mixing between mantle melts from two distinct reservoirs below Snfellsnes. The mantle melt would be more depleted in incompatible elements, but with^a higher³He/⁴He ratio (R/Ra≈ 20) beneath the centre of Iceland than at the tip of the Snfellsnes volcanic zone (R/Ra≈ 7.5).

From southwest to northeast along the SIVZ, the basalts change from alkali basalts to FeTi basalts and quartz-normative tholeiites. TheTh/U ratio of the Recent basalts increases and both (²³⁰Th/²³²Th) andδ¹⁸O values decrease in the same direction. This reflects an important crustal contamination of the FeTi-rich basalts and the quartz tholeiites. The two types of basalts could be produced through assimilation and fractional crystallization in which primary alkali basaltic and olivine tholeiitic melts ‘erode’ and assimilate the base of the crust. The increasingly tholeiitic character of the basalts towards the centre of Iceland, which reflects a higher degree of partial melting, is qualitatively consistent with increasing geothermal gradient and negative gravity anomaly.

The highest Sr isotope ratio in Recent basalts from Iceland is observed inÖrfajökull volcano, which has a³He/⁴He ratio (R/Ra≈ 7.8) close to the MORB value, and this might represent a mantle source similar to that of Mauna Loa in Hawaii.     

The metamorphic history of LL-group ordinary chondrites

A study of feldspar grain size and microprobe analysis of coexisting orthopyroxene, olivine, diopside, and metal in twenty equilibrated LL-group chondrites and clasts defines the conditions of metamorphism, indicates the relative position of samples within the parent body, and suggests an accretionary-metamorphic history.

Variation of mineral chemistry indicates a range of equilibrium temperature and parallels textural variation, providing two continuous relative temperature scales, based on mean feldspar grain size and the calcium content of orthopyroxene coexisting with diopside. Maximum equilibrium temperatures range from less than 600° up to 950°C. The total Tschermak's component in orthopyroxene indicates that maximum pressure was less than 1 kbar, and a pressure gradient correlates directly with temperature. A reduction sequence, evident from systematic variation of iron in ferromagnesian silicates and metal, correlates inversely with temperature.

The maximum permissible sampling depth for the formation of surface breccias (2/10 r) also limits pressure to less than 1 kbar. The LL parent body had a minimum radius of 200 km, determined using the permissible sampling depth and published cooling rates.

The data are broadly consistent with simple reheating of colder material within a completely accreted body. This model may have problems with volatile depletion, a required lack of impact sampling during metamorphism, and the timing of the event. The data are more consistent with accretion of colder material on a hot but cooling surface. Both models apparently require a recycling process during accretion to account for the reduction sequence.     

内蒙古中西部地区震源参数研究

选取内蒙古测震台网2009—2016年内蒙古中西部M_L≥2.8中小地震波形记录,采用多台联合反演方法,计算该区中小地震震源谱参数,获得地震矩M₀、矩震级M_W、应力降Δσ和震源尺度r,利用线性回归,分析近震震级M_L、地震矩M₀、矩震级M_W、应力降Δσ和震源尺度r之间的关系。结果表明,各参数之间存在一定线性关系,内蒙古中西部地区应力降模型为增加应力降模型（ISD）。     

Hydrological impact of Eucalyptus plantation in India

A brief review is given of the results of hydrological studies carried out in southern India on the effects of plantations of Eucalyptus and other fast growing exotic tree species on water resources, erosion and soil nutrients at sites of differing rainfall and soil depth in Karnataka. New results on the impacts of the plantations on raindrop-induced erosion and preliminary results from studies that are aimed at improving the water use efficiency of the plantations are also presented.

The erosion studies indicate that soil detachment resulting from net rainfall beneath Eucalyptus camaldulensis will be greater than beneath Pinus caribaea but less than that beneath the indigenous species, Tectona grandis.

The water use efficiency studies, which make use of a controlled environment facility, have confirmed that in the dry zone of southern India water availability is the principal limiting factor on growth although, when water is available, nutrient limitations become important. Removal of both water and nutrient stress results typically in a fivefold increase in volume growth for the first year of growth.     

Methyl mercury, total mercury, and selenium in snapper from two areas of the New South Wales coast, Australia

The muscle tissue of snapper (Chrysophrys auratus: Sparidae) from two areas of the New South Wales coast was analysed for methyl mercury, total mercury and selenium concentrations (134 samples from the Sydney area and 186 samples from the Nowra area, about 120 km to the south of Sydney).

The Sydney snapper ranged from 0.08 to 1.66 mg kg⁻¹ total mercury (mean 0.32 mg kg⁻¹), while the Nowra snapper ranged from 0.01 to 0.78 mg kg⁻¹ total Hg (mean 0.11 mg kg⁻¹). These latter results appear anomalously low relative to other studies in Australia and New Zealand.

Methyl mercury averaged 94.7% of the total mercury content of the Sydney snapper, and 91.3% for the Nowra snapper.

The concentrations of selenium in the muscle tissues of snapper from both areas were similar and did not increase with increasing mercury concentrations. The highest selenium concentration was 0.85 mg kg⁻¹.     

A metamorphic history from micron-scale Pb/Pb chronometry of Archean monazite

Ion microprobe measurements of Pb isotope ratios in monazites have been obtained, in situ, from thin sections using the Cambridge ISOLAB 120. Molecular interferences are sufficiently resolved at an RP of 6500 to allow ²⁰⁷Pb/²⁰⁶Pb dating of monazite with precisions as low as 4–5 Ma (2σ). The results presented here provide important information on the chronological history of the Late Archean metamorphism of the Wind River Range, Wyoming (USA).

Matrix monazites and monazite inclusions in garnets from a metapelite from the northern Wind River Range have been analysed by SIMS. In a previous study peak metamorphic conditions (T = 800°C; P = 8 ± 1 kb^*) were estimated using inclusion assemblages in garnets from this same sample. Isolated monazite inclusions in garnet yield ²⁰⁷Pb/²⁰⁶Pb age estimates of 2781 ± 6 to 2809 ± 10 Ma. Those along fractures yield lower ages (2603–2687 Ma) which are similar to TIMS and SIMS ages of matrix monazites. A single large (500 μm) monazite grain locally preserves growth zoning, but has a recrystallised core and a resorbed (recrystallised?) rim. Age estimates for these three regions are 2788 ± 9 Ma, 2663 ± 4 and 2523 ± 6 Ma, respectively. Thus the inclusion assemblages of Sharp and Essene^* may record peak metamorphic conditions at ca. 2.8 Ga, and indicate a phase of metamorphism that predates by over 100 Ma the emplacement of the Bridger Batholith, the major lithologic component of the northern Wind River Range.

The analysed monazite grains appear to preserve ca. 300 Ma history, even within a single grain. Monazite inclusions in garnet that are fully armoured may provide estimates for the time of garnet growth, even in high grade terranes where most chronometers are reset. The age pattern preserved by the large monazite grain cannot be simply related to diffusion controlled closure. Instead, a chronology is preserved which can be related to the petrographic setting of indicidual grains through in situ analysis.     

Transitions in thermal convection with strongly temperature-dependent viscosity in a wide box

Numerical models are systematically presented for time-dependent thermal convection of Newtonian fluid with strongly temperature-dependent viscosity in a two-dimensional rectangular box of aspect ratio 3 at various values of the Rayleigh number Ra_b defined with viscosity at the bottom boundary up to 1.6×10⁸ and the viscosity contrast across the box r_η up to 10⁸. We found that there are two different series of bifurcations that take place as r_η increases. One series of bifurcations causes changes in the behavior of the thermal boundary layer along the surface boundary from small-viscosity-contrast (SVC) mode, through transitional (TR) mode, to stagnant-lid (ST) mode, or from SVC mode directly to ST mode, depending on Ra_b. Another series of bifurcations causes changes in the aspect ratio of convection cells; convection with an elongated cell can take place at moderate r_η (10³–10^5.5 at Ra_b=6×10⁶), while only convection of aspect ratio close to 1 takes place at small r_η and large r_η. The parameter range of r_η and Ra_b for elongated-cell convection overlaps the parameter range for SVC and ST modes and include the entire parameter range for TR mode. In the elongated-ST regime, the lid of highly viscous fluid along the top boundary is not literally ‘stagnant’ but can horizontally move at a velocity high enough to induce a convection cell with aspect ratio much larger than 1.     

Anelastic structure of the upper mantle beneath the northern Philippine Sea

We investigated the upper mantle anelastic structure beneath the northern Philippine Sea region, including the Izu-Bonin subduction zone and the Shikoku Basin. We used regional waveform data from 69 events in the Pacific and the Philippine Sea slabs, recorded on F-net and J-array network broadband stations in western Japan. Using the S–P phase pair method, we obtained differential attenuation factors, δt^*, which represent the relative whole path Q. We conducted a tomographic inversion using 978 δt^* values to invert for a fine-scale (50–100 km) three-dimensional anelastic structure.

The results shows two high-Q regions (Q_P>1000) which are consistent with the locations of the Pacific and the Philippine Sea slabs. Also there is a low-Q (Q_P110) area extending to the deeper parts (350–400 km) of the model just beneath the old spreading center and the Kinan Seamount Chain in the Shikoku Basin. A small depth dependence of the laterally averaged Q_P was found, with values of 266 (0–250 km), 301 (250–400 km), and 413 (400–500 km).     

On the use of microtremor recordings in seismic microzonation

Experimental methods involving microtremor recordings are useful for determining site effects in regions of moderate seismic activity where ground motion records are few, and in urban or industrial contexts where the noise level is high. The aim of this study is to establish a microzonation by using the Nogoshi–Nakamura method,[1, 2] a simple experimental technique based on microtremor recordings. Since the physical phenomena underlying the method are only partially understood, the spectral responses obtained cannot be used alone. We, therefore, complete our experimental results by comparing them with the solutions of a one-dimensional numerical simulation (SHAKE91).[3, 4] The experimental programme was carried out on a plain near the Rhone Delta (South of France). H/V spectral ratios were calculated at 137 noise measurement points. In addition, we were able to compute the numerical transfer functions from soil columns defined by geotechnical characteristics inside the studied region. A comparison of the results obtained by the experimental and numerical methods showed that the fundamental frequencies are in good agreement, but that the amplitudes obtained by the two techniques are sometimes different. The analysis of H/V spectral ratios enabled us to establish maps to characterize the region: a resonance frequency map and maps of amplification levels as a function of frequency range, leading to a seismic microzonation for the whole of the region.     

基于槽波的TVSP超前探测方法:可行性研究

针对隧道工程中TVSP（Tunnel Vertical Seismic Profiling）方法不适合于煤矿巷道超前探测的问题,介绍了煤矿巷道中的槽波TVSP超前探测方法原理和实现过程.通过理论和数值模拟证明：（1）当采用垂向集中力源（z方向）激发地震波时,地震波场z分量受到巷道影响较小,且主要为Rayleigh型槽波,较为适合槽波超前探测;（2）采用传播相对稳定且易于识别的槽波埃里相峰值时间代替传统TVSP方法中体波的初至走时时间,能够较准确地确定断层在掌子面前方位置和断层角度;（3）当断层与巷道掘进方向夹角较大时,断层交点位置（断层线与巷道掘进方向的交点）具有很好的定位精度,但角度误差则会较大;当实际断层交点位置距离迎头较远时,断层交点位置定位精度会下降,但角度误差则影响较小;（4）基于槽波的TVSP方法对观测系统要求不高,数据处理过程简单,计算量小.     

Some magnetic properties of synthetic and natural haematite of different grain size

The behaviour of some magnetic properties of natural and synthetic haematite of different grain size is examined. The natural haematite was obtained from the hydrothermal deposit Kada (Czech and Slovak Federal Rep.). Six grain-size fractions ranging from 120 to 40 μm were prepared by means of sieving and two further fractions down to 5 μm by wet ultrasonic sieving. Since the behaviour of the fractions is similar, that of only four representative samples is reported. In addition, the behaviour of one submicron synthetic haematite fraction (0.5 μm) prepared by oxidation of ferrous sulphate (uniform in size and shape) was investigated.

The initial remanence value (J_r) seems to increase with decreasing grain size. During alternating field (AF) demagnetization, all fractions behaved similarly, except for the submicron fraction which is considerably softer than the others. Normalized (isothermal remanent magnetization) IRM acquisition curves were similar for all fractions.

Parameters of the anisotropy of magnetic susceptibility (AMS) display significant changes, mainly during IRM acquisition. During AF demagnetization, the anisotropy degree P exhibits a slight increase (some %), while the behaviour of the shape factor T is complicated. The anisotropy ellipsoid exhibits a tendency to rotate. Significant changes in the AMS parameters occur during IRM acquisition. Curves of P and T vs. IRM acquisition field, for various grain-size fractions, show no coherent pattern. For all the samples studied, the T vs. H curve exhibits a threshold value at which change in the type of arrangement of easy axis of magnetization occurs. For the IRM acquisition fields higher than some 320 kA m⁻¹, the minimum susceptibility axis parallels the direction of the IRM acquisition field.

Hysteresis curves of the fractions are similar to each other. The Preisach distribution function was determined and it indicates that the reversible part of the magnetization process plays an important role comparatively. Based on the coercivity data presented no unambiguous conclusion could be drawn from the single-domain (SD)-multidomain (MD) transition, associated with a coercivity maximum.     

Preliminary stratigraphic and structural architecture of Bhutan: Implications for the along strike architecture of the Himalayan system

High-precision in-situ ion microprobe (SIMS) oxygen isotope analysis of zircons from two diorite intrusions associated with the late Caledonian Lochnagar pluton in Scotland has revealed large differences in the degree of heterogeneity in zircon δ¹⁸O between the diorites. Zircon crystals from the Cul nan Gad diorite (CnG) show a unimodal distribution of oxygen isotope values (δ¹⁸O = 6.0 ± 0.6‰ (2σ)) and no or only minor grain-scale variation. Those from the Allt Darrarie diorite (AD1) show a large range in δ¹⁸O and an apparent bimodal distribution with modes of 6.6 ± 0.4‰ and 7.3 ± 0.4‰. Variations of up to 1.2‰ occur between and within grains; both an increase and decrease in δ¹⁸O with zircon growth has been observed. The δ¹⁸O composition of growing zircon can only change if open-system processes affect the magma composition, i.e. if material of contrasting δ¹⁸O composition is added to the magma. The variability in AD1 is interpreted to represent a cryptic record of magma mixing. A ‘deep crustal hot zone’ is a likely site for generation of the dioritic magmas which developed by mixing of residual melts and crustal partial melts or by melting of mafic lower crustal rocks. The overall small number of zircons with mantle-like δ¹⁸O values (5.3 ± 0.6‰ (2σ)) in the Lochnagar diorites is largely the product of crustal differentiation rather than crustal growth.

The δ¹⁸O of quartz from the CnG and AD1 diorites shows only minor variation (CnG: 10.9 ± 0.5‰ (2σ), AD1: 11.7 ± 0.6‰ (2σ)) within single populations, with no evidence of mixing. Quartz–zircon isotopic disequilibrium is consistent with later crystallisation of quartz from late magmatic fluids, and in case of the AD1 diorite after the inferred magma mixing from a homogenised, higher δ¹⁸O melt.

High-precision SIMS oxygen isotope analysis of zircon provides a new approach to identifying and resolving previously undetected early-stage magma mixing and constraining the compositions and origins of the component magmas. A combination of zircon, quartz and whole-rock data has proven to be a powerful tool in reconstructing the petrogenetic evolution of diorite from early crystallisation to late alteration.     

外核黏性对地磁场发电机数值模型的影响

地磁场发电机过程由数学方程组和输入参数所控制.为研究发电机参数对系统的影响,本文使用MoSST模型模拟了在外核黏性ν变化下发电机模型的输出.通过使用跨越近3个量级的ν值,着重研究了各物理场及其典型尺度随黏性的变化.发现黏性变化显著影响流场;磁场随ν增加而近乎单调减小,但变化幅度不超过30%;温度扰动随ν增加而小幅（6%）单调增加.经过拟合,得到外核流速u和黏性ν的比例关系：u~ν^0.49.流速随黏性增加的现象本质上是由于黏性增加打破了Taylor-Proudman约束,使得临界Rayleigh数减小,从而在相同的驱动力下带来了流速的增加.此外,作用力平衡分析发现,随着ν的变化,系统在几种平衡模式间切换.通过与之前比例关系的研究对比,本研究支持在一定范围内,磁场与黏性关系不大的结论;但反对流速与外核黏性无关的假设.