Earth and Planetary Science Letters

Transcription

1 Earth and Planetary Science Letters 278 (2009) Contents lists available at ScienceDirect Earth and Planetary Science Letters journal homepage: Mantle melting and magma supply to the Southeast Indian Ridge: The roles of lithology and melting conditions from U-series disequilibria Chris J. Russo a,b,, Ken H. Rubin b, David W. Graham a a College of Oceanic and Atmospheric Sciences, Oregon State University, Corvallis, OR USA b School of Ocean and Earth Science and Technology, University of Hawaii, Honolulu, HI USA article info abstract Article history: Received 19 May 2008 Received in revised form 14 November 2008 Accepted 14 November 2008 Available online 30 December 2008 Editor: R.W. Carlson Keywords: U-series disequilibria mid-ocean ridge basalt Southeast Indian Ridge mantle melting mantle lithology pyroxenite Intermediate-spreading Southeast Indian Ridge (SEIR) basalts display geographic gradients from 90 E to 118 E in Th/U, ( 230 Th/ 232 Th) and 238 U 230 Th 226 Ra disequilibria ( 230 Th-excesses of 1 to 21% and 226 Raexcesses of 0 to 160%). Highly correlated ( 238 U/ 232 Th) and ( 230 Th/ 232 Th) (r 2 =0.94) span nearly the entire global MORB range; basalts from three of four ridge morphologies form subparallel, vertically stacked arrays on an equiline diagram. ( 226 Ra/ 230 Th) is inversely correlated with ( 230 Th/ 238 U) within sub-regions of the study area. There is a convex ( 230 Th/ 238 U) along-axis pattern (lower 230 Th-excesses in the west and east, at shallowest and deepest depths), such that ( 230 Th/ 238 U) does not correlate with axial depth as described elsewhere, despite the ~2300 m west-to-east axial depth increase. κ and κ Pb also vary along axis from κbκ Pb in the west to κ κ Pb in the east. Variations in axial depth, crustal thickness, U-series disequilibria, κ and κ Pb are best explained by forward model scenarios in which the physical characteristics of melting (such as melting rate, degree of melting, porosity and melt initiation depth) vary along axis in response to an inferred long wavelength temperature gradient and small variations in the underlying source lithology (pyroxenite veins). This leads to decreased melt supply from west to east, a concomitant change in axial morphology and axial magma chamber depth, and systematic long wavelength geochemical gradients. Melting rate covaries with melt supply in the western and central regions, decreasing eastward along the SEIR, but supply and rate vary inversely from the central to the eastern regions, because cooler eastern mantle contains a few percent of fusible pyroxenite veins that melt productively to generate melts with high Th U Ra concentration but low 230 Th-excess and 226 Ra-excess. The veins contribute only a small proportion of the total melt volume from mantle that is N95% peridotite. Th, Pb and He isotopic gradients are consistent with a changing proportion of enriched pyroxenite veins along axis having a roughly 0.5 to 1 Gyr reservoir age Elsevier B.V. All rights reserved. 1. Introduction Corresponding author. School of Ocean Engineering Science and Technology, University of Hawaii, Honolulu, HI USA. Tel.: address: crusso@soest.hawaii.edu (C.J. Russo). Magmas produced beneath mid-ocean ridges (MORs) are indirect probes of upper mantle composition and melting conditions. Geographic variations in the conditions and extent of melting, and their linkages to mantle temperature, lithology and chemical composition are central to an understanding of mantle evolution and the origin of mid-ocean ridge basalts (MORBs). Important patterns in melt generation are revealed by the relations between MORB chemistry and physical attributes of MORs (e.g., Klein and Langmuir, 1987; Langmuir et al., 1992; Niu and Batiza, 1993; Shen and Forsyth, 1995; Bourdon et al., 1996b; Niu and Hekinian, 1997; Rubin and Sinton, 2007; Niu and O'Hara, 2008). Global studies of MORB composition must often make simplifying assumptions to relate compositional variations to changes in a particular physical characteristic, such as axial depth or spreading rate (Klein and Langmuir, 1987; Bourdon et al., 1996a; Lundstrom et al., 1998). This approach makes it difficult to decipher the competing effects of several independent parameters on MOR magmatism. Regional scale studies are a useful complement to global studies, particularly in regions with a limited variation in one or more key variables (e.g., spreading rate, axial depth), or where physical, morphologic and/or chemical gradients in ridge character result from a coherent geologic history. Here we present results of a regional scale U-series isotopic and geochemical study of melting conditions at the intermediate spreading rate Southeast Indian Ridge (SEIR) between 90 E and 118 E, an ~2100 km long region displaying significant along-axis physical and chemical variations that can be used to constrain mantle melting and source lithology. Uranium series disequilibria in MORBs record variations in U Th Ra elemental fractionation due to the type (e.g., fractional, batch), extent, rate and depth of melting beneath MORs (e.g., McKenzie, 1985, 2000; Williams and Gill, 1989; Spiegelman and Elliott, 1993; Bourdon et al., 1996a,b; Stracke et al., 1999; Spiegelman, 2000; Sims et al., X/$ see front matter 2008 Elsevier B.V. All rights reserved. doi: /j.epsl

2 56 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) ). U-series disequilibria also reflect MORB source mineralogy (e.g., garnet modal abundance, which dominates U Th fractionation, and clinopyroxene mode, which largely controls fertility and thus melting rate; Beattie, 1993; LaTourette et al., 1993; Pertermann and Hirschmann, 2003). There is a temporal aspect to observed U-series disequilibria, because parent-daughter activity ratios relax to a state of secular equilibrium at a rate governed by their half-lives following U Th Ra fractionation. We use these attributes in conjunction with chemical (Sours-Page, 2000) and radiogenic isotope compositions (Graham et al., 2001, 2006; Mahoney et al., 2002) of the same SEIR samples to test scenarios for physical variations in melting conditions of a heterogeneous mantle source. 2. Regional geologic background The modern SEIR (Fig. 1) marks the Australian-Antarctic plate boundary between the Rodrigues triple junction (25 S, 70 E) and the Macquarie triple junction (62 S, 151 E). This section of MOR contains no large transform offsets and exhibits a west-to-east gradient in axial depth, from ~2300 m to N4500 m (Cochran et al., 1997). The SEIR full spreading rate is nearly constant, from ~72 to 76 mm/yr across our study area (Cochran et al., 1997). The SEIR has a bathymetric high at the Amsterdam-St. Paul Plateau (ASP, E; Graham et al., 2001) and unusually deep, magma-starved ridges within the Australian- Antarctic Discordance (AAD; e.g., Klein et al., 1991; Sempéré et al., 1991). The west to east axial depth increase has been attributed to a long wavelength decrease in mantle temperature (Cochran et al., 1997; Sempéré et al., 1997), estimated to be ~ C using major element variations in N3,000 SEIR MORB glasses from the ASP (Douglas-Priebe, 1998) to the AAD (Klein et al., 1991; Pyle, 1994; Sours-Page, 2000) Along-axis morphologic variations Ridge morphology and crustal thickness variations indicate changes in melt supply along this portion of the SEIR. Three distinct axial morphological zones were discovered by a detailed geophysical survey in 1994/1995 (Cochran et al., 1997; Sempéré et al., 1997). Most western zone ridge segments (C17 C13; Fig. 1) have axial-high morphology, similar in dimension and shape to many segments of the fast-spreading East Pacific Rise (EPR). Segments C12 C4 have a transitional morphology typical of other intermediate-spreading ridges (e.g., Canales et al., 1997). Cochran et al. (1997) subdivided the transitional morphologies into: (a) rifted axial highs in the west (segments C12 C9), and (b) shallow axial valleys in the east (segments C8 C4). SEIR segments C3 and C2 have a deep axial valley, similar to much of the mid-atlantic ridge (MAR). We use the following abbreviations throughout this paper to refer Fig. 1. Satellite altimetry of the Indian Ocean basin between 40 E and 130 E and sample location map. Labeled features include: the Southwest Indian Ridge (SWIR), Central Indian Ridge (CIR), Rodrigues Triple Junction (RTJ), Amsterdam-St. Paul Plateau (ASP), Kerguelen Plateau (KP), Southeast Indian Ridge (SEIR), and the Australian-Antarctic Discordance (AAD). Enlarged panel defines ridge-axis trace of the SEIR from 40 S to 52 S and from 88 E to 118 E. Individual ridge segments are identified above the axial trace (C17 through C2), dredge locations for samples in this study are shown by black dots along the ridge-axis and are labeled below the axial trace (D71 through D145). Ridge-axis color in the online version indicates depth variation with red shades representing the shallowest segments and blue the deepest. We consider three geographical regions in this study (western, WG, central, CG, and eastern, EG, groups). The transform fault separating segments C13 and C12 forms a natural WG-CG divide. MORBs west of this transform tend to have higher 3 He/ 4 He, and different Sr, Nd, Pb isotope compositions, suggesting different mantle geochemical histories east and west of 100 E (Mahoney et al., 2002). The transition into deep-axial valley morphology between segments C3 and C4 forms the EG and CG divide.

