a self consistent model of melting magma migration and buoyancy driven circulation beneath mid ocean ridges numerical modeling and analysis are used to investigate the processes leading to the eruption of mantle derived magma at mid ocean ridges our model includes the following effects melting due to decompression magma migration by percolation and circulation of the mantle driven by both the oceanic plates and the distribution of buoyancy beneath the ridge the distribution of buoyancy is due to both the low density of the liquid and the difference in the density of the residual solids relative to unmelted mantle material the calculation of densities is based on a simple petrological model in which garnet spinel lherzolite melts to form a basaltic liquid and a harzburgite residue we find that the system spontaneously evolves to a state in which a rapid upwelling beneath the ridge axis faster than the plate velocity is confined laterally by stably stratified residual material beneath the newly formed plates this effect is exaggerated if a modest decrease in the shear viscosity of the solid upon melting is included our results provide a simple explanation for the narrowness of the zone of crustal formation at mid ocean ridges the model also predicts a transition from steady state to episodic crustal formation as the spreading velocity is reduced perhaps giving rise to along axis variations in the character of seafloor spreading the narrow rapid upwelling gives rise to substantial porosities at depths that are a large fraction of the depth to the solidus this may allow the liquid at depth to segregate into macroscopic channels which would account for the consensus from experimental petrology that the liquids parental to morb are derived from well below the base of the crust

transitions in axial morphology along the southeast indian ridge shipboard bathymetric and magnetic profiles across the southeast indian ridge seir were analyzed in order to examine the nature of along axis variations in axial morphology at this intermediate spreading rate ridge three types of axial morphology are observed along the seir an axial high a shallow 200 700 m deep axial valley and a deep 1000 m deep axial valley an axial high is found to the east of the australian antarctic discordance aad east of 128 e and between 82 e and 104 e a shallow rift valley is found from 104 e to 114 e and from 82 e westward past the amerstdam st paul hotspot asp to about 30 s 75 e deep rift valleys are found from 114 e to 128 e in the vicinity of the aad and from the indian ocean triple junction iotj at 25 s 70 e to about 30 s 75 e the transition near 30 s occurs in an area of constant zero age depth and does not appear to result from an increase in mantle temperature it could be the result of the rapid increase in spreading rate along that portion of the seir the most likely cause of the other transitions in axial morphology is variations in mantle temperature the transitions between the different types of axial morphology are well defined and occur over a limited distance transitions in axial morphology are accompanied by significant changes in ridge flank topographic roughness the transitions from axial valleys to axial highs are also accompanied by changes in the amplitude of the seafloor magnetic anomalies our observations suggest that there are distinct modes rather than a continuum of axial morphology on the seir and that there appears to be a threshold mechanism for a rapid change between different states of axial morphology the asp has only a limited influence on the seir the ridge axis is marked by an axial valley for the entire distance from the iotj up to and past the asp the ridge axis becomes shallower as the asp is approached from the northwest but only by about 300 m over a distance of 800 km in addition the ridge continues to become shallower away from amsterdam island toward the transition to an axial high at 82 e 350 km to the east of the asp the kerguelen hotspot appears to exert a major influence on the morphology of the seir by feeding asthenospheric material to the ridge axis a long narrow finger like gravity high extends ene away from the kerguelen plateau for a distance of 500 km shipboard data show that the gravity high results from a large volcanic ridge the ridge appears analogous to the rodriguez ridge extending from the reunion hotspot toward the central indian ridge a series of lower and broader lineated gravity highs extend from the volcanic ridge toward the seir in the ridge segment between the 81 e and 85 e transforms which is the westernmost segment with an axial high the only region of significant off ridge seismicity on the antarctic flank of the seir is a diffuse band of epicenters extending from kerguelen to the seir within the segment between the 81 e and 85 e fracture zones the along axis gradient in depth from 86 e to the aad and the transitions in axial morphology at 104 e and 114 e most likely reflect along axis variations in mantle temperature and melt production rate due to distance from the kerguelen hotspot and the influence of the aad

a mechanism for magmatic accretion under spreading centres the most active volcanic zone on earth occurs along the global mid ocean ridge system where the lithosphere is pulled apart and new oceanic crust is created at rates between 10 and 200 km myr 1 this volcanism is segmented and is concentrated along the axial valleys and highs of slow and fast spreading ridges respectively with transform fault zones separating and generally offsetting each active ridge segment volcanic activity may be nearly absent at the ridge transform intersection where a slow spreading ridge is offset and greatly reduced where a fast spreading ridge is offset implying thinner crust in these regions1 the source of magmas erupted along ocean ridges is the underlying mantle which undergoes decompression melting when rising to fill the gap between the spreading plates the resulting magma aggregates in the upper mantle and ascends to the crust due to its low density compared with the parent mantle rock the melt is believed to pool in crustal magma chambers whence it periodically erupts to the surface typically the formation ascent and aggregation of magmas along the mid ocean ridges is regarded two dimensionally and modelled in a section through the crust and mantle across the strike of the ridges see ref 2 here however we consider the problem with a three dimensional model for which we provide supporting geological evidence in which the magmas rise as a result of a gravitational instability modified rayleigh taylor in the underlying partial melt zone

