An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the Earth's mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.     

Approximate solutions for the formation of the lithosphere

The lithosphere is interpreted as a thermal boundary layer. Approximate solutions of the boundary layer cooling problem are developed which include mantle radioactivity, partial melt in the asthenosphere, a temperature gradient in the asthenosphere, and a non-zero lithospheric thickness at the ridge crests. The cooling history of oceanic lithosphere is found to be remarkably insensitive to assumptions about the amount of radioactivity in the upper mantle and the extent of melting in the asthenosphere. Determinations of the thickness of oceanic lithosphere and the depths of oceans as a function of age are in excellent agreement with boundary layer predictions which include a heat flux from the asthenosphere. However, the determinations do not resolve how much of the total asthenospheric heat flux might be caused by a temperature gradient in the asthenosphere. Simple thermal arguments indicate that the initial lithospheric thickness, L₀, at ridge crests should depend on the local half-spreading rate, V, as L₀ = 3 km/V(cm/year).     

Composition,temperature, and thickness of the lithosphere of the Archean Kaapvaal craton

A new model is proposed for the structure of the Kaapvaal craton lithosphere. Based on chemical thermodynamics methods, profiles of the chemical composition, temperature, density, and S wave velocities are constructed for depths of 100–300 km. A solid-state zone of lower velocities is discovered on the S velocity profile in the depth interval 150–260 km. The temperature profiles are obtained from absolute values of P and S velocities, taking into account phase transformations, anharmonicity, and anelastic effects. The examination of the sensitivity of seismic models to the chemical composition showed that relatively small variations in the composition of South African xenoliths result in lateral temperature variations of ～200°C. Inversion of some seismic profiles (including IASP91) with a fixed bulk composition of garnet peridotites (the primitive mantle material) leads to a temperature inversion at depths of 200–250 km, which is physically meaningless. It is supposed that the temperature inversion can be removed by gradual fertilization of the mantle with depth. In this case, the craton lithosphere should be stratified in chemical composition. The depleted lithosphere composed by garnet peridotites exists to depths of 175–200 km. The lithospheric material at depths of 200–250 km is enriched in basaltoid components (FeO, Al₂O₃, and CaO) as compared with the material of garnet peridotites but is depleted in the same components as compared with the fertile substance of the underlying primitive mantle. The material composing the craton root at a depth of ～275 km does not differ in its physical and chemical characteristics from the composition of the normal mantle, and this allows one to estimate the thickness of the lithosphere at 275 km. The results of this work are compared with data of seismology, thermal investigations, and thermobarometry.     

Siberian traps: Hypotheses and seismology data

Siberian traps are the result of huge basalt eruptions which took place about 250 Ma ago over a vast territory of Siberia. The genesis of Siberian traps is attributed to a mantle plume with a center in the region of Iceland or beneath the central Urals in terms of their present coordinates. The eruption mechanism is associated with delamination—replacement of the mantle lithosphere by the deep magma material. The receiver function analysis of the records from the Norilsk seismic station (NRIL) allows comparing these hypotheses with the factual data on the depth structure of the region of Siberian traps. The S-wave velocity section place the seismic lithosphere/asthenosphere boundary (LAB) at a depth of 155–190 km, commensurate with the data for the other cratons. The mantle lithosphere has a high S-wave velocity characteristic of cratons (4.6–4.8 km/s instead of the typical value 4.5 km/s). The seismic boundary, which is located at a depth around 410 km beneath the continents is depressed by ~10 km in the region of the NRIL station. The phase diagram of olivine/wadsleyite transformation accounts for this depression by a 50–100°С increase in temperature. At the depths of 350–400 km, the S-wave velocity drops due to partial melting. A new reduction in the S-wave velocities is observed at a depth of 460 km. The similar anomalies (deepening of the 410-km seismic boundary and low shear wave velocity at depths of 350–400 and 460–500 km, respectively) were previously revealed in the other regions of the Meso-Cenozoic volcanism. In the case of a differently directed drift of the Siberian lithosphere and underlying mantle at depths down to 500 km, these anomalies are barely accountable. In particular, if the mantle at a depth ranging from 200 to 500 km is fixed, the anomalies should be observed at the original locations where they emerged 250 Ma ago, i.e. thousands of km from the Siberian traps. Our seismic data suggest that despite the low viscosity of the asthenosphere, the mantle drift at depths ranging from 200 to 500 km is correlated with the drift of the Siberian lithospheric plate. Furthermore, the position of the mantle plume beneath the Urals is easier to reconcile with the seismic data than its position beneath Iceland because of the Siberian traps being less remote from the Urals.     

