Patterns and controls of surface sediment distribution: west-central Florida inner shelf

The west-central Florida inner shelf represents a transition between the quartz-dominated barrier-island system and the carbonate-dominated mid-outer shelf. Surface sediments exhibit a complex distribution pattern that can be attributed to multiple sediment sources and the ineffectiveness of physical processes for large-scale sediment redistribution. The west Florida shelf is the submerged extension of the Florida carbonate platform, consisting of a limestone karst surface veneered with a thin unconsolidated sediment cover. A total of 498 surface sediment samples were collected on the inner shelf and analyzed for texture and composition. Results show that sediment consists of a combination of fine quartz sand and coarse, biogenic carbonate sand and gravel, with variable but subordinate amounts of black, phosphorite-rich sand. The carbonate component consists primarily of molluskan fragments. The distribution is patchy and discontinuous with no discernible pattern, and the transition between sediment types is generally abrupt. Quartz-rich sediment dominates the inner 15 km north of the entrance into Tampa Bay, but south of the Bay is common only along the inner 3 km. Elsewhere, carbonate-rich sediment is the predominate sediment type, except where there is little sediment cover, in which cases black, phosphorite-rich sand dominates. Sediment sources are likely within, or around the periphery of the basin. Fine quartz sand is likely reworked from coastal units deposited during Pleistocene sea-level high stands. Carbonate sand and gravel is produced by marine organisms within the depositional basin. The black, phosphorite-rich sand likely originates from the bioerosion and reworking of the underlying strata that irregularly crop out within the study area. The distribution pattern contains elements of both storm- and tide-dominated siliciclastic shelves, but it is dictated primarily by the sediment source, similar to some carbonate systems. Other systems with similar sediment attributes include cool-water carbonate, sediment-starved, and mixed carbonate/siliciclastic systems. This study suggests a possible genetic link among the three systems.     

Environmental tracers as indicators of karst conduits in groundwater in South Dakota,USA

Environmental tracers sampled from the carbonate Madison aquifer on the eastern flank of the Black Hills, South Dakota, USA indicated the approximate locations of four major karst conduits. Contamination issues are a major concern because these conduits are characterized by direct connections to sinking streams, high groundwater velocities, and proximity to public water supplies. Objectives of the study were to estimate approximate conduit locations and assess possible anthropogenic influences associated with conduits. Anomalies of young groundwater based on chlorofluorocarbons (CFCs), tritium, and electrical conductivity (EC) indicated fast moving, focused flow and thus the likely presence of conduits. δ¹⁸O was useful for determining sources of recharge for each conduit, and nitrate was a useful tracer for assessing flow paths for anthropogenic influences. Two of the four conduits terminate at or near a large spring complex. CFC apparent ages ranged from 15 years near conduits to >50 years in other areas. Nitrate-N concentrations >0.4 mg/L in groundwater were associated with each of the four conduits compared with concentrations ranging from <0.1 to 0.4 mg/L in other areas. These higher nitrate-N concentrations probably do not result from sinking streams but rather from other areas of infiltration.     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

The elastic properties of single-crystal fayalite as determined by dynamical measurement techniques

We present new elasticity measurements on single-crystal fayalite and combine our results with other data from resonance, pulse superposition interferometry, and Brillouin scattering to provide a set of recommended values for the adiabatic elastic moduliC _ij and their temperature variations. We use a resonance method (RPR) with specimens that were previously investigated by pulse superposition experiments. The nineC _ij of fayalite are determined from three new sets of measurements. One set of our newC _ij data is over the range 300–500 K. We believe that the relatively large discrep ncies found in someC _ij are due in large part to specimen inhomogeneities (chemical and microstructural) coupled with differences in the way various techniques sample, rather than only systematic errors associated with experimental procedures or in the preparations of the specimens.Our recommendeaC _ij's (GPa) and (C _ij/T) _p (GPa/K) are: The resulting values for the isotropic bulk and shear moduli,K _s and , and their temperature derivatives are:K _s=134(4) GPa; =50.7(0.3) GPa; (K _s/T) _p =–0.024(0.005) GPa/K; and (/T) _p =–0.013(0.001) GPa/K. An important conclusion is thatK _s increases as the Fe/(Fe+Mg) ratio in olivine is increased.     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

