Rocket observation of atomic oxygen and night airglow: Measurement of concentration with an improved resonance fluorescence technique

An improved resonant fluorescence instrument for measuring atomic oxygen concentration was developed to avoid the Doppler effect and the aerodynamic shock effect due to the supersonic motion of a rocket. The shock effect is reduced by adopting a sharp wedge-shaped housing and by scanning of the detector field of view to change the distance between the scattering volume and the surface of the housing. The scanning enables us to determine absolute values of atomic oxygen concentration from relative variation of the scattered light signal due to the self-absorption. The instrument was calibrated in the laboratory, and the numerical simulation reproduced the calibration result. Using the instrument, the altitude profile of atomic oxygen concentration was observed by a rocket experiment at Uchinoura (31°N) on 28 January 1992. The data obtained from the rocket experiment were not perfectly free from the shock effect, but errors due to the effect were reduced by the data analysis procedure. The observed maximum concentration was 3.8× 10¹¹ cm⁻³ at altitudes around 94 km. The systematic error is estimated to be less than ±0.7×10¹¹ cm⁻³ and the relative random error is less than±0.07× 10¹¹ cm⁻³at the same altitudes. The altitude profile of the OI 557.7-nm airglow was also observed in the same rocket experiment. The maximum volume emission rate was found to be 150 photons cm⁻³ s⁻¹ at 94 km. The observed altitude profiles are compared with the MSIS model and other in situ observations.     

First observation of mesospheric wind shear as high as 330 m s−1 km−1

Mesospheric wind profiles with an altitude resolution of 25 m have been obtained by means of radar tracking of foil chaff clouds. Such experiments were performed during winter 1990 at Biscarrosse, France (44^°N, 1^°W). On one flight, a wind shear as high as 330 m s⁻¹ km⁻¹ at 87.4 km and a region of dynamical instability between 86 and 88 km was measured. This wind shear is believed to be the largest value ever measured in the mesosphere. The region of dynamical instability results from a superposition of two wave motions, and is found to link well with enhanced turbulence and small-scale wave activity.     

Simultaneous radio and optical observations of the mid-latitude atmospheric response to a major geomagnetic storm of 6–8 April 2000

Basic properties of the mid-latitude traveling ionospheric disturbances (TIDs) during the maximum phase of a major magnetic storm of 6–8 April 2000 are shown. Total electron content (TEC) variations were studied by using data from GPS receivers located in Russia and Central Asia. The nightglow response to this storm at mesopause and termospheric altitudes was also measured by optical instruments FENIX located at the observatory of the Institute of Solar-Terrestrial Physics (51.9°N,103.0°E), and MORTI located at the observatory of the Institute of Ionosphere (43.2°N, 77.0°E). Observations of the O (557.7 and 630.0 nm) emissions originating from atmospheric layers centered at altitudes of 90 and 250 km were carried out at Irkutsk and of the O₂(b¹∑_g⁺−X³∑_g⁻) (0-1) emission originating from an atmospheric layer centered at altitude of 94 km was carried out at Almaty. Our radio and optical measurement network observed a storm-induced solitary large-scale wave with duration of 1 h and a wave front width of no less than 5000 km, while it traveled equatorward with a velocity of 200 m/s from 62°N to 38°N geographic latitude. The TEC disturbance, basically displaying an electron content depression in the maximum of the F2 region, reveals a good correlation with growing nightglow emission, the temporal shift between the TEC and emission variation maxima being different for different altitudes. A comparison of the auroral oval parameters with dynamic spectra of TEC variations and optical 630 nm emissions in the frequency range 0.4–4 mHz (250–2500 s periods) showed that as the auroral oval expands into mid-latitudes, also does the region with a developed medium-sale and small-scale TEC structure.     

Recombination rates of small ions and their attachment to aerosol particles

It could be shown by measurements of the air conductivity and using a mean profile for the ionization rate that experimental and theoretical values of the recombination rate of small ions based on a three body recombination process (Thomson) are in very good agreement up to 20 km altitude. The divergency of the experimental and theoretical curves above 20 km can be interpreted by assuming that there exists in this altitude region a crossover from the three body recombination to a two body recombination process. The value of the recombination coefficient is about 4·10^–7 cm³ s^–1 in 25 km altitude, compared with 1.4·10^–6 cm³ s^–1 at ground level. Furtheron it was possible for the first time to get some experimental data of attachment coefficients up to 13 km from simultaneous measurements of the air conductivity and Aitken nuclei concentration. These values are in good agreement with those obtained by theoretical considerations.     