Stability of the boundary layer between the lithosphere and convecting mantle and the steady-state lithospheric geotherm

In the steady state, the convective boundary layer (CBL) (the transition from the lithosphere to the convecting mantle, the lithosphere-asthenosphere boundary) is on the verge of stability. This determines its depth, thickness, and the steady-state temperature distribution in the lithosphere. Had the mantle been homogeneous, the base of the lithosphere at the current potential temperature would lie globally at the same depth H _rh of 50 to 70 km. Actually, the regime of interaction of the mantle convection with the lithosphere is determined by the relationship between this depth and the thickness H _depl of the chemical boundary layer including the crust and the layer of the depleted rock. If the thickness of the chemical boundary layer is small H _depl < H _rh, as it is the case in the present-day oceanic mantle, the suboceanic regime is established with the mantle convection that does not reach the base of the chemical boundary layer. In this case, the top of CBL is located at depth H _rh, while the oceanic heat flow and the depth of the seafloor only depend on the potential temperature T _p and, within the areas where the crust is older than 60 to 70 Ma, are the same everywhere far from the disturbed territories (the hot points and the subduction zones). The absence of noticeable distinctions between the heat flow in the different oceanic basins suggests a global constancy of the potential temperature. If H _depl > H _rh, the subcontinental regime of the interaction of the mantle convection with the lithosphere is established. In this case, the CBL is immediately adjacent to the depleted lithosphere, its top is located at depth H _depl, and the surface heat flow q(T _p, H _depl) not only depends on the potential temperature T _p but also on the the thickness of the depleted lithosphere H _depl; it decreases with increasing H _depl and, therefore, with the age of the lithosphere. Given the potential temperature, the dependence q(T _p, H _depl) agrees well with the envelope of the results of kimberlite xenolith thermobarometry presented in the diagram of the deepest xenolith depth as a function of the heat flow. It is likely that in the lowest part of the continental lithosphere there is a zone of horizontal shear deformation, from where kimberlites entrain the strongly deformed and, at the same time, the deepest xenoliths. Besides, the azimuthal anisotropy of seismic velocities can be associated with this zone. The change in its direction with depth can be observed as the Lehmann discontinuity.     

Micro pulse lidar observations of mineral dust layer in the lower troposphere over the southwest coast of Peninsular India during the Asian summer monsoon season

A large aerosol plume with optical depth exceeding 0.7 engulfs most parts of the Arabian Sea during the Asian summer monsoon season. Based on Micro Pulse Lidar observations during the June–September period of 2008 and 2009, the present study depicts, for the first time, the existence of an elevated dust layer occurring very frequently in the altitude band of 1–3.5 km over the west coast of peninsular India with relatively large values of linear depolarization ratio (δ_L). Large values of δ_L indicate the dominance of significantly non-spherical aerosols. The aerosol optical depth of this layer (0.2) is comparable to that of the entire atmospheric column during dust-free days. Back-trajectory analysis clearly shows the advection of airmass from the arid regions of Arabia and the west Arabian Sea, through the altitude region centered around 3 km. This is in contrast to the airmass below 1 km originating from the pristine Indian Ocean region which contains relatively spherical aerosols of marine origin with δ_L generally <0.05.     

Formation time of the big mantle wedge beneath eastern China and a new lithospheric thinning mechanism of the North China craton—Geodynamic effects of deep recycled carbon

