scispace - formally typeset
Search or ask a question

Showing papers in "Journal of Geophysical Research in 1990"


Journal ArticleDOI
TL;DR: In this paper, the authors introduced the importance of subducted oceanic crustal age on arc petrogenesis and demonstrated that Archean TTD crustal generation processes are also present in selected high-Al Phanerozoic TTD terranes.
Abstract: The petrogenesis of trondhjemite-tonalite-dacite (TTD) involves all major petrologic models in various tectonic settings. A specific subtype of TTD, high-Al type, is the one most commonly associated with Archean gneiss terranes. During the Archean, continental crust formation was operating at an elevated rate relative to the Phanerozoic, and the generation of high-Al TTD played an integral role in its nucleation and growth. High heat flow, rapid convection, and subduction of hotter, smaller plates were unique tectonic elements to the Archean which optimized conditions required for transformation of subducted oceanic crust into sial via partial melting. Anatexis of Archean mid-ocean ridge basalt (MORB) under eclogitic to garnet amphibolitic conditions produced weakly peraluminous to metaluminous high-Al TTD with low heavy rare earth elements (HREE), Y, Nb, K/Rb, and Rb/Sr and high La/Yb and Sr/Y. This study demonstrates that Archean TTD crustal generation processes are also present in selected high-Al Phanerozoic TTD terranes. The Cenozpic high-Al TTD suites are commonly found in tectonic settings which are thought to recreate the elevated Archean thermal gradients, i.e., at sites of young, hot oceanic plate subduction. These relationships imply a petrologic continuity of TTD generation through time. A fertile zone of melting is envisioned at 23–26 kbar (75–85 km) and 700–775°C, where wet partial melting of the subducting slab occurs concurrently with dehydration reactions. At this depth, the converting mantle wedge continuously feeds hot mantle material to the wedge-slab interface, creating strong temperature gradients, intraslab fluid migration, and slab melting. In summary, in modern arc terranes where young ( 30 Ma) oceanic crust is subducted, mantle-derived magmas are dominant, giving rise to basaltandesite-dacite-rhyolite (BADR) fractionation suites. This study introduces the importance of subducted oceanic crustal age on arc petrogenesis.

1,351 citations


Journal ArticleDOI
TL;DR: In this article, the reflectance spectra of minerals are studied as a function of spectral resolution in the range from 0.2 to 3.0 microns, and selected absorption bands were studied at resolving powers as high as 2240.
Abstract: The reflectance spectra of minerals are studied as a function of spectral resolution in the range from 0.2 to 3.0 microns. Selected absorption bands were studied at resolving powers as high as 2240. At resolving powers of approximately 1000, many OH-bearing minerals show diagnostic sharp absorptions at the resolution limit. At low resolution, some minerals may not be distinguishable, but as the resolution is increased, most can be easily identified. As the resolution is increased, many minerals show fine structure, particularly in the OH-stretching overtone region near 1.4 micron. The fine structure can enhance the ability to discriminate between minerals, and in some cases the fine structure can be used to determine elemental composition.

1,349 citations


Journal ArticleDOI
TL;DR: In this article, a systematic study of rare earth and other trace elements in discrete diopsides from residual abyssal peridotites sampled from 5000 km of ocean ridge demonstrates that they are the residues of variable degrees of melting in the garnet and spinel peridotsite fields.
Abstract: A systematic study of rare earth and other trace elements in discrete diopsides from residual abyssal peridotites sampled from 5000 km of ocean ridge demonstrates that they are the residues of variable degrees of melting in the garnet and spinel peridotite fields. Further, the data clearly demonstrate that the peridotites are the residues of near-fractional melting, not batch melting, and that typical abyssal basalt can evolve from aggregated fractional melts. Ion microprobe analyses of diopsides in abyssal peridotites from fracture zones along the America-Antarctica and Southwest Indian ridges reveal ubiquitous extreme fractionation of rare earth elements (REE) ([Ce/Yb]n = 0.002–0.05); depletion of Ti (300–1600 ppm), Zr (0.1–10 ppm), and Sr (0.1–10 ppm); and fractionation of Zr relative to Ti (Ti/Zr = 250–4000). Ti and Zr in diopsides decrease with decreasing modal cpx in the peridotites, and samples dredged near hotspots are more depleted in incompatible elements than those dredged away from hotspots, consistent with higher degrees upper mantle melting in the former. All studied samples exhibit marked negative anomalies in Ti and Zr relative to REE. Incompatible element concentrations in peridotite clinopyroxenes are well modeled by repeated melting and segregation in ≤0.1% increments to a total of 5–25% melting, a process very close to Rayleigh (fractional) melting; batch melting of a LREE-depleted source cannot account for the observed trace element concentrations in abyssal peridotites. The shapes of some REE patterns are consistent with variable degrees of melting initiated within the garnet stability field. Trace element concentrations in calculated integrated fractional liquids approximate the composition of primitive ocean floor basalts, consistent with postsegregation aggregation of small increment melts produced over a depth and melting interval.

1,204 citations