3 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) to SEIR basalts from these morphologic groups: Axial High (AH), Rifted Axial High (RH), Shallow Axial Valley (SV), and Deep Axial Valley (DV). There are significant variations in upper crustal seismic structure along segments C12 C7 (Baran et al., 2005) and a crustal thickness gradient from ~6.3 km at 101 E to ~3.8 km at 115 E (Holmes et al., in press). The westernmost RH segment (C12) has the shallowest axial magma chamber (AMC), thinnest layer 2A (Baran et al., 2005) and greatest crustal thickness (~6.3 km between 100 and 116 E; Holmes et al., in press). Baran et al. (2005) divided segment C12 into a western AH and eastern RH segment, which we adopt for the dredge 96 and 100 morphological assignments (Table 1). There is currently no seismic data for the other AH segments (C17 C13). The AMC is ~1 km deeper and the crust is slightly thinner (~5.9 km) on RH segments. The SV segments do not have detectable AMCs. The DV segments have the thinnest crust (~4.6 km). They are outside the Baran et al. (2005) study area but likely lack an AMC (like most of the slow-spreading MAR; Detrick et al., 1990) Along-axis geochemical variations All SEIR lavas discussed here are Indian -type MORBs (Mahoney et al., 2002), with lower ε Nd,higher 87 Sr/ 86 Sr and elevated Δ 208 Pb/ 204 Pb (after Hart, 1984) compared to Pacific and North Atlantic MORBs. Eastward decreases in Fe 8 and increases in Na 8 suggest a gradient in the mean pressure and extent of melting (Sours-Page, 2000; Graham et al., 2001; Mahoney et al., 2002; Fe 8 and Na 8 are fractionation corrected oxide abundances, Klein and Langmuir, 1987). Broadly correlated 3 He/ 4 He and Fe 8 suggest that high 3 He/ 4 He mantle is preferentially sampled by deeper melting in the west (Graham et al., 2001). Conversely, the most strongly Indian -type mantle is associated with deeper axial depths and by inference, lower degrees of melting in the east (Mahoney et al., 2002). A heterogeneous upper mantle with lower 3 He/ 4 He and higher Δ 208 Pb/ 204 Pb material concentrated in the shallowest mantle and preferentially sampled by shallower melting at lower mantle temperature in the east, explains the Fe 8 3 He/ 4 He Δ 208 Pb/ 204 Pb relationship (Graham et al., 2001; Mahoney et al., 2002). Sr Nd Pb isotope data indicate little, if any, eastward mantle flow from the distant ASP or Kerguelen Heard hotspots into the study region (Mahoney et al., 2002; Nicolaysen et al., 2007). Loose correlations of Sr, Nd, Pb and Hf isotopic compositions with their respective parent/daughter ratios have been interpreted as either mixing trends or Myr reservoir ages (perhaps reflecting small amounts of trapped mantle melt from Gondwana breakup; Mahoney et al., 2002; Graham et al., 2006). 3. Samples and methods Sixteen basalt glass samples from 13 ridge segments, selected to represent the regional morphologic, geographic and geochemical variations within the study area (Fig. 1 and Table 1), were analyzed by thermal ionization mass spectrometry for Th U Ba concentrations and U and Th isotopic compositions; a 12 sample subset was analyzed for 226 Ra concentration. The MORB U-series analytical protocols used in this study are described in Rubin et al., 2005 (and references therein). Alteration-free glass samples were handpicked under a binocular microscope (see Table S6 for sample conditions and weights; the letter S indicates tables or figures that are in the supplement). Analytical blank and standards meta-data are reported in the Table 1 footnotes. Scenarios for producing observed U-series disequilibria were forward modeled, using both dynamic melting (McKenzie, 1985, Williams and Gill, 1989) and equilibrium porous flow (Spiegelman and Elliott, 1993; Spiegelman, 2000) models. Multiple model simulations were run to examine the relative roles of melt fraction (F), mantle porosity (ϕ), and melting rate (M) for three mantle lithologies, as described in Section 5.3 and the supplement. 4. Results Analytical results are presented in Table 1; activity ratios are denoted by parentheses in the table and throughout the text. ( 234 U/ 238 U) is effectively indistinguishable from unity in each of our samples (0.998±0.004 to 1.005±0.003), indicating that post-eruptive alteration/exchange with seawater has not affected U-series activities in this sample set Ba Th U concentrations Ba, Th, and U concentrations vary considerably in SEIR basalts (by factors of 8, 9 and 6, respectively) with much less variation in Ba/Th (2.3 times) and Th/U (1.6 times). Higher concentrations are more common in the east. Ba, Th, and U concentrations are each strongly correlated (r 2 N0.9) with (La/Sm) n over a threefold (La/Sm) n range ( 230 Th/ 232 Th) and ( 238 U/ 232 Th) SEIR basalts form a continuous array on an equiline diagram, from enriched (low) to depleted (high) ( 238 U/ 232 Th) values (Fig. 2) and span a larger range than any other ridge except the ultra-slow spreading Southwest Indian Ridge (SWIR), which has a similar range but noncontinuous data distribution (Standish, 2006). The SEIR ( 238 U/ 232 Th) range is greater than other intermediate spreading rate ridges such as the Gorda and Juan de Fuca (JDF) Ridges (Goldstein et al.,1992) and the eastern Galápagos Spreading Center (GSC; Kokfelt et al., 2005). Dividing the SEIR data into western, central and eastern geographic groups (WG, CG and EG; Table 1 and Fig. 2) reveals regional differences in our sample suite. SEIR basalts can also be divided into three slightly different relative source enrichment groups based on Th/U, forming distinct fields on an equiline diagram (Fig. 3A): enriched ( 3.4, n=5), intermediate (2.8bTh/Ub3.3, n=7), and depleted ( 2.4, n=4). Thorium concentrations in the intermediate and enriched groups overlap (0.17 to 0.58 and 0.25 to 1.16 µg/g, respectively), with implications for variations in melting style and mantle source composition (see section 5.4). 230 Th 238 U disequilibria vary systematically along the SEIR. ( 230 Th/ 232 Th) and ( 238 U/ 232 Th) are highly correlated in samples having 226 Raexcesses (r 2 =0.94; Fig. 3A). Because the data lie to the left of the equiline (the locus of secular equilibrium), the linear regression slope of 0.90 indicates slightly higher average 230 Th-excesses at the enriched end (~16%) than the depleted end (~6%) of this array. A similar relationship in other MORBs has been loosely interpreted to reflect a lithologic control on melting conditions (e.g., Goldstein et al., 1992; Rubin and Macdougall, 1992; Lundstrom et al., 1998) Post-eruptive decay effects on ( 230 Th/ 238 U) and ( 226 Ra/ 230 Th) Sample age is particularly important in U-series studies because initial disequilibria decay away at a rate governed by the shorter halflife of a parent-daughter pair (75.69 kyr and 1.6 kyr for 230 Th and 226 Ra, respectively). It is difficult to ensure that most MORBs are young enough to be unaffected by post-eruptive decay (especially those sampled by dredging, such as the SEIR basalts). Excess 226 Ra, which is commonly observed in young MORBs, occurs for ~8 kyr following eruption (e.g., Rubin and Macdougall, 1988), making it a good indicator of sample youth (i.e., ( 230 Th/ 238 U) will not have been changed by post-eruptive decay in lavas with ( 226 Ra/ 230 Th)N1). Assessing sample youth was the primary purpose of 226 Ra analyses in this study. Eleven of twelve SEIR basalts analyzed for 226 Ra have excess 226 Ra (ranging from 15±4% to 161±4%) and lie on or within a computed error envelope around the Fig. 3A linear regression. The

