three dimensional tomographic velocity structure of upper crust coaxial segment juan de fuca ridge implications for on axis evolution and hydrothermal circulation three dimensional models of compressional velocity and azimuthal anisotropy from tomographic inversions using 23 564 ocean bottom seismometer p wave arrivals define systematic lateral variations in seismic structure of the coaxial segment of the juan de fuca ridge jdfr over much of the segment the across axis structure is roughly axisymmetric characterized by a progressive increase in dike velocities moving away from the ridge axis this trend is most apparent in the basal dikes where on axis velocities are about 800 m s slower than those measured elsewhere within the rift valley the on axis sheeted dikes also exhibit ridge oriented azimuthal anisotropy with a peak to peak amplitude of about 600 m s outboard of the rift valley beneath ridge flanks with fault scarps velocities in the upper 1500 m of crust are reduced the maximum amplitude of this anomaly is about 700 m s located near the top of the sheeted dikes variations in the three dimensional velocity model are believed to reflect changes in crustal porosity from which we infer an axisymmetric porosity model for seismic layer 2 of the coaxial segment as the crust ages the evolution of layer 2 porosity could occur in the following way 1 the porosity of zero age on axis dikes is set at formation by the contraction of molten material 2 hydrothermal alteration fills pore spaces as the dikes move away from the center of the axial valley and 3 normal faulting on the ridge flank scarps opens fractures and increases porosity of the upper dikes as they move off axis at the north end of the segment dike velocities are several hundred meters per second slower on average and the across axis structure is lost the transition from a coherent aligned seismic structure to a less distinct pattern with reduced velocities may represent a shift from magmatic to amagmatic extension moving away from the cobb hotspot on the ridge axis the porosity structure we have derived for the coaxial segment suggests an alternative to the usual hydrothermal circulation model of cross axis convection cells a circulation model with along axis convection cells located entirely within the axial valley appears to be more compatible with our data

upper crustal seismic velocity structure and microearthquake depths at the endeavour segment juan de fuca ridge 1 we present the results of a study to invert microearthquake and explosive shot data from the endeavour segment of the intermediate spreading juan de fuca ridge the average isotropic p wave velocity structure derived from the shot data in the uppermost 1 5 km of the oceanic crust is characterized by an increase with age of 8 from the axis to at least 0 5 ma that is attributed to the sealing of layer 2a porosity by hydrothermal processes superimposed on this variation are axis parallel 2 km wide alternating bands of high and low velocity with a peak to peak variation of 5 12 high and low velocities away from the axis correspond to bathymetric trenches and ridges respectively and are likely due to variations in layer 2a thickness p wave azimuthal anisotropy is present in the data that is best fit with a model of 9 anisotropy at 750 m depth decreasing to 1 at 3 km depth and is likely due to the preferential alignment of vertical cracks and fissures in the along axis direction anisotropy and velocity heterogeneity are coupled anisotropy alone may explain the form but not the magnitude of the axis parallel bands there are strong trade offs between the hypocentral depths of microearthquakes and the p and s wave velocity structures changing the mean hypocentral depth by up to 0 5 km leads to only modest increases in the travel time rms but the resulting velocity models appear more feasible when the earthquakes are forced deeper than when they are forced shallower

vertical magnetic structure of ocean crust determined from near bottom magnetic field measurements consensus on the source of marine magnetic anomaly stripes has remained one of the most elusive aspects of the vine and matthews seafloor spreading model which has been of such fundamental importance to global plate tectonics the first attempt to measure the vertical magnetic structure of oceanic crust exposed at a scarp face using near bottom magnetic measurements is reported along with the theoretical basis for the reduction and analysis of these measurements magnetometers mounted to deep tow and submersible platforms were used to determine the horizontal and vertical variation of crustal magnetization of the upper 2 km of young oceanic crust exposed by a steep submarine escarpment at the blanco fracture zone in the northeast pacific ocean results show a large amplitude magnetic anomaly is consistently found at the extrusive to intrusive transition indicating that the extrusives are strongly magnetized compared to the weakly magnetized intrusive crust profiles also show the systematic presence of more than one polarity with depth indicating the possibility of dipping polarity boundaries within the extrusives forward models find that a dipping polarity boundary model correctly predicts the amplitude and shape of the overlying sea surface anomaly while a block model with vertical polarity boundaries overestimates the amplitude these in situ results are used to directly assess the contribution made by the extrusive section to the overlying magnetic anomaly signal in essence to define the magnetic source layer it is concluded that the extrusive lavas constitute the main magnetic source for magnetic anomalies in oceanic crust less than 2 m y in age

north east greenland shelf north of 79 n results of a reflection seismic experiment in sea ice abstract the first multichannel seismic reflection measurements of the north east greenland shelf north of 79 n were made as part of the r v polarstern cruise ark v 1988 this area is perennially covered by sea ice the sparse seismic data show a very deep rift basin filled with up to 10 km of high velocity rocks interpreted as carboniferous and permian sediments conspicuous structures interpreted as marginal reefs are associated with the western flank of the rift basin

upper crustal seismic structure of the slow spreading mid atlantic ridge 35 n constraints on volcanic emplacement processes 1 the upper crustal seismic structure of the slow spreading mid atlantic ridge is studied using a genetic algorithm based waveform inversion of multichannel streamer data four single ship multichannel profiles from 35 n are analyzed one in the rift valley and three in the rift mountains along 0 7 1 6 and 1 9 ma crust a layer 2a horizon is continuously imaged along three profiles and is associated with a sharp velocity increase from extrusives to dikes its depth and regularity in the rift valley indicate that most of the extrusive section is built on the inner valley floor through a pattern of deposition and fault bounded uplift into the rift mountains its variability along one line however shows that this process is disrupted during tectonically dominated periods a thickening of layer 2a toward the oceanographer fracture zone may be the result of along axis magma transport the interval of rapid velocity increase at the base of layer 2a thins with age a possible response to enhanced hydrothermal mineralization within the zone of mixed dikes and extrusives transition zone 200 m and off axis layer 2a thicknesses 350 600 m are similar to those at other spreading centers this indicates that equivalent extrusive volumes are produced at all spreading rates along a relatively narrow zone of dike emplacement however differences in on axis layer 2a thickness between this area and fast spreading ridges suggest that the exact pattern of thickening varies between spreading regimes relative to fast spreading ridges the moderate velocity increase with age recorded in the upper crust from 2 3 to 2 7 km s 1 within 2 myr may be due to a less active hydrothermal system and hence slower porosity reduction

detailed geomorphology and neotectonics of the endeavour segment juan de fuca ridge new results from seabeam swath mapping a 36 km wide corridor of seabeam bathymetry was collected on the juan de fuca ridge north of 48 00 n in may 1983 these data merged with previous seabeam bathymetric compilations for the cobb offset as well as with sea marc i and sea marc ii side scan sonar seismic reflection magnetometer and deep tow photography data have been used to identify morphological variations along the spreading center of the endeavour segment of the juan de fuca ridge between the cobb offset and sovanco fracture zone the endeavour segment consists of two ridge sections separated by the 13 km wide endeavour offset which formed when spreading jumped from middle valley to west valley 200 000 yr ago the locus of spreading occurs along the entire endeavour segment as a narrow 1 2 km wide shallow 10 30 m deep inner rift which is superimposed upon larger structures that range from broad 5 10 km wide deep as much as 3 000 m median valleys to narrow 5 km wide shallow 2 100 m volcanic ridges an abrupt discontinuity in the tectonic fabric occurs at about 48 05 n south of this latitude spreading occurs within an axial high endeavour ridge which gradually deepens to the south south endeavour valley as it approaches the zone of overlap with the propagating rift at the cobb offset north of 48 05 n the axis is characterized by broad deep median valleys north endeavour valley and west valley these large scale morphological variations reflect a complex interplay of the spreading center with seamount volcanism waning magma supply at the distal ends of ridge axes thermal contrasts across ridge offsets propagating rifts and an unstable triple junction the pronounced axial deep of west valley is interpreted as reflecting the youth of the structure 200 000 yr and may be primarily due to collapse with little extension and magmatism the rift axis jump to west valley fortuitously isolated endeavour seamount the youngest edifice in the heck seamount chain from the pacific plate and trapped it between a pair of overlapping spreading centers at the newly formed endeavour offset development of a remarkably simple overlapping conjugate rift pair at endeavour offset has resulted in spite of the complex and highly variable crustal thickness magma supply and pre existing structural grain rifting in south west valley has caused pre existing topography generated at the endeavour ridge to be destroyed by subsidence and burial thereby creating the apparent discontinuity in morphology north of 48 05 n

geology of a vigorous hydrothermal system on the endeavour segment juan de fuca ridge a high precision high resolution geologic map explicitly documents relationships between tectonic features and large steep sided sulfide sulfate silica deposits in the vigorously venting endeavour hydrothermal field near the northern end of the juan de fuca ridge water depth in the vent field varies from 2220 to 2200 m location of the most massive sulfide structures appears to be controlled by intersections of ridge parallel normal faults and other fracture fissure sets that trend oblique to and perpendicular to the overall structural fabric of the axial valley the fractured basaltic substrate is primarily composed of well weathered pillow and lobate flows as presently mapped the field is about 200 by 400 m on a side and contains at least 15 large 1000 m3 sulfide edifices and many tens of smaller commonly inactive sulfide structures the larger sulfide structures are also the most vigorously venting features in the field they are commonly more than 30 m in diameter and up to 20 m in height actively venting sulfide structures in the northern portion of the field stand higher and are more massive than active structures in the southern portion of the field which tend to be slightly to distinctly smaller maximum venting temperatures of 375 c are associated with the smaller structures in the southeastern portion of the field highest temperature venting fluids from the more massive structures in the northern portion of the field are consistently 20 30 c lower hydrothermal output from individual active sulfide features varies from no flow in the lower third of the edifice to vigorous output from fracture controlled black smoker activity near the top of the structures a different type of high temperature venting takes place from the upper sides of the structures in the form of overflow from fully exposed quiescent pools of buoyant 350 c vent water trapped beneath overhanging sulfide sulfate silica ledges or flanges these flanges are attached to the upper outer walls of the large sulfide edifices two types of diffuse venting in the endeavour field include a lower temperature 8 15 c output through colonies of large tubeworms and 25 50 c vent fluid that seems to percolate through the tops of overhanging flanges the large size and steep walled nature of the these structures evidently results from sustained venting in a mature hydrothermal system coupled with dual mineral depositional mechanisms involving vertical growth by accumulation of chimney sulfide debris and lateral growth by means of flange development