Rupture process of the 19 August 1992 Susamyr, Kyrgyzstan, earthquake

The Susamyr earthquake of August 19, 1992 in Kyrgyzstan is one of the largest events (M_s = 7.4, M_b = 6.8) of this century in this region of Central Asia. We used broadband and long period digital data from IRIS and GEOSCOPE networks to investigate the source parameters, and their space-time distribution by modeling both body and surface waves. The seismic moment (M₀ = 6.8 × 10¹⁹ N m) and the focal mechanism were determined from frequency-time analysis (FTAN) of the fundamental mode of long period surface waves (100–250 s). Then, the second order integral moments of the moment-rate release were estimated from the amplitude spectra of intermediate period surface waves(40–70 s). From these moments we determined a source duration of 11–13 s, major and minor axes of the source of 30 km and 10–22 km, respectively; and an instant centroid velocity of 1.2 km/s. Finally, we performed a waveform inversion of P and SH waves at periods from 5–60 s. We found a source duration of 18–20 s, longer than the integral estimate from surface wave amplitudes. All the other focal parameters inverted from body waves are similar to those obtained by surface waves ( = 87° ± 6°, = 49° ± 6°, = 105° ± 3°, h = 14 ± 2 km, and M₀ = 5.8 ± 0.7 × 10¹⁹ N m). The initial rupture of this shallow earthquake was located at the south-west border of Susamyr depression in the western part of northern Tien Shan. A finite source analysis along the strike suggests a westward propagation of the rupture. The main shock of this event was preceded 2 s earlier by small foreshock. The main event was almost immediately followed by a very strong series of aftershocks. Our surface and body wave inversion results agree with the general seismotectonic features of the region.     

Source Parameters of the ML 3.8 Earthquake on January 20, 2000 near Meckenheim, Germany

In the early morning hours on January 20,2000 at 03:03:17 (UTC) an M_L = 3.8earthquake occurred east of the city ofMeckenheim. It awakened many people in theMeckenheim-Bonn-Siegburg area. At least 200people called the police or fire brigade inBonn. The earthquake was felt as far asLimburg (Westerwald) and east Belgium. Themaximum intensity (EMS) in theMeckenheim-Adendorf area was 5. Theinstrumental epicenter was located at50.60^° N latitude and 7.08^° Elongitude. Source depth was 9.5 km. Localmagnitude M_L was 3.8 ± 0.3 andseismic moment was 1.86*10¹⁴ Nm,corresponding to a momentmagnitude M_W of 3.4± 0.2.Following the Brune model, the sourceradius was 0.45 km with an averagedislocation of 1.0 cm and a 0.8 MPa stressdrop. The source mechanism from a firstmotion fault plane solution shows a normalfault mechanism with a significant strikeslip contribution. Trend and plunge of themajor and minor principal axes P/T are168^°/54^° and275^°/12^°, respectively. Aninversion of the moment tensor fromamplitudes of direct waves reveals a doublecouple mechanism with a smaller strike slipcontribution than the first motion faultplane solution, but also shows asignificant 16% non-double couplecomponent of the total moment. Theintensity epicenter was determined by agrid search algorithm and was found to be14 km east from the instrumental epicenter.The macroseismic magnitude was 3.6.     

An anomalous subauroral red arc on 4 August, 1972: comparison of ISIS-2 satellite data with numerical calculations

