A New Model for Heat Flow in Extensional Basins: Estimating Radiogenic Heat Production

Radiogenic heat production (RHP) represents a significant fraction of surface heat flow, both on cratons and in sedimentary basins. RHP within continental crust—especially the upper crust—is high. RHP at any depth within the crust can be estimated as a function of crustal age. Mantle RHP, in contrast, is always low, contributing at most 1 to 2 mW/m² to total heat flow. Radiogenic heat from any noncrystalline basement that may be present also contributes to total heat flow. RHP from metamorphic rocks is similar to or slightly lower than that from their precursor sedimentary rocks. When extension of the lithosphere occurs—as for example during rifting—the radiogenic contribution of each layer of the lithosphere and noncrystalline basement diminishes in direct proportion to the degree of extension of that layer. Lithospheric RHP today is somewhat less than in the distant past, as a result of radioactive decay. In modeling, RHP can be varied through time by considering the half lives of uranium, thorium, and potassium, and the proportional contribution of each of those elements to total RHP from basement. RHP from sedimentary rocks ranges from low for most evaporites to high for some shales, especially those rich in organic matter. The contribution to total heat flow of radiogenic heat from sediments depends strongly on total sediment thickness, and thus differs through time as subsidence and basin filling occur. RHP can be high for thick clastic sections. RHP in sediments can be calculated using ordinary or spectral gamma-ray logs, or it can be estimated from the lithology.     

Global heat budget, plate tectonics and climatic change

For the past 2000 Ma, the temperature of the Earth's surface has fluctuated around a mean similar to that of today, although individual locations have undergone long-term changes of ∼30°C at different times in different places. Water bodies absorb at least five times as much solar radiation as land surfaces, and ocean currents transport the excess heat absorbed in the tropics towards the poles. Changes in the distribution of land and sea due to plate tectonics explain the major temperature fluctuations (>25°C) around the globe in the last 350 Ma, and are first-order controls. Large-scale changes in ocean currents and thermohaline circulations are probably second-order controls (15–25°C). The Milankovitch orbital cycles are third-order controls producing variations in air temperature of the order of 10°C, while massive volcanic eruptions and changes in carbon dioxide are amongst the fourth-order controls producing minor perturbations (<5°C). The major climatic fluctuations are continuous but regional in effect and not global. Extraterrestrial factors may not cause major changes in climate when viewed from a geological perspective.     

Characterisation of the 1997 Vulcanian explosions of Soufrière Hills Volcano,Montserrat, by video analysis

The activity of Convention at Montserrat Soufrière Hills Volcano, Montserrat, during the period 1995–1999 included numerous violent explosions. Two major cycles of Vulcanian explosions occurred in 1997: a first of 13 explosions between 4 and 12 August and a second of 75 between 22 September and 21 October. The explosions were short-lived events lasting a few tens of seconds during which partial fountain collapse generated pyroclastic surges and pyroclastic flows, and buoyant plumes ascended 3–15 km into the atmosphere. Each explosion discharged on average 3×10⁵ m³ (dense-rock equivalent, DRE) of magma, draining the conduit to depths of 1–2 km. The paper focuses on the first few seconds of three explosions of the 75 that occurred in September/October 1997: 6 October 1997 at 17:50, 7 October 1997 at 16:02 and 9 October 1997 at 12:32. Physical parameters such as exit velocities, magmatic water contents and magma pressures at fragmentation are estimated by following and modelling the ascent of individual momentum-dominated finger jets visible on videos during the initial stages of each explosion. The model treats each finger jet as an incompressible flow sustained by a steady flux of gas and particles during the few seconds of ascent, and produces results that compare favourably with those using a multiphase compressible code run using similar eruptive parameters. Each explosion reveals a progressive increase in eruptive intensity with time, jet exit velocities increasing from 40 m s^–1 at the beginning of the explosion up to 140 m s^–1 after a few seconds. Modelling suggests that the first magma to exit was largely degassed, whereas that discharged after a few seconds contained up to 2 wt% water. Magma overpressures up to ~10 MPa are estimated to have existed in the conduit immediately prior to each explosion. Progressive increases in jet exit velocity with time over the first few seconds of each explosion provide direct evidence for strong pre-eruptive gradients in water content and magma pressure in the upper reaches (probably 100–500 m) of the conduit. Fountain collapse occurred during the first 10–20 s of each explosion because the discharging jets had bulk densities up to 100 times that of the atmosphere and were unable to entrain enough air to become buoyant. Such high eruptive densities were due to the presence of partially degassed magma in the conduit.Editorial responsibility: A. Woods     