the accretion of oceanic crust by episodic sill intrusion seismic refraction data from mature oceanic regions show remarkable consistency in crustal thickness between sites yet the characteristics of different spreading centers show great variation this leads to the paradox that the similarity in mature oceanic crust suggests that the igneous processes in crustal accretion are the same across a range of spreading rates but that the differences in the spreading centers themselves imply the opposite we have developed a simple two dimensional model for crustal accretion at oceanic spreading centers by the episodic addition of sills at a high level within the crust as implied by recent geophysical studies of spreading centers using our model we have calculated velocity and temperature fields in the crust including the effects of hydrothermal cooling we show that oscillations in the temperature field due to episodicities of 200 20000 years have an effect localized to the region within 50 500 m of the intrusion by assuming that all latent heat and excess specific heat introduced by the injection of sills is convected away by hydrothermal circulation we estimate that the maximum heat transfer associated with hydrothermal cooling in the upper crust is of the order of 10 times greater than the conductive flux we have calculated seismic velocities associated with the model temperature field in the crust and we show that these are consistent with results from both seismic refraction and reflection experiments the model may be applied to both fast and slow spreading ridges and it accounts for the apparent increases in both the mean crustal temperature and the frequency of presence of melt with increasing spreading rate hence resolving the paradox we also show that the implications of this model are consistent with geological observations of layering in ophiolites

the distribution of geothermal fields on the juan de fuca ridge near bottom water temperatures were mapped along 400 km of the strike of the juan de fuca ridge as part of a combined sea marc seabeam experiment to image the variability of morphology and structure along a spreading center segment the water temperature data collected by a continuously towed thermistor chain in addition to salinity data indicate that there are four geothermal areas spaced at distances of 100 km from each other south of the cobb propagator and one field just to the north of the propagator on the endeavor ridge segment each thermal region is located above a morphological dome on the spreading center these domes are an average of 100 200 m shallower than the rest of the axis the structure of bottom water suggests that the geothermal regions are on average 20 km long and that the heat from these fields raises the temperature in the water column by a minimum of 0 06 sup 0 c up to 300 m above the bottom two simple models are used to estimate the heat flux associated with these features

shallow crustal structure beneath the juan de fuca ridge from 2 d seismic refraction tomography summary the formation of oceanic lithosphere along ocean ridges and the role that crustal magma chambers play in the accretionary process continue to be fundamental issues in plate tectonics to address these issues a multireceiver airgun oceanbottom seismograph refraction line was shot across the endeavour segment of the juan de fuca ridge near 48 n 129 w 2 d traveltime tomography applied to the data results in a three layer model for the upper crust layer 1 is 250 650m thick with u1 2 5 km s and vzvl 0 5 s layer 2 is 800 m thick u2 4 8 km s and vzuz 1 0 layer 1 and layer 2 probably represent the sequence of extrusives and the transition to layer 3 u3 5 8 km s l vzv3 0 5 s is associated with the extrusives dike complex transition an abrupt velocity transition between layer 1 and layer 2 may be a metamorphic front within the basalt pillows or it may be the depth to which surface fissures penetrate a low velocity anomaly velocities decreased by 0 45 km s exists beneath the ridge in layer 2 and upper layer 3 it is interpreted as a zone of increased fracture porosity and or permeability associated with axial hydrothermal circulation and correlates reasonably well with a sub axial reflector no evidence is found for the existence of a crustal magma chamber in the depth range of 1 5 3 5 km sub bottom however the velocity anomalies observed in layer 3 suggest that crustal temperatures in this layer are elevated by 150 250 c beneath and to the east of the ridge relative to temperatures west of the ridge

seismic attenuation structure beneath the juan de fuca ridge from tomographic inversion of amplitudes tomographic inversion of seismic amplitudes for attenuation structure has been applied to first arrival amplitudes from a seismic refraction profile across the juan de fuca ridge to further constrain a velocity structure previously obtained by travel time velocity tomography and to obtain a model for attenuation variations within the shallow crust amplitude analysis has identified an additional boundary defined by a decrease in vertical velocity gradient at 2 0 km subbottom which was not required by the travel time data this transition zone likely represents the oceanic layer 2 layer 3 boundary the analysis also identified zones of relatively high attenuation qe 20 100 at subbottom depths of 0 3 2 0 km beneath the ridge axis which correlate well with low velocity zones velocities reduced by 0 4 km s in the velocity model the subaxial zone is interpreted as a region of increased fracture porosity and or permeability associated with axial hydrothermal circulation within the extrusives and sheeted dike complex

volcanic episodicity and a non steady state rift valley along northeast pacific spreading centers evidence from sea marc i sea marc i side looking sonar images and sea beam bathymetry along a 400 km stretch of the juan de fuca ridge crest provide evidence that excessive extrusive volcanism periodically builds a crestal ridge along the axis of seafloor spreading an elongate summit depression esd or rift valley is commonly observed in the spine of this crestal ridge the crestal ridge volcanic landform has a distinctive shape that is recognized in bathymetric contours both along the spreading axis and at least up to 30 km away from the axis the landform has a plan form shape that resembles a side view of an archer s bow with the long dimension of the bow form parallel to the strike of the ridge in cross section the bow form is flat on top and has steep flanks these bow form shapes can be explained by magma that rises into the crust at a discrete center and flows laterally into belts of ridge parallel dikes similar to icelandic fissure eruptions both the variable dimensions of the esd along axis of the juan de fuca ridge and the relationship among volcanic flow morphologies within and beyond the esd suggest the four different segments of the juan de fuca ridge presented in detail here display different stages in a cycle of oceanic crust accretion this cycle includes episodes in which there is 1 extrusive volcanic construction which widen the crestal ridge prior to the collapse of the summit depression 2 collapse within the summit region of the crestal ridge to form an esd during a phase of volcanic inactivity and 3 renewed magmatism in the esd as its floor widens by extension and brittle fracture of the upper crust this episodic model implies that the width of the young seafloor affected by volcanic extrusion or dominated by tectonic stretching varies through time

morphology and genesis of slow spreading ridges seabed scattering and seismic imaging within the oceanic crust a grid of 32 across axis and five axis parallel multichannel seismic mcs reflection profiles were acquired at an axial volcanic ridge avr segment at 57 45 n 32 35 w on the slow spreading reykjanes ridge mid atlantic ridge to determine the along axis variation and geometry of the axial magmatic system and to investigate the relationship between magma chamber structure the along axis continuity and segmentation of melt supply to the crust the development of faulting and the thickness of oceanic layer 2a seismic reflection profiles acquired at mid ocean ridges are prone to being swamped by high amplitude seabed scattered noise which can either mask or be mistaken for intracrustal reflection events in this paper we present the results of two approaches to this problem which simulate seabed scatter and which can either be used to remove or simply predict events within processed mcs profiles the 37 mcs profiles show clear intracrustal seismic events which are related to the structure of oceanic layer 2 to the axial magmatic system and to the faults which dismember each avr as it ages through its tectono magmatic life cycle and which form the median valley walls the layer 2a event can be mapped around the entirety of the survey area between 0 1 and 0 5 s two way traveltime below the seabed being thickest at avr centres and thinning both off axis and along axis towards avr tips both avr parallel and ridge parallel trends are observed with the pattern of on axis layer 2a thickness variation preserved beneath relict avrs which are rafted off axis largely intact each active avr is underlain by a mid crustal melt lens reflection extending almost along its entire length similar reflection events are observed beneath the offset basins between adjacent avrs these are interpreted as new avrs at the start of their life cycle developing centrally within the median valley the east west spacings of relict avrs and offset basins is 5 7 km corresponding to a life span of the order of 0 5 0 7 myr during which avrs appear to undergo multiple 20 60 kyr tectono magmatic cycles

an ocean bottom gravity study of the southern juan de fuca ridge we use seafloor and sea surface gravity data to model density structure along the southern juan de fuca ridge we determine the average density of the shallow oceanic crust at the ridge using seafloor gravity measurements and then use these data in conjunction with sea surface gravity observations to construct density structure models of 63 seafloor gravity measurements obtained 42 observations were distributed over the 3 km left lateral overlapping rift zone orz located at latitude 45 03 n and separating the juan de fuca ridge into the cleft and vance segments a 21 measurement seafloor gravity transect was made perpendicular to the ridge strike at latitude 44 52 n a region of linear ridge geometry on the cleft segment these seafloor gravity measurements which are sensitive to shallow crustal density variations were augmented by roughly 800 km of sea surface gravity measurements using a bathymetry bouguer anomaly analysis of the seafloor gravity data we determined the average density of the shallow roughly upper 2 km oceanic crust at these two locations to be 2630 kg m 3 50 kg m 3 within the uncertainties of the density determinations there is no difference between the average shallow oceanic crustal density at the linear cleft segment and at the cleft vance orz using the seafloor measurements we modeled the study area s underlying density to provide constraints on its fine scale structure localized porosity of up to 17 within the upper 500 m of oceanic crust layer 2a provide a possible explanation for the observed gravity anomalies at the cleft segment two dimensional density models of the crust underlying the linear cleft segment 44 52 n show that no deep source is required to explain the topographical asymmetry observed between the east and west sides of the ridge axis for the cleft vance orz three dimensional modeling suggests low density material between the rifts in agreement with thickened layer 2a from seismic observations in contrast layer 2a variations do not explain the low density inferred at the northern end of the cleft segment since seismic observations suggest a thin layer 2a the magmatic activity recently observed in this region may have a distinct deeper source