This study compares the Isis II satellite measurements of the electron density and temperature, the integral airglow intensity and volume emission rate at 630 nm in the SAR arc region, observed at dusk on 4 August, 1972, in the Southern Hemisphere, during the main phase of the geomagnetic storm. The model results were obtained using the time dependent one-dimensional mathematical model of the Earth’s ionosphere and plasmasphere (the IZMIRAN model). The major enhancement to the IZMIRAN model developed in this study to explain the two component 630 nm emission observed is the analytical yield spectrum approach to calculate the fluxes of precipitating electrons and the additional production rates of N⁺₂, O⁺₂, O⁺(⁴S), O⁺(²D), O⁻(²P), and O⁺(²P) ions, and O(¹D) in the SAR arc regions in the Northern and Southern Hemispheres. In order to bring the measured and modelled electron temperatures into agreement, the additional heating electron rate of 1.05 eV cm⁻³ s⁻¹ was added in the energy balance equation of electrons at altitudes above 5000 km during the main phase of the geomagnetic storm. This additional heating electron rate determines the thermally excited 630 nm emission observed. The IZMIRAN model calculates a 630 nm integral intensity above 350 km of 4.1 kR and a total 630 nm integral intensity of 8.1 kR, values which are slightly lower compared to the observed 4.7 kR and 10.6 kR. We conclude that the 630 nm emission observed can be explained considering both the soft energy electron excited component and the thermally excited component. It is found that the inclusion of N₂(v > 0) and O₂(v > 0) in the calculations of the O⁺(⁴S) loss rate improves the agreement between the calculated N_e and the data on 4 August, 1972. The N₂(v > 0) and O₂(v > 0) effects are enough to explain the electron density depression in the SAR arc F-region and above F2 peak altitude. Our calculations show that the increase in the O⁺ + N₂ rate factor due to the vibrationally excited nitrogen produces the 5–19% reductions in the calculated quiet daytime peak density and the 16–24% decrease in NmF2 in the SAR arc region. The increase in the O⁺ + N₂ loss rate due to vibrationally excited O₂ produces the 7–26% decrease in the calculated quiet daytime peak density and the 12–26% decrease in NmF2 in the SAR arc region. We evaluated the role of the electron cooling rates by low-lying electronic excitation of O₂(a¹δ_g) and O₂(b¹σ_g⁺), and rotational excitation of O₂, and found that the effect of these cooling rates on T_e can be considered negligible during the quiet and geomagnetic storm period 3–4 August, 1972. The energy exchange between electron and ion gases, the cooling rate in collisions of O(³P) with thermal electrons with excitation of O(¹D), and the electron cooling rates by vibrational excitation of O₂ and N₂ are the largest cooling rates above 200 km in the SAR arc region on 4 August, 1972. The enhanced IZMIRAN model calculates also number densities of N₂(B³π_g⁺), N₂(C³π_u), and N₂(A³σ_u⁺) at several vibrational levels, O(¹S), and the volume emission rate and integral intensity at 557.7 nm in the region between 120 and 1000 km. We found from the model that the integral integral intensity at 557.7 nm is much less than the integral intensity at 630 nm.     

Hydrothermal vent waters at 13°N on the East Pacific Rise: isotopic composition and gas concentration

Gas concentrations and isotopic compositions of water have been measured in hydrothermal waters from 13°N on the East Pacific Rise. In the most Mg-depleted samples ( 5 × 10⁻³ moles/kg) the gas concentrations are: 3–4.5 × 10⁻⁵ cm³ STP/kg helium, 0.62–1.24 cm³ STP/kg CH₄, 10.80–16.71 × 10⁻³ moles/kg CO₂. The samples contain large quantities (95–126 cm³/kg) of H₂ and some carbon monoxide (0.26–0.36 cm³/kg) which result from reaction with the titanium sampling bottles. δ¹³C in methane and CO₂ (−16.6 to −19.5 and −4.1 to −5.5 respectively) indicate temperatures between 475 and 550°C, whereas δ¹³C_CO is compatible with formation by reduction of CO₂ on Ti at 350°C close to the sampling temperature.³He/⁴He are very homogeneous at (7.5 ± 0.1)R_A(³He/⁴He = 1.0 × 10⁻⁵) and very similar to already published data as well as CH₄/³He ratios between 1.4 and 2.1 × 10⁶.¹⁸O and D in water show enrichments from 0.39 to 0.69‰ and from 0.62 to 1.49‰ respectively. These values correspond to W/R ratios of 0.4–7. The distinct¹⁸O enrichments indicate that the isotopic composition of the oceans is not completely buffered by the hydrothermal circulations. The³He-enthalpy relationship is discussed in terms of both hydrothermal heat flux and³He mantle flux.     

Heterogeneous loss of OH on NaCl and NH4 NO 3 at tropospheric temperatures

The uptake coefficients () for OH radicals on some dry salts of tropospheric interest (NaCl and NH₄NO₃) have been investigated as a function of temperature using the flow tube technique combined with an EPR spectrometer as a detection method. The temperature dependence of -values measured over the temperature range 245–340 K can be expressed in Arrhenius form: ^OH_NaCl=(1.2±0.7)×10⁻⁵exp[(1750±200)/T] and ^OH_NH4NO3=(1.4±0.5)×10⁻⁴exp[(1000±100)/T]. These Arrhenius expressions lead to very similar -values (\approx4×10⁻³) for both salts studied at 300 K. It is shown that the heterogeneous OH sinks on solids aerosol play a very minor role in tropospheric chemistry in comparison with the homogeneous sinks.     