4 Table 1 Southeast Indian Ridge U-series disequilibria and Th, U, Ra and Ba concentrations measured by thermal ionization mass spectrometry (TIMS) Sample (WW10-XX-x) Lat (S) Long (E) Depth (m) EorN MORB Ridge morphology [Th] ng/g ± [U] ng/g ± ( 234 U/ 238 U) ± 232 Th/ 230 Th ± ( 230 Th/ 232 Th) ± ( 238 U/ 232 Th) ± ( 230 Th/ 238 U) ± Ra fg/g ± ( 226 Ra/ 230 Th) ± Ba µg/g Western Group (WG) N AH N AH N AH N AH E RH WG averages Central group (CG) E AH E RH E RH N SV E RH N SV E SV E SV E SV CG averages Eastern group (EG) N DV E DV EG averages Table notes: 1. Isotope dilution concentrations and isotopic compositions were measured on the Univ. of Hawaii VG Sector 54 thermal ionization mass spectrometer equipped with a high abundance sensitivity (WARP) filter and ion-counting Daly detector (see also Rubin et al., 2005). 2. All measurements were made on splits from the same sample dissolution. Reported errors are 2σ. 3. Activities are denoted by ( ); Decay constants used to calculate activities are: λ 226 = yr 1, λ 230 = yr 1, λ 232 = yr 1, λ 234 = yr 1, λ 238 = yr Typical procedural blanks on concentration measurements were [Ba]= g; [Th]= g; [U]= g. 5. Duplicate 232 Th/ 230 Th measurements on the same dissolutions for samples denoted by are given in the supplementary material. 6. Ridge morphologies are AH=axial high, RH=rifted high, SV=shallow valley, DV=deep valley; see section 2 of the text. 7. Sea water alteration of samples can be detected with ( 234 U/ 238 U), which is 1 in fresh, unaltered lavas and 1.14 in seawater (e.g., Thurber, 1967). ( 234 U/ 238 U) is effectively indistinguishable from unity in each of our samples. 8. Regional Group averages are given only for parameters that are discussed in the text. The WG ( 226 Ra/ 230 Th) average does not include the secular equilibrium sample. 9. Mean values for replicate analyses of Th (UCSC-A) and U (CRM112A) isotope standards are; 232 Th/ 230 Th= ±2844 (2σ, n=31) and 234 U/ 238 U= ± (2σ, n=10), respectively. Analyses of the in-house rock standard K1919 conducted during this study are reported in Rubin et al. (2005). 58 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) 55 66