High-resolution P wave tomography shows that the subducting Pacific slab is stagnant in the mantle transition zone and forms a big mantle wedge beneath eastern China. The Mg isotopic investigation of large numbers of mantle-derived volcanic rocks from eastern China has revealed that carbonates carried by the subducted slab have been recycled into the upper mantle and formed carbonated peridotite overlying the mantle transition zone, which becomes the sources of various basalts. These basalts display light Mg isotopic compositions(δ26 Mg = –0.60‰ to –0.30‰) and relatively low87 Sr/86 Sr ratios(0.70314–0.70564) with ages ranging from 106 Ma to Quaternary, suggesting that their mantle source had been hybridized by recycled magnesite with minor dolomite and their initial melting occurred at 300-360 km in depth. Therefore, the carbonate metasomatism of their mantle source should have occurred at the depth larger than 360 km, which means that the subducted slab should be stagnant in the mantle transition zone forming the big mantle wedge before 106 Ma. This timing supports the rollback model of subducting slab to form the big mantle wedge. Based on high P-T experiment results, when carbonated silicate melts produced by partial melting of carbonated peridotite was raising and reached the bottom(180–120 km in depth) of cratonic lithosphere in North China, the carbonated silicate melts should have 25–18 wt% CO2 contents, with lower Si O2 and Al2 O3 contents, and higher Ca O/Al2 O3 values, similar to those of nephelinites and basanites, and have higher εNdvalues(2 to 6). The carbonatited silicate melts migrated upward and metasomatized the overlying lithospheric mantle, resulting in carbonated peridotite in the bottom of continental lithosphere beneath eastern China. As the craton lithospheric geotherm intersects the solidus of carbonated peridotite at 130 km in depth, the carbonated peridotite in the bottom of cratonic lithosphere should be partially melted, thus its physical characters are similar to the asthenosphere and it could be easily replaced by convective mantle. The newly formed carbonated silicate melts will migrate upward and metasomatize the overlying lithospheric mantle. Similarly, such metasomatism and partial melting processes repeat, and as a result the cratonic lithosphere in North China would be thinning and the carbonated silicate partial melts will be transformed to high-Si O2 alkali basalts with lower εNdvalues(to-2). As the lithospheric thinning goes on,initial melting depth of carbonated peridotite must decrease from 130 km to close 70 km, because the craton geotherm changed to approach oceanic lithosphere geotherm along with lithospheric thinning of the North China craton. Consequently, the interaction between carbonated silicate melt and cratonic lithosphere is a possible mechanism for lithosphere thinning of the North China craton during the late Cretaceous and Cenozoic. Based on the age statistics of low δ26 Mg basalts in eastern China, the lithospheric thinning processes caused by carbonated metasomatism and partial melting in eastern China are limited in a timespan from 106 to25 Ma, but increased quickly after 25 Ma. Therefore, there are two peak times for the lithospheric thinning of the North China craton: the first peak in 135-115 Ma simultaneously with the cratonic destruction, and the second peak caused by interaction between carbonated silicate melt and lithosphere mainly after 25 Ma. The later decreased the lithospheric thickness to about70 km in the eastern part of North China craton.     

A worldwide comparison of predicted S-wave delays from a three-dimensional upper mantle model with P-wave station corrections

This paper is a comparison, on a worldwide scale, between P-station corrections deduced from ISC residuals (Dziewonski and Anderson) and synthetic S-station corrections computed using a three-dimensional upper mantle model obtained from mantle wave data (Woodhouse and Dziewonski). The upper mantle S-velocity model is described by a spherical harmonic expansion up to degree 8; the P-station corrections are smoothed using a similar expansion, in order that the two data sets can be compared.Correlations between P-station corrections, δt_P, and synthetic S-station corrections, δt_S relevant to various depths of integration indicate that a station correction contains information about structures down to at least 670 km. For this depth of integration, the correlation coefficient of the two data sets is 0.59; the slope ‘a’ of the relation δt_S = aδt_P + b, obtained for the worldwide distribution of stations, is in good agreement with results of previous regional studies using direct readings of P and S arrival times (a = 3.61 ± 0.13).An analysis of regional variations of the relation δt_S = aδt_P + b is carried out on the basis of two published global tectonic patterns (Okal; Jodan). Results for oceanic regions are not reliable, due to the lack of data. On continental areas, a significant difference appears between mountains (a = 2.7 ± 0.3) and shields (a = 4.5 ± 0.4 for Okal's pattern, a = 5.7 ± 1.5 for Jordan's pattern). The largest a-value for shields rules out an explanation by partial melting, as proposed in previous studies. Thermal heterogeneities lead to low a-values; undulations of the lithosphere-asthenosphere boundary appear to be the most feasible explanation of the high slope beneath shields; they are also able to explain the range of variation of the station corrections; the lowest values of the station corrections correspond to a total vanishing of the low velocity zone beneath the oldest shields. For mountains, the mean values of the station corrections as well as the low a-value can be accounted for by a slight increase of Poisson's ratio together with a significant density increase.     

