Regional variations of heat flow differences with depth in Alberta, Canada

Summary. Differences between estimated average heat flow values for the Mesozoic and Cenozoic formations ( Q ₁) and estimated average heat flow values for the Palaeozoic formations below the erosional unconformity ( Q ₂) are calculated for the Alberta part of the western Canadian sedimentary basin. Significant heat flow differences exist for these two intervals and the map of Δ Q = Q ₁– Q ₂ shows that Q ₂ is generally greater than Q ₁ in the western and south-western part of Alberta, while in the northern part of the province Q ₂ is generally less than Q ₁. The regional variations of Δ Q are large, with standard deviation of 26 mW m⁻² and average value –13.5 mW m⁻². A regional trend of Δ Q correlates with topographic relief and the hydraulic head variations in the basin. It is shown that there is a heat flow increase with depth in water recharge areas and a decrease in heat flow with depth in the low topographic elevation water discharge areas when comparing the average heat flow in Mesozoic + Cenozoic and Palaeozoic formations.     

Heat refraction and heat production in and around granite plutons in north-east England

Summary. Existing surface heat flow determinations in north-east England indicate a difference in thermal structure between the Alston block of the North Pennines and the coastal regions. New heat flow determinations in the Askrigg block and on the coast support the original indications of higher heat flow in the positive blocks, which are underlain by 400 Myr old granites, than in their marginal sedimentary troughs along the present coastline. Although the two blocks are geophysically and geologically similar, surface heat flow in the Alston block is 30mW m^-2 higher than it is in the Askrigg block 50km to its south. Heat flow refraction around high conductivity plutons with high heat production may significantly affect the interpretation of heat flow measurements made on a regional scale but in the present case its effect, if any, is to emphasize the heat flow difference between the two otherwise similar blocks. The chemical differences between the Wensleydale and the Weardale granites – cupolas of the granites underlying the Askrigg and Alston may be representative of the granites as a whole and the blocks respectively – difference in surface heat flows may result from a contrast in radiogenic heat production between the two granites extending to mid-crustal depths.     

TERRESTRIAL HEAT FLOW IN THE SWISS ALPS

i
Terrestrial heat flow has been measured in three Alpine railroad tunnels. The geothermal gradients were calculated from temperatures measured during the construction of the tunnels, and corrections for topographic irregularities were made. The thermal conductivity of 113 rock specimens from the vicinity of the tunnels was measured. The heat flow in the Gotthard tunnel was found to be 1.6 10^-6 cal/cm² sec, in the Simplon 2.2 10^-6 cal/cm² sec, and in the Loetschberg 1.9 10^-6 cal/cm² sec. Most of the flux at the surface can be attributed to radioactive decay in a thickened crust, but a non-uniform distribution of radioactive elements may be required to explain the relatively high heat flow in the Simplon and Loetschberg tunnels.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

Long-period (30 days1 year) electromagnetic sounding and the electrical conductivity of the lower mantle beneath Europe

The C -response connects the magnetic vertical component and the horizontal gradient of the horizontal components of electromagnetic variations and forms the basis for deriving the conductivitydepth profile of the Earth. Time-series of daily mean values at 42 observatories typically with 50 years of data are used to estimate C -responses for periods between 1 month and 1  yr. The Z : Y method is applied, which means that the vertical component is taken locally whereas the horizontal components are used globally by expansion in a series of spherical harmonics.
In combination with results from previous analyses, the method yields consistent results for European observatories in the entire period range from a few hours to 1  yr, corresponding to penetration depths between 300 and 1800  km.
1-D conductivity models derived from these results show an increase in conductivity with depth z to about 2  S  m^-1 at z =800  km, and almost constant conductivity between z =800 and z =2000  km with values of 310  S  m^-1, in good agreement with laboratory measurements of mantle material. Below 2000  km the conductivity is poorly resolved. However, the best-fitting models indicate a further increase in conductivity to values between 50 and 150  S  m^-1.     

Preliminary Observations of 1–16 Hz Natural Background Infrasound and Signals from Apollo 14 and Aircraft