geologic framework of the mississippian barnett shale barnett paleozoic total petroleum system bend arch fort worth basin texas this article describes the primary geologic characteristics and criteria of the barnett shale and barnett paleozoic total petroleum system tps of the fort worth basin used to define two geographic areas of the barnett shale for petroleum resource assessment from these two areas referred to as assessment units the u s geological survey estimated a mean volume of about 26 tcf of undiscovered technically recoverable hydrocarbon gas in the barnett shale the mississippian barnett shale is the primary source rock for oil and gas produced from paleozoic reservoir rocks in the bend arch fort worth basin area and is also one of the most significant gas producing formations in texas subsurface mapping from well logs and commercial databases and petroleum geochemistry demonstrate that the barnett shale is organic rich and thermally mature for hydrocarbon generation over most of the bend arch fort worth basin area in the northeastern and structurally deepest part of the fort worth basin adjacent to the muenster arch the formation is more than 1000 ft 305 m thick and interbedded with thick limestone units westward it thins rapidly over the mississippian chappel shelf to only a few tens of feet the barnett paleozoic tps is identified where thermally mature barnett shale has generated large volumes of hydrocarbons and is 1 contained within the barnett shale unconventional continuous accumulation and 2 expelled and distributed among numerous conventional clastic and carbonate rock reservoirs of paleozoic age vitrinite reflectance ro measurements show little correlation with present day burial depth contours of equal ro values measured from barnett shale and typing of produced hydrocarbons indicate significant uplift and erosion furthermore the thermal history of the formation was enhanced by hydrothermal events along the ouachita thrust front and mineral wells newark east fault system stratigraphy and thermal maturity define two gas producing assessment units for the barnett shale 1 a greater newark east fracture barrier continuous barnett shale gas assessment unit encompassing an area of optimal gas production where dense impermeable limestones enclose thick 300 ft 91 m barnett shale that is within the gas generation window ro 1 1 and 2 an extended continuous barnett shale gas assessment unit covering an area where the barnett shale is within the gas generation window but is less than 300 ft 91 m thick and either one or both of the overlying and underlying limestone barriers are absent

migration of mid ocean ridge volcanic segments the propagation of small offset volcanic spreading centre segments along mid ocean ridge crests may reflect absolute motion of the plate boundary relative to the underlying mesospheric frame such a relationship could be caused by a purely vertical flow of the mantle under spreading centres and would have value in constraining past plate motions from non transform trends generated during along ridge propagation and in linking the major element variability of oceanic crust and upper mantle to the bulk composition and temperatures of mantle ascending under mid ocean ridges

three dimensional gravity and bathymetry analysis of the mid atlantic ridge between 20 n and 24 n flow geometry and temporal evolution of the segmentation bathymetry and gravity data covering 10 m y on both flanks of the mid atlantic ridge south of the kane fracture zone reveal a complex segmentation pattern well correlated with short scale temporal variations in crustal production segment and paleosegment centers are associated with crustal thickening while their extremities appear as elongated bands of thin crust subsidence rates vary between centers and extremities of a segment with higher than world average values at segmert midpoints and closer to average values at the edges this is consistent with segment centers being the site of focused mantle upwellings with associated thin lithosphere and with thick crust a pronounced circular minimum in the long wavelength of the mantle bouguer anomaly associated to a maximum in the long wavelength of the bathymetry marks the site of a mantle upwelling associated with a particularly robust segment located north of a small fracture zone 21 20 n an offset is often observed between extreme values of topography and crustal thickness resulting in a lack of isostatic equilibrium for off axis short wavelength features an obliquity between maxima of crustal thickening associated with active segments and the topographic axial domain suggests that the shift between topography and crustal structure may arise partly from a decoupling between crust and topography formation processes

a systematic analysis of the mid atlantic ridge morphology and gravity between 15 n and 40 n constraints of the thermal structure multibeam bathymetry data obtained along a 2400 km long section of the mid atlantic ridge mar from 15 n to the azores platform 40 n and satellite derived gravity data were used to calculate the mantle bouguer anomaly mba along this portion of the mar both data sets were used to determine the relations between gravity anomalies and topographic variations and discuss these in terms of thermal difference a long wavelength influence of the azores hot spot is characterized by a gentle continuous slope of the average ridge axial depth and a general gradient in the along axis mba profile this thermal influence of the azores hot spot controls a systematic southward propagation of the spreading segments at least to 26 30 n south of 26 30 n the direction of the segment propagation is controlled by the local difference in thermal state between adjacent segments except on the azores platform the systematic along axis 11 90 km long wavelength segmentation is independent of the long wavelength influence of the azores at the segment center the axial morphology is linked to the thermal state of the segments between 1 hotter segments characterized by a smooth axial morphology a well defined shallow inner valley high mba and a long length 2 colder segments which present a rough axial morphology with a deep wide and well defined rift valley a low mba and a small length for hotter segments the formation of the abyssal hills is mainly due to a magmato tectonic cycle over periods of 0 3 to 1 myr whereas on colder segments the axial morphology is mainly controlled by a tectonic rift valley formation we propose that these different segment types correspond to a temporal evolution of the rift valley morphology over periods of several million years