Journal ArticleDOI
TL;DR: In this article, the authors examined the available data, mainly topography, geoid, and heat flow, describing hotspots worldwide to constrain the mechanisms for swell uplift and to obtain fluxes and excess temperatures of mantle plumes.
Abstract: The available data, mainly topography, geoid, and heat flow, describing hotspots worldwide are examined to constrain the mechanisms for swell uplift and to obtain fluxes and excess temperatures of mantle plumes. Swell uplift is caused mainly by excess temperatures that move with the lithosphere plate and to a lesser extent hot asthenosphere near the hotspot. The volume, heat, and buoyancy fluxes of hotspots are computed from the cross-sectional areas of swells, the shapes of noses of swells, and, for on ridge hotspots, the amount of ascending material needed to supply the length of ridge axis which has abnormally high elevation and thick crust. The buoyancy fluxes range over a factor of 20 with Hawaii, 8.7 Mg s -1, the largest. The buoyancy flux for Iceland is 1.4 Mg s -1 which is similar to the flux of Cape Verde. The excess temperature of both on-ridge and off-ridge hotspots is around the 200oC value inferred from petrology but is not tightly constrained by geophysical considerations. This observation, the similarity of the fluxes of on-ridge and offridge plumes, and the tendency for hotspots to cross the ridge indicate that similar plumes are likely to cause both types of hotspots. The buoyancy fluxes of 37 hotspots are estimated; the global buoyancy flux is 50 Mg s -1, which is equivalent to a globally averaged surface heat flow of 4 mWm -2 from core sources and would cool the core at a rate of 50 o C b.y. -1. Based on a thermal model and the assumption that the likelihood of subduction is independent of age, most of the heat from hotspots is implaced in the lower lithosphere and later subducted. I.NTRODUCWION ridge plumes using Iceland as an example. The geometry of flow implied by the assumed existence of a low viscosity Linear seamount chains, such as the Hawaiian Islands, are asthenospheric channel is illustrated by this exercise. Then the frequently attributed to mantle plumes which ascend from deep methods for obtaining the flux of plumes on a rapidly moving in the Earth, perhaps the core-mantle boundary. The excessive plate are discussed with Hawaii as an example. These methods volcanism of on-ridge hotspots, such as Iceland, is also often involve determining the flux from the plume from the crossattributed to plumes. If on-ridge and midplate hotspots are sectional area of the swell and taking advantage of the kinematreally manifestations of the same phenomenon, one would ics of the interaction of asthenospheric flow away from the expect that the temperature and flux of the upwelling material plume and asthenospheric flow induced by the drag of the would be similar under both features. In particular, the core- lithospheric plate. The methods for extending this approach to mantle boundary is expected to be nearly isothermal so that the hotspots on slowly moving plates are then discussed which Cape temperature of plumes ascending from the basal boundary layer Verde as an example. An estimate of the global mass and heat should be the same globally provided that cooling by entrain- transfer by plumes is then obtained by applying the methods to ment of nearby material and thermal conduction are minor. 34 additional hotspots. The magnitude of this total estimated Finally, the global heat loss from plumes should imply a reason- flux is compatible with the heat flux expected from cooling the

1,087 citations


Journal ArticleDOI
TL;DR: The Regional Acid Deposition Model (RADM2) as discussed by the authors is a state-of-the-art gas phase chemical mechanism for modeling atmospheric chemistry on a regional scale.
Abstract: A state-of-the-art gas phase chemical mechanism for modeling atmospheric chemistry on a regional scale is presented. The second generation Regional Acid Deposition Model (RADM2) gas phase chemical mechanism, like its predecessor RADM1, is highly nonlinear, since predicted ozone, sulfate, nitric acid and hydrogen peroxide concentrations are complicated functions of NO{sub x} and nonmethane hydrocarbon concentrations. The RADM2 chemical mechanism is an upgrade of RADM1 in that (1) three classes of higher alkanes are used instead of one, (2) a more detailed treatment of aromatic chemistry is used, (3) the two higher alkene classes now represent internal and terminal alkenes, (4) ketones and dicarbonyl species are treated as classes distinct from aldehydes, (5) isoprene is now included as an explicit species, and (6) there is a more detailed treatment of peroxy radical-peroxy radical reactions. As a result of these improvements the RADM2 mechanism simulates the concentrations of peroxyacetyl nitrate, HNO3, and H{sub 2}O{sub 2} under a wide variety of environmental conditions. Comparisons of RADM2 mechanism with the RADM1 mechanism predictions and selected environmental chamber experimental results indicate that for typical atmospheric conditions, both mechanisms reliably predict O{sub 3}, sulfate and nitric acid concentrations. The RADM2 mechanism gives lower and presumably moremore » realistic predictions of H{sub 2}O{sub 2} because of its more detailed treatment of peroxy radical-peroxy radical reactions.« less

1,086 citations


Journal ArticleDOI
TL;DR: In this paper, the authors provided an intercomparison and interpretation of climate feedback processes in 19 atmospheric general circulation models and found that a roughly threefold variation in one measure of global climate sensitivity was found among the 19 models.
Abstract: The present study provides an intercomparison and interpretation of climate feedback processes in 19 atmospheric general circulation models. This intercomparison uses sea surface temperature change as a surrogate for climate change. The interpretation of cloud-climate interactions is given special attention. A roughly threefold variation in one measure of global climate sensitivity is found among the 19 models. The important conclusion is that most of this variation is attributable to differences in the models' depiction of cloud feedback, a result that emphasizes the need for improvements in the treatment of clouds in these models if they are ultimately to be used as reliable climate predictors. It is further emphazied that cloud feedback is the consequence of all interacting physical and dynamical processes in a general circulation model. The result of these processes is to produce changes in temperature, moisture distribution, and clouds which are integrated into the radiative response termed cloud feedback.

863 citations


Journal ArticleDOI
TL;DR: In this article, a least-squares program was developed to fit magnetic field data within a cloud, while estimating such cloud properties as its size, maximum field strength, and axis inclination.
Abstract: Interplanetary magnetic clouds emerge as a feature of the solar wind at 1 AU, exhibiting enhanced field strength and lower plasma temperature and density than the surrounding plasma. A least-squares program has been developed which fits magnetic field data within a cloud, while estimating such cloud properties as its size, maximum field strength, and axis inclination. The results obtained from a study of 12 clouds observed at 1 AU point to a probable cloud axis direction within 15 deg of the ecliptic plane and about 100 deg from the sun's direction, when projected into the ecliptic plane. A wide variety of orientations is observed; some extend to 80 deg from the ecliptic.