Heat and mass transfer in a thermochemical plume under an oceanic plate far from the mid-ocean ridge axis

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit. The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v. The main equations of heat and mass transfer are obtained for a thermochemical plume interacting with a horizontal convective mantle flow. The joint multiparameter problem of heat and mass transfer is solved for a thermochemical plume located far from an MOR axis. The dope concentration at the base of the plume is found as a function of the Lewis number. The Lewis numbers and, accordingly, the diffusion coefficients of the chemical dope in the plume conduit far from the MOR axis are determined.     

Velocity and turbulence variations at the edge of saltmarshes

Saltmarsh vegetation significantly influences tidal currents and sediment deposition by decelerating the water velocity in the canopy. In order to complement previous field results, detailed profiles of velocity and turbulence were measured in a laboratory flume. Natural Spartina anglica plants were installed in a 3 m length test section in a straight, recirculating flume. Different vegetation densities, water depths and surface velocities were investigated. The logarithmic velocity profile, which existed in front of the vegetation, was altered gradually to a skimming-flow profile, typical for submerged saltmarsh vegetation. The flow reduction in the denser part of the canopy also induced an upward flow (the current was partially deflected by the canopy). The skimming flow was accompanied by a zone of high turbulence co-located with the strongest velocity gradient. This gradient moved upward and the turbulence increased with distance from the edge of the vegetation. Below the skimming flow, the velocity and the turbulence were low. The structure of the flow in the canopy was relatively stable 2 m into the vegetation. The roughness length (z₀) of the vegetation depends only on the vegetation characteristics, and is not sensitive to the current velocity or the water depth. Both the reduced turbulence in the dense canopy and the high turbulence at the top of the canopy should increase sediment deposition. On the other hand, the high turbulence zone just beyond the vegetation edge and the oblique upward flow may produce reduced sedimentation; a phenomenon that was observed near the vegetation edge in the field.     

Velocity and density contrasts across the Moho interface of Guangdong province in south China constrained by receiver functions*

The P receiver function includes P-to-SV converted phases and multiple reverberations of the discontinuities in the crust and mantle. The time of these phases is related to the crustal thickness and v_P/v_S ratio, and the amplitude of these phases is mainly controlled by the velocity and density contrast of interfaces. By using H-κ stacking method, this work estimated the crustal thickness and v_P/v_S ratio beneath the stations in the Guangdong province of South China. The velocity and density contrast (δβ-δρ) scanning stacking algorithm of the receiver function is applied to constrain the velocity and density contrast of the Moho in Guangdong province. This work analyzed the results of the crustal thickness, v_P/v_S ratio, and the velocity and density contrasts of Moho. The results indicate that the velocity contrast is higher beneath Yangjiang area in western Guangdong province and Nanao area in eastern Guangdong, which has a strong correlation with the distribution of geothermal springs in local areas and the characteristics of high heat flow. The velocity contrast of Moho has also a good correlation with the v_P/v_S ratio and the crustal thickness, which indicates that there is a strong material composition contrasts of the Moho in the study area. Velocity and density contrasts of Moho in some local area (such as western Guangdong) are somewhat consistent with the seismic activities.     

The thermodynamic properties of the earth's lower mantle