the spacing of rift axis highs dependence upon diapiric processes in the underlying asthenosphere abstract detailed bathymetric mapping has revealed a near harmonic distribution of intermediate wavelength morphologic highs on slow to fast spreading centers one model suggests that the highs may result from the thermal instabilities generated by the interaction of cold fracture zone walls with the axial asthenosphere this may force a deepening of the isotherms near the fracture zone and an uplift of the isotherms midway between the fracture zones a preferred model based on newly accumulated data suggests that the highs may be a result of upwelling within the asthenosphere underlying the ridge axis where a linear region of instability develops at the base of a partial melting zone less dense material at depths of 40 80 km may penetrate the overlying material at discrete and ideally evenly spaced points harmonically spaced upwellings should generate excess magma at discrete points on the ridge axis creating volcanic highs the spacing between highs may depend upon the continuity and thickness of the deep axial instability layer the temperature and viscosity ratios within the magma and the distance between major transform faults instability of the diapirs may stimulate the migration of the highs through time

resolution of the scripps noaa marine gravity field from satellite altimetry the july 1995 declassification of the entire geosat gm satellite altimeter data set enabled a joint scripps noaa effort to compute a new version 7 2 marine gravity field on a 2 minute grid this gravity field covers the world oceans between 72 n and 72 s and is derived from a combination of ers 1 and geosat gm and erm data an earlier noaa geosat only gravity field solution was confined to the southern latitudes because the 1992 declassification was limited to gm data south of 30 s a simple coherence analysis between accurately navigated ship gravity profiles and comparable gravity profiles obtained from the gravity grids reveals that the scripps noaa gravity field is coherent with ship gravity down to 23 30 km this slight increase in resolution over the previous noaa geosat only gravity field short wavelength resolution of 26 30 km implies that the increased spatial coverage provided by the ers i altimeter when combined with geosat improves the solution coherence analyses between satellite gravity and ship topography and ship gravity and ship topography show that even shorter wavelength gravity anomalies 13 km are present in sea surface measurements made by ship even so the scripps noaa marine gravity field does an excellent job of resolving most of the short wavelength gravity anomalies covering the world oceans

asymmetric sea floor spreading caused by ridge plume interactions crustal accretion at mid ocean ridges is generally modelled as a symmetric process regional analyses however often show either small scale asymmetries which vary rapidly between individual spreading corridors or large scale asymmetries represented by consistent excess accretion on one of the two separating plates over geological time spans1 2 3 4 5 6 in neither case is the origin of the asymmetry well understood here we present a comprehensive analysis of the asymmetry of crustal accretion over the past 83 myr based on a set of self consistent digital isochrons7 and models of absolute plate motion8 9 we find that deficits in crustal accretion occur mainly on ridge flanks overlying one or several hotspots we therefore propose that asymmetric accretion is caused by ridge propagation towards mantle plumes or minor ridge jumps sustained by asthenospheric flow10 11 between ridges and plumes quantifying the asymmetry of crustal accretion provides a complementary approach to that based on geochemical12 and other geophysical data13 14 in helping to unravel how mantle plumes and mid ocean ridges are linked through mantle convection processes

characteristics and evolution of the segmentation of the mid atlantic ridge between 20 n and 24 n during the last 10 million years abstract high resolution bathymetry and geophysical data collected along the slow spreading axis and flanks of the mid atlantic ridge between 20 n and 24 n reveal the characteristics and history of different wavelengths of segmentation during the last 10 m y the bathymetric data exhibit a morphotectonic pattern dominated by ridge normal and oblique bathymetric lows that partition the ridge flanks into rhomb shaped areas of relatively high elevation at least four different types of oblique bathymetric lows have been identified which represent the off axis traces of axial discontinuities and suggest a complex and ongoing evolution of ridge axis segmentation one group of oblique structures is represented by two deep ridge normal depressions with typical fracture zone characteristics that are connected to the present active transform by oblique depressions near the ridge axis these oblique traces correspond to the southward shift of axial discontinuities associated with the propagation of the ridge axis while maintaining a constant offset of the latter two other types of oblique structures correspond to elongate bathymetric lows and oblique alignments of ridge parallel bathymetric lows symmetric about the ridge axis both types of oblique structures frequently change their orientation from normal to subparallel to the ridge axis and appear to merge and diverge off axis these oblique depressions are characterized by positive filtered mantle bouguer anomalies high magnetizations complex magnetic anomaly patterns and possible exposure of mantle lithologies the ridge segments defined by these oblique depressions lengthen or shorten along the ridge axis with propagation rates varying from 0 to 25 km m y 1 the last and smallest discontinuities observed in this area correspond to small ridge axis offsets and off axis traces identified by alignments of the terminations of abutting abyssal hills the ridge flank morphotectonic patterns produced by the evolution of these elementary segments of accretion may represent temporally variable upwelling volumes of melt the centres of the rhomb shaped areas correspond to maximum crust production and thin lithosphere and the discontinuities correspond to a thick lithosphere with very thin crust and possible outcrops of peridotites we propose a model which accounts for the punctuated injection of magma and the evolution of elementary segments of accretion over periods of several million years