Summary Infrasound with frequencies of 1–16 Hz, detected by an array of four thermistor flow-meter microphones in Sterling Forest, New York, was observed to have a continuous background with peak energy distributed near 16 Hz in frequency, with amplitudes of about 1 dyn cm^-2, and arriving from the south-west and south-east at slightly above the speed of sound in air at ground level. The same array of microphones detected 5 dyn cm^-2 signals from the Apollo 14. The earlier part of the 10-min signal arrived from the first stage re-entry, the later from the launch site vicinity. It is shown that aircraft beyond the visible and audible range can be detected and tracked by monitoring the infrasound emitted throughout most of the 1–16 Hz frequency band.     

Estimating the crust permeability from fluid-injection-induced seismic emission at the KTB site

During the hydraulic-fracturing experiment in the German Continental Deep Drilling Borehole (KTB) in December 1994, microseismic activity was induced. Here we develop a technique for estimating permeability using the spatio-temporal distribution of the fluid-injection-induced seismic emission. The values we have obtained for the KTB experiment (0.25times10^-16 to 1.0times10^-16 1.0times10^-16 m²) are in a very good agreement with the previous hydraulic-type permeability estimates from KTB deep-observatory studies. In addition, our estimates of the hydraulic diffusivity support the previously calculated value for the upper crust, which is of the order of 1 m² s^-1. However, this estimate now relates to the depth range 7.5-9 km.     

Explosive volcanic eruptions — V. Observations of plume dynamics during the 1979 Soufrière eruption, St Vincent

Summary During the 1979 eruption of the Soufrière of St Vincent, West Indies, a 14 min period of explosive activity occurred on April 22. The first three minutes of this activity were filmed. Five individual explosions formed distinct plumes which fed an eruption column which eventually ascended to a height of over 18 km. The ascent velocities of the fronts of these plumes ranged from 8.5 to 61.7 ms^-1. The half-angle of spreading of the plume fronts ranged from 21.5 to 24°. One of the plumes was observed to 8 km height and decelerated steadily from 51.5 to 23 ms^-1. The main column fed by individual explosions was, however, observed to accelerate from 8.5 to 58 m s^-1 in the same height interval.
A theoretical treatment of volcanic plume motion is presented. Measurements of dimensions, velocities and spreading rates from the film are used to estimate plume parameters such as temperature, particle content and volume discharge rate of magma from the theoretical relationships. These calculations show that the individual plumes became hotter and richer in juvenile ash with time. The acceleration of the main eruption column was the result of being fed by increasingly hotter and more ash-rich explosions. An average volume discharge rate of 12600m³ s^-1 is estimated from the analysis of the plume motions. This value agrees closely with an estimate of discharge rate based on the heat flux required to form an 18 km high column. This agreement suggests that theories of convective motions in plumes can be successfully adapted to the volcanic case, as long as the effect of hot entrained ash particles is considered.     

Calibration of the Yellowknife Seismic Array with First Zone Explosions

Summary Recordings from a crustal seismic experiment, which was conducted in the Yellowknife area in 1966, were used for calibration of the Yellow-knife seismic array. In the immediate vicinity of the array the crust is found to be very uniform. A superficial layer with an intercept time of 0–172 ± 0–012s and unknown velocity is underlain by a crust with a P wave velocity of 6.04 ± 0–01 km s^-1 near the top: assuming this velocity constant throughout the second layer, the total thickness of the crust is about 34 ± 2 km. The Mohorovicic discontinuity is horizontal under the array within the resolution of this experiment and the apparent P_n velocity is 8.15 km s^-1. At a distance of a few tens of kilometres the crustal uniformity breaks down. The distances are such that, for most teleseismic signals, the effect of these in homogeneities should be negligible.     

Crust and upper mantle shear velocities from controlled sources