Conductivity as a predictor of a total cations and salinity in Ethiopian lakes and rivers: revisiting earlier models

We used regression analyses of water samples from 18 lakes, nine rivers, and one spring in Ethiopia to (a) test the hypothesis that water bodies of relatively higher salinity (K₂₅>1000 μS cm⁻¹) have a different conductivity to salinity relationship than waters of lower salinity (K₂₅ < 1000 μS cm⁻¹), and (b) develop models to predict total cations and salinity from conductivity that can be used for Ethiopian waters and other African aquatic systems of similar chemical composition. We found no statistical difference in the bilogarithmic relationships (total cations vs. conductivity; salinity vs. conductivity) for waters of higher salinity (K₂₅ > 1000 μS cm⁻¹) and waters of lower salinity (K₂₅ < 1000 μS cm⁻¹). However, comparison among our models and models from the literature suggests that developing separate equations for low and high salinity water bodies has some merit. We believe that the equations developed in this study can be used for Ethiopian waters and other African waters within the range of conductivity in this study.     

The mean ozone distribution from several series of rocket soundings to 52 km at latitudes from 58°S to 64°N

The results of 21 rocket flights of Arcas optical ozonesondes have been combined to produce estimates of the mean ozone distribution and its variability which apply to a broad range of latitudes. The flights were launched at sites from near the equator to 58°S and to 64°N in the years from 1965 to 1971. The local-noon mean ozone densities in molecules/cubic centimeter are 7.0×10¹⁰ at 50 kw, 6.7×10¹¹ at 40 km, 3.1×10¹² at 30 km, and 3.1×10¹² at 20 km. The maximum density is 4.5×10¹² at 24 km. The range of observed densities is about ±30% of the mean value at 50 km, ±40% at 40 km, ±40% at 30 km and +200%, –66% at 20 km. The variabilities at the higher altitudes in this set of observations are much less than that indicated from previous measurements.     

Observations of the response time of high-latitude ionospheric convection to variations in the interplanetary magnetic field using EISCAT and IMP-8 data

We have combined ∼300 h of tristatic measurements of the field-perpendicular F region ionospheric flow measured overhead at Tromsø by the EISCAT UHF radar, with simultaneous IMP-8 measurements of the solar wind and interplanetary magnetic field (IMF) upstream of the Earth’s magnetosphere, in order to examine the response time of the ionospheric flow to changes in the north-south component of the IMF (B_z). In calculating the flow response delay, the time taken by field changes observed by the spacecraft to first effect the ionosphere has been carefully estimated and subtracted from the response time. Two analysis methods have been employed. In the first, the flow data were divided into 2 h-intervals of magnetic local time (MLT) and cross-correlated with the “half-wave rectifier” function V²B_s, where V is the solar wind speed, and B_s is equal to IMF B_z if the latter is negative, and is zero otherwise. Response delays, determined from the time lag of the peak value of the cross-correlation coefficient, were computed versus MLT for both the east-west and north-south components of flow. The combined data set suggests minimum delays at ∼1400 MLT, with increased response times on the nightside. For the 12-h sector centred on 1400 MLT, the weighted average response delay was found to be 1.3 ± 0.8 min, while for the 12-h sector centred on 0200 MLT the weighted average delay was found to increase to 8.8 ± 1.7 min. In the second method we first inspected the IMF data for sharp and enduring (at least ∼5 min) changes in polarity of the north-south component, and then examined concurrent EISCAT flow data to determine the onset time of the corresponding enhancement or decay of the flow. For the case in which the flow response was timed from whichever of the flow components responded first, minimum response delays were again found at ∼1400 MLT, with average delays of 4.8 ± 0.5 min for the 12-h sector centred on 1400 MLT, increasing to 9.2 ± 0.8 min on the nightside. The response delay is thus found to be reasonably small at all local times, but typically ∼6 min longer on the nightside compared with the dayside. In order to make an estimate of the ionospheric information propagation speed implied by these results, we have fitted a simple theoretical curve to the delay data which assumes that information concerning the excitation and decay of flow propagates with constant speed away from some point on the equatorward edge of the dayside open-closed field line boundary, taken to lie at 77° magnetic latitude. For the combined cross-correlation results the best-fit epicentre of information propagation was found to be at 1400 MLT, with an information propagation phase speed of 9.0 km s⁻¹. For the combined event analysis, the best-fit epicentre was also found to be located at 1400 MLT, with a phase speed of 6.8 km s⁻¹.     