5 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) MORBs (~4). Instead, this basalt is likely substantially younger than 1600 years. ( 226 Ra/ 230 Th) is much lower in the other WG/EG basalts. It is difficult to tell a priori whether this is because these samples are from older lava flows or if lower ( 226 Ra/ 230 Th) reflects different melt generation conditions beneath these ridge segments. This latter possibility is considered below in conjunction with ( 230 Th/ 238 U) in the lavas and melt model predictions, but we stress that primary 226 Ra-excesses are not well-constrained in our data set, so we rely more heavily on 230 Th-excesses to quantify melting conditions. 5. Discussion 5.1. U-series variations along the SEIR U-series signatures vary along the SEIR, both within and between the aforementioned WG, CG and EG geographic groups (Figs. 2, 3A, 4 6A). Fig. 2. Equiline diagram. ( 230 Th/ 232 Th) ( 238 U/ 232 Th) diagram with SEIR basalts plotted by geographic groupings. The diagonal line in the center of the plot is the locus of secular equilibrium. Data that sit to the left of this line display excess 230 Th. Global MORB data are shown for comparison from Elliott and Spiegelman (2003) with additions from Kokfelt et al. (2005), Standish (2006) and Tepley et al. (2004). three CG and one WG samples not analyzed for 226 Ra, plus the one sample (WW ) with ( 226 Ra/ 230 Th)=1, lie within or above this envelope, indicating that post-eruptive decay has not significantly lowered their 230 Th excesses compared to those samples with ( 226 Ra/ 230 Th)N1. Thus we believe that our measured 230 Th 238 U disequilibria record magmatic values at the time of eruption. One should also evaluate the potential of post-eruptive decay on ( 226 Ra/ 230 Th) to use it with ( 230 Th/ 238 U) to constrain magmatic processes and timing, as we attempt to do here. A few studies have used MORBs of known eruption age or geologic context from submersible collection to largely avoid post-eruptive 226 Ra decay effects (e.g., Sims et al., 2002; Rubin et al., 2005). Most MORB studies, including this one, use dredged samples, for which age constraints and field relationships are limited, yet it is still possible that ( 226 Ra/ 230 Th) disequilibria retain aspects of disequilibria relationships at eruption, particularly when U-series patterns resemble those in age-constrained, submersible collected samples. There is an inverse relationship of ( 226 Ra/ 230 Th) and ( 230 Th/ 238 U) in age-constrained EPR basalts from the 9 10 N (Lundstrom et al., 1999; Sims et al., 2002); our SEIR basalts exhibit the same inverse relationship, albeit at lower ( 226 Ra/ 230 Th) (Fig. 4). SEIR CG samples are correlated on this diagram and fall close to a line that represents 1 half life of decay (1600 years) from the EPR trajectory (although one datum is displaced an additional half life below this line). Such a wellcorrelated trajectory is unlikely to result from randomly dispersed initial ratios experiencing different amounts of decay. One can therefore speculate that this CG data array mimics the data dispersion as it was upon eruption, implying that the conditions that arose in the EPR trajectory likewise exist along this part of the SEIR. However, it is not possible with the present data set to determine whether these CG basalts initially followed the EPR array and have since decayed, or if zero-age SEIR lavas are offset to lower ( 226 Ra/ 230 Th) than the EPR. A similar explanation probably does not apply to WG and EG lavas, which form a separate, steeper array at generally lower 226 Ra-excess for a given 230 Th-excess (Fig. 4). Post-eruptive aging of a simple, EPRlike array could only produce this if sample aging increased smoothly from basalts having low ( 230 Th/ 238 U) to those having high ( 230 Th/ 238 U). One WG sample (84 7) has 160% 226 Ra-excess; if this sample has decayed by even one 226 Ra half life, it would have erupted with ( 226 Ra/ 230 Th)=4.2, making it one of the highest values observed in Fig. 3. A) SEIR Equiline diagram. SEIR basalts plotted by regional groupings; symbols are the same as in Fig. 2; black dots additionally indicate samples with measured 226 Ra and thus agecontrol. Regional groupsare furtherdivided into enrichmentgroupsdesignatedby dashed vertical lines. The solid black line is a least-squares regression through the data having measured 226 Ra, following Lundstrom et al. (1998). Curved grey lines are the ±2 standard error envelope on the linear regression. B) Equiline diagram and axial morphology. SEIR basalts grouped by axial morphology as described in Cochran et al. (1997) and Baran et al. (2005). Lines are least squares regressions through the data (with radium control) for each morphologic group; 4 samples that may have suffered 230 Thdecay are excluded. Note the generally lower 230 Th-excess of AH lavas and the progression to lower 230 Th-excess for RH, SV and DV groups. AH morphologies are EPR-like, with m axial highs (segments C17-C14), RH segments have m axial highs with ~2 km wide by m deep rifts (segments C13-C9), SV segments have km wide by m deep rifts (segments C8-C4), and DV morphologies having well-defined (N1000 m deep) axial valleys.