The thermodynamic properties of the lower mantle are determined from the seismic profile, where the primary thermodynamic variables are the bulk modulus K and density ρ. It is shown that the Bullen law (K ∝ P) holds in the lower mantle with a high correlation coefficient for the seismic parametric Earth model (PEM). Using this law produces no ambiguity or trade-off between ρ₀ and K₀, since both K₀ and K′₀ are exactly determined by applying a linear K?ρ relationship to the data. On the other hand, extrapolating the velocity data to zero pressure using a Birch-Murnaghan equation of state (EOS) results in an ambiguous answer because there are three unknown adjustable parameters (ρ₀, K₀, K′₀) in the EOS.From the PEM data, K = 232.4 + 3.19 P (GPa). The PEM yields a hot uncompressed density of 3.999 ± 0.0026 g cm^?3 for material decompressed from all parts of the lower mantle. Even if the hot uncompressed density were uniform for all depths in the lower mantle, the cold uncompressed mantle would be inhomogeneous because the decompression given by the Bullen law crosses isotherms; for example, the temperature is different at different depths. To calculate the density distribution correctly, an isothermal EOS must be used along an isotherm, and temperature corrections must be placed in the thermal pressure P_TH.The thermodynamic parameters of the lower mantle are found by iteration. Values of the three uncompressed anharmonic parameters are first arbitrarily selected: α₀ (hot), the coefficient of thermal expansion; γ₀, the Grüneisen parameter; and δ, the second Grüneisen parameter. Using γ₀ and the measured ρ₀ (hot) and K₀ (hot), the values of θ₀ (Debye temperature) and q = dlnγ/dlnρ are found from the measured seismic velocities. Then from (αK_T)₀ and q the thermal pressure P_TH at all high temperatures is found. Correlating P_TH against T to the geotherm for the lower mantle, P_TH is found at all depths Z. The isothermal pressure, along the 0 K isotherm, at every Z is found by subtracting P_TH from the measured P given by the seismic model. Using the isothermal pressure at depth Z, the solution for the cold uncompressed density ρ_0C and the cold uncompressed bulk modulus, K_T0 is found as a trace in the K_T0?_ρ0C plane. A narrow band of solutions is then found for ρ_0C and K_T0 at all depths.The thermal expansion at all T is found from [ρ_0C ? ρ₀ (hot)/ρ_0C. From Suzuki's formula, the best fit to the thermal expansion determines γ₀ and α₀ (hot). When the values of these two parameters do not agree with the original assumptions, the calculation is repeated until they do agree. In this way all the important thermodynamic parameters are found as a self-consistent set subject only to the assumptions behind the equations used.     

Generation of komatiite magma and gravitational differentiation in the deep upper mantle

The densities of silicate liquids with basic, picritic, and ultrabasic compositions have been estimated from the melting curves of minerals at high pressures. Silicate liquids generated by partial melting of the upper mantle are denser than olivine and pyroxenes at pressures higher than 70 kbar, and garnet is the only phase which is denser than the liquid at pressures from 70 kbar to at least 170 kbar. In this pressure range, garnet and some fraction of liquid separate from ascending partially molten diapirs. It is therefore suggested that aluminium-depleted komatiite with a high Ca/OAl₂O₃ ratio may be derived from diapirs which originated in the deep upper mantle at pressures from 70 kbar to at least 140 kbar (200–400 km in depth), where selective separation of pyropic garnet occurs effectively. On the other hand, aluminium-undepleted komatiite is probably derived from diapirs originating at shallower depths (< 200 km). Enrichment of pyropic garnet is expected at depths greater than 200 km by selective separation of garnet from ascending diapirs. The 200-km discontinuity in the seismic wave velocity profile may be explained by a relatively high concentration of pyropic garnet at depths greater than 200 km.     

Regional induction studies: A review of methods and results

The geomagnetic skin-effect is specified by setting three length scales in relation to each other: L₁ for the overhead source. L₂ for the lateral non-uniformity of the subsurface conductor, L₃ for the depth of penetration of a quasi-uniform transient field into this conductor. Relations for the skin-effect of a quasi-uniform source in layered conductors are generalized to include sources of any given geometry by introducing response kernels as functions of frequency and distance. They show that only those non-uniformities of the source which occur within a distance comparable to L₃ from the point of observation are significant. The skin-effect of a quasi-uniform source in a laterally non-uniform earth is expressed by linear transfer functions for the surface impedance and the surface ratio of vertical/horizontal magnetic variations. In the case of elongated structures and E-polarisation of the source, a modified apparent resistivity is defined which as a function of depth and distance gives a first orientation about the internal distribution of conductivity. The skin-effect of a non-uniform source in a non-uniform earth is considered for stationary and “running” sources. Recent observations on the sea floor and on islands indicate a deep-seated change of conductivity at the continent—ocean transition, bringing high conductivity close to the surface, a feature which may not prevail, however, over the full width of the ocean. There is increasingly reliable evidence for high conductivities (0.02 to 0.1 micro ^?1 m^?1) at subcrustal or even at crustal depth beneath certain parts of the continents, in some cases without obvious correlation to geological structure.     

QLM9: A new radial quality factor (Qμ) model for the lower mantle

We employ a niching genetic algorithm to invert ∼30,000 differential ScS/S attenuation values for a new spherically symmetric radial model of shear quality factor (Q_μ) with high sensitivity to the lower mantle. The new radial Q_μ model, QLM9, possesses greater sensitivity to Q_μ at large mantle depths than previous studies. On average, lower mantle Q_μ increases with depth, which supports models of increasing viscosity with depth [B.M. Steinberger, A.R. Calderwood. Mineral physics constraints on viscous flow models of mantle flow, J. Conf. Abs., 6, 2001., 2001.]. There are two higher-Q_μ regions at ∼1000 and ∼2500 km depth, which roughly correspond to high-viscosity regions observed by Forte and Mitrovica [A.M. Forte and J.X. Mitrovica, Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data, Nature 410, 1049–1056, 2001.]. There is a lower-Q_μ layer at the core–mantle boundary and a relatively low-Q_μ region in the mid-lower mantle. With several caveats, we infer a divergence of the solidus and geotherm in the lower mantle and a convergence within Dʺ by relating Q_μ to homologous temperature.     

Temperature dependence of the elastic moduli of ringwoodite

The elastic moduli of polycrystalline ringwoodite, (Mg_0.91Fe_0.09)₂SiO₄, were measured up to 470 K by means of the resonant sphere technique. The adiabatic bulk (K_S) and shear (μ) moduli were found to be 185.1(2) and 118.22(6) GPa at room temperature, and the average slopes of dK_S/dT and dμ/dT in the temperature range of the study were determined to be −0.0193(9) and −0.0148(3) GPa/K, respectively. Using these results, we estimate seismic wave velocity jumps for a pure olivine mantle model at 520 km depth. We find that the jump for the S-wave velocity is about 1.5 times larger than that for the P-wave velocity at this depth. This suggests that velocity jumps at the 520 km discontinuity are easier to detect using S-waves than P-waves.     

Should average shear-wave velocity in the top 30 m of soil be used to describe seismic amplification?

The average velocity of shear waves in the top 30 m of soil, ν_L, has become the parameter used by many engineering design codes and most recently by published empirical-scaling equations to estimate the amplitudes of strong ground motion. Yet there are few studies to determine whether this is a meaningful parameter to use—and whether estimates that do use it are reliable. In 1995, the authors studied this problem and concluded that ν_L should not be used. We reported then that an older site characterization in terms of the soil site parameter proposed by Seed et al. [1], s_L, worked better because it included a measure of the thickness of the soil layers to considerably greater depths. Our report, however, made no difference; numerous papers continued to be published based only on scaling in terms of ν_L, and worse, they also ignored the geological site conditions. The purpose of this paper is to emphasize that the average shear-wave velocity in the top 30 m of soil should not be the only site parameter used to scale strong-motion amplitudes. While the search continues for the more meaningful site parameters to use in empirical scaling of strong earthquake ground motion, it is better to use s_L to describe the amplification of seismic waves by soil deposits near the surface.     

Seismicity and fluid geochemistry at Lassen Volcanic National Park,California: Evidence for two circulation cells in the hydrothermal system

Seismic analysis and geochemical interpretations provide evidence that two separate hydrothermal cells circulate within the greater Lassen hydrothermal system. One cell originates south to SW of Lassen Peak and within the Brokeoff Volcano depression where it forms a reservoir of hot fluid (235–270 °C) that boils to feed steam to the high-temperature fumarolic areas, and has a plume of degassed reservoir liquid that flows southward to emerge at Growler and Morgan Hot Springs. The second cell originates SSE to SE of Lassen Peak and flows southeastward along inferred faults of the Walker Lane belt (WLB) where it forms a reservoir of hot fluid (220–240 °C) that boils beneath Devils Kitchen and Boiling Springs Lake, and has an outflow plume of degassed liquid that boils again beneath Terminal Geyser. Three distinct seismogenic zones (identified as the West, Middle, and East seismic clusters) occur at shallow depths (< 6 km) in Lassen Volcanic National Park, SW to SSE of Lassen Peak and adjacent to areas of high-temperature (≤ 161 °C) fumarolic activity (Sulphur Works, Pilot Pinnacle, Little Hot Springs Valley, and Bumpass Hell) and an area of cold, weak gas emissions (Cold Boiling Lake). The three zones are located within the inferred Rockland caldera in response to interactions between deeply circulating meteoric water and hot brittle rock that overlies residual magma associated with the Lassen Volcanic Center. Earthquake focal mechanisms and stress inversions indicate primarily N–S oriented normal faulting and E–W extension, with some oblique faulting and right lateral shear in the East cluster. The different focal mechanisms as well as spatial and temporal earthquake patterns for the East cluster indicate a greater influence by regional tectonics and inferred faults within the WLB. A fourth, deeper (5–10 km) seismogenic zone (the Devils Kitchen seismic cluster) occurs SE of the East cluster and trends NNW from Sifford Mountain toward the Devils Kitchen thermal area where fumarolic temperatures are ≤ 123 °C. Lassen fumaroles discharge geothermal gases that indicate mixing between a N₂-rich, arc-type component and gases derived from air-saturated meteoric recharge water. Most gases have relatively weak isotopic indicators of upper mantle or volcanic components, except for gas from Sulphur Works where δ¹³C–CO₂, δ³⁴S–H₂S, and δ¹⁵N–N₂ values indicate a contribution from the mantle and a subducted sediment source in an arc volcanic setting.     

Models of density and mechanical Q-factor distributions according to new data on the nutations and free oscillations of the Earth: 1. Ambiguity in the inverse problem

Ambiguity in the inverse problem of retrieval of the mechanical parameters of the Earth’s shell and core from the set of data on the velocities V _p and V _S , of longitudinal and transverse seismic body waves, the frequencies f _i and quality factors Q _i , of free oscillations, and the amplitudes and phases of forced nutation is considered. The numerical experiments show that the inverse problem of simultaneous retrieval of the density profile ρ in the mantle-liquid core system and the mechanical quality factor Q _μ of the mantle (if the total mass M and the total mean moment of inertia I of the Earth, and V _p and V _S are constant at all depths) has most unstable solutions. An example of depth distributions of ρ and Q _μ which are alternative to the well-known PREM model is given. In these distributions, the values of M and I and the velocities V _p and V _S at all depths for the period of oscillations T = 1 s exactly coincide with their counterparts yielded by PREM model (T = 1 s); however, the maximum deviations of the ρ and Q _μ profiles from those in the PREM model are about 3% and 40%, respectively; the mass and the moment of inertia of the liquid core are smaller than those for the PREM model by 0.75% and 0.63%, respectively. In this model, the root mean square (rms) deviations of all the measured values of f _i and Q _i from their values predicted by theory are half to third the corresponding values in the PREM model; the values of Δ for natural frequencies of the fundamental tone and overtones of radial oscillations, the fundamental tones of torsional oscillations, and the fundamental tones of spheroidal oscillations, which are measured with the highest relative accuracy, are smaller by a factor of 30, 6.6, and 2 than those in the PREM model, respectively. Such a large ambiguity in the solution of the inverse problem indicates that the current models of the depth distribution of density have relatively low accuracy, and the models of the depth distribution of the mechanical Q in the mantle are extremely unreliable. It is shown that the ambiguity in the models of depth distribution of density considerably decreases after the new data on the amplitudes and phases of the forced nutation of the Earth are taken into account. Using the same data, one may also refine by several times the recent estimates of the creep function for the lower mantle within a wide interval of periods ranging from a second to a day.     

Heat transfer between a thermochemical plume channel and the surrounding mantle in the presence of horizontal mantle flow

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.