Fabry-Perot interferometer measurements of neutral winds and F2 layer variations at the magnetic equator

This letter presents some night-time observations of neutral wind variations at F2 layer levels near the dip equator, measured by the Fabry-Perot interferometer set up in 1994 at Korhogo (Ivory Coast, geographic latitude 9.25°N, longitude 355°E, dip latitude -2.5°). Our instrument uses the 630 nm (O¹D) line to determine radial Doppler velocities of the oxygen atoms between 200 and 400 km altitude. First results for November 1994 to March 1995 reveal persistent eastward flows, and frequent intervals of southward winds of larger than 50 ms^–1 velocity. Compared with the simultaneous ionospheric patterns deduced from the three West African equatorial ionosondes at Korhogo, Ouagadougou (Burkina-Faso, dip latitude +1.5°) and Dakar (Sénégal, dip latitude + 5°), they illustrate various impacts of the thermospheric winds on F2 layer density: (1) on the mesoscale evolution (a few 10³ km and a few 100 minutes scales) and (2) on local fluctuations (hundreds of km and tens of minutes characteristic times). We report on these fluctuations and discuss the opportunity to improve the time-resolution of the Fabry-Perot interferometer at Korhogo.     

Crustal structure of the Peninsular India

Summary Earthquake parameters for the forty aftershocks of the main Koyna earthquake of 10 December, 1967, have been determined. Depths of the foci of the earthquakes have been found to vary between 2 to 17 km. The velocities for the phasesP _g ,P ^*,P _n have been observed to be 5.78±0.00, 6.58±0.04, 8.19±0.02 km/sec, and forS _g ,S ^*,S _n to be 3.42±0.00, 3.92±0.01 and 4.62±0.01 km/sec respectively. A two-layered crustal model has been interpreted for the Peninsular shield with the average thickness of the granitic layer as 20 km and that for the basaltic layer as 18.7 km. A plot of the epicenters suggests a NNE to SSW orientation of the fault.     

Determination of the Geopotential Scale Factor from TOPEX/POSEIDON Satellite Altimetry

The geopotential scale factor R _o = GM/W _o (the GM geocentric gravitational constant adopted) and/or geoidal potential Wo have been determined on the basis of the first year's (Oct 92 – Dec 93) ERS-1/TOPEX/POSEIDON altimeter data and of the POCM 4B sea surface topography model: R _o °=(6 363 672.58°±0.05) m, W _o °=(62 636 855.8°±0.05)m ² s ^–2 . The 2°–°3 cm uncertainty in the altimeter calibration limits the actual accuracy of the solution. Monitoring dW _o /dt has been projected.     

A fracture mechanics study of subcritical tensile cracking of quartz in wet environments

Load relaxation and cross-head displacement rate-change experiments have been used to establish log₁₀ stress intensity factor (K) versus log₁₀ crack velocity (v) diagrams for double torsion specimens, of synthetic quartz cracked on thea plane in liquid water and moist air.For crack propagation normal toz and normal tor at 20°C,K _Ic (the critical stress intensity factor) was found to be 0.852±0.045 MN·m^–3/2 and 1.002±0.048 MN·m^–3/2, respectively.Subcritical crack growth at velocities from 10^–3 m·s^–1 to 10^–9 m·s^–1 at temperatures from 20°C to 80°C is believed to be facilitated by chemical reaction between the siloxane bonds of the quartz and the water or water vapour of the environment (stress corrosion). The slopes, of isotherms in theK-v diagrams are dependent upon crystallographic orientation. The isotherms have a slope of 12±0.6 for cracking normal tor and 19.9±1.7 for cracking normal toz. The activation enthalpy for crack propagation in the former orientation in liquid water at temperatures from 20°C to 80°C is 52.5±3.8 kJ·mole^–1.A discussion is presented of the characteristics of theK-v diagrams for quartz.     