6 60 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) sample 84 7). Additionally, the six CG samples have the highest ( 226 Ra/ 230 Th) as a group, the WG samples have intermediate values, and the two EG samples have the lowest mean value (Table 1). Some of these differences may be due to post-eruptive 226 Ra decay (section 4.3), but we also consider the implications of these regional Ra variations assuming some or all are magmatic Axial depth and 230 Th/ 238 U disequilibria Fig. 4. ( 226 Ra/ 230 Th) versus ( 230 Th/ 238 U). (A) SEIR basalts separated into regional groups. Compare the inverse relationships in SEIR basalts to submersible collected, ageconstrained samples from the East Pacific Rise, plotted as black circles (Lundstrom et al., 1999; Sims et al., 2002). SEIR regional groupings fall along offset trends, with central group (CG) lavas lying along a trend at relatively higher values of ( 226 Ra/ 230 Th) than western group (WG) and eastern group (WG) lavas. Some of this difference is probably due to post-eruptive decay of 226 Ra, although decay is not expected to move randomly distributed zero-age samples onto a well correlated line (such as the CG samples). The solid black line is a least squares regression through the EPR data. The dashed line is the slope of the EPR trend after one half-life of decay. (B) CG measured data, and data corrected for potential 226 Ra decay. Much of the CG data follows this same relationship as the EPR if corrected for one half-life of decay (open symbols). Correction for two halflives would raise the lowest CG sample to the same trend as the others. (C) EG and WG measured data and data corrected for potential 226 Ra decay. The WG and EG samples have generally lower 226 Ra-excess, which may be due either to more extensive decay than in CG samples or to different melting conditions. All basalts have 230 Th-excesses, ranging from 21±1.2% to 1.4±1.0%. This 20% 230 Th-excess range is smaller than in young (i.e., with excess 226 Ra), intermediate spreading rate MORBs from the JDF-Gorda Ridges ( 230 Thexcess of 4% to 40%; Goldstein et al., 1992; Volpe and Goldstein, 1993), but slightly greater than GSC MORBs (5% to 21%; Kokfelt et al., 2005). The along-axis SEIR ( 230 Th/ 238 U) variation forms a slightly convex pattern, with lower values associated with the western and eastern ends of the study area (Fig. 5). The CG has much larger 230 Th-excesses (16% on average) compared to lower WG and EG values (8% and 11%, respectively; Table 1). SEIR ( 230 Th/ 232 Th) shows a different geographic pattern, generally decreasing from WG to CG to EG segments and primarily reflecting mantle source composition before melting. Aspects of ( 226 Ra/ 230 Th) variations in the sample set likely reflect initial conditions at the time of eruption (e.g., the CG ( 230 Th/ 238 U) vs. ( 226 Ra/ 230 Th) inverse correlation and the high ( 226 Ra/ 230 Th) in WG Bourdon et al. (1996a) interpreted a weak negative correlation between ridge segment averaged axial depth and ( 230 Th/ 238 U) in a substantially filtered early global MORB data set to result from mantle temperature control on depth/degree of melting. Their model involves melting of hotter mantle beneath shallower ridges producing larger ( 230 Th/ 238 U) via a deeper peridotite solidus, higher overall F, thicker oceanic crust, and more melt generated in the garnet stability field (where Th is more incompatible then U; Beattie, 1993; LaTourrette et al., 1993; Blundy and Wood, 2003). The ( 230 Th/ 238 U) vs. axial depth relationship is less evident when newer data from the deep southern MAR and the shallow Reykjanes Ridge are included in the global compilation (Lundstrom et al., 1998; Peate et al., 2001; Lundstrom, 2003). MORBs from intermediate axial depths (~ m) are highly variable and span nearly the entire MORB range (Fig. 6). The global data set still displays a weak negative covariation but the statistical significance of the global trend has been questioned, and it is largely defined by end-member samples from the shallow Azores platform (MAR) and the unusually deep AAD (Elliott and Spiegelman, 2003). The two deepest samples from the AAD have not been analyzed for 226 Ra, and since their ( 230 Th/ 238 U) is substantially different from our nearest EG samples, there are legitimate questions concerning whether they have experienced significant post-eruptive 230 Th decay. Axial depth and ( 230 Th/ 238 U) are not negatively correlated along the SEIR, although SEIR samples do largely fit within the published MORB swath (Fig. 6). Yet ( 230 Th/ 238 U) displays the entire SEIR range in the narrow axial depth range from 2300 to 2700 m. In addition, many WG basalts have lower than expected 230 Th-excesses for their axial depths and the two EG basalts erupted at the greatest depth have the same ( 230 Th/ 238 U) as CG and WG basalts erupted much shallower. If there truly is a weak global trend, the lack of a ( 230 Th/ 238 U)-axial depth correlation along the SEIR suggests that melt initiation depth does not respond to regional mantle temperature differences in the way envisioned to explain the global data set. Other variables, including M, ϕ and mantle composition likely vary along the SEIR in a way that obscures any simple potential ( 230 Th/ 238 U)-axial depth relationship Melt supply, axial morphology and U-series disequilibria General observations and modeling Melt supply (i.e., the volumetric flux of melt from mantle to crust) affects ridge axis morphology, AMC depth and lava chemistry. Thermal models predict that AMC depth and melt supply are inversely related (Phipps Morgan and Chen, 1993; Chen 1996); ridge-axis morphology reflects interplays between melt supply, mantle temperature, and hydrothermal cooling (e.g., Chen, 2000); and strong correlations between lava chemistry, AMC depth and spreading rate demonstrate how such thermal conditions control magma accumulation conditions and MORB compositions (Rubin and Sinton, 2007). Coincident changes in AMC depth and ridge morphology at some ridges also support these thermal models (e.g., at EPR segments north and south of the Orozco transform; Chen, 2000, Carbotte et al., 1998). Intermediate spreading ridges such as the SEIR are generally more sensitive to magma supply variations than ridges with faster or slower spreading rates (e.g., Niu and Batiza 1993; Carbotte et al., 2006; Detrick et al., 2002), and small magma supply variations can produce large variations in ridge axis morphology (Chen, 1996).

7 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) Fig. 5. Longitudinal variations within the SEIR study area. (A) Axial-depth profile and smoothed axial depth curve (grey line) from Sempéré et al. (1997). (B) κ is the measured 232 Th/ 238 U ratio for SEIR basalt glasses. (C) κ Pb calculated from the Pb isotope data of Mahoney et al. (2002); colored data points are samples with Th-U isotope dilution data from this study. κ Pb is the 232 Th/ 238 U ratio assuming closed system terrestrial evolution over 4.56 Ga, with Canyon Diablo troilite 206 Pb/ 204 Pb and 208 Pb/ 204 Pb as initial values (Tatsumoto et al., 1973). (D) ( 230 Th/ 238 U) for SEIR basalt glasses. Grey curve is the ( 230 Th/ 238 U) for the SEIR basalts, predicted from the smoothed axial depth in panel (A) using the global ( 230 Th/ 238 U)-axial depth relationship of Bourdon et al. (1996a). (E) Measured ( 230 Th/ 232 Th). (F) Measured ( 226 Ra/ 230 Th). Shallower axial depth, axial high morphology and greater crustal thickness all point to greater melt supply in the western portion of our SEIR study area. Systematic variations in these parameters as well as in U-series disequilibria (at nearly constant, intermediate spreading rate) provide a unique opportunity to quantify and investigate the underlying melt supply variations. U-series disequilibria in volcanic rocks are often modeled to examine melting conditions and melt supply effects beneath MORs and mid-plate hotspot volcanoes at passive and active mantle upwelling conditions, respectively (e.g., Speigelman and Elliott, 1993, Sims et al., 1995, McKenzie, 2000; Stracke et al., 2006), and at hotspot-influenced ridges (e.g., Bourdon et al., 1996b; Kokfelt et al.,

8 62 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) of melting parameters have the biggest effect on 230 Th-excess in MORB, with excesses being reduced by: 1) greater porosity, ϕ; 2) greater melt fraction, F; 3) higher melting rate, M; and in some instances, 4) shallower melt initiation depth. A wide range of conditions result in excess 226 Ra, but the magnitude of the excess is highly dependent on melting and melt transport rates because the 226 Ra half life is so short. Fig. 6. ( 230 Th/ 238 U) versus axial depth. (A) SEIR basalts are divided into regional groups. Small grey symbols are the global dataset from Bourdon et al. (1996a), as updated by Lundstrom (2003). (B) SEIR basalts grouped by ridge morphology (Symbols are the same as in Fig. 3B). Small grey circles as in panel A. Other data include Galápagos Spreading Center samples with measured 226 Ra-excesses (open circles; Kokfelt et al., 2005). 2005). Two general types of models have been developed that constrain the rate, duration and extent of mantle melting (e.g., Elliott, 1997); both are one-dimensional, consider variations with depth beneath the ridge axis, and produce radioactive disequilibria primarily by variations in M, ϕ and F. The main differences between the models are in specific melting conditions and whether or not melt and solid maintain chemical equilibrium throughout the melting column. Dynamic melting (e.g., McKenzie, 1985; Williams and Gill, 1989) predicts compositions of melts that remain with the solid residue until a porosity threshold, followed by melt extraction everywhere from the melt column at a rate balanced by continued melt production, resulting in constant ϕ. This model behaves like accumulated fractional melting for non-radioactive elements (e.g., Williams and Gill, 1989). The equilibrium porous flow (EPF) model predicts the combined effects of chemical fractionation, radioactive decay, and melt transport time on U-series disequilibria (Spiegelman and Elliott, 1993). M is not an input parameter but can be calculated from EPF model inputs (i.e., depths of melt initiation and cessation, plus the amount of melt present at each depth). Melt simulations described in this paper were run at the same conditions of M, ϕ and F for both types of models (see the Supplementary data). In general these melt models do not uniquely and quantitatively predict melting conditions that lead to specific disequilibria values (e.g., for 230 Th 238 U) because (a) model outputs are highly dependent on poorly constrained input parameters, (b) multiple combinations of parameters can generate the same result, and (c) such parameterizations have greater degrees of freedom when they are used to match just one U-series element pair. Here we use the results of model simulations as qualitative indicators of conditions that may result in disequilibria variations between groups of SEIR samples without emphasizing the absolute values of a particular model simulation. Many previous authors have used these same models to constrain mantle melt column parameters, some more literally than others (e.g., Williams and Gill, 1989; Spiegelman and Elliott, 1993; Lundstrom et al., 1995; Bourdon et al., 1996a,b, 1998; Sims et al., 1995, 2002; Stracke et al., 1999, 2003). These previous studies agree that a handful Melt supply along the western portion of the study area The SEIR samples were subdivided into four groups (Table 1; Section 2.1) that reflect their axial morphology, AMC depth and crustal thickness. These morphologically grouped samples cluster discretely on the ( 230 Th/ 238 U) vs. axial depth diagram (Fig. 6B). The RH, SV and DV groups form an array that resembles the global trend of Bourdon et al. (1996a). However, the AH group (i.e., with highest melt supply) plots well below that trend at lower 230 Th-excess for the same axial depth as RH samples. Shallower AMC depth beneath segment C12 (Baran et al., 2005) is one indication that higher melt supply may be the underlying cause of this RH-to-AH morphological and 230 Thexcess shift. Observations at the intermediate spreading rate Galápagos Spreading Center provide context for these observations. The GSC has a strong axial depth gradient from interaction of the ridge and the Galápagos hotspot, and exhibits AH morphology in a central bulge region nearest the hotspot due to excess magma supply, RH morphology near 94 W, and SV morphology west of 95 W (e.g., Detrick et al., 2002; Christie et al., 2005). MORBs erupted in the central bulge (Fig. 6B) have lower 230 Th-excesses than nearby non- AH segments (Kokfelt et al., 2005). The nearest hotspots to our study area are 700 to 1600 km away and do not directly influence magma supply, so something else must be increasing magma supply to account for the AH morphology, shallower ridge depth and lower 230 Th-excesses than RH segments. In principle, all could result from some combination of larger F, greater ϕ, and faster M. Increased melt productivity (or melting rate) will produce larger F at constant melt column height if the mantle is more fertile or at higher temperature. Nothing from the trace element or isotope geochemistry implies that SEIR AH basalts were produced by more fertile mantle, but a mantle temperature difference with RH segments cannot be ruled out. Producing extra melt by only deepening the melt initiation depth beneath AH segments would lead to a greater proportion of melt generated in the garnet stability field, thereby increasing rather than lowering melt 230 Th-excesses, which is inconsistent with the SEIR observations. Likewise, it is possible to increase melt supply at constant melt column properties if a greater proportion of melt were sampled from off-axis portions of the triangular melting region beneath the ridge, but smaller degree melts produced there would have greater 230 Th-excess. Producing a lower proportion of melt within the garnet stability field would generate lower melt 230 Th-excesses (i.e., if melt initiation depth were shallower beneath AH versus RH segments). For example, 230 Th-excess is reduced by ~25% when initiation depth is decreased from 40 kb to 25 kb (Fig. S4). Garnet no longer influences the melt composition if initial melting depth is reduced further, dramatically reducing melt 230 Th-excess (e.g., compare dynamic and EPF melting simulations of spinel lherzolite and garnet lherzolite in Figs. S1 and S2). Some spinel lherzolite model simulations reproduce the 230 Thexcesses observed in AH and RH samples, but only at unrealistically low values of M, ϕ and F (b kg m 3 yr 1, b0.1% and b0.08, respectively). These conditions are also inconsistent with the relatively low 226 Ra excess in some WG samples (i.e., except WW ). Furthermore the very small degrees of melting in these simulations are inconsistent with 6.3 km thick crust in the western region (see calculations in the supplement). Thus shallow melt initiation depths are unlikely to have caused the low 230 Th-excesses in the WG samples. However, the unusual characteristics of sample