844 citations


Journal ArticleDOI
TL;DR: In this paper, the authors analyzed more than 270,000 ion moments and magnetic field measurements with respect to the occurrence rates and typical characteristics of high-speed ion flows with velocities in excess of 400 km/s.
Abstract: Using 8 months of tail data obtained with the AMPTE/IRM satellite, more than 270,000 ion moments and magnetic field measurements were analyzed with respect to the occurrence rates and typical characteristics of high-speed ion flows with velocities in excess of 400 km/s. The occurrence rates in the plasma sheet boundary layer, the outer central plasma sheet and the neutral sheet neighborhood have a 4:1:2 ratio for flows of 400-600 km/s. For flows in excess of 800 km/s, there is only a minimal chance to detect them in the outer central plasma sheet but equal chances in the two other regions. For high AE the chances to detect high-speed flows in the inner central plasma sheet are greater than to find them in the plasma sheet boundary layer. In the outer central plasma sheet the high-speed flow occurrence rate is small and independent of AE. In all three regions the largest occurrence rates are found near the midnight meridian at the largest radial distances accessible to IRM. High-speed flow occurrence rates and ion densities are anticorrelated. The high-speed slows are bursty with the majority of the flows lasting less than 10 s. The occurrence of the high-speed flows is strongly peakedmore » in the sunward direction. Virtually no tailward high-speed ion flow could be detected. About 60-70% of all high-speed flows near the neutral sheet have a dominant component perpendicular to the magnetic field and are associated with comparatively large northward and duskward magnetic field directions. At times, also appreciable duskward flow components appear. Overall, the results indicate that both the plasma sheet boundary layer and the inner central plasma sheet are important regions for the dynamics of the Earth's plasma sheet.« less

724 citations


Journal ArticleDOI
TL;DR: In this paper, the resolving power of teleseismic P waveforms for receiver structure was studied using a time domain waveform inversion scheme with a range of initial models to estimate the range of acceptable velocity structures.
Abstract: To study the resolving power of teleseismic P waveforms for receiver structure, the authors model synthetic waveforms using a time domain waveform inversion scheme beginning with a range of initial models to estimate the range of acceptable velocity structures. To speed up the waveform inversions, they implement Randall's (1989) efficient algorithms for calculating differential seismograms and include a smoothness constraint on all the resulting velocity models utilizing the jumping inversion technique of Shaw and Orcutt (1985). They present the results of more than 235 waveform inversions for one-dimensional velocity structures that indicate that the primary sensitivity of a receiver function is to high wavenumber velocity changes, and a depth-velocity product, not simply velocity. The range of slownesses in a typical receiver function study does not appear to be broad enough to remove the depth-velocity ambiguity; the inclusion of a priori information is necessary. They also present inversion results for station RSCP, located in the Cumberland Plateau, Tennessee. The results are similar to those from a previous study by Owens et al. (1984) and demonstrate the uncertainties in the resulting velocity estimate more clearly.

712 citations


Journal ArticleDOI
TL;DR: In this paper, a formal analysis for profile retrieval is developed which is independent of the nature of the retrieval method, provided that the measurement process can be characterized adequately, and the relationship between the retrieved and true profiles is expressed in terms of a smoothing function which can be straightforwardly calculated.
Abstract: The characterization and error analysis of profiles retrieved from remote measurements present conceptual problems, particularly concerning interlevel correlations between errors, the smoothing effect of remote sounding and the contribution of a priori information to profile. A formal analysis for profile retrieval is developed which is independent of the nature of the retrieval method, provided that the measurement process can be characterized adequately. The relationship between the retrieved and true profiles is expressed in terms of a smoothing function which can be straightforwardly calculated. The retrieval error separates naturally into three components, (1) random error due to measurement noise, (2) systematic error due to uncertain model parameters and inverse model bias, and (3) null-space error due to the inherent finite vertical resolution of the observing system. A recipe is given for evaluating each of the components in any particular case. Most of the error terms appear as covariance matrices, rather than simple error variances. These matrices can be interpreted in terms of “error patterns”, which are statistically independent contributions to the total error. They are the multidimensional equivalent of “error bars”. An approach is described which clarifies the relation of a priori data to the retrieved profile, and identifies a priori in cases where it is not an explicit part of the retrieval.

708 citations


Journal ArticleDOI
TL;DR: In this paper, a simple eddy kinetic energy parameterization of the oceanic vertical mixing is presented, which is designed to simulate vertical mixing at all depths, from the upper boundary layer down to the abyss.
Abstract: A simple eddy kinetic energy parameterization of the oceanic vertical mixing is presented. The parameterization scheme is based on recent works on atmospheric turbulence modeling. It is designed to simulate vertical mixing at all depths, from the upper boundary layer down to the abyss. This scheme includes a single prognostic equation for the turbulent kinetic energy. The computation of the turbulent length scales is diagnostic, rather than prognostic. In weakly turbulent regions the simulated vertical diffusivity is inversely proportional to the Brunt-Vaisala frequency. In the first validation experiments presented here, the vertical mixing scheme is embedded into a simple one-dimensional model and used for upper ocean simulations at two very different test sites: the station Papa in the Gulf of Alaska and the Long-Term Upper Ocean Study (LOTUS) mooring in the Sargasso Sea. At station Papa the model successfully simulates the seasonal evolution of the upper ocean temperature field. At LOTUS the focus is on a short 2-week period. A detailed analysis of the oceanic heat budget during that period reveals a large bias in the bulk-derived surface heat fluxes. After correction of the fluxes the model does well in simulating the evolution of the temperature and wind-driven current. In particular, the large observed diurnal cycles of the sea surface temperature are well reproduced. During the second (windy) week of the selected period the model accounts for about two thirds of the kinetic energy of the observed upper ocean currents at periods larger than 6 hours. The local wind forcing thus appears to be the dominant generation mechanism for the near-inertial motions, which are the most energetic. The velocity simulation is especially good at the low frequencies. During the second simulated week the model accounts for as much as 78% of the kinetic energy at subinertial frequencies. The simulated mean velocity profile is reminiscent of an Ekman spiral, in agreement with the observations.