Charged dust in the outer planetary magnetospheres

We assume that typical interplanetary grains are fragile, aggregates of the Brownlee type, and discuss the physical and dynamical processes associated with their entry into the Jovian magnetosphere. Limiting ourselves to the equatorial plane of the planet, we show that grains traversing the outer edge of the co-rotating magnetodisc (r35R _J ) are rapidly charged up to large negative potentials on both the day and the night sides. A parent grain of radiusR _g 20 is electrostatically disrupted when it attains a potential of about –220 V. While the eventual potential achieved by the smallest fragments (R _g 0.1 ) are controlled by the rapid field emission of electrons, those of the larger fragments (R _g 1 ) are set by the plasma and photoemission currents.All the negatively charged fragments are strongly attracted towards the planet by the (radial) corotational electric field and some are stably trapped. We suggest that the sudden enhancement by about 2 orders of magnitude of the interplanetary dust flux measured by Pioneer 10, at about 30R _J from Jupiter result from the combination of these two effects.The different brightness asymmetries between the leading and the trailing sides of the outer and inner Galilean satellites appear to be a natural consequence of the way the trajectories of these charged dust grains intersect these satellite surfaces. Finally, the similarity in the brightness asymmetries between the Jovian and Saturnian satellites, and our belief that they have a similar cause, leads us to the expectation that Saturn's magnetic momentM and spin , are parallel as in the case of Jupiter, with the limit of plasma co-rotation lying between the satellites Rhea and Iapetus.     

Numerical simulation of the final stages of terrestrial planet formation

Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are presented. The mass and relative-velocity distributions of the bodies in these simulations are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies, therefore, arise in a natural way and the assumption of expressions for the relative velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and the Earth is possible, at least for the case of a two-dimensional system. Simulations assuming an initial uniform distribution of orbital eccentricities on the interval from 0 to e_max are found to produce final states containing too many bodies with masses which are too small when e_max < 0.10, while simulations with e_max > 0.20 result in too many catastrophic collisions between bodies thus preventing rapid accretion of planetary-size bodies. The e_max = 0.15 simulation ends with a state surprisingly similar to that of the present terrestrial planets and, therefore, provides a rough estimate of the range of radial sampling to be expected for the terrestrial planets.     

Global tectonics of a despun planet

