Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

Effects of Coal Gangue Content on Water Movement and Solute Transport in a China Loess Plateau Soil

The mining industry has grown strongly in China in recent decades, resulting in large amounts of coal gangues, which cause water and soil pollution, soil erosion, and various other environmental problems. They are often used in reclamation projects in attempts to restore land damaged by mining, hence they are frequently present (in widely varying proportions) in the topsoil in areas around mines. Their presence can strongly affect key soil variables, including its bulk density, structure, water retention, water movement, and solute transport rates. In the study presented here, the effects of gangue contents on infiltration, saturated hydraulic conductivity, and solute transport parameters of a Chinese Loess plateau soil were examined. The results show that infiltration rates and saturated hydraulic conductivity decreased with increasing gangue content. The Peck–Watson equation modeled these relationships well, but Bouwer–Rice equations provided poorer matches with the acquired data. Cumulative infiltration over time was described well by both the Philip equation and Kostiakov equation. Both the simplified convection–dispersion equation and a two‐region model described the solute transport processes well. In addition, the dispersion increased, while both the Peclet number and mobile water fraction decreased, with increases in gangue contents.     

Evaluating the quality of hydraulic conductivity estimates from piezometer slug tests in peat

Although widely used in wetland hydrological studies, hydraulic conductivity (K) estimates from piezometer slug tests are often of questionable validity. Frequently, this is because insufficient attention is paid to the details of the test procedure. Further, in a potentially heterogeneous and anisotropic medium such as peat, the use of slug tests is prone to error. In this paper we address some of the methodological issues surrounding piezometer slug tests in peat. We compare slug test data with laboratory determinations of vertical and horizontal K obtained using a new method. Piezometers were installed at three depths in a floodplain fen peat in Norfolk, UK. Slug tests were initiated by both slug insertion and slug withdrawal, and repeat tests were conducted to examine the robustness of our K estimates. Most of the tests displayed departures from the log‐linear model of Hvorslev, the form of departure being consistent with compressible soil behaviour. The results suggest that insertion tests gave similar results to those initiated by withdrawal. Repeat testing showed that withdrawal data, in particular, gave highly reproducible normalized responses that were independent of the initial head. Values for K estimated using the slug tests were in the range 1 × 10⁻⁴ to 1·6 × 10⁻³ cm s⁻¹, which is towards the upper end of the range reported for peats generally. Laboratory tests yielded similar values of K to those obtained from the slug tests. Although the laboratory tests showed that the peat was anisotropic, the K values generated by slug testing proved relatively good estimates of both vertical and horizontal K. Copyright © 2005 John Wiley & Sons, Ltd.     

A thermo‐mechanical model for the catastrophic collapse of large landslides

In this work, a new thermo‐mechanical model is developed, applicable to large‐scale, deep‐seated landslides consisting of a coherent mass sliding on a thin clayey layer. The considered time window is that of catastrophic acceleration, starting at incipient failure and ending when the acquired displacement and velocity are such that the sliding material begins to break up into pieces. The model accounts for temperature rise in the slip zone due to the heat produced by friction, leading to water expansion, thermoplastic collapse of the soil skeleton, and subsequent increase of pore water pressure. The model incorporates the processes of heat production and diffusion, pore pressure generation and diffusion, and an advanced constitutive law for the thermo‐mechanical behavior of soil. An analysis of the Vajont landslide is presented as an example. A sensitivity analysis shows that friction softening is the mechanism most affecting the timescale of the final collapse of a slide, but also that the mechanism of thermal pressurization alone can cause a comparably catastrophic dynamic evolution. It is also shown that, all other factors being equal, thermo‐mechanical collapse will cause thicker slides to accelerate faster than shallow ones. Copyright © 2010 John Wiley & Sons, Ltd.     

Gamma-ray burst precursors as the remnant of the thermal radiation initially trapped in the fireball

In the standard fireball model of gamma-ray bursts (GRBs), the fireball starts with an optically thick phase. As it expands, the fireball becomes optically thin at some stage. The thermal radiation trapped in the originally opaque fireball then leaks out, producing a transient event. The appearance of the event is investigated in the framework of a homogeneous, spherically symmetric and freely expanding fireball produced instantly by an explosive process without continuous injection of mass and energy. We find that, generally, the event has a time duration shorter than that of the main burst, which is presumably produced by the internal shock after the fireball becomes optically thin. The event is separated from the main burst by a quiescent time interval, and is weaker than the main burst at least in a high-energy band. Hence, the event corresponds to a GRB precursor. The precursor event predicted by our model has a smooth and Fast Rise and Exponential Decay (FRED) shaped light curve, and a quasi-thermal spectrum. Typically, the characteristic blackbody photon energy is in the X-ray band. However, if the distortion of the blackbody spectrum by electron scattering is considered, the characteristic photon energy could be boosted to the gamma-ray band. Our model may explain a class of observed GRB precursors – those having smooth and FRED-shaped light curves and quasi-thermal spectra.     