Journal ArticleDOI
TL;DR: In this article, the authors proposed a model of the solar wind as a fluid which contains both classical transverse Alfvenic fluctuations and a population of quasi-transverse fluctuations.
Abstract: Assuming that the slab and isotropic models of solar wind turbulence need modification (largely due to the observed anisotropy of the interplanetary fluctuations and the results of laboratory plasma experiments), this paper proposes a model of the solar wind. The solar wind is seen as a fluid which contains both classical transverse Alfvenic fluctuations and a population of quasi-transverse fluctuations. In quasi-two-dimensional turbulence, the pitch angle scattering by resonant wave-particle interactions is suppressed, and the direction of minimum variance of interplanetary fluctuations is parallel to the mean magnetic field. The assumed incompressibility is consistent with the fact that the density fluctuations are small and anticorrelated, and that the total pressure at small scales is nearly constant.

Journal ArticleDOI
TL;DR: In this paper, the impact of clouds on the earth's radiation balance is assessed in terms of longwave, shortwave, and net cloud forcing by using monthly averaged clear-sky and cloudy-sky flux data derived from the NASA Earth Radiation Budget Experiment (ERBE).
Abstract: The impact of clouds on the earth's radiation balance is assessed in terms of longwave, shortwave, and net cloud forcing by using monthly averaged clear-sky and cloudy-sky flux data derived from the NASA Earth Radiation Budget Experiment (ERBE). Emphasis is placed on regional measurements, regional cloud forcing, zonal cloud forcing, and snow and ice contributions. It is shown that the global mean cooling varied from 14 to 21 W/sq m between April 1985 and January 1986; hemispherically, the longwave and shortwave cloud forcing nearly cancel each other in the winter hemisphere, while in the summer the negative shortwave cloud forcing is significantly lower than the longwave cloud forcing, producing a strong cooling. The ERBE data reveal that globally, hemispherically, and zonally, clouds have a significant effect on the radiative heating gradients.

Journal ArticleDOI
TL;DR: In this paper, the authors derived the complete solutions for seismic wave reflection, conversion, and transmission across a displacement and velocity discontinuity between two half-spaces with different densities and elastic properties for all angles of the incident wave.
Abstract: Fractures and other nonwelded contacts are important mechanical and hydrological features of rock masses. Their effects on seismic wave propagation can be modeled as a boundary condition in the seismic wave equation. Seismic stress is continuous across such a boundary, but seismic particle displacement and seismic particle velocity are not. The complete solutions for seismic wave reflection, conversion, and transmission across a displacement and velocity discontinuity between two half-spaces with different densities and elastic properties are derived for all angles of the incident wave. The ratio between the seismic stress across this boundary and the seismic particle displacement and velocity are described by a specific stiffness and a specific viscosity, respectively. A displacement discontinuity results in frequency-dependent reflection and transmission coefficients and a frequency-dependent group time delay. The velocity discontinuity results in frequency-independent coefficients and zero delay. Results of laboratory experiments on compressional and shear wave transmission across three different natural fractures in a quartz monzonite are described. Measurements were made at different effective stresses under dry and saturated conditions at room temperature. It is shown that the effect of these fractures on the spectral amplitudes for compressional and shear pulses transmitted across these fractures are described well by a displacement discontinuity for compressional pulses under dry and saturated conditions and by a combined displacement and velocity discontinuity for shear wave pulses under dry and saturated conditions. Values of specific stiffness and specific viscosity vary between fractures and increase with increasing effective stress, as does the static specific stiffness of these fractures. Changes in the spectral amplitudes of transmitted pulses are also analyzed in terms of attenuation using the seismic quality factor Q, which is found to be a function of frequency.

Journal ArticleDOI
TL;DR: In this article, the authors investigated the factional properties and stability of frictional sliding for simulated fault gouge and found that the sliding occurs at approximately constant shear stress and net compaction from one load/unload cycle to another ceases.
Abstract: This paper presents an investigation of the factional properties and stability of frictional sliding for simulated fault gouge. In these experiments we sheared gouge layers (quartz sand) under saturated drained conditions and at constant normal stress (50–190 MPa) between either rough steel surfaces or Westerly granite surfaces in a triaxial apparatus. Surface roughness (60 to 320 grit) and gouge layer thickness (0–4.0 mm) were varied in the experiments with granite samples. Porosity ϕ was monitored continuously during shear. Our measurements indicate that granular gouge exhibits strain hardening and net compaction for shear strains γ less than 0.5–1.0. For γ > 0.5–1.0, sliding occurs at approximately constant shear stress and net compaction from one load/unload cycle to the next ceases. Dilatancy occurs at 1/3 to 1/2 the shear stress required for sliding and d2ϕ/dγ2 becomes negative at about the peak stress in a given loading cycle, indicating the onset of shear localization. Oblique shear bands appear in the layers at γ = 1.3–1.5. Experiments with an initial gouge layer exhibit velocity strengthening (the coefficient of friction increases with slip velocity), and initially bare granite surfaces exhibit velocity weakening. The magnitude of velocity strengthening varies inversely with normal stress and directly with gouge thickness and surface roughness. In the gouge experiments the dilatancy rate dϕ/dγ also varies with slip rate. Using a simple energy balance to relate volume change and frictional resistance, we find quantitative agreement between the measured change in dilatancy rate and friction following changes in slip rate. This indicates that velocity strengthening within granular gouge is the result of dilatancy. The slip rate dependence of dϕ/dγ increases with gouge thickness and surface roughness, in agreement with the friction data. Our data therefore suggest that slip within unconsolidated granular material, such as some natural fault gouges, is inherently stable. The results thus provide an explanation for (1) the tendency of gouge accumulation to stabilize slip in laboratory samples, and (2) the tendency for aseismic slip within shallow (< 3–5 km) unconsolidated fault gouge and within unconsolidated sediments such as shallow alluvium and accretionary prisms.