Mercury, the Moon, and many large satellites of the major planets have been tidally despun from an initially faster rotation. These bodies probably possessed equatorial bulges which relaxed as they lost their spin. An analysis of the stresses induced in an elastic shell by the relaxation of an equatorial bulge indicates that differential stresses may reach a few kilobars and that the tectonic pattern developed depends mainly upon the shell thickness. In every model studied the azimuthal stress σ_?? is larger (more compressive) than the meridional stress σθθ. For a thin elastic shell (thickness less than one-twentieth of the planet's radius) the zone from the equator to 48° latitude is characterized by strike-slip faulting. Poleward of this, normal faults and graben trending east-west are expected. Thicker elastic shells acquire an equatorial belt of thrust faults with east-west throw and rough north-south trends. These tectonic styles may be modified by a small (0.05-0.1%) radial expansion or contraction. Expansion shifts the polar normal faulting province toward the equator, while contraction shifts the equatorial provinces poleward. These patterns are not substantially altered by plastic yielding of the shell, although the equatorial thrust fault province is suppressed by strike-slip faulting until strike-slip faults occur poleward of 64.8° latitude. We conclude that there are many tectonic patterns consistent with despinning and radial contraction or expansion, but they must all be consistent with σ_?? > σθθ. These results also indicate that the polar regions of a despun planet are of particular interest in deciding whether a given lineament system is due to stresses induced by the relaxation of the planet's equatorial bulge.     

The influence of amphibole fractionation on the evolution of calc-alkaline andesite and dacite tephra from the central Aleutians,Alaska

Petrologic studies of tephra from Kanaga, Adak, and Great Sitkin Islands indicate that amphibole fractionation and magma mixing are important processes controlling the composition of calc-alkaline andesite and dacite magmas in the central Aleutians. Amphibole is ubiquitous in tephra from Kanaga and Adak Islands, whereas it is present only in a basaltic-andesite pumice from Great Sitkin. Dacitic tephra from Great Sitkin do not contain amphibole. Hornblende dacite tephra contain HB+PLAG+OX±OPX±CPX phenocrysts with simple zoning patterns, suggesting that the dacites evolved in isolated magma chambers. Andesitic tephra from Adak contain two pyroxene and hornbelende populations, and reversely zoned plagioclase, indicating a more complex history involving mixing and fractional crystallization. Mass balance calculations suggest that the andesitic tephra may represent the complements of amphibole-bearing cumulate xenoliths, both formed during the evolution of high-Al basalts. The presence of amphibole in andesitic and dacitic tephra implies that Aleutian cale-alkaline magmas evolve in the mid to lower crust under hydrous (>4 wt.% H₂O) and oxidizing (Ni–NiO) conditions. Amphibole-bearing andesites and pyroxene-bearing dacites from Great Sitkin indicates fractionation at several levels within the arc crust. Despite its absence in many calc-alkaline andesite and dacite lavas, open system behavior involving amphibole fractionation can explain the trace element characteristies of lavas found on Adak Island. Neither open nor closed system fractionation involving a pyroxene-bearing assemblage is capable of explaining the trace element concentrations or ratios found in the Adak suite. We envision a scenario where amphibole was initially a liquidus phase in many calc-alkaline magmas, but was later replaced by pyroxenes as the magmas rose to shallow levels within the crust. The mineral assemblage in these evolved lavas reflects shallow level equilibration of the magma, whereas the trace element chemistry provides evidence for a earlier, amphibole-bearing, mineral assemblage.     

Magmatism in the eastern Aleutian Arc: temporal characteristic of igneous activity on Akutan Island

Lavas from Akutan Island, located in the eastern Aleutian arc at the transition between continental and oceanic crust, show a gradual change in their petrologic and chemical characteristics over the last 4 million years. The oldest lavas exposed on the island, the Hot Springs Bay Volcanics (HSBV), range from magnesian basalt to dacite (45%–62% SiO₂). The most mafic basalts contain salitic clinopyroxene, Cr- and Al-rich spinel, and pargasitic amphibole suggesting that they were derived from relatively hydrous magmas at greater pressures than lavas from the younger Akutan Volcanics (AKV) and the modern volcano (MOD). AKV lavas also range between basalt and dacite (46%–63% SiO₂), but contain no hydrous phenocrysts and seem to have fractionated within a shallow level magma chamber. Lavas from the modern volcano are andesitic (52%–57% SiO₂) and have a mineral assemblage similar to that of AKV lavas of similar composition. With the exception of clinopyroxene and spinel in the most mafic lavas, the compositions of plagioclase (An_92?45), olivine (Fo_88?51), orthopyroxene (En_69?56), and titanomagnetite (15%–21% TiO₂) phenocrysts found in these lavas are within the range observed in lavas from other Aleutian volcanoes. Variations in the major element chemistry of the older lavas can be reproduced by fractional crystallization of the observed mineral assemblages, however closed system crystal fractionation models are inadequate to explain the trace element variations. During the last 4 million years, La/Yb ratios have decreased (6.5–3.3 for HSBV lavas and 2.9–1.9 for MOD lavas) whereas Ba/La ratios appear to have increased slightly (37–43 for HSBV and AKV, and 41–45 of MOD). The lower La/Yb ratios of MOD lavas correspond with lower total abundances of the REE and slightly higher Sr and Pb isotopic ratios. The increased⁸⁷Sr/⁸⁶Sr ratios and Pb isotopic ratios in the MOD lavas, the less enriched LREE, and the higher Ba/La ratios may result from partial melting of an arc source which has experienced previous melting events but has continued to be contaminated by a component from the subducting slab. It may also indicate a change in the degree of partial melting of the underlying mantle, which corresponds to a different percentage of a slab derived component being incorporated into the overlying mantle.