seismic structure and anisotropy of the juan de fuca ridge at 45 n a seismic refraction experiment was conducted with air guns and ocean bottom seismometers on the juan de fuca ridge at 45 n at the northern cleft segment and at the overlapping rift zone between the cleft and vance segments these data determine the average velocity structure of the upper crust and map the thickness variability of the shallow low velocity layer which we interpret as the extrusive volcanic layer the experiment is unique because a large number of travel times were measured along ray paths oriented at all azimuths within a small 20 km by 35 km area these travel times provide evidence for compressional velocity anisotropy in the upper several hundred meters of oceanic crust presumed to be caused by ridge parallel fracturing compressional velocities are 3 35 km s in the ridge strike direction and 2 25 km s across strike travel time residuals are simultaneously inverted for anisotropy as well as lateral thickness variations in the low velocity layer extrusive layer thickness ranges from approximately 200 m to 550 m with an average of 350 m the zone of the thinnest low velocity layer is within the northern cleft segment axial valley in a region of significant hydrothermal activity layer thickness variability is greatest near the cleft vance overlapping rift zone where changes of 300 m occur over as little as several kilometers laterally these low velocity layer thickness changes may correspond to fault block rotations in an episodic spreading system where the low side of each fault block accumulates more extrusive volcanics

east pacific rise from siqueiros to orozco fracture zones along strike continuity of axial neovolcanic zone and structure and evolution of overlapping spreading centers we have used the sea beam multibeam echo sounding system to survey the east pacific rise epr from 8 20 n to 18 30 n obtaining complete coverage of the epr axial neovolcanic zone and all intervening transform faults here we focus on the epr neovolcanic zone and the definition of overlapping spreading centers ocs s between the orozco and siqueiros transform faults the neovolcanic zone is narrow 0 5 2 0 km and continuous along strike and occurs within and near a strike continuous axial summit graben the neovolcanic zone and summit graben reside along the crest of a volcanic axial high which is 2 10 km wide and which continues uninterrupted along strike for 40 140 km within 20 50 km of an intersection with a transform fault the axial neovolcanic zone narrows and deepens and within only 3 6 km of the intersection the neovolcanic zone turns sharply into the transform valley at seven locations between the siqueiros and orozco transform faults the axial neovolcanic zone is discontinuous and is laterally offset a short distance 1 15 km in contrast with a classic ridge transform ridge plate boundary however offset rise terminations overlap each other by a distance approximately 3 times greater than their offset they curve toward each other and often one merges into the other along strike separating the osc s is an elliptical overlap basin up to 600 m deep whose long axis is subparallel to the adjacent spreading centers the overlap basin is characterized by volcanic constructional edifices a key difference between these offsets and small transform faults in the atlantic is that there are no transform fault structures in the overlap basin which link the offset rise axes a continuous axial depth profile reveals a long wavelength 20 60 km undulation of the rise axis which may be associated with variations in the magmatic budget along strike the lowest points in the axial depth profile occur near transform faults and at osc s suggesting that these features are associated with along strike minima in the magmatic budget we suggest that major phases of volcanic activity propagate episodically along the axis away from centers of magmatic inflation and that osc s develop where magmatic pulses fail to meet due to misalignment of fracture systems sinuous bends in the rise axis or other heterogeneities transform faults fail to develop on fast spreading centers where the lateral offsets are small 15 km because the lithosphere is too thin and weak to maintain a classic rigid plate ridge transform fault pattern the osc geometry is unstable and evolves rapidly and the overlapping rises may propagate at rates faster than the local spreading rate one of the two osc s prevails while the other is abandoned while the locations of osc s may not be fixed relative to the rise axis they may tend to recur near the same place creating scars in older lithosphere which resemble fracture zones but whose fine scale structure is quite different there is excellent agreement between the observed structure of osc s and laboratory studies of osc s using wax models photoelastic studies of overlapping cracks and boundary element modeling of overlapping cracks seismic refraction and reflection near two osc s suggests that large offset osc s have separate magma chambers beneath the offset rise axes while small offset osc s may share a single axial magma chamber three dimensional inversion of magnetic anomalies suggests that significant petrologic anomalies may occur near some osc tips waxing and waning of overlapping magma chambers may give rise to the eruption of highly fractionated feti basalts and associated magnetic anomalies

relationship between axial morphology crustal thickness and mantle temperature along the juan de fuca and gorda ridges we analyze gravity data over the juan de fuca and northern gorda ridges to understand the lithospheric structure of two ridges with contrasting axial morphologies spreading at the same intermediate rate 28 mm yr half rate the cleft segment at the south end of the juan de fuca ridge has an axial high morphology while the northern segment of the gorda ridge has a rift valley residual mantle bouguer anomalies rmba on the northern gorda ridge are elevated relative to the cleft segment by 10 20 mgal indicating thinner crust and or a colder mantle the minimum value 50 mgal of the rmba along the juan de fuca ridge is over axial seamount and gradually increases south toward the blanco transform the observed rmba are interpreted to result from along axis variations in crustal thickness and mantle density both of which are controlled by temperatures in the upper mantle where decompression melting occurs we estimate that mantle temperatures are elevated by 30 40 c beneath axial seamount resulting in an excess crustal thickness of 1 7 km the cleft segment is associated with crust that is estimated to be only 300 700 m thicker and mantle temperatures are only 10 15 c higher than beneath the northern gorda ridge however even these small differences in crustal thickness and mantle temperature appear to be sufficient to produce a major change in lithospheric strength and axial morphology these results are consistent with the predicted sensitivity of recent thermo mechanical models for rift valley formation to small changes in crustal thickness and mantle temperature at intermediate spreading rates we attribute the systematic differences in axial morphology crustal thickness mantle temperature and lithospheric strength along the juan de fuca gorda ridge system to the presence of the cobb thermal anomaly at axial seamount the hotter mantle beneath the juan de fuca ridge results in greater amounts of decompression melting and the formation of a thicker crust and a thinner weaker lithosphere than along the gorda ridge

in situ electrical resistivity and bulk porosity of the oceanic crust costa rica rift in situ electrical resistivity was measured in dsdp hole 504b to a depth of 1 013 m into oceanic basement apparent resistivities are about 10 m in about 600 m of pillow lavas sharply increasing to nearly 1 000 m in the underlying dykes bulk porosities calculated from resistivities show a distinct layering apparently corresponding to geophysical layers 2a 2b and 2c

seismic attenuation in upper ocean crust at hole 504b seismic attenuation and its relationship to borehole stratigraphy in the upper 1 8 km of ocean basement at hole 504b are determined from analysis of vertical seismic profile vsp data vsp data provide unambiguous measurements of seismic amplitude decay along a vertical propagation path through the crust and ancillary borehole measurements enable detailed modeling of the relative contributions from geometrical spreading scattering and intrinsic loss mechanisms to this decay about 60 of the total observed amplitude decay occurs in the pillow basalt section and is due mostly to geometrical spreading and scattering from impedance contrasts the remaining amplitude decrease is concentrated in two layers at 500 650 and 800 900 meters below seafloor mbsf 225 375 and 525 625 m below basement across which amplitude rapidly decays and the frequency characteristics of the downgoing wave field are significantly and permanently modified attenuation in these layers is not due to scattering but rather to an intrinsic mechanism that can be characterized by q of 10 and 8 respectively it is likely that the q structure of both of these intervals is formed with the crust near the ridge and thus related to fundamental ocean crust forming processes the shallow interval coincides with a change in alteration mineralogy deposited by late stage fluid flow and may separate lower lavas that were emplaced within the rift zone from upper lavas that were emplaced by off axis flow through large lava tubes intrinsic attenuation in the deeper horizon is probably due to an increase in porosity and cracking associated with either intracrustal deformation or subhorizontal faulting there is negligible attenuation of seismic frequencies in the dikes below 1000 mbsf 725 m subbasement

crustal structure of endeavour ridge segment juan de fuca ridge from a detailed seismic refraction survey a detailed seismic refraction survey was carried out over the endeavour segment of the juan de fuca ridge a medium rate spreading center which lies off western north america to investigate the creation and evolution of oceanic crust a bathymetric high and the presence of hydrothermal vents suggested that the study area was the most recent locus of spreading travel time and amplitude data from 15 in line air gun ocean bottom seismometer profiles were forward modeled using asymptotic ray theory to obtain two dimensional velocity models consisting of four primary layers which correlate well with classic models of oceanic crust significant lateral variations in thicknesses and velocities on the scale of a few to 10 km are superimposed on this basic velocity structure but they appear to be random rather than distributed symmetrically about the ridge we attribute them to fracturing which causes porosity changes hydrothermal circulation which fills voids and fractures with alteration products and variations in magmatic and or deformational processes at the spreading center layer 2a is found to have low 2 6 2 8 km s velocities to average 0 4 km in thickness with variations up to 0 2 km and to be bounded at its base by a sharp velocity increase to 4 8 km s along the axial ridge velocities 0 4 0 6 km s higher than average are interpreted for layers 2b and 2c but these values are confined to a 2 km wide zone centered below the ridge velocities along ridge parallel lines offset 10 km are normal indicating that maturation to off ridge structure has occurred within at most 0 3 ma layer 3 velocities decrease by 0 1 0 2 km s for arrivals traveling along and under the axial ridge perhaps caused by higher temperatures however we find no anomalously low velocities beneath the ridge indicating that no large crustal magma chamber exists on the basis of this study we conclude that magmatic accretion is a fully three dimensional process within ridge segments such as endeavour ridge

the east pacific rise near 21 n 13 n and 20 s inferences for along strike variability of axial processes of the mid ocean ridge abstract data obtained along various segments of the mid ocean ridge mor are used to construct an idealized model for crustal accretion the model seeks to predict the topographic volcanic tectonic and hydrothermal characteristics of any given spreading segment of the mor as a function of distance away from the bounding transform faults this model is based on a series of detailed mapping efforts carried out on segments of the mor having a broad range of spreading rates between 2 and 16 cm yr this paper includes the results of the french searise program carried out during the summer of 1980 aboard the n o jean charcot and two american cruises conducted in 1981 aboard the r v melville