Journal ArticleDOI
TL;DR: In this paper, the authors examined the sensitivity of ozone concentrations in rural areas of the US to emissions of NO{sub x} and hydrocarbons using a regional photochemical model.
Abstract: The authors examine the sensitivity of ozone concentrations in rural areas of the US to emissions of NO{sub x} and hydrocarbons using a regional photochemical model. Ozone production in rural areas appears to be limited by the availability of NO{sub x}. Rural ozone is strongly dependent on emission rates for NO{sub x} but is almost independent of hydrocarbons. This relationship is quite different from that in urban air, where ozone levels depend on both NO{sub x} and hydrocarbons. The predicted relationship between ozone and nitrogen oxides appears to be consistent with observations in rural air. For the low NO{sub x} regime (< 2 ppb) in rural areas, increases in NO{sub x} lead to increases in OH and to corresponding increases in the oxidation rate of hydrocarbons and in levels of ozone. Ozone concentrations in urban plumes appear to be related to regional scale production in addition to production within the plume.

Journal ArticleDOI
TL;DR: In this paper, a general class of continuous multiplicative processes is constructed to provide a theoretical framework for these observations, namely, the log-log linearity of moments with spatial scale and the concavity of corresponding slopes.
Abstract: Two common properties of empirical moments shared by spatial rainfall, river flows, and turbulent velocities are identified: namely, the log-log linearity of moments with spatial scale and the concavity of corresponding slopes with respect to the order of the moments. A general class of continuous multiplicative processes is constructed to provide a theoretical framework for these observations. Specifically, the class of log-Levy-stable processes, which includes the lognormal as a special case, is analyzed. This analysis builds on some mathematical results for simple scaling processes. The general class of multiplicative processes is shown to be characterized by an invariance property of their probability distributions with respect to rescaling by a positive random function of the scale parameter. It is referred to as (strict sense) multiscaling. This theory provides a foundation for studying spatial variability in a variety of hydrologic processes across a broad range of scales.

Journal ArticleDOI
TL;DR: For example, in this article, three Sm-Nd isochrons on tholeiitic lavas were obtained at 2.229±0.042 Ga and initial eNd = 3.6±1.1 Ga.
Abstract: Birimian terranes from West Africa (Mauritania, Senegal, Ivory Coast, Burkina Faso, Niger) comprise two major units: a dominantly mafic bimodal volcanic unit and a volcano-detrital unit with mostly felsic to intermediate protolith. Stratigraphic relationships of these units are still a matter of debate but current work suggest that they both formed in a short time interval around 2.1 Ga. Widespread basaltic magmas from the bimodal unit have been analyzed for REE distributions and Sr-Nd isotopes. Three Sm-Nd isochrons on tholeiitic lavas were obtained at 2.229±0.042 Ga and initial eNd = 3.6±1.0 for Mauritania, 2.126±0.024 Ga and initial eNd = 2.9±0.7 for Burkina Faso, 2.063±0.041 Ga and initial eNd = 3.1± .0 for Eastern Senegal, data which compare with the age of 2.11±0.09 Ga and initial eNd = 2.1±1.8 obtained in Guyana by Gruau et al. (1985). Samples from other localities (Ivory Coast, Niger) give generally similar results. Although the variations of Sm/Nd ratios and the scatter of eNd(T) values from +1.2 to +4.3 preclude a single origin for these magmas, initial isotopic heterogeneities are unlikely to bias significantly the ages given by the isochrons which are in good agreement with U-Pb zircon ages (Boher et al., 1989; unpublished data, 1990). Presence of lavas with frequent pillow structures and sediments virtually free of older recycled components suggests that Birimian terranes formed in ocean basins far from continental influence. The isotopic heterogeneities are not consistent with a MORB-like mantle source. Most lavas are slightly depleted in LREE and inversion of the data through a melting model suggests 5–15 percent melting of a slightly depleted Iherzolite. Strong depletion (Burkina Faso) and slight enrichment (Senegal) are occasionally observed. With a noticeable trend of Ti enrichment with differentiation intermediate between that of MORB and IAT, the geochemical signature of Birimian basalts does not fit the best known geodynamic environments. Back-arc or low-Ti continental flood basalts provide a marginally good agreement but still face some difficulties. Oceanic flood basalts similar to those which form oceanic plateaus (e.g. in the Nauru basin) and later accreted to continents as allochtonous terranes represent the most acceptable modern analogue of many Proterozoic basalts. It is suggested that deep plumes piercing young lithosphere can generate huge amounts of tholeiites in a short time. Birimian basalts, like many Early Proterozoic basalts, may also be viewed as recent equivalents of the Archean greenstone belts. The modern komatiite of Gorgona Island is suggested to fit this model of intraplate volcanism. Although the 2.1 Ga magmatic event in West Africa has gone virtually unnoticed in the literature, it extends over several thousand kilometers and compares with the distribution of mantle-derived magmatic activity in other major orogenic provinces (e.g. Superior). It shows that the growth rate of continents cannot be extrapolated from the data obtained solely from the best studied continents (North America, Europe, Australia). If such large crustal segments were overlooked, a spurious pattern of episodic activity of the mantle could arise.