Velocity structure of uppermost mantle beneath North China from Pn tomography and its implications

20301 Pn arrival time data are collected from the seismological bulletins of both national and regional seismic networks. Pn travel time residuals are tomographically inverted for the Pn velocity structure of uppermost mantle beneath North China. The result indicates that the average Pn velocity in North China is 7.92 km/s, and the velocity varies laterally from ?0.21 to +0.29 km/s around the average. The approximately NNE trending high and low velocity regions arrange alternatively west-eastward. From west to east we can see high velocity in the middle Ordos region, the Shanxi graben low, the Jizhong depression high, the west Shandong uplift and Bohai Sea low, and the high velocity region to the east of the Tanlu fault. In the southern boundary zone of the North China block, except for the high velocity in the Qingling Mountains region, the velocity is generally lower than the average. Obvious velocity anisotropy is seen in the Datong Cenozoic volcanic region, with the fast velocity direction in NNE-SSW. Notable velocity anisotropy is also seen around the Bay of Bohai Sea, and the fast velocity directions seem to show a rotation pattern, possibly indicating a flow-like deformation in the uppermost mantle there. The Pn velocity variations show a reversed correlation with the Earth's heat flow. The low Pn velocity regions generally show high heat flow, e.g., the Shanxi graben and Bohai Sea region. While the high Pn velocity regions usually manifest low heat flow, e.g., the region of Jizhong depression. This indicates that the Pn velocity variation in the study region is mainly aroused by the regional temperature difference in the uppermost mantle. Strong earthquakes in the crust tend to occur in the region with the abnormal low Pn velocity, or in the transition zone between high and low Pn velocity regions. The earthquakes in the low velocity region are shallower, while that in the transition zone are deeper.     

Simulation of marine boundary layer characteristics using a 1-D PBL model over the Bay of Bengal during BOBMEX-99

The characteristic features of the marine boundary layer (MBL) over the Bay of Bengal during the southwest monsoon and the factors influencing it are investigated. The Bay of Bengal and Monsoon Experiment (BOBMEX) carried out during July–August 1999 is the first observational experiment under the Indian Climate Research Programme (ICRP). A very high-resolution data in the vertical was obtained during this experiment, which was used to study the MBL characteristics off the east coast of India in the north and south Bay of Bengal. Spells of active and suppressed convection over the Bay were observed, of which, three representative convective episodes were considered for the study. For this purpose a one-dimensional multi-level PBL model with a TKE-ε closure scheme was used. The soundings, viz., the vertical profiles of temperature, humidity, zonal and meridional component of wind, obtained onboard ORV Sagar Kanya and from coastal stations along the east coast are used for the study. The temporal evolution of turbulent kinetic energy, marine boundary layer height (MBLH), sensible and latent heat fluxes and drag coefficient of momentum are simulated for different epochs of monsoon and monsoon depressions during BOBMEX-99.The model also generates the vertical profiles of potential temperature, specific humidity, zonal and meridional wind. These simulated values compared reasonably well with the observations available from BOBMEX.     

Magnetotelluric Investigation of the Hydrothermal System and Heat Source in the Muine-Toyoha Geothermal Area, Hokkaido, Japan

