Differences between Mean Sea Levels for the Pacific,Atlantic and Indian Oceans From Topex/Poseidon Altimetry

Geopotential values W of the mean equipotential surfaces representing the mean ocean topography were computed on the basis of four years (1993 - 1996) TOPEX/POSEIDON altimeter data: W = 62 636 854.10m ² s ^–2 for the Pacific (P), W = 62 636 858.20m ² s ^–2 for the Atlantic (A), W = 62 636 856.28m ²s^–2 for the Indian (I) Oceans. The corresponding mean separations between the ocean levels were obtained as follows: A – P = – 42 cm, I– P = – 22 cm, I – A = 20 cm, the rms errors came out at about 0.3 cm. No sea surface topography model was used in the solution.     

Mantle heat flow and geotherms for major tectonic units in central Europe

The results of seismic measurements along the deep seismic sounding profile VII and terrestrial heat flow measurements used for construction of heat generation models for the crust in the Paleozoic Platform region, the Sudetic Mountains (Variscan Internides) and the European Precambrian Platform show considerable differences in mantle heat flow and temperatures. At the base of the crust variations from 440–510°C in the models of Precambrian Platform to 700–820°C for the Paleozoic Platform and the Variscan Internides (Sudets) are found. These differences are associated with considerable mantle heat flow variations.The calculated models show mantle heat flow of about 8.4–12.6 mW m^–2 for the Precambrian Platform and 31 mW m^–2 to 40.2 mW m^–2 for Paleozoic orogenic areas. The heat flow contribution originating from crustal radioactivity is almost the same for the different tectonic units (from 33.5 mW m^–2 to 37.6 mW m^–2). Considerable physical differences in the lower crust and upper mantle between the Precambrian Platform and the adjacent areas, produced by lateral temperature variations, could be expected. On the basis of carbon ratio data it can be concluded that the Carboniferous paleogeothermal gradient was much lower in the Precambrian Platform area than in the Paleozoic Platform region.     

The oceans as source of cloud-forming nuclei

Summary The number and size of salt particles produced by the bursting of air bubbles in sea water has been measured. Bubbles of diameters varying between 1/2 and 2 mm each produced about 300 nuclei the sizes of which, under the electron microscope, were mainly between 0.1 and 0.5 diameter. They appeared to consist mainly of sodium chloride, the smallest ones containing only 10^–15 g of salt.These results, together with measurements of the size distribution of salt nuclei collected over the oceans in areas of spray formation, indicate that the total concentrations of salt nuclei over the oceans in winds of up to 15 m sec^–1 probably do not exceed 100 cm^–3. The corresponding rate of production of salt nuclei at the sea surface is estimated to be 1000 cm^–2 sec^–1. It is therefore inferred that sea spray contributes perhaps only one-fifth of the nuclei involved in cloud formation, the majority being the products of combustion, either natural or man-made.     

On the heat flux related to stretching in the NW-Mediterranean continental margins

Summary The surface thermal flux of the continental margins of the northwestern Mediterranean Sea is interpreted on the basis of a 1-D instantaneous pure shear stretching model of the lithosphere in terms of three components: the background heat flowing out from the asthenosphere (38 mW m^–2), the transient contribution depending on the rift age and extension amount (35 mW m^–2 at the most), and the contribution due to the radiogenic elements of the lithosphere. The radiogenic component is estimated at the continental margins of the Ligurian-Provençal basin and Valencia trough, and in the surrounding mainland areas by means of available data of surface heat generation from Variscan Corsica, Maures-Estérel and the Central Massif along with a geophysical-petrological relationship between heat production and seismic velocity. The lithosphere radiogenic heat contribution q_l decreases with the thinning factor according to the exponential law: q_l() = a exp(-b), in which factor b is greater for that part of the lithosphere below the uppermost 10 km. Considering also the heat generated by radioactive isotopes in sediments, the stable Variscan lithosphere produces an average thermal flux of 30 mW m^–2 which decreases by about one half where the lithosphere is thinned by one third. Although the surface heat generation is 2·1 – 3·3 µW m^–3 in the Maures-Estérel massif — excepting small outcrops of dioritic rocks with lower heat production — and 1·8 µW m^–3 for most of Corsica, the radiogenic heating within the lithosphere for such areas is nearly the same and does not explain the higher heat flux of the Corsica margin. This asymmetric thermal pattern with surface heat flux which is 10 – 15 mW m^–2 higher than predictions is probably of upper mantle origin, or can be ascribed to penetrative magmatism.     