Journal ArticleDOI
TL;DR: In this article, an analysis of more than 22,000 ozone profiles from Stratospheric Aerosol and Gas Experiment I (SAGE I) and SAGE II between 50 deg N and 50 deg S is used in conjunction with 9 years (1979-1987) of daily global depictions of total ozone from the TOMS instrument aboard Nimbus 7 to investigate the spatial distribution and seasonal cycle of the integrated amount of ozone in the troposphere.
Abstract: An analysis of more than 22,000 ozone profiles from Stratospheric Aerosol and Gas Experiment I (SAGE I) (1979-1981) and SAGE II (1984-1987) between 50 deg N and 50 deg S is used in conjunction with 9 years (1979-1987) of daily global depictions of total ozone from the TOMS instrument aboard Nimbus 7 to investigate the spatial distribution and seasonal cycle of the integrated amount of ozone in the troposphere. In the tropics, highest concentrations are found in the eastern Atlantic Ocean downwind (west) of Africa and maximize during the time when biomass burning is most prevalent, between July and October. A different seasonal cycle in the tropics is also observed over Indonesia, where a relative maximum is present in the March-April time frame, likewise consistent with when biomass burning is most prevalent. At mid-latitudes, highest concentrations are found downwind of Asia and maximize in the summer. Relatively higher amounts of tropospheric ozone are similarly observed downwind of North America and Europe. At mid-latitudes, the ratio between the amount of tropospheric ozone in the Northern Hemisphere and the amount in the Southern Hemisphere is 1.4, in good agreement with in situ measurements.

Journal ArticleDOI
TL;DR: The absorption spectrum of O2 and O2-O2 collision pairs were measured over the wavelength range from 330 to 1140 nm using pressures of O 2 from 1 to 55 atm at 298 K as mentioned in this paper.
Abstract: The absorption spectrum of O2 and O2-O2 collision pairs were measured over the wavelength range from 330 to 1140 nm using pressures of O2 from 1 to 55 atm at 298 K. Absorption cross sections, pressure dependences, band centers, and full widths at half maximum of the observed absorption bands centered at 343.4, 360.5, 380.2, 446.7, 477.3, 532.2, 577.2, 630.0, 688, 762, and 1065.2 nm are reported. The absorption bands centered at 360.5, 380.2, and 477.3 nm were also measured at 196 K and their temperature dependences were characterized.

Journal ArticleDOI
TL;DR: In this paper, a suite of porous sandstones (Boise, Kayenta, St Peter, Berea, and Weber) were selected because of their wide range of porosity and grain size.
Abstract: The hydrostatic compaction behavior of a suite of porous sandstones was investigated at confining pressures up to 600 MPa and constant pore pressures ranging up to 50 MPa. These five sandstones (Boise, Kayenta, St. Peter, Berea, and Weber) were selected because of their wide range of porosity (5–35%) and grain size (60–460 μm). We tested the law of effective stress for the porosity change as a function of pressure. Except for Weber sandstone (which has the lowest porosity and smallest grain size), the hydrostat of each sandstone shows an inflection point corresponding to a critical effective pressure beyond which an accelerated, irrecoverable compaction occurs. Our microstructural observations show that brittle grain crushing initiates at this critical pressure. We also observed distributed cleavage cracking in calcite and intensive kinking in mica. The critical pressures for grain crushing in our sandstones range from 75 to 380 MPa. In general, a sandstone with higher porosity and larger grain size has a critical pressure which is lower than that of a sandstone with lower porosity and smaller grain size. We formulate a Hertzian fracture model to analyze the micromechanics of grain crushing. Assuming that the solid grains have preexisting microcracks with dimensions which scale with grain size, we derive an expression for the critical pressure which depends on the porosity, grain size, and fracture toughness of the solid matrix. The theoretical prediction is in reasonable agreement with our experimental data as well as other data from soil and rock mechanics studies for which the critical pressures range over 3 orders of magnitude.

Journal ArticleDOI
TL;DR: In this article, the spatial and temporal variability of nearshore sand bar morphology is quantified using a unique data set spanning 2 years, which consists of daily time exposure images of incident wave breaking on an open coast sandy beach which may be used to infer bar morphology.
Abstract: The spatial and temporal variability of nearshore sand bar morphology is quantified using a unique data set spanning 2 years. The data consist of daily time exposure images of incident wave breaking on an open coast sandy beach which may be used to infer bar morphology (Lippmann and Holman, 1989). The morphology in each image is classified into an eight state morphologic scheme in which bars are uniquely defined by four independent criteria. The most frequently observed morphologies are the longshore-periodic (rhythmic) bars, observed in 68% of the data. Linear bars occur under highest wave conditions (H¯s=1.78 m) and are unstable (mean residence time ≈ 2 days). Shore-attached rhythmic bars are the most stable (mean residence time ≈ 11 days) and generally form 5–16 days following peak wave events. Non-rhythmic, three-dimensional bars are very transient (mean residence time ≈ 3 days). Eighty-seven percent of transitions to lower bar types (defined in text) occurred one state at a time, supporting our selection of the ordering of states, and suggesting the suitability of a sequential morphology model. Transitions to higher states occurred under rising wave energy and were evenly spread among the possible higher states, with more substantial changes in morphology resulting from larger wave increases. This suggests that up-state, erosional transitions (based on offshore bar migration) are better described by an equilibrium model where response is better correlated with incident wave energy than with preceding morphological state. Time exposure images were also digitized to yield quantitative estimates of bar crest location as a function of longshore distance. Principal component analysis was used to decompose bar position into two-dimensional (linear) and three-dimensional (longshore variable) components. Cross-shore (linear) bar position ranges ±50 m about the 2-year mean (27 m standard deviation) and dominates bar variability (74.6%). Three-dimensional bar structure accounts for ∼14% of the variance (12 m standard deviation). Changes in incident wave height precede cross-shore bar migration by less than 1 day. Changes in longshore variability are inversely correlated to changing wave conditions, with bar morphology becoming linear rapidly during storms (on time scales of less than 1 day). Evolution to significantly three-dimensional structure typically occurs over 5–7 days following peak wave events.