errors in daily ablation measurements in northern greenland 1993 94 and their implications for glacier climate studies ablation climate studies were made at two locations in northern green land in the summers of 1993 and 1994 respec tively daily ablation ii as measured at ten stakes within a small area and the data were compared with each other to detect gross errors for example high deviations for data taken on the same day or low cor rela tions between data series at differ ent stakes indicate erroneous data a fter cl iscard i ng data for one stake in 1993 and t wo stakes in 1994 random errors in daily ablation data for individ ual stakes arc 5 kg m 2 d i which is further reduced to only about 2 kg m 2 d 1 by averag ing over eight or nine stakes random errors in cal culated energy bala nces using the pres ent abla tion data are much imover than found in earlier studies in vest greenland where ablation was only measured on three stakes without any attem pt to detect gross errors aside from day to day errors there arc 10 diffe rences in mca n ablation at different sta kes which are probably caused by small s cal e variations in surface albedo such inters take diffe rences gi e 10 uncertainty in positive degree day faetors which arc 9 8 0 9 and 5 9 0 6 kg m 2 d 1 deg 1 for the two sites

rivera ocean seismic experiment rose overview the rivera ocean seismic experiment rose was designed as a combined sea and land seismic program to utilize both explosive and earthquakes to study a number of features of the structure and evolution of a mid ocean ridge a major oceanic fracture zone and the transition region between ocean and continent the primary region selected for the experiment included the rivera fracture zone the crest and eastern flank of the east pacific north of the rivera and adjacent areas of baja california and mainland mexico the experiment included 1 study of the east pacific rise south of the orozco fracture zone primarily using ocean bottom recording and explosive sources 2 a seismicity program at the orozco and 3 a land based program of recording natural events along the coastal region of mexico a considerable amount of useful data was obtained in each of the three subprograms in the marine parts of the experiment we were able to address a variety of problems including structure and evolution of young oceanic crust and mantle structure and dynamics of the east pacific rise seismicity of the orozco fracture zone and partitioning of energy transmission between the ocean volume and the crust lithosphere on land the fortuitousmore occurrence of the petatlan m7 6 earthquake of march 14 1979 permitted the acquisition of an excellent data set of foreshocks and aftershocks of this large event which provide new insight into the filling of a major seismic gap in the region this overview describes the scientific rationale and the design of the experiments along with some general results less

internal structure of uppermost oceanic crust along the western blanco transform scarp implications for subaxial accretion and deformation at the juan de fuca ridge 1 the uppermost 2 km of oceanic crust created near the southern end of the juan de fuca ridge 60 mm yr full spreading rate is exposed in a tectonic window created along the north wall of the blanco transform fault a total of 27 submersible dives transect this escarpment over a distance of 30 km corresponding to 1 m y of spreading history in this exposure weakly deformed basaltic lavas grade downward into much more intensely fractured lavas cut by basaltic dikes lava flow surfaces dip persistently to the northwest toward the spreading center and the dikes dip southeast away from the spreading center the underlying sheeted dike complex is composed of subparallel dikes that dip moderately to the southeast discrete faults and more distributed cataclastic deformation zones separate the panels of sheeted dikes from one another disruption of the lava and dike units are interpreted as a result of dramatic subaxial vertical mass transport that created the accommodation space for the accumulation of a thickness of 1000 1500 m of basaltic lava the observed patterns of tilting of lavas and dikes imply subsidence beneath the region of active volcanism 3 km wide that is most rapid near the axis and diminishes quickly off axis deformation that accommodates this subsidence has generated extensive fracture porosity in the subaxial region that could strongly influence upper crustal hydrothermal pathways alteration history seismic properties and the distribution of subsurface biological communities

linking rift propagation barriers to excess magmatism at volcanic rifted margins break up related extrusive magmatism imaged in reflection seismic data as seaward dipping reflectors sdrs extends symmetrically along the volcanic margins of the atlantic ocean recent research found distinct along margin variations in the distribution of sdrs and abundance of volcanic material was found to be spatially linked to transfer fault systems these segmented the propagating rift that later developed into the ocean and are interpreted as rift propagation barriers based on these observations we develop a numerical model which shows that rift parallel mantle flow and locally enhanced rates of volcanism are the result of delays in rift propagation and segmented opening our model suggests that segmentation is one of the major factors in the distribution and localization of riftrelated extrusive magmatism we conclude that in addition to mantle temperature and inherited crustal structures e g weaknesses from previous rift episodes rift propagation delay plays an important role in the distribution of extrusive volcanism at volcanic passive margins by controlling the mantle flow beneath the rift axis

a finite element computation of seismic diffraction about the rise axis magma chamber in a previous study models of the crust and magma chamber at 12 50 n on the east pacific rise were derived with a two dimensional ray tracing analysis of rise normal refraction lines the best model featured a 4 km wide axial magma chamber but the data contained diffraction effects that were not modeled by geometric ray theory this raised questions as to whether this model would reproduce the waveform data and whether diffraction might be obscuring a much larger magma chamber we have tested the model with a finite element calculation of synthetic seismograms this method incorporates full wave effects including diffraction the synthetics accurately reproduce the data and exhibit diffraction phenomena similar to those in the data the previous travel time analysis was largely successful and was only slightly compromised by diffraction diffraction in rise normal refraction experiments is therefore not obscuring a large low velocity magma chamber

new observations of the shallow seismic structure of young oceanic crust the results are presented of three new experiments carried out on the mid atlantic ridge mar near latitude 23 n using a unique method of studying the seismic structure of the uppermost few hundred meters of the oceanic crust the data were collected using a fixed ocean floor hydrophone receiver and a controllable explosive source that was towed within a few tens of meters of the rugged bottom topography these 1 to 2 km long refraction lines produced for the first time direct observations of the compressional wave velocity structure of the uppermost 200 300 m of the young igneous crust one experiment was carried out over the site of hole 648b of the ocean drilling program on a small volcano within the median valley of the mar the seafloor velocity was observed to be 2 1 km s 1 underlain by an approximately linear velocity gradient of 4 s 1 given that we know the crust at this location consists of fresh basalt lavas with laboratory measured velocities in excess of 5 8 km s 1 porosities at the seafloor of as high as 30 50 are inferred the two other experiments were located over 7 m y old crust near deep sea drilling project site 395 west of the mar separated by only 14 km laterally but by over 1400 m in water depth the experiment positioned in the 100 m thick sediment pond and that located 14 km to the south atop a prominent topographic high produced results that within the data resolution were indistinguishable the velocity of the uppermost basaltic basement was 4 1 km s 1 and the velocity gradient was less than 0 5 s 1 this 2 km s 1 difference in velocity between the two sites is assumed to be a consequence of age related modifications to the physical properties of the crust if the primary change is in only the total porosity then a 15 20 reduction is required to explain these observations this seems unlikely because any such pervasive process would have been previously observed in geophysical and sampling data e g off axis volcanism secondary mineral deposition tectonic compression a combination of many processes acting in unison would be necessary to produce such a decrease in porosity if the geometry of the cracks and voids were modified in just the right way then no change in total porosity is required to explain the velocity difference these on bottom refraction profiles are proven to be practical experiments that can provide precise information on a scale appropriate to many of the ocean floor s geologic features they produce results that may be correlated in a meaningful way with results of ocean drilling and downhole measurements

microearthquakes on the endeavour segment of the juan de fuca ridge 1 we report the results of a 55 day microearthquake experiment on the endeavour segment of juan de fuca ridge the network covered a 5 km section of the ridge axis centered on the main vent field and extended 15 km off axis on the west flank the ridge axis and flanks were seismically active and 1750 earthquakes were located with a minimum of five travel time picks including at least one s wave over half the earthquakes occurred in swarms and a waveform cross correlation technique was used to obtain relative locations on the western flank the hypocenters for four swarms at midcrustal depths are compatible with steeply dipping fault planes that strike at 035 050 n and oblique to the abyssal hills focal mechanisms determined from p wave first motions and p s amplitude ratios are predominantly strike slip with north south compression and appear to be affected by the reorganization of the explorer plate earthquakes beneath the ridge axis are concentrated in a band of intense seismicity at 1 5 3 5 km depth to the north of high rise vent field the seismicity defines a plane striking parallel to the ridge axis and dipping east at 70 and the earthquakes appear to extend beneath an axial reflector previously imaged at 2 3 km depth farther south the hypocenters are not compatible with a single fault plane focal mechanisms are characterized by subhorizontal tension axes oriented in all directions except parallel to the ridge and suggest a stress field that is about equally influenced by ridge spreading and hydrothermal cooling

dynamic behaviour analysis of glacier de saint sorlin france from 40 years of observations 1957 97 mass balance and dynamic measurements carried out on glacier de saint sorlin since 1957 provide a good opportunity to study the dynamics of this glacier ice flow analysis shows that dynamic changes have been important over the last 40 years and that these changes are not consistent with the concepts usually used in glacier modelling present velocities are larger than the 1960 velocities although the thickness decreased everywhere 10 30 m in the ablation zone a simple numerical ice flow model which does not include longitudinal stress gradients has been used to investigate these phenomena this model allows us to infer the sliding velocity from observed surface and calculated deformation velocities we conclude that 1 the sliding velocity cannot be described by weertman analysis or empirical relations which link the sliding to the thickness and surface slope 2 the inferred sliding velocity is uniform over at least half of the glacier and 3 there is no clear link between the sliding process and the quantity of water coming from surface ablation furthermore it may not be reasonable to calibrate model flow parameters from geometry changes because the surface geometry is relatively insensitive to velocity changes over some decades