Vertical velocities associated with gravity waves measured in the mesosphere and lower thermosphere with the EISCAT VHF radar

The EISCAT VHF radar (69.4°N, 19.1°E) has been used to record vertical winds at mesopause heights on a total of 31 days between June 1990 and January 1993. The data reveal a motion field dominated by quasi-monochromatic gravity waves with representative apparent periods of 30–40 min, amplitudes of up to 2.5 m s^–1 and large vertical wavelength. In some instances waves appear to be ducted. Vertical profiles of the vertical-velocity variance display a variety of forms, with little indication of systematic wave growth with height. Daily mean variance profiles evaluated for consecutive days of recording show that the general shape of the variance profiles persists over several days. The mean variance evaluated over a 10 km height range has values from 1.2 m²s^–2 to 6.5 m²s^–2 and suggests a semi-annual seasonal cycle with equinoctial minima and solsticial maxima. The mean vertical wavenumber spectrum evaluated at heights up to 86 km has a slope (spectral index) of -1.36 ± 0.2, consistent with observations at lower heights but disagreeing with the predictions of a number of saturation theories advanced to explain gravity-wave spectra. The spectral slopes evaluated for individual days have a range of values, and steeper slopes are observed in summer than in winter. The spectra also appear to be generally steeper on days with lower mean vertical-velocity variance.     

Comparison of model electron densities and temperatures with Millstone Hill observations during undisturbed periods and the geomagnetic storms of 16–23 March and 6–12 April 1990

Measurements of F-region electron density and temperature at Millstone Hill are compared with results from the IZMIRAN time-dependent mathematical model of the Earths ionosphere and plasmasphere during the periods 16–23 March and 6–12 April 1990. Each of these two periods included geomagnetically quiet intervals followed by major storms. Satisfactory agreement between the model and the data is obtained during the quiet intervals, provided that the recombination rate of O⁺(⁴S) ions was decreased by a factor of 1.5 at all altitudes during the nighttime periods 17–18 March, 19–20 March, 6–8 April and 8–9 April in order to increase the NmF2 at night better to match observations. Good model/data agreement is also obtained during the storm periods when vibrationally excited N₂ brings about factor-of-2-4 reductions in daytime NmF2. Model calculations are carried out using different expressions for the O⁺ – O collision frequency for momentum transfer, and the best agreement between the electron-density measurements and the model results is obtained when the CEDAR interim standard formula for the O⁺ – O collision frequency is used. Deviations from the Boltzmann distribution for the first five vibrational levels of NI were calculated. The calculated distribution is highly non-Boltzmann at vibrational levels j > 2, and the Boltzmann distribution assumption results in the increase of 10–30% in calculated NmF2 during the storm-time periods. During the March storm at solar maximum the model results obtained using the EUVAC solar flux model agree a little better with the observations in comparison with the EUV94 solar flux model. For the April storm period of moderate solar activity the EUV94X model results agree better with the observations in comparison to the EUVAC model.     

Geodynamic implications of deep mantle upwelling in the source of Tertiary volcanics from the Veneto region (South-Eastern Alps)

Major and trace element and Sr–Nd–Hf–Pb isotopic data for the most primitive Tertiary lavas from the Veneto region (South-Eastern Alps, Italy) show the typical features of HIMU hotspot volcanism, variably diluted by a depleted asthenospheric mantle component (⁸⁷Sr/⁸⁶Sr_i=0.70306–0.70378; _Ndi=+3.9 to +6.8; _Hfi=+6.4 to +8.1, ²⁰⁶Pb/²⁰⁴Pb_i=18.786–19.574). P-wave seismic tomography of the mantle below the Veneto region shows the presence of low-velocity anomalies at depth, which is consistent with possible upwellings of plume material. Between the depths of 100–250 km the velocity anomalies are approximately 2–2.5% slower than average, implying a temperature excess of about 220–280 K, in agreement with estimates for other mantle plumes in the world. In this context, the Veneto volcanics may represent the shallow expression of a mantle upflow. The presence of a HIMU-DM component in a collision environment has significant geodynamic implications. Slab detachment and ensuing rise of deep mantle material into the lithospheric gap is proposed to be a viable mechanism of hotspot magmatism in a subduction zone setting.     