Journal ArticleDOI
TL;DR: In this paper, a one-dimensional radiative-convective equilibrium model of surface-temperature changes induced by changes in the vertical distribution of O3 in the atmosphere is parameterized and applied to measurement data for the northern midlatitudes during the 1970s.
Abstract: Results from a one-dimensional radiative-convective equilibrium model (similar to that of Lacis et al., 1981) of surface-temperature changes induced by changes in the vertical distribution of O3 in the atmosphere are parameterized and applied to measurement data for the northern midlatitudes during the 1970s. The construction of the model is outlined, and the results are presented in tables and graphs and discussed in detail. It is found that decreases in O3 in the lower stratosphere in the 1970s caused atmospheric cooling greater than that due to CO2 (thus possibly obscuring the CO2 greenhouse effect) and surface cooling which outweighed the warming effect of increased O3 in the troposphere. The net surface cooling is shown to be equal to about half of the estimated CO2 warming effect for the period.

Journal ArticleDOI
TL;DR: In this paper, high spatial resolution column atmospheric water vapor amounts were derived from spectral data collected by the airborne visible-infrared imaging spectrometer (AVIRIS), which covers the spectral region from 0.4 to 2.5 μm in 10-nm bands and has a ground instantaneous field of view of 20×20 m from an altitude of 20 km.
Abstract: High spatial resolution column atmospheric water vapor amounts were derived from spectral data collected by the airborne visible-infrared imaging spectrometer (AVIRIS), which covers the spectral region from 0.4 to 2.5 μm in 10-nm bands and has a ground instantaneous field of view of 20×20 m from an altitude of 20 km. The quantitative derivation is made by curve fitting observed spectra with calculated spectra in the 1.14-μm and 0.94-μm water vapor band absorption regions using an atmospheric model, a narrowband spectral model, and a nonlinear least squares fitting technique. The derivation makes use of the facts that (1) the reflectances of many ground targets vary approximately linearly with wavelength in the 0.94- and 1.14-μm water vapor band absorption regions, (2) the scattered radiation near 1 μm is small compared with the directly reflected radiation when the atmospheric aerosol concentrations are low, and (3) the scattered radiation in the lower part of the atmosphere is subjected to the water vapor absorption. The technique is directly applicable for retrieving column water vapor amounts from AVIRIS spectra measured on clear days with visibilities 20 km or greater. The precision of the retrieved column water vapor amounts from several data sets is 5% or better. Based on the analyses of an AVIRIS data set that was acquired within an hour of radiosonde launch, it appears that the accuracy approaches the precision. The derived column water vapor amounts are independent of the absolute surface reflectances. It now appears feasible to derive high spatial resolution column water vapor amounts over land areas from satellite altitude with the proposed high resolution imaging spectrometer (HIRIS). Curve fitting of spectra near 1 μm from areas covered with vegetation, using an atmospheric model and a simplified vegetation reflectance model, indicates that both the amount of atmospheric water vapor and the moisture content of vegetation can be retrieved simultaneously because the band centers of liquid water in vegetation and the atmospheric water vapor are offset by approximately 0.05 μm.

Journal ArticleDOI
TL;DR: In this article, measurements and analyses of the aerosol size distributions and optical properties found in the marine boundary layer (MBL) during the 1983 USNS Lynch cruise from Charleston, South Carolina, to Scotland via Canary Islands are presented.
Abstract: Measurements and analyses of the aerosol size distributions and optical properties found in the marine boundary layer (MBL) during the 1983 USNS Lynch cruise from Charleston, South Carolina, to Scotland via Canary Islands are presented. The data given are the most extensive and accurate measurements of the submicron marine aerosol size distribution to date and are supplemented by extensive meteorological observations. Eight detailed case studies of the evolution of the size distribution that occurred under different meteorological conditions are presented and discussed. The data indicate that repeated cycling of MBL air through nonprecipitating clouds at the top of the MBL is a major factor in shaping the size distribution and that new particle formation by heteromolecular, homogeneous nucleation is the most likely mechanism for sustaining the particle concentration below 0.04-μm radius. Calculations of the scattering and extinction coefficients and optical depth of the MBL as a function of wavelength directly from the measured size distribution and MBL vertical structure are compared to measured values of the scattering coefficient and optical depth. These measured and calculated optical properties correlate well throughout the cruise and the results give a relatively consistent picture of the relationship between the aerosol size distribution and electromagnetic properties in the MBL.

Journal ArticleDOI
TL;DR: In this paper, the diffraction patterns obtained by XRD indicate that pure iron remains in the hcp structure to 304 GPa and that Fe(0.8)Ni( 0.2) is stable to at least 255 GPa.
Abstract: Results are reported on a room temperature static compression study of iron and the Fe(0.8)Ni(0.2) alloy to above 260 GPa, which provide direct pressure-volume measurements on geophysically important materials at the conditions close to those at the inner core boundary. The diffraction patterns obtained by XRD indicate that pure iron remains in the hcp structure to 304 GPa and that Fe(0.8)Ni(0.2) is stable to at least 255 GPa. The results of this study, in conjunction with work at higher temperatures, will make it possible to address directly the question of the composition of the inner core with a level of certainty that were previously applicable only to the mantle.