tomographic inversion for three dimensional velocity structure at mount st helens using earthquake data tomographic inversion is applied to 17 659 p phase observations at 21 stations from 2023 earthquakes in the vicinity of mount st helens to study the three dimensional velocity structure block size for the inversion is 2 km horizontally and 2 km or more vertically locations of hypocenters are assumed known and are based on a reference one dimensional layered velocity structure a conjugate gradient technique lsqr is used to invert the large sparse system of equations augmented by regularization with a laplacian roughening matrix resolution is estimated by computing the impulse response of the inversion for various critical locations and uncertainties of the estimates are determined by a jackknife approach the results of the inversion show a remarkable correlation with known geological and geophysical features the spirit lake and spud mt plutons are characterized by high velocity regions extending to approximately 9 km depth the st helens seismic zone a band of diffuse seismicity extending nnw from the volcano is evident as a prominent low velocity lineation the change in character of the velocity anomalies south of st helens corresponds well with the near cessation of seismic activity there a low velocity anomaly beneath the crater from 6 to 16 km depths may represent modern magma accumulations

recent movements of the juan de fuca plate system analysis of the magnetic anomalies of the juan de fuca plate system allows instantaneous poles of rotation relative to the pacific plate to be calculated from 7 ma to the present by combining these with global solutions for pacific america and absolute relative to hot spot motions a plate motion sequence can be constructed this sequence shows that both absolute motions and motions relative to america are characterized by slower velocities where younger and more buoyant material enters the convergence zone pivoting subduction the resistance provided by the youngest portion of the juan de fuca plate apparently resulted in its detachment at 4 ma as the independent explorer plate in relation to the hot spot framework this plate almost immediately began to rotate clockwise around a pole close to itself such that its translational movement into the mantle virtually ceased after 4 ma the remainder of the juan de fuca plate adjusted its motion in response to the fact that the youngest material entering the subduction zone was now to the south differences in seismicity and recent uplift between northern and southern vancouver island may reflect a distinction in tectonic style between the normal subduction of the juan de fuca plate to the south and a complex underplating occurring as the explorer plate is overridden by the continent the history of the explorer plate may exemplify the conditions under which the self driving forces of small subducting plates are overcome by the influence of larger adjacent plates the recent rapid migration of the absolute pole of rotation of the juan de fuca plate toward the plate suggests that it too may be nearing this condition

structure of the earth oceanic crust and uppermost mantle the past four years have witnessed the introduction of a variety of new instruments and methods for the study of the seismic structure of the oceanic crust and lithosphere the application of these and existing tools has led to the discovery of a number of new phenomena and to a fuller understanding of the genesis and evolution of the oceanic ithosphere borehole seismic instrumentation has become more widely employed ocean bottom seismographs while generally decreasing in number have become significantly more reliable and useful and multichannel seismic systems have been employed in innovative experiments ranging from studies of fracture zones to the regular detection of magma chambers beneath rise axes the techniques available for the analysis of seismic data have become more sophisticated waveforms collected in seismic experiments can now be used directly in constructing and evaluating seismic velocity models and travel time data are regularly inverted directly for structure trial and error modeling has become increasingly unimportant marine seismologists are becoming increasingly involved in understanding the coupling between the ocean and the underlying oceanic lithosphere this has led to a more complete understanding of seafloor noise processes and the partitioning of energy between acoustic and elastic waves

seismic amplitude modeling and the shallow crustal structure of the east pacific rise at 12 n velocity structures of the upper 2 km of the oceanic crust were used to model seismic refraction data collected on the eastern flank of the east pacific rise at 12 n during the rivera ocean seismic experiment rose data generated with a bolt 1500c air gun source show multiple narrow 1 km in range peaks in the amplitude of the first refracted arrival record sections separated by as much as 150 km indicate similar amplitude versus range patterns suggesting that the amplitude peaks are caused by consistent structural features and not by random crustal inhomogeneities these amplitude highs do not consistently correlate with seafloor bathymetry and hence cannot be explained solely by topographic focussing of seismic energy through travel time analysis and ray tracing we find a crustal velocity model with two vertical zones of concave gradients which focus rays at 6 and 8 km range synthetic seismograms and power versus range patterns computed for this model by the reflectivity algorithm are consistent with the two peak amplitude pattern observed on nearly all air gun record sections these results suggest that the vertical gradient zone structure of the uppermost 2 3 km of ocean crust is continuous to first order over the rose area

wavelet analysis of deep tow magnetic profiles modeling the magnetic layer thickness over oceanic ridges 1 the interpretation of marine magnetic anomalies usually consists either in determining the magnetization distribution assuming the source geometry and magnetization direction or in determining the magnetic layer thickness assuming the magnetization direction and intensity in this paper we introduce a new technique that allows modeling of the thickness of the magnetic source layer with very few a priori assumptions about the magnetization the magnetic layer is assumed to be made of a series of bodies each having a constant unknown magnetization and an unknown size this technique is based upon the application of the continuous wavelet transform recently introduced for the interpretation of potential field data as a multipole decomposition we present applications to synthetic data to one deep tow magnetic profile recorded across the juan de fuca ridge jdf and to three deep tow magnetic profiles recorded across the central indian ridge cir our results confirm that despite significant source thickness variations 100 1200 m across the cir measured magnetic anomalies mostly reflect past geomagnetic field intensity fluctuations however we show that within the axial region of high magnetization thickness variations have a significant contribution to short wavelength variations of deep tow magnetic signals 100 nt

seismic constraints on shallow crustal emplacement processes at the fast spreading east pacific rise we present the results of nine on bottom seismic refraction experiments carried out over young east pacific rise crust the experiments are unusual in that both the source and receiver are located within a few meters of the seafloor allowing high resolution determinations of shallow crustal structure three experiments were located within the axial summit caldera asc over zero age crust the seismic structure at these three locations is fundamentally the same with a thin 60 m surficial low velocity 2 5 km s layer a 100 to 150 m thick transition zone with velocities increasing by 2 5 km s and a layer with velocities of 5 km s at a depth beneath the seafloor of 130 190 m the surficial low velocity layer and transition zone are defined as seismic layer 2a and the 5 km s layer is defined as the top of layer 2b both the surficial low velocity layer and the transition zone double in thickness within 1 km of the rise axis we model layer 2a as the extrusive sequence and transition zone and the 2a 2b boundary as the top of the sheeted dikes the primary implication of this interpretation is that the depth to the top of the sheeted dikes deepens from 150 m to 300 m within 1 km of the asc the thickening of the extrusive layer is interpreted to be due to lava that either overflows the asc walls is emplaced through eruptions outside of the asc or travels laterally from the asc through conduits the most probable cause for the thickening of the transition zone is sill emplacement outside of the asc either from magma that does not reach the surface in an off axis eruption or magma that is transported laterally during the drainage process creating the asc we suggest that the mechanism controlling the magnitude and rate of the dike subsidence is the mechanism that determines the thickness of the extrusive section and the total thickness of layer 2a

ne atlantic continental rifting and volcanic margin formation the ne atlantic rift system which has been under extensional deformation since the caledonian orogeny experienced a series of rift episodes in late paleozoic to early cenozoic times fig 1 late paleozoic rifting is poorly constrained particularly with respect to timing however rifted basin geometries inferred to be of this age are observed at depth in seismic data from the flanks of the younger rift structures

seismic structure of the endeavour segment juan de fuca ridge correlations with seismicity and hydrothermal activity 1 multichannel seismic reflection data collected in july 2002 at the endeavour segment juan de fuca ridge show a midcrustal reflector underlying all of the known high temperature hydrothermal vent fields in this area on the basis of the character and geometry of this reflection its similarity to events at other spreading centers and its polarity we identify this as a reflection from one or more crustal magma bodies rather than from a hydrothermal cracking front interface the endeavour magma chamber reflector is found under the central topographically shallow section of the segment at two way traveltime twtt values of 0 9 1 4 s 2 1 3 3 km below the seafloor it extends approximately 24 km along axis and is shallowest beneath the center of the segment and deepens toward the segment ends on cross axis lines the axial magma chamber amc reflector is only 0 4 1 2 km wide and appears to dip 8 36 to the east while a magma chamber underlies all known endeavour high temperature hydrothermal vent fields amc depth is not a dominant factor in determining vent fluid properties the stacked and migrated seismic lines also show a strong layer 2a event at twtt values of 0 30 0 09 s 380 120 m below the seafloor on the along axis line and 0 38 0 09 s 500 110 m on the cross axis lines a weak moho reflection is observed in a few locations at twtt values of 1 9 2 4 s below the seafloor by projecting hypocenters of well located microseismicity in this region onto the seismic sections we find that most axial earthquakes are concentrated just above the magma chamber and distributed diffusely within this zone indicating thermal related cracking the presence of a partially molten crustal magma chamber argues against prior hypotheses that hydrothermal heat extraction at this intermediate spreading ridge is primarily driven by propagation of a cracking front down into a frozen magma chamber and indicates that magmatic heat plays a significant role in the hydrothermal system morphological and hydrothermal differences between the intermediate spreading endeavour and fast spreading ridges are attributable to the greater depth of the endeavour amc and the corresponding possibility of axial faulting

the east pacific rise in cross section a seismic model in 1982 we undertook a seismic refraction experiment known as the magma expedition to examine the detailed structure of the east pacific rise near 12 50 n this segment of the rise where the full spreading rate is about 110 mm yr is near the projected trace of the o gorman fracture zone and is the site of an overlapping spreading center as well as black smoker hydrothermal activity in this paper we describe the analysis of a subset of the travel time data collected during the magma expedition namely the data from profiles which were oriented normal to the rise axis these profiles provide a data set roughly equivalent to those collected on other experiments and sample the cross sectional structure of the rise we have modeled these data using a two dimensional ray tracing program we have found that the seismic velocities in the young oceanic crust are rather high with velocity gradients of 4 0 5 5 s 1 in the uppermost crust the highest velocities at the seafloor occur beneath the rise axis itself and seem to decrease as the crust ages to 0 1 ma this decrease in velocity must result from an increase in porosity in the upper crust and may coincide with the development of abundant surface fissures as the crust spreads the decrease in velocity does not appear to penetrate deeper than about 0 5 km and may reverse itself as hydrothermal alteration fills the pores and cracks the fact that the highest velocities occur under the rise axis suggests that the hydrothermal circulation responsible for the black smokers is confined to seismically unresolved conduits a result consistent with the high temperatures and discrete nature of the vents our best model includes a magma chamber some 4 km wide and extending from the moho to within about 1 1 km of the seafloor this magma chamber is far smaller than many models for the rise axis have predicted but larger than those inferred from seismic refraction experiments at other sites on the east pacific rise these discrepancies probably arise because the magma chamber under the east pacific rise is not a steady state feature but changes with time because of hydrothermal cooling and perhaps because of an episodic supply of magma from the mantle or along the rise axis