b
A two ship refraction profile was undertaken on the Australian continental shelf during the Banda Sea geophysical program, carried out by the Woods Hole Oceanographic Institution, the Scripps Institution of Oceanography and the Geological Survey of Indonesia. S waves originating close to the sea bottom were observed to distances of up to 1150 km at an array of stations in northern Australia.
These observations are interpreted as implying S mantle velocities of 4.60 km s^-1 from a depth of 45 km to a depth of 76 km and 4.72 km s^-1 below a depth of 76 km.
Ratios of the P and S travel times (V_p/V_s) have been determined to be 1.74 in the crust rising to a value of greater than 1.79 below a velocity discontinuity at a depth of 200 km. It is inferred that this high value arises because the effect of temperature is greater for S than for P .
Using the data from this and other studies in the shield region of Northern Australia it has been found that the _S travel times are significantly less than predicted by the Jeffreys—Bullen tables.     

The outflow of polar water through the Arctic Archipelago and the oceanographic conditions in Baffin Bay

The transports of mass, heat, and salt through the Arctic Archipelago are estimated by two different approaches. First: A baroclinic, rotationally controlled flow is assumed through the Arctic Archipelago and the Davis Strait, with all passages transporting at maximum capacity. Second: Mass, salt, and energy balances are considered for the upper 'cold core' in Baffin Bay. The mass transport in the first case decreases and in the second case increases with increasing salinity in the Beaufort Sea. It is thus possible to determine a salinity at which both approaches give the same mass transport. The outflow in the upper layers is estimated at 0.7 · 10⁹ kg s⁻¹ with a salinity of 32.9. The amount of water from the Polar Ocean entering the deeper (Atlantic) layer in Baffin Bay is more difficult to assess. A tentative value of 0.3 10·kg s⁻¹ with a salinity cf 34.3 is proposed.     

Axial Intrusion Zone beneath the Median Valley of the Mid-Atlantic Ridge at 37°N Detected by Explosion Seismology

Summary. Four seismic refraction lines, three of which had shots every 250 m, were shot across, along and parallel to the median valley of the Mid-Atlantic Ridge at 37° N. A method has been developed for calculating the effect on the travel times of the rough sea-floor relief beneath the profiles and has been used to correct all the travel times for this effect. Most arrivals were from a main refractor of apparent velocity 5·4 to 6·3 km s⁻¹; only beyond 35 km were faster arrivals observed from an 8·09 ± 36 km s⁻¹ refractor. The main refractor corresponds in depth, at least approximately, to the top of Layer 3 of the ocean basins but its velocity is significantly less than normal for Layer 3, perhaps due to dip. A study of time residuals along two profiles across the median valley indicates the presence of a 2 to 3 km wide low velocity zone (about 3·2 km s⁻¹) beneath the median valley floor. This zone extends over the upper 2·5 km of the crust and is believed to represent a zone of intrusion through which magma passes on its way to the sea floor.     

Heat Flow and Geothermal Potential in the South-Central United States

Geothermal exploration is typically limited to high-grade hydrothermal reservoirs that are usually found in the western United States, yet large areas with subsurface temperatures above 150°C at economic drilling depths can be found east of the Rocky Mountains. The object of this paper is to present new heat flow data and to evaluate the geothermal potential of Texas and adjacent areas. The new data show that, west of the Ouachita Thrust Belt, the heat flow values are lower than east of the fault zone. Basement heat flow values for the Palo Duro and Fort Worth Basins are below 50 mW/m² while, in the frontal zone of the belt, they can exceed 60 mW/m². Further east, along the Balcones fault system the heat flow is in general higher than 55 mW/m². The eastern most heat flow sites are in Louisiana and they show very high heat flow (over 80 mW/m²), which is associated with the apparently highly radioactive basement of the Sabine uplift. The geothermal resource in this area is large and diverse, and can be divided in high grade (temperature above 150°C) convective systems, conductive based enhanced geothermal systems and geothermal/geopressured systems. One of the most attractive areas east of the cordillera extends from eastern Texas across Louisiana and Arkansas to western Mississippi. Here temperatures reach exploitation range at depths below 4 km, and tapping such a resource from shut in hydrocarbon fields is relatively easy. The initial costs of the development can be greatly reduced if existing hydrocarbon infrastructure is used, and therefore using shut-in hydrocarbon fields for geothermal purposes should not be neglected.     

Crustal structure of the Faeroe–Shetiand Channel