mantle downwelling beneath the australian antarctic discordance zone evidence from geoid height versus topography the australian antar ctic discordance zone aad is an anomalously deep and rough segment of the southeast indian ridge between 120 and 128 e a large negative deeper than predicted depth anomaly is centered on the discordance and a geoid low is evident upon removal of a low order geoid model and the geoid height age relation we investigate two models that may explain these anomalies a deficiency in ridge axis magma supply that produces thin oceanic crust i e shallow airy compensation and a downwelling and or cooler mantle beneath the aad that results in deeper convective type compensation to distinguish between these models we have calculated the ratio of geoid height to topography from the slope of a best line fit by functional analysis i e non biased linear regression a method that minimizes both geoid height and topography residuals geoid topography ratios of 2 1 0 9 m km for the entire study area 380 60 s 105 140 e 2 3 1 8 m km for a subset comprising crust 25 ma and 2 7 2 0 m km for a smaller area centered on the aad were obtained these ratios are significantly larger than predicted for thin oceanic crust 0 4 m km and 2 7 m kin is consistent with downwelling convection beneath young lithosphere average compensation depths of 27 29 and 34 km respectively estimated from these ratios suggest a mantle structure that deepens towards the aad the deepest compensation 34 km of the aad is below the average depth of the base of the young lithosphere 30 kin and a downwelling of asthenospheric material is implied the observed geoid height age slope over the discordance is unusually gradual at 0 133 m m y we calculate that an upper mantle 170 c cooler and 0 02 g cm 3 denser than normal can explain the shallow slope unusually fast shear velocities in the upper 200 km of mantle beneath the discordance and major element geochemical trends consistent with small amounts of melting at shallow depths provide strong evidence for cooler temperatures beneath the aad

three dimensional mantle convection beneath a segmented spreading center implications for along axis variations in crustal thickness and gravity segmentation and along axis variations within individual segments indicate the inherently three dimensional nature of mantle up welling and melting beneath oceanic spreading centers numerical convection experiments are used to explore the effects of local buoyancy forces on upwelling and melt production beneath a segmented spreading center the experiments are conducted in a region consisting of a thermally defined rigid lithosphere and a uniform viscosity asthenosphere overlying a higher viscosity mantle half space a periodic plate boundary geometry is imposed consisting of spreading segments and transform offsets buoyancy forces are caused by thermal expansion and the compositional density reduction due to the extraction of partial melt the relative magnitudes of the buoyant and plate driven components of mantle flow are controlled by the spreading rate and mantle viscosity with buoyant flow more important at lower spreading rates and viscosities buoyant flow beneath the spreading axis amplifies along axis variations in upwelling near a ridge transform intersection and distributes the variations along the entire spreading axis buoyant flow may thus be responsible for the more three dimensional character of slow spreading centers away from the spreading axis thermal buoyancy drives convective rolls that align with the direction of plate motion and which have an along axis wavelength controlled by the prescribed thickness of the asthenosphere however the position and stability of rolls are influenced by the segmentation geometry in cases where the spreading center geometry does not allow a stable configuration of rolls the flow is time dependent along axis variations in upwelling cause variations in melt production which imply large variations in crustal thickness that dominate the surface gravity signal the crustal thickness distributions implied by these numerical experiments produce bulls eye shaped negative mantle bouguer anomalies centered over spreading segments as observed at several spreading centers the amplitude of the anomaly increases with decreasing spreading rate

segmentation and crustal structure of the western mid atlantic ridge flank 25 25 27 10 n and 0 29 m y we conducted a detailed geological geophysical survey of the west flank of the mid atlantic ridge between 25 25 n and 27 10 n and from the ridge axis out to 29 ma crust acquiring hydrosweep multibeam bathymetry hawaii mr1 sidescan sonar imagery gravity magnetics and single channel seismic reflection profiles the survey covered all or part of nine spreading segments bounded by mostly nontransform right stepping discontinuities which are subparallel to flow lines but which migrated independently of one another some discontinuities alternated between small right and left stepping offsets or exhibited zero offset for up to 3 4 m y despite these changes the spreading segments have been long lived and extend 20 m y or more across isochrons a large shift 9 in relative plate motion about 24 22 ma caused significant changes in segmentation pattern the nature of this plate boundary response together with the persistence of segments through periods of zero offset at their bounding discontinuities suggest that the position and longevity of segments are controlled primarily by the subaxial position of buoyant mantle diapirs or focused zones of rising melt within segments there are distinct differences in seafloor depth morphology residual mantle bouguer gravity anomaly and apparent crustal thickness between inside corner and outside corner crust this demands fundamentally asymmetric crustal accretion and extension across the ridge axis which we attribute to low angle detachment faulting near segment ends cyclic variations in residual gravity over the crossisochron run of segments also suggest crustal thickness changes of at least 1 2 km every 2 3 m y these are interpreted to be caused by episodes of magmatic versus relatively amagmatic extension controlled by retention and quasiperiodic release of melt from the upwelling mantle detachment faulting appears to be especially effective in exhuming lower crust to upper mantle at inside corners during relatively amagmatic episodes creating crustal domes analogous to turtleback metamorphic core complexes that are formed by low angle detachment faulting in subaerial extensional environments