Journal ArticleDOI
TL;DR: In this article, single crystals from San Carlos in Arizona have been annealed at temperatures between 800° and 1000°C under hydrothermal conditions at a confining pressure of 300 MPa.
Abstract: To investigate the kinetics of diffusion of hydrogen in olivine, single crystals from San Carlos in Arizona have been annealed at temperatures between 800° and 1000°C under hydrothermal conditions at a confining pressure of 300 MPa. The hydrogen diffusivities were determined for the [100], [010], and [001] directions from concentration profiles for hydroxyl in the samples. These profiles were obtained from infrared spectra taken at 100-μm intervals across a thin slice which was cut from the central portion of each annealed crystal. The rate of diffusion is anisotropic, with fastest transport along the [100] axis and slowest along the [010] axis. The fit of the data to an Arrhenius law for diffusion parallel to [100] yields an activation enthalpy of 130±30 kJ/mol with a preexponential term of (6±3)×10−5 m2 s−1. For diffusion parallel to [001], as there are insufficient data to calculate the activation enthalpy for diffusion, we used the same value as that for diffusion parallel to [100] and determined a preexponential term of (5±4)×10−6 m2 s−1. The diffusion rate parallel to [010] is about 1 order of magnitude slower than along [001]. The measured diffusivities are large enough that the hydrogen content of olivine grains which are millimeters in diameter will adjust to changing environmental conditions in time scales of hours at temperatures as low as 800°C. As xenoliths ascending from the mantle remain at high temperatures (i.e., >1000°C) but experience a rapid decrease in pressure, and hence hydrogen fugacity, olivine grains may dehydrate during ascent. By comparison, slow rates of carbon diffusion (Tingle et al., 1988) suggest that carbon will not be lost from olivine during ascent. Thus, low hydrogen contents within olivine and within fluid inclusions in olivine cannot be taken as support for low water contents in the mantle.

Journal ArticleDOI
TL;DR: An improved approach to spectral deconvolution is presented in this article that accurately represents absorption bands as discrete mathematical distributions and resolves composite absorption features into individual absorption bands, and a modified Gaussian model is derived using a power law relationship of energy to average bond length.
Abstract: Although visible and near IR reflectance spectra contain absorption bands that are characteristic of the composition and structure of the absorbing species, deconvolving a complex spectrum is nontrivial. An improved approach to spectral deconvolution is presented that accurately represents absorption bands as discrete mathematical distributions and resolves composite absorption features into individual absorption bands. The frequently used Gaussian model of absorption bands is shown to be inappropriate for the Fe(2+) electronic transition absorptions in pyroxene spectra. A modified Gaussian model is derived using a power law relationship of energy to average bond length. The modified Gaussian model is shown to provide an objective and consistent tool for deconvolving individual absorption bands in the more complex orthopyroxene, clinopyroxene, pyroxene mixtures, and olivine spectra.

Journal ArticleDOI
TL;DR: In this paper, the authors reviewed experiments and models describing polynyas and leads in the polar oceans and showed that large sensible heat fluxes, together with evaporation and longwave radiation from a very small percentage of open water and thin ice dominate regional heat budgets.
Abstract: Polynyas and leads are openings in pack ice due to divergences in ice drift and to local melting. They are the vents and windows to the polar oceans. In winter they are a major source of brine during freezing and a locus for gas exchange. Large sensible heat fluxes, together with evaporation and longwave radiation from a very small percentage of open water and thin ice, dominate regional heat budgets. In summer, solar radiation is absorbed by open water but is reflected from snow-covered pack ice. Experiments and models describing these processes are reviewed.

Journal ArticleDOI
TL;DR: In this article, the authors analyzed four expanded spread profiles acquired at distances of 0, 2.1, 3.1 and 10 km from the axis of the East Pacific Rise between 9° and 10°N.
Abstract: We analyze four expanded spread profiles acquired at distances of 0, 2.1, 3.1, and 10 km (0–0.2 m.y.) from the axis of the East Pacific Rise between 9° and 10°N. Velocity-depth models for these profiles have been obtained by travel time inversion in the τ-p domain, and by x-t forward modeling using the WKBJ and the reflectivity methods. We observe refracted arrivals that allow us to determine directly the uppermost crustal velocity structure (layer 2A). At the seafloor we find very low Vp and VS/Vp values around 2.2 km/s and ≤ 0.43. In the topmost 100–200 m of the crust, Vp remains low (≤ 2.5 km/s) then rapidly increases to 5 km/s at ∼500 m below the seafloor. High attenuation values (Qp < 100) are suggested in the topmost ∼500 m of the crust. The layer 2–3 transition probably occurs within the dike unit, a few hundred meters above the dike-gabbro transition. This transition may mark the maximum depth of penetration by a cracking front and associated hydrothermal circulation in the axial region above the axial magma chamber (AMC). The on-axis profile shows arrivals that correspond to the bright AMC event seen in reflection lines within 2 km of the rise axis. The top of the AMC lies 1.6 km below the seafloor and consists of molten material where Vp ≈ 3 km/s and VS = 0. Immediately above the AMC, there is a zone of large negative velocity gradients where, on the average, Vp decreases from ∼6.3 to 3 km/s over a depth of approximately 250 m. Associated with the AMC there is a low velocity zone (LVZ) that extends to a distance no greater than 10 km away from the rise axis. At the top of the LVZ, sharp velocity contrasts are confined to within 2 km of the rise axis and are associated with molten material or material with a high percentage of melt which would be concentrated only in a thin zone at the apex of the LVZ, in the axial region where the AMC event is seen in reflection lines. Away from the axis, the transition to the LVZ is smoother, the top of the LVZ is deeper, and the LVZ is less pronounced. The bottom of the LVZ is probably located near the bottom of the crust and above the Moho. Moho arrivals are observed in the profiles at zero and at 10 km from the rise axis. Rather than a single discontinuity, these arrivals indicate an approximately 1-km-thick Moho transition zone.