First in situ measurement of electric field fluctuations during strong spread F in the Indian zone

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip 14°N) along with other experiments, as a part of equatorial spread F (ESF) campaign, to study the nature of irregularities in electric field and electron density. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of vertical and horizontal electric field fluctuations are presented here. Scale sizes of electric field fluctuations were measured in the vertical direction only. Strong ESF irregularities were observed in three regions, viz., 160/190 km, 210/257 km and 290/330 km. Some of the valley region vertical electric field irregularities (at 165 km and 168 km), in the intermediate-scale size range, observed during this flight, show spectral peak at kilometer scales and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of a new type. Scale sizes of vertical electric field fluctuations showed a decrease with increasing altitude. The most prominent scales were of the order of a few kilometers around 170 km and a few hundred meters around 310 km. Spectra of intermediate-scale vertical electric field fluctuations below the base of the F region (210/257 km) showed a tendency to become slightly flatter (spectral index n = –2.1 ± 0.7) as compared to the valley region (n = –3.6 ± 0.8) and the region below the F peak (n = –2.8 ± 0.5). Correlation analysis of the electron density and vertical electric field fluctuations suggests the presence of a sheared flow of current in 160/330 km region.     

Average crustal structure of Sweden

Summary Records obtained at the permanent stations of the Swedish seismograph network from explosions carried out in Scandinavian waters in June 1969 are evaluated. The study includes determination of velocities for all crustal phases observed, furthermore of layer thicknesses, Poisson ratios and amplitude ratios. The purpose of the study is partly to provide a first approximation to the crustal structure in Sweden, partly to provide regional data for location of earthquakes and explosions in the area in the future. Average velocities (km/sec) are forPn 7.88±0.05,Pg1 6.25±0.08,Pg2 5.70,Sn 4.58±0.04,S ^* 3.70±0.04,Sg1 (Lg1) 3.58±0.03,Sg2 (Sg) 3.40±0.03,Rg 3.02±0.07. The average thickness is 12 km for the granitic layer, and 23 km for the basaltic layer, thus making the average crustal thickness equal to 35 km. Relative amplitudes plotted versus distance complete the dynamical side of the study and they are useful for identification of waves. A regional travel-time table is presented for the distance range 0°–10° with entries for each 0.1° and including all crustal phases read.     

Dissolved carbon dioxide in basaltic glasses: concentrations and speciation

Carbon dioxide dissolved in both synthetic Ca±Mg-bearing silicate glasses and natural basaltic glasses has been characterized using infrared spectroscopy. CO₂ is inferred to be dissolved in these glasses as distorted Ca or Mg carbonate ionic complexes that result in unique infrared absorption bands at 1515 cm⁻¹ and 1435 cm⁻¹. This speciation contrasts with the case of CO₂-bearing sodium aluminosilicate glasses, which contain both dissolved molecular CO₂ and dissolved Na-carbonate ionic-complexes. The difference in speciation in Ca±Mg-bearing melts may result in part from a higher activity of oxygens that react with CO₂ molecules to produce carbonate.Dissolved CO₂ contents of natural basaltic glasses can be determined from the intensities of the carbonate absorption bands at 1515 cm⁻¹ and 1435 cm⁻¹. The uncertainty of the method is estimated to be ± 15% of the amount present. The infrared technique is a powerful tool for the measurement of dissolved CO₂ contents in natural basaltic glasses since it is non-destructive, can be aimed at regions of glass a few tens of microns in size, and can discriminate between dissolved carbonate and carbon present as carbonate alteration, contained in fluid inclusions, or adsorbed on the glass.A set of submarine basaltic glasses dredged from a variety of locations contain 0–400 ppm dissolved CO₂, measured using the infrared technique. These concentrations are lower than most previous reports for similar basaltic glasses. No general relationship is observed between dissolved CO₂ content and depth of magmatic eruption, although some correlation might be present in restricted geographic locales.     

The global-average production rate ofBe

Precipitation collected in continuously open containers for about a year at seven sites around the United States was analyzed for¹⁰Be,⁹⁰Sr,²¹⁰Pb and²³⁸U. Based on these data and long-term precipitation,⁹⁰Sr and²¹⁰Pb delivery patterns, the stratospheric, tropospheric and recycled¹⁰Be components in the collections were estimated and the global¹⁰Be production rate was assessed. Single station production rate estimates range from 0.52 × 10⁶ atoms cm⁻² yr⁻¹ to 2.64 × 10⁶ atoms cm⁻² yr⁻¹. The mean value is 1.21 × 10⁶ atoms cm⁻² yr⁻¹ with a standard error of 0.26 × 10⁶ atoms cm⁻² yr⁻¹.