Natural recolonization of a productive tropical pond: Day to day variations in the photosynthetic parameters

Chlorophyll pigments (CHL), primary productivity (PP) and particulate nitrogen (N_p) in relation to several environmental factors were monitored during planktonic colonization of an aquaculture pond (Layo, Côte d'Ivoire). How interactions between the organisms are established in an initially azoic environment were investigated. From March, 15 (D1) to March, 31 (D16), the system transformation went through three stages. First, a precolonization by heterotrophic microbial community from D1 to D2 (N_p < 1 m maximum at D2: 243 mg m^–2; CHL around 0). Then, a pioneer microalgal community developped from D3 to D7 (maximum CHL on D6: 19 mg m^–2; PP: 1.0 g C m^–2 d^–1) with a significant contribution of picoplankton (CHL and PP < 3 m: 33 and 23% of the total, respectively). Finally, a second microalgal colonization was noticed from D9 to D12 (maximum CHL: 55 mg m^–2, PP: 2.8 g C m^–2 d^–1), largely dominated by nanoplankton (CHL and PP > 3 m: 95 and 99% of the total, respectively). Overall, photosynthetic activity appeared to be closely linked to algal biomass. The study of autotrophic biomass and activity in different size classes in relation to the other parameters allowed us to precise the origin of the biomass fluctuations. The first bloom appeared to be controlled by selective grazing on small algae. The second algal development ended when N requirement represented at least 69% of N supply (in the N — NH₄ form). This control was enhanced by the appearance of rotifers, leading to a more complex equilibrium.     

A study of the upward continuation of gravity data from a plane surface

Summary A summation method of upward continuation of gravity data has been considered under the assumption that observations are available at regular intervals. The upward continued value has been obtained as the sum of products of individual gravity values and corresponding theoretical coefficients. Besides the usual parameter involving horizontal and vertical distances, the theoretical coefficients have been generalized to be dependent also on i) the order of a low order polynomial assumed to represent the gravity variation around a grid point and ii) the weights assigned to the gravity values at the nearest four grid points used for least-squares determination of the polynomial. Since the observations in practical cases are available over a finite area only, the effect of truncation of the area of the integration has also been discussed separately. The method has been programmed and tested on a three-dimensional model, whose true gravity effects were computed at various levels over a finite area. Upward continued values have been computed under various assumptions about the gravity field in the outside region. Comparisons of these results with the true values indicate that the truncation effect becomes increasingly important for larger values of the ratio of elevation to grid separation and/or when the gravity field is not a local one. It has also been found that the greater is the above ratio, the less important is the effect of weights on the theoretical coefficients and practically vanishes(<10 ^–4 ) when the ratio is greater than5.0.     

Terrestrial heat flow in Hungary

Summary The geothermal gradient in the Carpathian Basin lies between 40–70 C/km. According to careful measurements in shafts the value of terrestrial heat flow in the southern part of Hungary is (2.055–3.066)·10^–6 cal/cm² sec. These measurements are believed the first ever attempted in continental Europe. Systematic heat flow measurement are being extended to other part of this country.     

Fracture surface energy of olivine

A relatively simple indentation technique for the rapid measurement of fracture surface energy, , of small samples is described. The reliability of this technique is assessed by testing soda-lime glass for which there are good independent fracture mechanics determinations of fracture surface energy. The indentation technique gives a value for of 4.33 J m^–2 which compares favourably with the accepted value of 3.8 J m^–2. Fracture surface energies of the {010} and {001} cleavage planes of single crystal olivine (modal composition Fo₈₈Fa₁₂) are then determined and compared with theoretical estimates of the thermodynamic surface energy, , calculated from atomistic parameters ( is equal to in the absence of dissipative processes during crack extension). The experimental values for _{010} and _{001} are respectively 0.98 J m^–2 and 1.26 J m^–2. The calculated values of _{010} and _{001} are respectively in the range from 0.37 J m^–2 to 8.63 J m^–2 and 12.06 J m^–2. The particular advantages of the indentation technique for the study of the fracture surface energies of geological materials are outlined.     

Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