mass balance measurements problems and two new methods of determining variations les deux nouvelles methodes d analyse des bilans sont basees sur le modele lineaire de variation des bilans propose par lliboutry 1974 la premiere n utilise qu un systeme de mesures limitees a quelques sites fixes la seconde basee sur l equation de continuite permet d obtenir un bilan sur un secteur limite par deux profils transversaux ou ont ete mesures les altitudes et vitesses moyennes en surface

processes associated with ocean crustal formation the juan de fuca ridge recent geophysical geological and geochemical studies on the juan de fuca ridge have been gathered into a special collection of manuscripts that demonstrate the interdisciplinary nature of the crustal formation and evolutionary processes on the ridge axis these new studies include long term sea floor and water column measurements that indicate real time changes due to volcanic activity and related hydrothermal flux sea floor measurements of seismic activity detecting swarms of earthquakes near the axial zone and submersible observations of a 2 km along strike spacing of large scale hydrothermal fields detailed sampling of axial basalts show that in spite of abundant evidence of excess volcanic activity the juan de fuca ridge is not a traditional hot spot on a larger scale geophysical surveys of the ridge flanks show that the upper crust can vary in thickness by a factor of two that crustal porosity systematically decreases with age and the topographic forcing is the dominant control of off axis hydrothermal circulation finally the juan de fuca ridge has long been a testing ground for innovative survey techniques and this issue includes reports on the development of acoustic backscatter as a quantitative measure of basement roughness and as an indicator of the presence of lithified sediments overlying zones of fluid upwelling

crystalline lithology across the kapuskasing uplift determined using in situ poisson s ratio from seismic tomography seismic refraction tomography has been used to image lateral velocity variations both vp and vs within the shallow crystalline basement using travel times from a high resolution reflection profile across the kapuskasing uplift structure in northern ontario comparison of these in situ velocities and the outcrop geology indicates a correlation of velocity with the major lithologic units anorthosite vp 6150 6350 m s and vs 2700 3000 m s mafic gneiss vp 6000 6250 m s and vs 2775 3000 m s tonalite vp 5850 6000 m s and vs 2950 3025 m s vp measurements on core samples from the main lithologic units confirm that tonalite can be be clearly distinguished from mafic gneiss and anorthosite based on vp the observed vp values are attributed to mineralogy based on modal analyses of core samples from the same region poisson s ratio calculated from in situ vp and vs helps to discriminate between velocity variations caused by porosity effects from those due to lithology values of poisson s ratio obtained for the three lithologic units within the upper 100 m of the basement are anorthosite 0 35 0 39 mafic gneiss 0 34 0 38 and tonalite 0 32 0 34 these relatively high values of poisson s ratio indicate that the shallow basement is porous and water saturated poisson s ratio is a particularly effective discriminator in this instance because of its extremely low value in quartz 0 08 relative to other silicate minerals 0 21 0 32 and the greater quartz content of tonalite 27 36 relative to anorthosite and mafic gneiss 0 23 reflection coefficients calculated for the contrasting lithologies range up to 0 065 for tonalite mafic gneiss contacts and up to 0 044 for tonalite anorthosite contacts confirming that an interlayered sequence of tonalitic anorthositic and mafic rocks mapped at the surface may be the cause of the shallow reflections observed in the coincident seismic reflection data

aqueous volatiles in hydrothermal fluids from the main endeavour field northern juan de fuca ridge temporal variability following earthquake activity the main endeavour field northern juan de fuca ridge experienced intense seismic activity in june 1999 hydrothermal vent fluids were collected from sulfide structures in september 1999 and july 2000 and analyzed for the abundance of h2 h2s ch4 co2 nh3 mg and cl to document temporal and spatial changes following the earthquakes dissolved concentrations of co2 h2 and h2s increased dramatically in the september 1999 samples relative to pre earthquake abundances and subsequently decreased during the following year in contrast dissolved nh3 and ch4 concentrations in 1999 and 2000 were similar to or less than pre earthquake values aqueous cl abundances showed large decreases immediately following the earthquakes followed by increases to near pre earthquake values the abundances of volatile species at the main endeavour field were characterized by strong inverse correlations with chlorinity phase separation can account for 20 50 enrichments of co2 ch4 and nh3 in low chlorinity fluids while temperature and pressure dependent fluid mineral equilibria at near critical conditions are responsible for order of magnitude greater enrichments in dissolved h2s and h2 the systematic variation of dissolved gas concentrations with chlorinity likely reflects mixing of a low chlorinity volatile enriched vapor generated by supercritical phase separation with a cooler gas poor hydrothermal fluid of seawater chlorinity decreased abundances of sediment derived nh3 and ch4 in 1999 indicate an earthquake induced change in subsurface hydrology elevated co2 abundances in vent fluids collected in september 1999 provide evidence that supports a magmatic origin for the earthquakes temperature salinity relationships are consistent with intrusion of a shallow dike and suggest that the earthquakes were associated with movement of magma beneath the ridge crest these data demonstrate the large and rapid response of chemical fluxes at mid ocean ridges to magmatic activity and associated changes in subsurface temperature and pressure

structure and variability of oceanic crust on the flanks of the east pacific rise between 11 and 13 n we use over 500 seismograms collected by five ocean bottom seismometers to examine the structure and variability of 0 5 m y old crust in the rivera ocean seismic experiment rose area of the east pacific rise between 11 and 13 n there are no significant differences among the first arrival travel times for 10 independent refraction line segments along the isochron amplitude patterns recorded on the individual lines however suggest that structural boundaries within the upper crust vary by hundreds of meters in depth over lateral distances of tens of kilometers within this 200 km segment of lithosphere on the basis of the modeling of amplitude patterns at ranges of 5 10 km and of comparisons with the upper crust sampled at deep sea drilling project hole 504b and in ophiolite complexes layer 2 can be divided into three units 2a corresponds to extrusive volcanics with a steep velocity gradient 2b to extrusives with a near zero velocity gradient and 2c to a complex region of transition from extrusives to sheeted dikes the total thickness of layer 2 varies by about 0 7 km over the length of the line another amplitude peak which varies in range between 14 and 21 km may be related to the downward transition from sheeted dikes to isotropic gabbros within layer 3 the variability of these two amplitude patterns provides direct evidence that the along strike structure of oceanic crust in the rose area is definably heterogeneous and that the process of crustal accretion is likewise variable on the scale resolved by amplitude and travel time information

dsdp hole 504b the first reference section over 1 km through layer 2 of the oceanic crust hole 504b provides a reference section for correlation of in situ petrological and geophysical studies with regional crustal models based on dredged samples and remote surveys and at least for the medium spreading rate costa rica rift hole 504b confirms that the main features of the upper part of an idealized ophiolite sequence occurs in the oceanic crust

the laptev sea rift abstract in north eastern siberia the active mid ocean gakkel ridge interacts with the continental shelf of the laptev sea extension has affected the shelf since at least the early tertiary and has resulted in the formation of a complex horst and graben system we present new seismic data from the laptev sea including deep seismic soundings the most prominent rift basin is the ust lena rift with a minimum e w width of 300 km at latitude 75 n and a cenozoic infill up to 13 km in thickness the asymmetric shape of the basin and conclusive evidence for a detachment imply a simple shear geometry the suggested rift model combines a ramp and flat geometry for the detachment with ductile stretching beneath the detachment a major west dipping hingeline listric fault separates the ust lena rift from the laptev horst the 100 150 km wide laptev horst is subdivided into three units by narrow rift grabens another prominent rift graben is the anisin basin which is located in the northern shelf area though the laptev sea rift formed in interaction with an active mid oceanic ridge there are indications that the laptev sea rift is of the passive rift type the rift was developed east of a sw ne trending transfer zone which links the gakkel ridge to the laptev sea rift

mid cenozoic post breakup deformation in the passive margins bordering the norwegian greenland sea abstract since the initial breakup of the norwegian greenland sea in early eocene time chron 24b 53 7 ma the conjugate passive margins have been subject to compression extension magmatism and broad uplifts and subsidence mid to late cenozoic compressional doming basin inversion and reverse faulting occurred predominantly offshore within the rheologically weak deep cretaceous depocenters on the norwegian margin the domes seem to have formed primarily in middle eocene to early oligocene and in early miocene while doming in the faroes rockall region developed in late paleocene to early eocene oligocene and mid to late miocene we suggest that the main driving force was plume enhanced ridge push a testable working hypothesis is that discrete phases of compression were associated with periods of high plume flux in addition the oligocene chron 13 35 ma change in plate motion appears to have generated local shear adjustments along pre existing basement weaknesses forming compressional structures in the overlying sedimentary section this change in plate motion was first manifest by significant extension between the sw barents sea and the conjugate nne greenland margin and led to spreading along the knipovitch ridge mid cenozoic extension and magmatism also preceded the development of the kolbeinsey ridge the associated separation of the jan mayen microcontinent from east greenland must have been strongly influenced by the position of the iceland plume center in addition to the change in extension direction with respect to initial breakup of the norwegian greenland sea in early eocene time the mid cenozoic extension and magmatism is a post breakup phenomena widely spaced areas of neogene uplift occurred along all bordering margins these enigmatic uplifted areas acted as nucleation sites for growth of ice caps and initiated large glacial deposits that built out on the adjacent shelves