Abstract. Magnetotelluric (MT) surveys were carried out around the Muine volcano, Hokkaido, Japan, where it is expected that the heat and metal source forming the polymetallic Ag-Pb-Zn-Cu-In Toyoha deposit is present at depth. Measurements were performed at 20 sites, 18 of which were located along a WSW-ENE profile traversing the north ridge of Mt. Muine. A resistivity model obtained from 2D inversion of the MT data shows subsurface specific conductive and resistive features. Conductive layers are present at the surface of Mt. Muine. The low resistivity is probably due to the clay-rich rocks associated with the hydrothermal alteration. A high resistivity layer, which corresponds to the pre-Tertiary Usubetsu Formation, crops out east of Mt. Muine and dips westward. At the west foot of Mt. Muine, relatively high resistive layers are widely exposed. The resistivity increases with depth and exceeds 1000 ohm-m. This fact indicates that this region is not influenced by the recent hydrothermal activity. An extremely conductive zone about 3–6 km wide and 6–9 km thick exists at a depth of 2 km below Mt. Muine. This zone mostly corresponds to an elastic wave attenuation zone detected by a seismic survey. It is interpreted as a large hydrothermal reservoir or melted magma, which is a heat source of the hydrothermal system in this area.     

粘粒含量对南京粉细砂液化影响的试验研究

南京粉细砂是一种以片状颗粒成分为主的粉细砂，与通常的圆形颗粒石英砂有一定区别，片状颗粒成分使得南京粉细砂具有各向异性的性质。通过对不同粘粒含量的南京粉细砂进行液化试验，分析其试验结果发现：对于粘粒含量ρc＝5％、10％和15％三种情况，在均压固结情况下，当南京粉细砂达到初始液化时，采用变形标准和孔压标准是一致的；而偏压固结时，只能采用变形标准作为南京粉细砂达到初始液化的标准，此时振动孔压仅达到围压的50％-70％。粘粒含量对南京粉细砂的抗液化强度影响很大，随粘粒含量增加，在不同固结比时，其抗液化强度并不是单调变化，而是在粘粒含量ρc＝10％左右处于一个低谷，其抗液化强度最低。     

Determination of Roughness Lengths for Heat and Momentum Over Boreal Forests

The roughness length for momentum (z_0m), zero-plane displacementheight (d), and roughness length for heat (z_0h) are importantparameters used to estimate land-atmosphere energy exchange. Although many different approaches have been developed to parameterizemomentum and heat transfer, existing parameterizations generally utilizehighly simplified representations of vegetation structure. Further, a mismatch exists between the treatments used for momentum and heat exchange and those used for radiative energy exchanges. In this paper, parameterizations are developed to estimate z_0m, d, and z_0h for forested regimes using information related to tree crown density and structure. The parameterizations provide realistic representationfor the vertical distribution of foliage within canopies, and include explicit treatment for the effects of the canopy roughness sublayer and leaf drag on momentum exchange. The proposed parameterizationsare able to realistically account for site-to-site differences in roughness lengths that arise from canopy structural properties.Comparisons between model predictions and field measurements show good agreement, suggesting that the proposed parameterizations capture the most important factors influencing turbulent exchange of momentumand heat over forests.     

A Technique for Direct Retrieval of Surface Specific Humidity Over Open Oceans from Irs-P4/Msmr Satellite Data

The Multi frequency Scanning Microwave Radiometer (MSMR) onboard Oceasat-1 was used to develop a retrieval method fornear-surface specific humidity by means of multivariate regressiontechnique. The MSMR measures the microwaveradiances in 8 channels at the frequencies of 6.6, 10.7, 18 and 21 GHzfor both vertical and horizontal polarizations. Regression coefficients were derived using the ship reports of the Comprehensive Ocean-Atmosphere Data Set (COADS) for the months of July, October and December, in 1999. Daily near-surface specific humidity data from COADS in 2° × 2° latitude/longitude bins and collocated brightness temperature data from MSMR were used to derive the coefficients. The derived coefficients werevalidated with humidity given in COADS.A linear relationship is established to determine the near-surface specifichumidity from MSMR brightness temperature (Tb) with an rms error of 1.2 g kg^-1 for individual situations and an rms errorof 0.84 g kg^-1 for monthly time scales over global oceans.The retrieval algorithm is validonly for the open sea regions.