Multi‐depth monitoring of electrical conductivity and temperature of groundwater at a multilayered coastal aquifer: Jeju Island,Korea

To supplement conventional geophysical log data, this study presents temporal variations in electrical conductivity (EC) and temperature with depth in a multilayered coastal aquifer, on the eastern part of Jeju Island, Korea. One‐month time‐series data obtained at eight points from a multi‐depth monitoring system showed that semidiurnal and semimonthly tidal variations induced dynamic fluctuations in EC and temperature. At some depths, EC ranged from 1483 to 26 822 µS cm^?1, while some points showed no significant variations. The results of EC log and time‐series data revealed that a sharp fresh‐saltwater interface occurred at low tide, but the diffusion zone broadened to 20 m at high tide. EC, temperature, and tide level data were used for the cross‐correlation analysis. The response time of EC and temperature to tide appears to range from less than 30 min to 11 h. Using end‐member mixing analysis (EMMA), the fraction of variations of chloride concentration in the multilayered aquifer was explained, and a conceptual model was developed which subdivided the coastal aquifer into four vertical zones. The percentage of water derived from seawater varied from 2 to 48 at specific depth, owing to tidal fluctuations. Continuous observations of EC and temperature at multiple depths are powerful tools for quantifying the transport of saline water by tidal variations in multilayered coastal aquifers. Copyright © 2008 John Wiley & Sons, Ltd.     

Estimating finite difference block equivalent hydraulic conductivity for numerically solving the Richards' equation

A method is proposed for calculating the equivalent hydraulic conductivity (EHC) within a finite difference block (FDB). Application of the constant‐flux assumption of Darcy's Law, the EHC equals to the integration of effective hydraulic conductivity (K_w) as a function of pressure head (h_w) divided by the head difference at the ends of the FDB. Error analysis show that the constant‐flux (CF) EHC estimates are better than those computed by the commonly used arithmetic‐mean (AM), geometric‐mean (GM), and harmonic‐mean (HM) techniques. CF EHC results are even more superior at larger interblock head difference situations. Simulations of water infiltration experiments show that simulations using the CF EHC or AM or GM weighting technique have only slight difference while applying the Neumann type boundary condition at the ground surface. In case of the Dirichlet type boundary condition, however, the CF EHC is superior to the other two in correctly estimating the depth of infiltration while enlarging the grid size. Therefore, it is recommended to adopt the CF EHC with a larger grid size to the more stable and more efficient results. Copyright © 2007 John Wiley & Sons, Ltd.     

Metamorphic P–T conditions and thermal structure of Chinese Continental Scientific Drilling main hole eclogites: Fe–Mg partitioning thermometer vs. Zr-in-rutile thermometer

Core rocks recovered from the main hole (5158 m deep) of the Chinese Continental Scientific Drilling (CCSD‐MH) project, southern Sulu UHP terrane, east‐central China, consist of eclogites, various gneisses and minor metaperidotite cumulates; this lithological section underwent subduction‐zone UHP metamorphism. Coesite‐bearing eclogites are mainly present between the depths of 100–2000 m, but below 2000 m, mafic eclogites are rare. Selected elements (Zr, Nb, Cr, Fe, Si, Mg, Al & Ti) in rutile from 39 eclogite cores from 100 to 2774 m, and major elements of minerals from representative eclogites were analysed by electron microprobe. Zirconium and Nb concentrations of rutile cluster ～100–400 and 200–700 ppm respectively. However, Zr and Nb contents in rutile from strongly retrograded eclogites show larger variations than those of fresh or less retrograded eclogites, implying that somehow fluid infiltration affected rutile chemistry during retrograde metamorphism. Zr contents in rutile inclusions in garnet and omphacite are slightly lower than those of the matrix rutile, suggesting that the rutile inclusions formed before or close to the peak temperature. The P–T conditions of the CCSD‐MH eclogites were estimated by both Fe–Mg exchange and Zr‐in‐rutile thermometers, as well as by the Grt–Cpx–Phn–Ky geothermobarometer. The maximum temperature range of 700–811 °C calculated at 40 kbar using the Zr‐in‐rutile thermometer is comparable with temperature estimates by the Fe–Mg exchange thermometer. The temperature estimates of eclogites in a ～3000 m thick section define a continuous gradient, and do not show a distinct temperature gap, suggesting that the rocks from 100 to 3000 m depth might belong to a single, large‐scale UHP slab. These data combined with P–T calculations for CCSD‐MH peridotites yield a low geotherm (～5 °C km^?1) for the Triassic subduction zone between the Sino‐Korean and Yangtze cratons; it lies ～30–35 mW m^?2 conductive model geotherm.