symmetric conjugation of continent ocean boundary structures along the norwegian and east greenland margins abstract a characteristic structural feature of the norwegian continental margin that has been described from multichannel seismic data as well as dsdp and odp drilling is a major buildup of layered extrusive basaltic rocks forming wedges that dip seaward along the entire margin models for their origin describe them as either products of a late stage intra continental rift phenomenon with the volcanic rocks overlying attenuated continental crust or as an integral component of the oceanic crust formed during the first few million years of seafloor spreading in the norwegian greenland sea from studies of the norwegian margin alone it was not possible to unequivocally distinguish between the two proposed models so an investigation of the conjugate portions of the east greenland margin was undertaken which used two ship multichannel seismic acquisition as its focus we found that the wedges were not only present on the east greenland margin but that details of their along strike variability are mirrored on both margins together with their spatial relationship to the oldest magnetic lineation pattern the wedges are relatively small in the northern part of the margin where they occur almost entirely landward of lineation 24 moving south we see that the wedges broaden encroach into the pattern of lineation 24 then develop into two separate wedges the seaward of which is associated with lineation 23 lineations 24 and 23 on both margins are distorted when associated with wedges we believe that these data lend strength to a tectonic association of the wedges with nascent plate boundary processes their observation in many margins throughout the world suggest that a prolific volcanic phase may be commonly associated with the inception of a divergent plate boundary and we believe that the most likely mechanism to generate the excess volcanism is the convective partial melting scheme recently proposed by mutter and buck 1986 and mutter et al in press

using remote sensing data to determine equilibrium line altitude and mass balance time series validation on three french glaciers 1994 2002 alpine glaciers are very sensitive to climate fluctuations and their mass balance can be used as an indicator of regional scale climate change here we present a method to calculate glacier mass balance using remote sensing data snowline measurements from remotely sensed images recorded at the end of the hydrological year provide an effective proxy of the equilibrium line mass balance can be deduced from the equilibrium line altitude ela variations three well documented glaciers in the french alps where the mass balance is measured at ground level with a stake network were selected to assess the accuracy of the method over the 1994 2002 period eight mass balance cycles results obtained by ground measurements and remote sensing are compared and show excellent correlation r 2 0 89 both for the ela and for the mass balance indicating that the remote sensing method can be applied to glaciers where no ground data exist on the scale of a mountain range or a given climatic area the main differences can be attributed to discrepancies between the dates of image acquisition and field measurements cloud cover and recent snowfalls constitute the main restrictions of the image based method

crustal and uppermost mantle velocity structure of northern eurasia along the profile quartz the velocity structure of the crust and uppermost mantle beneath the dss profile quartz in northwestern russia is determined using a linearized seismic tomography scheme the analysis is performed on the northern 2200 km of the profile spanning from the kola peninsula the timan pechora basin and the uralian fold belt to the west siberian basin refracted and reflected phases recorded from 28 chemical and two nuclear shot points are analyzed the inversion technique allows determination of the p wave velocity structure and estimation of the depth of the crust mantle discontinuity resolution analysis shows that the upper crustal structure and sections of the lower crust and uppermost mantle are well resolved the modeled velocity and crustal thickness anomalies are confirmed by the resulting structure of the resolution matrix poorer ray coverage in parts of the lower crust and most sections of the uppermost mantle limits their resolvable structures a thick over 10 km high velocity block vp 7 0 km sec in the lower crust of the precambrian beneath the kola peninsula crustal thinning and a wide range of compressional velocities in the sedimentary basin structures and a significant crustal root beneath the uralian fold belt 42 to 45 km moho depth with high lower crustal velocities exceeding 7 0 km sec were identified in the tectonic structure of northwestern eurasia

photogrammetric reconstruction of glacier mass balance using a kinematic ice flow model a 20 year time series on grubengletscher swiss alps the kinematic boundary condition at the glacier surface can be used to provide glacier mass balance at individual points if changes in surface elevation horizontal and vertical surface velocities and surface slope are known vertical ice velocity can in turn be estimated from basal slope basal ice velocity and surface strain this relation is applied to reconstruct a 20 year mass balance curve of grubengletscher swiss alps largely using repeated aerial photogrammetry but only a minimum of field work for individual years the mass balance distribution on the glacier tongue was modelled with an accuracy of about 0 9 m a 1 ice mechanical assumptions and errors in glacier bed geometry markedly affect discrete mass balance patterns but are eliminated to a major extent when calculating year to year mass balance changes the resulting curve for the grubengletscher tongue 1973 1992 shows reasonable consistency with meteorological data and other glaciologically derived mass balance series large changes in measured ice speed on the glacier tongue 50 significantly governed the long term variability of ice thickness over the observational period

fine scale magnetic anomaly field over the southern juan de fuca ridge axial magnetization low and implications for crustal structure sea surface deep tow and submersible magnetic surveys show that the spreading axis of the cleft segment of the southern juan de fuca ridge sjdf is overlain by a distinctive magnetic anomaly low in contrast to the more typical magnetic anomaly high found over the adjacent vance segment and most mid ocean ridge spreading centers forward iterative thickness modeling based on the near bottom magnetic data shows that this axial magnetic low can be produced by a thin magnetic source layer beneath the cleft axial valley thick crust is required to produce the magnetic anomaly high over the vance axial valley this magnetic modeling demonstrates that variations in the base of the magnetic source layer which is assumed to be equivalent to the base of the extrusive lavas can produce a significant proportion of the observed short wavelength magnetic anomaly signal these magnetic layer thickness models also correlate quite well with nearby seismic results which show that seismic layer 2a thins beneath the axial valley of the cleft and is thicker beneath the vance axial valley recent seismic results over other fast and medium spreading rifts also show rapid changes in the thickness of the uppermost crust near the axis of spreading which could produce significant magnetic contrasts the results of this study and the correlation with independent seismic measurements suggest that while magnetization may vary due to geomagnetic behavior or because of petrology or alteration first order change in upper crustal thickness is one of the most important variables to consider in the production of fine scale magnetic anomalies in young ocean crust

jasper seamount structure seafloor seismic refraction tomography the velocity structure of jasper seamount was modeled using one and three dimensional inversions of p wave travel times the results represent the first detailed seismic images of a submerged intraplate volcano two seismic refraction experiments were completed on jasper seamount incorporating ocean bottom seismometers and navigated seafloor shots the p wave travel times were first used to compute a one dimensional velocity profile which served as a starting model for a three dimensional tomographic inversion the seamount p velocities are significantly slower than those observed in typical oceanic crust at equivalent subbasement depths this suggests that jasper seamount is constructed predominantly of extrusive lavas with high average porosity the velocity models confirm morphological predictions jasper seamount is a shield volcano with rift zone development high seismic velocities were detected beneath the large radial ridges while low velocities characterize the shallow summit and flanks comparisons between p velocity models of jasper seamount and the island of hawaii reveal that these two shield volcanoes are not structurally proportional jasper seamount is far smaller than hawaii yet both volcanoes exhibit an outer extrusive layer of similar thickness this suggests that seamount size influences the intrusive extrusive proportions density equilibrium between melt and country rock may explain this behavior

hydrothermal alteration of a 1 km section through the upper oceanic crust deep sea drilling project hole 504b mineralogy chemistry and evolution of seawater basalt interactions deep sea drilling project hole 504b penetrates 1076 m into oceanic layer 2 and is the first hole to pass through the transition from pillow basalts altered at low temperatures into hydrothermally metamorphosed sheeted dikes alteration of the crust at site 504 occurred in four stages related to the movement of the crust away from the spreading axis 1 dikes reacted with seawater 200 300 c in the upwelling zone of an axial convection cell at the spreading axis resulting in the formation of greenschist facies parageneses in veins and host rocks mixing of the upwelling hydrothermal fluids with seawater circulating in the overlying more permeable pillow section occurred in the upper part of the lithologic transition zone causing a steep temperature gradient at the base of the pillow section secondary minerals formed in the lower pillow section from the resultant reducing mixed fluids at temperatures of around 100 c at the same time the initial effects of seafloor weathering began in the upper 320 m of the pillow section at low temperatures 50 c and under conditions of open circulation of oxidizing seawater 2 following refracturing of the dikes a second stage of axial upwelling occurred hydrothermal fluids 200 380 c were probably similar to those presently sampled from spreading ridges on the seafloor mixing of these fluids with seawater in the lithologic transition zone caused deposition of quartz epidote and sulfides in veins and formation of a sulfide rich stockworklike zone within this transition zone alteration of the pillow section proceeded under the prior conditions with the effects of seafloor weathering extending progressively downward into the crust 3 seawater recharge penetrated to depths of at least 1075 5 m subbasement and deposited anhydrite locally in veins 4 off axis alteration of the dikes was characterized by formation of zeolites calcite and prehnite in veins and host rocks from more highly evolved and lower temperature 100 250 c fluids alteration in the upper pillow section evolved from seafloor weathering conditions to more reducing and rock dominated as cracks in the basalt were sealed with secondary minerals and the basement was covered with a layer of sediment zeolites and calcite were the last phases to form throughout the pillow section sealing of the crust to convective cooling also allowed conductive reheating of the crust from below calculations of seawater crustal chemical fluxes from whole rock data are complicated by the low recovery and heterogeneity of alteration effects but indicate that basalt seawater interactions are a sink for seawater mg and k and a source for si and ca

active hydrothermal vents and sulfide deposits on the southern juan de fuca ridge massive sulfide deposits rich in zinc and silver were recovered from the juan de fuca ridge 500 km west of oregon in september 1981 the samples recovered are composed largely of zinc sulfide with lesser amounts of iron lead and copper sulfide most of the deposits occur at a series of hydrothermal vents within a relatively continuous depression in the center of a smooth 1 km wide valley along the ridge axis the depression appears to be formed by collapse of a lava lake possibly modified by extensional faulting the axial valley floor outside the depression is underlain by fresh glassy ferrobasalt sheet and lobate flows two types of sulfide mineral deposits were dredged from one of the hydrothermal vents 1 angular slabs of dark gray zinc sulfide and 2 subrounded fragments of porous light gray zinc sulfide the samples contain fresh sphalerite zoned wurtzite pyrite and minor amounts of marcasite galena and chalcopyrite cubanite the spreading rate of the juan de fuca ridge and the composition of the sulfide samples are generally similar to the east pacific rise lat 21 n site however the texture and geologic setting of the sulfide deposits are significantly different