segmentation of mid ocean ridges with an axial valley induced by small scale mantle convection the small scale segmentation of mid ocean ridges with an axial rise has been modelled by considering each ridge segment as a giant crack in the lithosphere with a tip propagating along the ridge axis1 2 for ridges with an axial valley however this type of model fails because the lithosphere is too thick to tear3 yet such ridges are clearly segmented as defined by morphology gravity and structure at the 50 100 km length scale the ridge offsets are large4 5 and vary dramatically with time this type of segmentation is commonly related to a three dimensional small scale mantle flow occurring in the partially molten asthenosphere below the ridge6 8 here we propose a model for segmentation in such ridges in which the convective flow below the ridge favours an asymmetrical breaking of the axial valley lithosphere this leads to the development and separation of ridge segments in a pattern that mimics the observed geometry and temporal evolution of the segmentation of most ridges with an axial valley if our model is correct it implies that coupling of oceanic lithosphere to small scale convection controls the dynamics of mid ocean ridges

three dimensional gravity study of the mid atlantic ridge evolution of the segmentation between 28 and 29 n during the last 10 m y during the high resolution survey sara segmentation ancienne de la ride atlantique sea beam bathymetry magnetic gravity and seismic reflection data were collected on the flanks of the mid atlantic ridge between 28 and 29 n this survey was designed to provide off axis information up to approximately 10 m y and to complement a detailed on axis survey carried out in 1988 and 1989 sempere et al 1990 between the atlantis and kane fracture zones a previous gravity study had revealed the existence of bulls eye shaped gravity lows centered on the axial segments and gravity highs centered on the non transform discontinuities lin et al 1990 we carried out a three dimensional calculation of the mantle bouguer anomaly in order to investigate if the axial pattern of circular anomaly lows can be followed on the flanks of the spreading center the off axis gravity anomalies are characterized by anomaly lows over the centers of the segments and anomaly highs over the discontinuities after correcting for the gravity effect of lithospheric cooling away from the ridge the segmentation configuration determined from gravity data appears to be very similar to that deduced from bathymetry off axis bathymetry is characterized by southward pointing v shaped basins which indicate that the traces of the nontransform discontinuities do not follow plate motion flow lines and that each ridge segment advances and retreats continuously however the gravity trace of the discontinuities is always slightly offset with respect to the bathymetric lows northward on the western flank and southward on the eastern flank the sense of offset between the gravity and bathymetric traces appears to be related to the right stepping nature of the axial discontinuities

high resolution inversion for south atlantic plate kinematics using joint altimeter and magnetic anomaly data we present an inversion for plate kinematics that solves for finite rotation parameters using fracture zone fz and magnetic anomaly location data jointly we define misfit functions that incorporate properties unique to each data type in particular the fz misfit function does not depend upon alignment of conjugate fz traces in the same way as magnetic lineations under the finite rotations this property is useful for fz locations in which the signals on conjugate sides of the ridge may include systematic differences or where data from one side of the ridge are sparse or missing formal error bounds estimated for the pole parameters show that the magnetic and fz data are complementary in their information content error bounds computed for the joint inversion are substantially smaller than for either the fz or magnetics data used separately indicating that simultaneous use of the data in an inversion is crucial we apply this method to seasat altimeter data and magnetic anomaly picks in the south atlantic we solve for finite poles corresponding to magnetic anomalies 5 6 8 13 21 22 25 30 33 33r and 34 these define a smooth path over the past 84 m y indicating that the plate motions have not been as erratic as found previously

the paradox of the axial profile isostatic compensation along the axis of the mid atlantic ridge mantle bouguer gravity anomalies mba and bathymetry on three profiles covering more than 1000 km along the axis of the mid atlantic ridge mar are highly correlated suggesting that along axis topographic relief is locally compensated by variations in crustal thickness and or mantle density structure the quantitative relationship between topography and gravity on these profiles could be explained by the flexure of a thin narrow elastic strip representing the response to isostatic loads of an inner rift valley isolated from the rest of the plate by weak bounding normal faults the paradox is that across axis profiles show that the median valley is an uncompensated feature apparently created by a dynamic mechanism new extensive off axis coverage of the mar from 31 to 36 s shows that the high correlation does not persist outside the axial zone we suggest that the on axis correlation exists because the mechanism creating the median valley is controlled by the mantle thermal structure and along axis variations in crustal thickness both of which contribute to the mba if the mechanism is extension of a brittle ductile lithosphere the critical parameter controlling topographic relief is the thickness of the relatively stiff mantle layer immediately beneath the crust at the ridge axis a three dimensional 3 d thermal model incorporating passive mantle flow hydrothermal circulation the plate boundary geometry and variable magmatic heating associated with observed variations in crustal thickness predicts variations in the thickness of the stiff mantle layer that correlate with the observed axial topography we model the expected topography using cross axis sections of the 3 d thermal model and a 2 d finite element model of an extending lithosphere that incorporates temperature and strain dependent rheology as well as the flexural response of a thickening plate the predicted topographic signal produced by the combined dynamic and isostatic effects matches the amplitude of the observed axial bathymetry at the 18 mm yr 1 spreading half rate in our south atlantic survey the 2 d models also reproduce the observed rapid transition along the axis between rift valley no rift valley morphology we conclude that extensional forces acting on segmented oceanic lithosphere with varying rates of crustal production produce the highly variable morphology of the southern mar and the along axis correlation between mba and bathymetry