Extensions of Leovy's linear mesospheric circulation model

Linear mesospheric circulation models can quite realistically simulate mid-latitude mesospheric circulation and structure. In doing so they provide estimates for dissipation (Rayleigh friction 2×10^–6s^–1) and vertical mixing (30 m²s^–1) appropriate near the mesopause.     

Determination of characteristics of steam reservoirs by Radon-222 measurements in geothermal fluids

The authors conducted a Rn²²² survey in wells of the Larderello geothermal field (Italy) and observed considerable variations in concentrations. Simple models show that flow-rate plays an important part in the Rn²²² content of each well, as it directly affects the fluid transit time in the reservoirs. Rn²²² has been sampled from two wells of the Serrazzano area during flow-rate drawdown tests. The apparent volume of the steam reservoir of each of these two wells has been estimated from the Rn²²² concentration versus flow-rate curves.List of symbols Q Flow-rate (kg h^–1) - Decay constant of Rn²²² (=7.553×10^–3 h^–1) - Porosity of the reservoir (volume of fluid/volume of rock) - ₁ Density of the fluid in the reservoir (kg m^–3) - ₂ Density of the rock in the reservoir (kg m^–3) - M Stationary mass of fluid filling the reservoir (kg). - E Emanating power of the rock in the reservoir (nCi kg _rock ^–1 h^–1). - P Production rate of Rn²²² in the reservoir: number of atoms of Rn²²² (divided by 1.764×10⁷) transferred by the rock to the mass unit of fluid per unit time (nCi kg _fluid ^–1 h^–1). - N Specific concentration of Rn²²² in the fluid (nCi kg^–1) - Characteristic time of the steam reservoir at maximum flow-rate (=M/Q)     

The large volume multi-anvil press as a highP-T deformation apparatus

The rheological properties of mantle materials are being investigated up to pressures of 16 GPa and temperatures of 1600°C for times up to 24 h, using a new sample assembly for the 6–8 multi-anvil apparatus. Al₂O₃ pistons, together with a liquid confining medium, are used to generate deviatoric stress in the specimen. Strain rates are estimated by monitoring the relative displacement of the guide blocks of the multi-anvil apparatus, scaled to the total axial strain of the sample. The applied stress on the sample is estimated using grain size piezometry. Strain rates and flow stresses of approximately 10^–4 to 10^–6 s^–1 and 50 to 250 MPa respectively, are presently attainable.Preliminary results on San Carlos olivine single crystals, partially dynamically recrystallized to a grain size of 10 to 300 m, indicate that the effective viscosity of polycrystalline olivine is consistent with values obtained from olivine single crystal creep laws. Assuming a dislocation creep mechanism (n3.5) with (010)[001] as the dominant slip system, the data are best fit using a creep activation volume of 5 to 10×10^–6 m³ mol^–1.     

Multivariate statistical analysis of geochemical data as indicative of the hydrogeochemical evolution of groundwater in a sedimentary rock aquifer system