lava flow layer east pacific rise a multichannel seismic reflection profile across the east pacific rise near 9 n reveals an approximately 400 m thick layer at the top of the oceanic crust with seismic velocity of about 2 8 km s at the rise crest this velocity increases to about 3 8 km s 210 km east of the rise crest these low velocities suggest that this layer is probably a porous lava flow layer layer 2a of seismic refraction studies generated by volcanic activity at the spreading center reflections from the base of the layer imply a relatively sharp transition between the lava flows and the underlying layer presumably a sheeted dike complex the underlying layer is a zone about 1 4 km thick with velocity about 5 7 km s whose base may coincide with the top of a low velocity magma chamber from which the crust develops

asymmetric deep crustal structure across the juan de fuca ridge a multichannel line was shot across the endeavour segment of the juan de fuca ridge careful processing of 25 km of the line centered about the rise axis shows that young crustal structure is asymmetric an intermittent reflector from the top 1 km of crust follows topography on both the pacific and the juan de fuca plates it could be either the pillow basalt diabase dike contact or a metamorphic front within the pillow basalts a reflection exists 2 3 km below the rift valley floor diffractions from its edges indicate that the rift valley9s velocity structure is markedly different from the structure of the flanks the deep crustal structure is highly asymmetric moho exists within 4 km of the rift valley on the pacific plate in contrast to the juan de fuca plate which shows a reflector at intermediate depths dipping away from the rise to reach moho depth 12 km from the rift valley this reflector9s origin is enigmatic we speculate that the thermal history and structure of the two plates differ because of local differences in mantle structure caused by either ridge migration or a melting anomaly

pore pressure and oceanic crustal seismic structure summary marine refraction studies during the past decade have found considerable lateral variability in the seismic properties of the upper basaltic regions of the oceanic crust in many localities compressional and shear wave velocities are quite low at the top of the basalt section and velocities increase rapidly with depth it is concluded that pore pressure may be at least in part responsible for the low upper crustal velocities and contribute significantly to the lateral variability this is supported by compressional and shear wave velocity measurements as functions of confining pressure and pore pressure for basalt from the juan de fuca ridge and dolerite from the samail ophiolite oman within the oceanic crust regions of overpressure and underpressure will possess anomalous velocities the magnitude of which will depend upon the porosity and the deviation of the pore pressure from hydrostatic the influence of pore pressure on velocities is expected to diminish with depth and is unlikely to be significant at lower crustal depths where porosity is extremely low of significance poisson s ratio is shown to be dependent on pore pressure as well as confining pressure at constant confining pressure poisson s ratio increases with increasing pore pressure thus overpressured regions within the upper oceanic crust are likely to have relatively high poisson s ratios as well as low compressional and shear wave velocities

crustal accretion at the reykjanes ridge 61 62 n we report results of a seismic gravity and magnetic survey of the reykjanes ridge spreading center at 61 62 n about 600 km from the center of the iceland mantle plume anomalously shallow water on the ridge crest enabled us to record seismic refractions on a 2 4 km hydrophone streamer the velocity within layer 2a is 2 4 0 3 km s 1 and its mean thickness is 400 100 m the velocity at the base of layer 2a is 3 3 0 3 km s 1 on the ridge axis increasing with crustal age to 4 0 km s 1 at 1 5 ma and 4 5 km s 1 at 5 ma assuming that seismic layer 2a on the ridge axis is also the extrusive layer i e the magnetic source layer we have successfully modeled the variations in amplitude of the magnetic field the best magnetic model includes enhanced magnetization within layer 2a at the sites of recent volcanic activity as independently recognized in side scan sonar data we also present a full crustal seismic model based on wide angle seismic recordings on digital ocean bottom hydrophones and disposable sonobuoys the seismic model is complemented by gravity modeling which further suggests that the ridge crest is in isostatic equilibrium the zero age crust is 10 0 km thick while crust of age 5 ma is 7 8 km thick these crustal thicknesses are greater than those of normal oceanic crust which we attribute to the presence of anomalously hot asthenospheric mantle beneath the spreading center we suggest that the variation in thickness between 0 ma and 5 ma crust is caused by temporal variation in the plume fed asthenospheric temperature beneath the reykjanes ridge

upper crustal densities derived from sea floor gravity measurements northern juan de fuca ridge a transect of sea floor gravity stations has been analyzed to determine upper crustal densities on the endeavour segment of the northern juan de fuca ridge data were obtained using alvin along a corridor perpendicular to the axis of spreading over crustal ages from 0 to 800 000 years calculated elevation factors from the gravity data show an abrupt increase in density with age distance for the upper 200 m of crust this density change is interpreted as a systematic reduction in bulk porosity of the upper crustal section from 23 for the axial ridge to 10 for the off axis flanking ridges the porosity decrease is attributed to the collapse and filling of large scale voids as the abyssal hills move out of the crustal formation zone forward modeling of a plausible density structure for the near axis region agrees with the observed anomaly data only if the model includes narrow along strike low density regions adjacent to both inner and outer flanks of the abyssal hills the required low density zones could be regions of systematic upper crustal fracturing and faulting that were mapped by submersible observers and side scan sonar images and whose presence was suggested by the distribution of heat flow data in the same area

relative contribution of surface mass balance and ice flux changes to the accelerated thinning of the mer de glace alps over 1979 2008 by subtracting surface topographies from 1979 1994 2000 and 2008 we measured ice thinning rates increasing from 1 m a 1 1979 1994 to more than 4 m a 1 2000 2008 on the tongue of the mer de glace the relative contributions of changes in surface mass balance and ice fluxes to this acceleration in the thinning are estimated using field and remote sensing measurements between 1979 1994 and 2000 2008 surface mass balance diminished by 1 2 m w e a 1 mainly because of atmospheric warming mass balance changes induced by the growing debris covered area and the evolving glacier hypsometry compensated each other meanwhile the mer de glace slowed down and the ice fluxes through two cross sections at 2200 m a s l and 2050 m a s l decreased by 60 between 1979 1994 and 2000 2008 two thirds of the increase in the thinning rates was caused by reduced ice fluxes and one third by rising surface ablation however those numbers need to be interpreted cautiously given our inability to respect mass conservation below our upper cross section an important implication is that large errors would occur if one term of the continuity equation e g surface mass balance was deduced from the two others e g elevation and ice flux changes

upper crustal velocity structure in the rose area of the east pacific rise we describe a method to constrain the seismic velocity structure of the upper 500 800 m of young oceanic crust conventional interpretations of seismic refraction profiles provide little or no information about this shallowmost crust the method described here requires the assumption of a linear velocity gradient and utilizes measured values of range travel time and slowness of the observed refracted waves results are reported from several refraction profiles that sample crust between 0 and 4 m y in age on the flanks of the east pacific rise at 12 n these data give an average gradient in the uppermost 800 m of the crust of 3 5 s 1 and a velocity at the seafloor of 2 5 km s ray path and travel time modeling show that the refracted waves sampling the uppermost part of the crust are effectively masked by reflections from the seafloor

variations in oceanic crustal structure and implications for the fine scale magnetic anomaly signal marine magnetic anomalies are traditionally modelled using a constant thickness source layer usually defined as the extrusive lavas associated with seismic layer 2a recent seismic studies reveal that the thickness of layer 2a can vary by more than a factor of two an observation which has important implications for the magnetic anomaly source layer we model deep tow magnetometer data from the endeavour segment of the northern juan de fuca ridge in combination with nearby seismic data and find a positive correlation between the thickness of the seismically defined layer 2a and the magnetic anomaly source layer this suggests that short wavelength magnetic anomalies observed within individual magnetic polarity units can arise from variations in the thickness of the extrusives rather than geomagnetic field variations or secondary magnetization

parsimonious wave equation travel time inversion for refraction waves we present a parsimonious wave equation travel time inversion technique for refraction waves a dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line with n geophones evenly deployed along the line these two reciprocal shots contain approximately 2n refraction travel times which can be spawned into o n2 refraction travel times by an interferometric transformation then these virtual refraction travel times are used with a source wavelet to create n virtual refraction shot gathers which are the input data for wave equation travel time inversion numerical results show that the parsimonious wave equation travel time tomogram has about the same accuracy as the tomogram computed by standard wave equation travel time inversion the most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location

structure of young oceanic crust at 13 n on the east pacific rise from expanding spread profiles we present the results of the analysis of expanding spread profiles esps collected on and near the axis of the east pacific rise at 13 n these profiles were collected at 0 1 1 2 1 3 6 and 9 5 km from the rise axis and all but the most distant profile show a distinct low velocity zone lvz located within layer 3 of the oceanic crust at the ridge crest the top of the magma chamber is at the base of layer 2 while 3 6 km off axis the roof of the lvz is 1 1 km below the top of layer 3 the profile farthest from the ridge could possibly have a residual lvz confined to the lower 1 5 km of the crust the total width of the lvz as determined from the esp data is at least 6 km and possibly much greater this wide lvz apparently contradicts multichannel seismic data which show cross axis reflections from the magma chamber with a width of 5 km we suggest that a resolution of this apparent contradiction lies in a model of the rise axis with a small and transient central magma chamber of high partial melt fraction surrounded by a much larger and permanent region of hot rock with only isolated pockets of partial melt the esp data are sensitive to this larger region while the reflection data accurately map the presence or absence of the central magma chamber with its high impedance contrast we identify the presence of a layer at the top of the oceanic crust with initial p wave velocities between 2 35 and 2 6 km s while the s wave velocity is estimated as being 0 8 km s the layer thickness lies between 100 and 200 m these velocities are consistent with previous estimates for the pacific and recent results for the atlantic the thickness of this layer is consistent with that of layer 2a determined from geophysical measurements at deep sea drilling project hole 504b

mass balance of rhonegletscher during 1882 83 1986 87 les bilans de masse effectues sur le glacier des sources du rhone alpes suisses pendant differentes periodes 1884 85 1908 09 1979 80 et 1981 82 sont presentes