high resolution mean sea surface computed with altimeter data of ers 1 geodetic mission and topex poseidon the remote sensing satellite ers 1 launched in 1991 to study the earth s environment was placed on a geodetic 168 day repeat orbit between 1994 april and 1995 march to map through altimetric measurements the gravity field over the whole oceanic domain with a resolution of 8 km at the equator in both along track and cross track directions we have analysed the precise altimeter data of the geodetic mission and by also using one year of topex poseidon altimeter data we have computed a global high resolution mean sea surface the various steps involved in pre processing the ers 1 data consisted of correcting the data for environmental factors editing and reducing through crossover analyses the radial orbit error which directly affects sea surface height measurements for this purpose we adjusted sinusoids at 1 and 2 cycle rev 1 along the ers 1 profiles in order to minimize crossover differences between ers 1 and yearly averaged topex poseidon profiles in effect the orbit of topex poseidon is very accurately determined within 2 3 cm for the radial component so topex poseidon altimeter profiles can serve as a reference to reduce the ers 1 radial orbit error the ers 1 residual orbit error was further reduced through a second crossover analysis between all ascending and descending profiles of the geodetic mission the along track ers 1 and topex poseidon data were then interpolated over the whole oceanic domain on a regular grid of 1 16 1 16 size the mapping of the gridded sea surface heights reveals the very fine structure of the marine geoid up until now unknown at a global scale this new data set will be most useful for marine geophysical and tectonic investigations

thermomechanical evolution of oceanic fracture zones a fracture zone fz model is constructed from existing models of the thermal and mechanical evolution of the oceanic lithosphere as the lithosphere cools by conduction thermal and mechanical boundary layers develop and increase in thickness as age 1 2 surface expressions of this development such as topography and deflection of the vertical i e gravity field are most apparent along major fz s because of the sharp age contrast a simple model including the effects of lateral heat transport but with no elastic layer predicts that variations in seafloor depth and deflection of the vertical will become increasingly smooth and ultimately disappear as the fz ages observations however show that both the fz topography and the deflection of the vertical remain sharp as the fz evolves these two observations as well as the observed asymmetry in deflection of the vertical profiles across the udintsev romanche and mendocino fz s are explained by including a continuous elastic layer in the model the asymmetry in deflection of the vertical is a consequence of elastic thickness variations across the fz modeling also shows that the evolution of the fz topography is extremely sensitive to the initial thermal structure near the ridge transform intersection model geoid steps and their development with age are used to access techniques for measuring geoid offsets across fz s reasonable step estimation techniques will underestimate the overall step amplitude by up to 50 this implies that abnormally thin thermal boundary layers derived from studies of geoid height versus age are not required by the data

earthquake slip vectors and estimates of present day plate motions although horizontal slip directions determined for earthquakes along transform faults and shallow subduction thrust faults are commonly assumed to parallel the direction of motion between two plates systematic biases of earthquake slip vectors relative to the local long term average plate direction can be introduced by a variety of causes including deformation of forearcs above subduction zones unmodeled lateral velocity heterogeneities near earthquake sources and unrecognized changes in plate velocities within the averaging interval of the kinematic model given that earthquake slip vectors comprise nearly two thirds of the data used to derive the nuvel 1 model for present day global plate velocities it is important to determine the extent to which potential systematic biases of earthquake slip vectors could degrade estimates of plate velocities given by nuvel 1 here two alternative models for global plate velocities are derived and compared to nuvel 1 the first model nuvel sz is derived from a data set that omits the 240 subduction zone sup vectors in the nuvel 1 data set but includes all remaining 882 kinematic data the velocities predicted by nuvel sz and nuvel 1 differ little with maximum differences of only 1 mm yr 1 and 2 along all plate boundaries except those surrounding the caribbean plate the second model nuvel g is derived from a data set that omits all 724 earthquake sup vectors but includes the remaining 277 spreading rates and 121 transform fault azimuths velocities predicted by nuvel g and nuvel 1 also differ little with maximum differences of 2 mm yr 1 and 4 except for the caribbean plate boundaries these results indicate that even in the unlikely event that all earthquake slip vectors in the nuvel 1 data set are unreliable recorders of long term plate directions the slip vectors do not significantly degrade the model it thus appears that the decision to use or not to use earthquake slip vectors to derive a global plate motion model has little practical effect on the ability to derive an accurate description of present day plate velocities little emphasis is placed on the changes in estimates of caribbean plate velocities that occur upon exclusion of earthquake slip vectors nearly all of the data from caribbean plate boundaries are suspect and it is unlikely that the nuvel sz or nuvel g models represent a significantly improved description of caribbean plate velocities comparison of 677 slip vectors from transform faults to directions predicted by nuvel g shows statistically insignificant differences along 12 of 15 spreading centers the good agreement between the nearly instantaneous estimates of slip directions provided by slip vectors and the longer term average directions given by nuvel g suggest that in general transform fault slip vectors parallel the longer term average directions along transform faults interestingly a statistically significant difference between slip vectors from right slipping and left slipping transform faults is noted for nearly all spreading centers with slip vectors along right slipping and left slipping faults rotated clockwise and counterclockwise respectively from the predicted direction the cause of this bias is unknown but may be related to biases introduced by unmodeled lateral heterogeneities in the mantle near transform faults