9 C.J. Russo et al. / Earth and Planetary Science Letters 278 (2009) WW (extremely low 1.4% 230 Th-excess and very high 160% 226 Ra-excess) make it the one exception, in that spinel peridotite melting simulations can produce these characteristics and garnet peridotite simulations cannot. This sample was taken from a small dying ridge segment between a westward propagating segment and the transform near 96 E, a likely place for a shallow, low productivity melt column. It is unlikely that melt initiation depth varies significantly beneath most AH and RH segments since changing it cannot produce the observed decrease in average 230 Th-excess at AH segments (6%) relative to nearby RH segments (19%). Assuming constant initial melting depth, one can forward model our U-series observations to qualitatively determine how much melt supply must increase from RH to AH segments to change the axial morphology and reduce the 230 Thexcess. An increase in both M and ϕ will reduce 230 Th-excess, whereas only an increase in M will also increase the amount of melt generated. A seven fold increase in M (from to kg m 3 yr 1 )at constant ϕ (0.25%) and bulk K d, will produce this 230 Th-excess shift by dynamic melting with essentially no effect on 226 Ra-excess (Table S3). This change in M can be reduced slightly to a five-fold increase (from 1 to kg m 3 yr 1 )ifϕ is also varied from ~0.25% (RH) to ~1.0% (AH). In the latter case, 226 Ra-excess would decrease along with 230 Thexcess (Fig. S1), consistent with lower 226 Ra-excess observed in all but one AH sample. The ( 230 Th/ 238 U) difference could also be made by increasing ϕ at constant M, although M would be relatively low ( kg m 3 yr 1 ) and the ϕ shift would be extreme (i.e., 12-fold, from ~0.2%, RH to ~2.5%, AH). It is difficult to imagine a scenario that would allow such large ϕ increase at constant upwelling and melting rate, so we consider model simulations requiring some change in both M and ϕ to be more realistic. We speculate that higher ϕ might increase melt transport efficiency and thus melt supply, but such melt transport mechanics are beyond the scope of this work Melt supply along the eastern portion of the study area It is conceptually straightforward to envision how an increase in M could supply the extra melt needed to produce the U-series, AMC depth and morphology of AH segments on the SEIR. However, this rationale does not explain the lower 230 Th-excesses (and very low 226 Ra-excesses) in EG segments with morphology and crustal thickness indicative of decreased melt supply. SEIR basalts from lowest melt supply DV segments have 230 Th-excesses between those of high melt supply AH segments and moderate melt supply RH and SV segments (i.e., the average 230 Th-excess shift between SV-RH (15%) and DV (11%) segments is smaller than the AH (6%) to RH (19%) shift). Low melt supply morphology and thinner crust both imply that melting conditions beneath DV segments must suppress melt 230 Thexcess without increasing the net melt flux from the mantle. If increased ϕ and/or M were the cause of reduced 230 Th-excess, DV segments would require a shorter melt column than SV segments. For example, reducing the initial melting depth from 40 kb to 25 kb shifts ( 230 Th/ 238 U) from 1.16 to 1.11 with the EPF model at constant M ( ) and ϕ (0.5%) (Fig. S4), because less melt is produced within the garnet stability field. However, a shallower melt initiation depth alone cannot explain the mean 230 Th-excess shift from SV to DV samples because the lavas also require a change to more enriched mantle source composition (e.g., Th/U, radiogenic isotope ratios). The difficulty with generating the combined variations in 230 Thexcess, 226 Ra-excess, Th/U and melt supply for the EG basalts with only physical changes to the melt column leads us to consider compositional changes that can influence melt productivity, Th and U bulk distribution coefficients, and source enrichment Mantle source lithology and 230 Th 238 U disequilibria Mantle composition might affect the magnitude of 230 Th 238 U disequilibrium in MORBs (e.g., Goldstein et al., 1992; Rubin and Macdougall, 1992; Lundstrom et al., 1998). E-MORBs (which often have elevated Th/U) can have larger 230 Th-excesses than N-MORBs from the same region, producing equiline diagram data arrays with slopes b1 (Lundstrom et al., 1998). Lundstrom et al. (1998) postulated that increasing the modal proportion of garnet (to lower the relative Th bulk distribution coefficient) and clinopyroxene (to increase fertility) in the E-MORB source increases 230 Th 238 U fractionation during melting. Melts of such a source would have larger 230 Thexcesses than melts of more depleted mantle, and their mixtures would produce the observed equiline diagram arrays. Variations in mantle upwelling rate can be superimposed on this mantle composition effect to flatten the array slope at slower compared to faster spreading ridges because slower upwelling leads to relatively greater melt 230 Th-excesses by 230 Th production in the melting column. The slope of the full SEIR array (0.90) is significantly steeper than slopes for other intermediate spreading ridges (JDF slope=0.42 and Gorda Ridge slope=0.32; Lundstrom et al., 1998). The SEIR slope is, in fact, steeper than the steepest slope (0.82) for the superfast spreading EPR at S, where the high upwelling rate apparently counteracts any source lithologic effect on melt ( 230 Th/ 238 U) (Lundstrom et al., 1998). Given the nearly constant intermediate spreading rate and large Th/U range (plus other isotopic variations), our Fig. 3A data array slope indicates that different mantle compositions beneath the SEIR do not melt in the manner predicted by Lundstrom et al. (1998). Furthermore, generally higher Fe 8.0 in SEIR WG lavas (which also have lower Th/U) are inconsistent with requirements of the Lundstrom et al. model (i.e., shallower melt initiation depth of a less productive mantle source to produce lower 230 Th-excess), as are the inflated axial morphology and shallow ridge depth in the area. N-MORBs and E-MORBs (defined as having (La/Sm) N N0.80; Mahoney et al., 2002) occur throughout the SEIR study area with no systematic offset in ( 230 Th/ 238 U) or ( 230 Th/ 232 Th) (Figs. S5a and S5b). Thus the origin of the large Th/U (and Ba/Th and La/Sm) range along the SEIR (e.g., Fig. 3A) must be different from that which manifests itself as Th/U and ( 230 Th/ 238 U) variations in N-MORBs and E-MORBs elsewhere. Although sample numbers are small in this study, there are systematic variations in source composition, melt supply and ridge morphology that suggest a link between variations in mantle lithology, mantle temperature and melting conditions along the SEIR axis. Average Th/U of each of our morphological groups indicate different mean mantle compositions beneath ridge sections of some morphologies (Th/U is 2.5±0.2, 3.2±0.3, 3.3±0.3, 3.6±0.3 along AH, RH, SV, DV segments, respectively). There are also 230 Th 238 U disequilibrium differences between our various morphologic groups. The AH segment has the lowest Th/U and 230 Th-excess (6%). The other three morphologic groups (RH, SV and DV) form vertically stacked arrays on an equiline diagram, with decreasing average 230 Th-excesses (i.e., 19%, 15% and 11%, respectively; Fig. 3B). Because it is difficult to fractionate Th from U during MORB mantle melting (e.g., Elliott, 1997), decreasing 230 Th-excesses at increasing Th/U in these latter three groups implies a change in mantle composition (e.g., lithology, which can affect Th/U) and how that composition melts (which affects 230 Thexcess). The trace element and isotopic gradients along this part of the SEIR also imply a gradient in mean mantle composition, such as would be produced by variable proportions of depleted and enriched lithologies. Least squares regression slopes for the individual morphological data arrays on Fig. 3B are 1 for the moderate to high melt supply groups (SV, RH, AH) and N1 for the DV group. The more enriched samples in the SV, RH and AH groups have slightly larger 230 Th-excess than the depleted samples, whereas for the DV group the more enriched sample has smaller 230 Th-excess. The stacked RH, SV and DV data arrays rotate toward steeper slopes as the average ( 238 U/ 232 Th) decreases and as the SEIR deepens to the east, so that increasing the proportion of melts from mantle with higher (i.e., less depleted)