Summary. The deep structure of the Faeroe–Shetland Channel has been investigated as part of the North Atlantic Seismic Project. Shot lines were fired along and across the axis of the Channel, with recording stations both at sea and on adjacent land areas. At 61°N, 1.7 km of Tertiary sediments overlies a 3.9–4.5 km s^-1 basement interpreted as the top of early Tertiary volcanics. A main 6.0–6.6 km s^-1 crustal refractor interpreted as old oceanic crust occurs at about 9 km depth. The Moho (8.0 ° 0.2 km s^-1) is at about 15–17 km depth. There is evidence that P _n may be anisotropic beneath the Faeroe–Shetland Channel. Arrivals recorded at land stations show characteristics best explained by scattering at an intervening boundary which may be the continent–ocean crustal contact or the edge of the volcanics.
The Moho delay times at the shot points, determined by time-term analysis, show considerable variation along the axis of the Channel. They correlate with the basement topography, and the greatest delays occur over the buried extension of the Faeroe Ridge at about 60° 15'N, where they are nearly 1 s more than the delays at 61°N after correction for the sediments. The large delays are attributed to thickening of the early Tertiary volcanic layer with isostatic downsagging of the underlying crust and uppermost mantle in response to the load, rather than to thickening of the main crustal ayer.
The new evidence is consistent with deeply buried oceanic crust beneath the Faeroe–Shetland Channel, forming a northern extension of Rockall Trough. The seabed morphology has been grossly modified by the thick and laterally variable pile of early Tertiary volcanic rocks which swamped the region, accounting for the anomalous shallow bathymetry, the transverse ridges and the present narrowness of the Channel.     

Crustal and upper mantle structure of the northwestern North Island, New Zealand, from seismic refraction data

The crustal and upper mantle structure of the northwestern North Island of New Zealand is derived from the results of a seismic refraction experiment; shots were fired at the ends and middle of a 575 km-long line extending from Lake Taupo to Cape Reinga. The principal finding from the experiment is that the crust is 25 ± 2 km thick, and is underlain by what is interpreted to be an upper mantle of seismic velocity 7.6 ± 0.1 km s⁻¹, that increases to 7.9 km s⁻¹ at a depth of about 45 km. Crustal seismic velocities vary between 5.3 and 6.36 km s⁻¹ with an average value of 6.04 km s⁻¹. There are close geophysical and geological similarities between the north-western North Island of New Zealand and the Basin and Range province of the western United States. In particular, the conditions of low upper-mantle seismic velocities, thin crust with respect to surface elevation, and high heat-flow (70–100 mW m⁻²) observed in these two areas can be ascribed to their respective positions behind an active convergent margin for about the past 20 Myr.     

The Residual Water Flow through the Solent, South England

Summary. A reduced equation of motion is used to compute the residual velocity and the residual transport through the West Solent from the water levels recorded over an eight-month period at tide gauges at either end of the channel. A coefficient of bottom friction of 5·0 × 10⁻³ is assumed. There was a spring-neap variation and a significant correlation of fluctuations in the residual velocity with meteorological conditions. Westward residual velocities occurred at spring tides with low barometric pressure and south-westerly winds. Eastward residual velocities occurred at neap tides with high pressure and north-easterly winds. Because of the progressive nature of the tidal wave the long term residual transport appeared to be towards the west and the flushing time for the Solent system was long for considerable periods. The maximum velocities experienced during a tidal cycle half way along the channel are towards the west with a probability of values exceeding 160 cm s⁻¹ for 10 min in 5·4 days.     

Anelasticity of the Crust and Upper Mantle of South America From the Inversion of Observed Surface Wave Attenuation

Fundamental-mode Rayleigh and Love waves generated by several earthquakes situated along great-circle paths between pairs of seismograph stations have been analysed to obtain coefficients of attenuation, group velocities, phase velocities, and specific quality factors in the period range 18–80s in two regions of the South American continent. One set of paths crosses the shield region which lies on the eastern coast and another set traverses the mountainous region inland. the average attenuation coefficient values are clearly higher in the tectonically active western region throughout the entire period range than in the eastern or shield region.
Inversion of the attenuation data yielded shear wave internal friction ( Q ^-1_β) models as a function of depth in the crust and upper mantle in both regions. A low- Q zone below the lithosphere is prominent in both regions. the results show that substantial variations of Q _β occur in the two regions of South America. the Q_β values were found to be inversely related to the heat flow values or to the temperature.     