the relationship between buoyant mantle flow melt migration and gravity bull s eyes at the mid atlantic ridge between 33 n and 35 n abstract crustal thickness and mantle bouguer anomaly mba patterns along the mid atlantic ridge mar between the oceanographer and hayes transforms 33 35 n are predicted from a three dimensional mantle flow model driven by plate spreading and mantle density variations assuming a two or three dimensional model of melt extraction our two dimensional melt extraction model focuses melt to the nearest ridge segment by assuming that melt moves only vertically and parallel to the spreading direction in contrast the three dimensional model assumes melt migrates vertically to the top of the melting region and then up the local slope of this surface the shape of this surface which is controlled by enhanced conductive cooling at ridge offsets directs melt toward the center of segments our results suggest that buoyant mantle flow driven by temperature variations and a small retained melt fraction

temporal variations of the segmentation of slow to intermediate spreading mid ocean ridges 1 synoptic observations based on satellite altimetry data 1 the high resolution geoid and gravity maps derived from ers 1 and geosat satellite geodetic missions reveal a set of small scale lineations on the flanks of slow to intermediate spreading mid ocean ridges assuming that these lineations reflect the variations in crustal structure induced by mid ocean ridge axial discontinuities we use them to investigate how the discontinuities and the segments they bound appear migrate and disappear we provide a synoptic description of the main characteristics of the crustal structure variations as well as their evolution in time over the flanks of the mid atlantic indian and pacific antarctic ridges the second order segment length does not appear to vary with the spreading rate for the slow to intermediate spreading ridges investigated here the amplitude of the gravity signal associated with off axis discontinuity traces increases with the obliquity of the ridge to spreading and decreases with spreading rate and with the proximity of a ridge section to a hot spot the patterns of the gravity lineations appear to be very homogeneous over 500 to 1000 km large corridors bounded by large fracture zones far from hot spots corridors are characterized either by segments bounded by discontinuities migrating back and forth along the axis implying a lifetime of 10 30 myr for the segments or by segments and discontinuities very stable in space and time surviving for 40 50 myr closer to hot spots the segmentation is affected in two ways first segments tend to migrate along axis away from hot spots or toward cold spots second asymmetric spreading tends to keep sections of ridges closer to hot spots than normal spreading would these observations support the hypothesis that ridge segmentation and its evolution are controlled by mantle dynamics our analysis provides observational constraints for further models of crustal production along ridges which are presented in the companion paper by rabinowicz and briais 2002

models of ocean ridge lithospheric deformation dependence on crustal thickness spreading rate and segmentation we use three dimensional 3 d temperature and rheology models to investigate the effects of crustal thickness and ridge segmentation on mid ocean ridge lithospheric structure we find that crustal thickness variations associated with focused magma accretion at a slow spreading ridge segment can cause significant along axis gradients in crustal temperature and a pita pocket shaped weak zone in the lower crust that decouples brittle upper crust from upper mantle in contrast fast spreading ridge segments with little crustal thickness variation are found to be uniformly weak along axis the overthickened crust produced by ridge hotspot interaction alters the heat balance between magma emplacement and hydrothermal cooling creating an extremely weak lithosphere in a steady state model we apply a simple two dimensional 2 d cyclic faulting model to across axis sections of temperature and rheological structure slow spreading segments are predicted to have relatively large along axis variations in fault height and spacing in contrast fast spreading and hotspot affected segments are predicted to have much smaller variations in faulting styles except in the immediate vicinity of a major transform fault the presence or lack of a steady state axial rift valley is also predicted to depend on crustal thickness spreading rate and the ability of rift bounding normal faults to propagate along axis for fast spreading and hotspot affected segments the sizes of axial rift valleys are predicted to be small and may not be distinguishable from those of neovolcanic or magma chamber supported isostatic features this is in contrast to most slow spreading segments where axial rift valleys are predicted to be a dominant morphological feature for a slow spreading segment of great length however the axial rift valley size is predicted to diminish toward the segment center because large faults formed at segment distal ends cannot propagate through a large decoupling zone created by locally thickened crust at the segment center