The study of groundwater hydrogeochemistry of the Paleozoic Basses-Laurentides sedimentary rock aquifer system in Québec produced a large geochemical dataset. Groundwater samples were collected at 153 sites over a 1500 km² study area and analyzed for major and minor ions. The large number of data can lead to difficulties in the integration, interpretation and representation of the results. Two multivariate statistical methods, hierarchical cluster analysis (HCA) and principal components analysis (PCA), were applied to a subgroup of the dataset to evaluate their usefulness to classify the groundwater samples, and to identify geochemical processes controlling groundwater geochemistry. This subgroup consisted of 144 samples and 14 parameters (Ca²⁺, Mg²⁺, Na⁺, K⁺, , Cl⁻, , Fe²⁺, Mn²⁺, Br⁻, Sr²⁺, F⁻, Ba²⁺, HS⁻). Seven geochemically distinct clusters, C1–C7, resulted from the HCA. Samples from clusters C3, C4, C6 and C7 are mostly located in preferential recharge areas. The majority of these samples have Ca–Mg–HCO₃ recharge groundwater (C3, C6, C7) and Na–HCO₃ evolved groundwater (C4). Samples from the other three clusters (C1, C2, C5) are characteristic of an aquifer system under confined conditions. The majority of these samples have Na–HCO₃ evolved groundwater (C1, C5) and Na–Cl ancient groundwater that exhibits elevated concentrations in Br⁻ (C2). In addition to recognizing the importance of hydrogeological conditions on groundwater geochemistry, the distribution of clusters also showed the importance of the geological formations on minor and trace elements, such as Fe²⁺, Mn²⁺, Sr²⁺, F⁻ and Ba²⁺. The first five components of the PCA account for 78.3% of the total variance in the dataset. Component 1 is defined by highly positive loadings in Na⁺, Cl⁻ and Br⁻ and is related to groundwater mixing with Champlain Sea water and solute diffusion from the marine clay aquitard. The high positive loadings in Ca²⁺ and Mg²⁺ of component 2 suggest the importance of dissolution of carbonate rocks in this aquifer system. From their characteristic loadings, the first two components are defined as the “salinity” and “hardness” components, respectively. Components 3–5 are related to more local and geological effects. The integration of the HCA and the PCA, with conventional classification of groundwater types, as well as with the hydrogeological and geological contexts, allowed the division of the region into four main geochemical areas, providing an improved regional picture of the aquifer system dynamics and hydrogeochemical evolution of groundwater. The following factors were recognized as influencing the evolution of groundwater identified in every geochemical area: (1) geological characteristics including sedimentary rock type and till mineralogy; (2) hydrogeological characteristics represented by the level of confinement and the hydraulic gradient; and (3) the geological history including the latest glaciation and the Champlain Sea invasion. With its integrated approach, this hydrogeochemical study contributes to the characterization and understanding of complex groundwater flow systems, and provides an example of the long-term geochemical evolution of hydrogeological systems after a major perturbation, in this case seawater invasion.     

Anatomy of 1986 Augustine volcano eruptions as recorded by multispectral image processing of digital AVHRR weather satellite data

Eighteen digital AVHRR (advanced very high resolution radiometer) data sets from NOAA-6 and NOAA-9 polar-orbiting satellites recorded between 27 March and 7 April 1986 depict the eruptive activity of Augustine volcano, located 280 km SW of Anchorage, Alaska. The synoptic view (resolution of either 1.1 or 4.4 km), frequent coverage (often twice a day), and multispectral coverage (five bands: 0.58–0.68; 0.72–1.1; 3.55–3.93; 10.5–11.3; and 11.5–12.5 m) makes the AVHRR broadly applicable to analyzing explosive eruption clouds. The small scale of the Augustine activity (column heights of 2–13 km and eruption rates of 2x10⁶–8x10⁷ metric tonnes/day) facilitated intensive multispectral study because the plumes generally covered areas within the 550x550 km area of one easily manipulated image field. Hourly ground weather data and twice-daily radiosonde measurements from stations surrounding the volcano plus numerous volcanological observations were made throughout the eruption, providing important ground truth with which to calibrate the satellite data. The total erupted volume is estimated to be at least 0.102 km³. The pattern of changing eruption rates determined by satellite observations generally correlate with more detailed estimates of explosion magnitudes. Multispectral processing techniques were used to distinguish eruption clouds from meteorological clouds. Variable weather during the Augustine eruption offered an opportunity to test various trial algorithms. A ratio between thermal IR channels four and five, served to delineate the ashbearing eruption plumes from ordinary clouds. Future work is needed to determine whether the successful multispectral discrimination is caused by wavelength-dependent variable emission of silicate ash or reflects a spectral role of sulfuric acid aerosol in the plume.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

Two-dimensional kinematic and rheological modeling of the 1912 pyroclastic flow,Katmai, Alaska

Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s^–1. The flow progressively decelerated and was travelling only 2–8 m s^–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s^–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 10³ P. If the flow was Newtonian, its viscosity was 4.2 × 10³ P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.