The crustal structure of the eastern Marmara region, Turkey by teleseismic receiver functions

We study the crustal structure of eastern Marmara region by applying the receiver function method to the data obtained from the 11 broad-band stations that have been in operation since the 1999 İzmit earthquake. The stacked single-event receiver functions were modelled by an inversion algorithm based on a five-layered crustal velocity model to reveal the first-order shear-velocity discontinuities with a minimum degree of trade-off. We observe crustal thickening from west (29–32 km) to east (34–35 km) along the North Anatolian Fault Zone (NAFZ), but we observe no obvious crustal thickness variation from north to south while crossing the NAFZ. The crust is thinnest beneath station TER (29 km), located near the Black Sea coast in the west and thickest beneath station TAR (35 km), located inland in the southeast. The average crustal thickness and S -wave velocity for the whole regions are  31 ± 2  km and  3.64 ± 0.15 km s⁻¹  , respectively. The eastern Marmara region with its average crustal thickness, high heat flow value (101 ± 11 mW m⁻²) and with its remarkable extensional features seems to have a Basin and Range type characteristics, but the higher average shear velocities (∼3.64 km s⁻¹) and crustal thickening from 29 to 35 km towards the easternmost stations indicate that the crustal structure shows a transitional tectonic regime. Therefore, we conclude that the eastern Marmara region seems to be a transition zone between the Marmara Sea extensional domain and the continental Anatolian inland region.     

North Atlantic Water in the Barents Sea Opening, 1997 to 1999

North Atlantic Water (NAW) is an important source of heat and salt to the Nordic seas and the Arctic Ocean. To measure the transport and variability of one branch of NAW entering the Arctic, a transect across the entrance to the Barents Sea was occupied 13 times between July 1997 and November 1999, and hydrography and currents were measured. There is large variability between the cruises, but the mean currents and the hydrography show that the main inflow takes place in Bjørnøyrenna, with a transport of 1.6 Sv of NAW into the Barents Sea. Combining the flow field with measurements of temperature and salinity, this results in mean heat and salt transports by NAW into the Barents Sea of 3.9×10¹³ W and 5.7×10⁷ kg s⁻¹, respectively. The NAW core increased in temperature and salinity by 0.7 °C yr⁻¹ and 0.04 yr⁻¹, respectively, over the observation period. Variations in the transports of heat and salt are, however, dominated by the flow field, which did not exhibit a significant change.     

Uplift and heat flow following the injection of magmas into the lithosphere

Summary. The thermal effect of a rapid injection of hot magmas into the lower part of the lithosphere is modelled as an increase in heat production through the invaded region. The change in surface heat flow and the uplift resulting from the thermal expansion are determined in three-dimensional axially symmetric geometry: they are expressed as the space time convolutions of a Green's function with the anomalous heat production.
The anomalies with shorter wavelength (compared to the lithospheric thickness) are attenuated. This filtering affects the surface uplift more than the heat flow anomaly; the attenuation effect is larger when only the lower part of the lithosphere is invaded.
The uplift time constant is of the same order as the heat conduction time if the lower lithosphere is invaded by magmas at a moderate rate (i.e. the rate of injection does not exceed the equivalent of 0.1 per cent of the lithospheric volume in 10⁶yr). Fifty per cent of the total uplift takes place in about 80 × 10⁶yr for a lithosphere 100 km thick. The uplift is slightly faster when the whole lithosphere is invaded. The heat flow anomaly is delayed when the lower part of the lithosphere is invaded.
The spatial extent and the timing of the uplift and heat flow anomalies are critical in determining the mechanism's feasibility. Magma injections explain rapid uplifts [> 100 m (10⁶ yr)⁻¹] only if the magma is supplied at a very high rate (i.e. at least 10 per cent of the lithosphere volume per 10⁶yr). It is a feasible mechanism for uplifts that occur over longer periods of time (≊ 30 × 10⁶yr) such as those that seem to have occurred when the African plate came to rest with respect to the mantle.