evidence for variable upper mantle temperature and crustal thickness in and near the australian antarctic discordance abstract the southeast indian ridge seir in and near the australian antarctic discordance aad exhibits at a constant spreading rate almost the full range of the many geophysical and geochemical parameters characteristic of the slow mid atlantic ridge and fast east pacific rise we used satellite derived gravity data in combination with seamarc ii bathymetry in and near the aad to examine regional density variations in the upper mantle beneath the aad through three dimensional gravity analysis we found that at least two end member models satisfy the gravity observations regional crustal thickness variations of at least 3 km along the seir near the aad or a temperature anomaly of the order of 150 c in the upper mantle beneath the seir these new observations combined with other geophysical and geochemical characteristics of the australian antarctic discordance provide further evidence that the temperature structure of a mid ocean ridge is a controlling factor in addition to spreading rate in the crustal accretionary process numerical models of mantle flow beneath mid ocean ridges offer one means of investigating the dynamic effect of a variable upper mantle temperature on the accretionary process our results indicate that temperature is important especially at intermediate and slower spreading rates where thermal effects can dominate mantle flow beneath a mid ocean ridge and result in increasing crustal production with decreasing spreading rate at the constant intermediate spreading rate of 37 mm yr characteristic of the seir in and near the aad our numerical models show that significant crustal thinning 2 4 km can occur with relatively small variations in upper mantle temperature all else being equal thus combined with our end member gravity models these observations and results suggest that both anomalously cool upper mantle and thin crust exist beneath the aad

rift valley no rift valley transition at mid ocean ridges a simple mechanical model for the axial topographic and gravity variations with spreading rate is developed the model emphasizes the rapid transition from one mechanism to another with changing spreading rate our model has a brittle layer over a ductile layer which are being pulled by forces away from the ridge axis there exists a potential failure zone of width lf where the stresses within the brittle plate exceeds its yield strength if the yield strength of the brittle plate is nearly constant this width is largely determined by the strain rate and the viscosity of the strong mantle the high ductility of oceanic crust at high temperatures relative to that of mantle at the same temperatures leads to the concept of a decoupling region a critical feature in this model the size of the decoupling region ld is determined by the thickness of the oceanic crust and temperatures within the crust if the decoupling region is small ld lf as shown to be the case for fast spreading ridges then each of the two layers responds separately to the horizontal stretching the brittle plate will be broken only at the axis and any buoyant force of the region would lift it up to form the rise crest high as a result of local isostasy the changing size of the decoupling region gives rise to the transition from atlantic type axial topography to pacific type as well as the transition from a high amplitude highly variable gravity signature to a low amplitude very uniform gravity signature parameters can be chosen so that such a transition occurs at a half rate of 35 mm yr which is the critical half rate suggested by both topographic and gravity data

long lived mid ocean ridge segmentation of the pacific antarctic ridge and the southeast indian ridge the pacific antarctic and the southeast indian ridges are examined for evidence of long lived tens of millions of years tectonic segmentation of the mid ocean ridge by examining a variety of geophysical parameters similar to those previously documented in the south atlantic kane and hayes 1992 these parameters include subsidence rate of the oceanic crust individual flanks and asymmetry zero age depth and geoid height decrease with age geoid rate the variability in these parameters along strike of the mid ocean ridge is systematic and serves to define a large scale ridge segmentation of the order of hundreds of kilometers the southeast indian ridge exhibits long lived segmentation in all parameters examined and is subdivided into five flow line corridors with boundaries between corridors occurring at or near major fracture zones while the pacific antarctic ridge exhibits conspicuous asymmetric crustal subsidence evidence of segmentation exists primarily in the along strike geoid rate the observations in these two areas are found to be internally inconsistent with simple thermal conduction models for oceanic crust strongly suggesting that factors other than cooling and isostasy are influencing the creation and subsequent modification of the oceanic crust and lithosphere

crenulated seafloor evidence for spreading rate dependent structure of mantle upwelling and melting beneath a mid oceanic spreading center abstract the availability of dense geosat gm satellite altimetry data south of 30 s 1 provides the first opportunity to compare at 25 km resolution seafloor generated at ridges spreading with a wide range of spreading rates the axis of the slow spreading southwest indian ridge as well as the seafloor that it has generated has a distinctly crenulated appearance this crenulated seafloor shows regularly spaced gravity lineations in the spreading direction crenulations that record stationary centers of mantle upwelling over the life of the southwest indian ridge a similar pattern of crenulation is seen along sections of the slow spreading northern mid atlantic ridge but along the faster intermediate spreading southern mid atlantic ridge while regularly spaced highs are apparent along sections of the spreading axis stationary crenulation does not persist off the axis this pattern implies that the faster spreading southern mid atlantic ridge has variable centers of time dependent upwelling that migrate along the ridge crenulations are not visible near the axis of the fast spreading pacific antarctic ridge crenulation amplitudes decrease strongly with increasing spreading rate the absence of crenulation at high spreading rate could be due to a 3d 2d transition in flow structure or the effect of mechanisms to redistribute crust along axis while crenulations are not ubiquitous even at slowly spreading ridges where they exist they provide a time record to observe the characteristics of mid ocean ridge upwelling melting and melt migration not only beneath active spreading centers but through a much longer spreading history