Magmatically assisted off-rift extension—The case for broadly distributed strain accommodation

As a continental rift progresses toward an oceanic spreading center, the mechanism of strain accommodation must transition from faulting and thinning of the lithosphere to focused magmatic intrusion (Buck, 2004; Buck, 2006; Casey et al., 2006; Corti, 2009, 2012; Mazzarini et al., 2013; Rooney, 2010; Rooney et al., 2014; Wolfenden et al., 2005). Specifically, strain is thought to migrate from high-angle border faults to zones of focused magmatic activity within the rift that are adjacent to the rift-border faults (Ebinger and Casey, 2001). Though the specific mechanisms remain controversial, these zones of focused magmatic intrusion must eventually migrate toward the rift axis as possible precursors to oceanic spreading centers (Ebinger, 2005). While elegant in its simplicity, observations from rifting environments reveal magmatic features that are unexplained by this model. For example, magmatism outside of the rift proper is frequently observed to be contemporaneous with episodes of both mechanical extension and focused magmatic intrusion within the rift (Abebe et al., 1998; Keranen and Klemperer, 2008). Seemingly related to the rifting process, this magmatism is not adequately explained by the currently accepted models of continental rifting.

The East African Rift System (EARS) preserves within it a broad range of rifting morphologies from incipient rifting in the south to the transition from continental to oceanic crust in the north (Corti, 2009; Ebinger, 2005; Ebinger and Casey, 2001; Gregory et al., 1896; Gregory, 1920; Sueß, 1891; Woldegabriel et al., 1990). The existing models of strain accommodation within continental rifting settings were borne of studies of the EARS. As studies of this system progressed, the initial models of rift-centered strain accommodation have evolved into models where the extensional strain is accommodated both by the rift-border faults and the magmatic belts within them (Ebinger, 2005; Wolfenden et al., 2005). The amount of extensional strain that is accommodated by the rift-border faults and magmatic belts, and the timing of this accommodation, remain controversial (e.g., Agostini et al., 2011a; Bendick et al., 2006; Bilham et al., 1999; Casey et al., 2006; Molin and Corti, 2015; Pizzi et al., 2006). The Main Ethiopian Rift (MER), located directly south of Afar, is a transitional region whereby strain is accommodated both by rift-border faults and focused magmatic intrusion (Casey et al., 2006; Corti, 2009; Ebinger and Casey, 2001; Pizzi et al., 2006) (Fig. 1A). Within the Main Ethiopian Rift, recent magmatic activity has been focused along linear belts known as the Silti-Debre Zeyit Fault Zone (SDFZ) and the Wonji Fault Belt (WFB) (Acocella et al., 2003; Bonini et al., 2005; Casey et al., 2006; Corti, 2009; Ebinger and Casey, 2001; Kurz et al., 2007; Mohr, 1962; Morton et al., 1979; Rooney et al., 2007; Rooney et al., 2014; Rooney et al., 2005; Woldegabriel et al., 1990; Wolfenden et al., 2004) (Fig. 1B). While much attention has been given to the magmatic activity along these linear belts within the rift proper (Medynski et al., 2015; Rooney et al., 2007; Rooney, 2010; Rooney et al., 2011; Rooney et al., 2014), there is a growing realization that magmatic activity associated with rifting is more diffuse (e.g., Abebe et al., 1998; Meshesha and Shinjo, 2007; Rooney et al., 2016; Trestrail et al., 2016).

Here we undertake a petrographic and geochemical analysis of an area of focused magmatic intrusion on the Southeastern Ethiopian Plateau, which is interpreted to be contemporaneous with rift development but occurs outside of the Ethiopian Rift valley (Kennan et al., 1990; Mohr and Potter, 1976). The Galema range is a NNE-trending series of en echelon dikes, lavas, and eroded cinder cones, located ∼35 km east of the Asela-Sire Border Fault (Fig. 1C). We focus upon constraining the magmatic plumbing system of the Galema range through the use of geochemical and petrographic analyses, offering insight into the architecture of the magmatic plumbing system in which the magma ascended and the depths of crystal phase fractionation. The depths at which a magma commonly stalls on ascent are, in part, determined by the extensional stresses acting on the lithosphere at the time, and these volumes at depth represent presumed zones of magmatic strain accommodation (Corti et al., 2003). In this manner, a profile of the lithospheric strain at an off-rift location may be constructed. Our results demonstrate that the Galema range accommodated extensional strain after the initial formation of the border faults that define the rift, but before the inception of zones of focused magmatic intrusion within the rift (i.e., WFB and SDFZ). This suggests that strain accommodation by focused magmatic intrusion can occur over a broad area in continental extension settings. Furthermore, the emplacement of the Galema range also shows that migration of focused magmatic activity toward the rift axis is not always a straightforward progression as some models of continental rifting evolution assert.

Cenozoic magmatic activity in East Africa commenced in the Eocene with basaltic volcanism in the broadly rifted zone in the Turkana area of northern Kenya and southern Ethiopia at ca. 45 Ma (Corti, 2009; Ebinger et al., 1993; George and Rogers, 2002) (Fig. 1A). Large-scale flood basalts were emplaced during the Oligocene as a result of the impingement of the Afar plume (Baker et al., 1996; Hofmann et al., 1997; Kieffer, 2004; Pik et al., 1999; Rochette et al., 1998). In the early Miocene, a resurgence of volcanism occurred throughout the region, with increased basaltic activity in Turkana and both shield volcanism and fissure-fed basalts on the Ethiopian Plateau (Rooney, 2017). Subsequent magmatism was largely localized to the developing MER and took the form of large silicic centers and smaller basaltic eruptions (Rooney et al., 2011).

The MER runs from Afar in northeastern Ethiopia to Turkana in Kenya and changes from a ∼NE-SW trend in the north to a ∼N-S trend in the south. The MER is divided into three sectors: the Northern Main Ethiopian Rift (NMER), the Central Main Ethiopian Rift (CMER), and the Southern Main Ethiopian Rift (SMER) (Corti, 2009; Woldegabriel et al., 1990) (Fig. 1A). These sectors evolved at different times and into a nonlinear trend (Abebe et al., 2010; Balestrieri et al., 2016; Bonini et al., 2005; Corti, 2009; Keranen and Klemperer, 2008; Woldegabriel et al., 1990; Wolfenden et al., 2004).

The CMER is bounded by the Boru-Toru structural high and the Goba-Bonga structural lineament (Bonini et al., 2005) (Fig.1A). The main border fault in this area of the CMER, known as the Asela-Sire Border Fault, is a segmented system of high-angle, normal faults (>60°) that define the MER. These faults trend at approximately N30°E–N40°E, separating the CMER from the Southeastern Ethiopian Plateau (Abebe et al., 2010; Boccaletti et al., 1998; Woldegabriel et al., 1990) and are mid-Miocene (Ebinger and Casey, 2001), no older than ca. 8 Ma (Abebe et al., 2010; Bonini et al., 2005; Woldegabriel et al., 1990).

The Asela-Sire Border Faults formed during the southward propagation of rifting from Afar (Abebe et al., 2007; Bonini et al., 2005; Keranen and Klemperer, 2008). The N-S trend of the Asela-Sire Border Faults led to the interpretation that the orientation of extensional stress had rotated to roughly E-W from the previous NW-SE direction (Boccaletti et al., 1998; Bonini et al., 2005). It should be noted that this difference in trend has been attributed, not to rotation of the stress field, but to the interference of preexisting lithospheric weaknesses interacting with an extensional stress field oblique to the trend of the rift (Agostini et al., 2011b; Corti, 2008; Corti et al., 2013). Postdating the initiation of the Asela-Sire Border Faults, the magmatic activity within the CMER was focused on the Quaternary Wonji Fault Belt (WFB) and the Silti-Debre Zeyit Fault Zone (SDFZ) (Mohr, 1962; Rooney et al., 2014; Woldegabriel et al., 1990) (Fig. 1B). While the WFB accommodates some of the extensional strain in the CMER, the Asela Sire rift-border faults also currently continue to accommodate a portion of this strain (Agostini et al., 2011a; Corti, 2009; Pizzi et al., 2006). However, the amount of extensional strain accommodated by the rift-border faults and the timing of this strain accommodation by both mechanisms remain controversial.

Analysis of the river network patterns in the CMER indicates that a portion of the extensional strain is accommodated by the Asela-Sire Border Faults (Molin and Corti, 2015). GPS velocities show that the recent extension rate across the MER is ∼4.0 ± 0.9 mm/yr, accommodated by both mechanical stretching and faulting of the lithosphere as well as aseismic dike injection (Bendick et al., 2006).

Structural analysis of the CMER indicates that the Aslea-Sire border faults have recently accommodated extensional strain—as recorded by fault-slip data recorded across the rift (Agostini et al., 2011a). This analysis also indicates that extensional strain accommodation of the WFB, via faulting of the crust, is subordinate to that accommodated at the eastern rift-border fault. However, forces exerted by the magma overpressure of an intruding dike can reduce the amount of observable slip along nearby faults (Rubin and Pollard, 1988). The compressive force exerted by the intruding dike in the direction of σ₃, as a result of the buoyancy of an ascending magma, can change the local stress field in the host rock by increasing compression to each side of the dike and reducing compression above (Rubin and Pollard, 1988). One consequence of this is that preexisting faults that are located adjacent to the intruding dike (i.e., those not generated by intrusion of the dike) can become locked (Rubin and Pollard, 1988). Exploitation of preexisting faults by magmatic intrusion on ascent adds volume to the host rock, accommodating extensional strain (Rubin and Pollard, 1988). In areas experiencing both faulting and intrusion (such as the WFB), the total amount of deformation is partitioned between fault-slip and volumetric increase, but the observable amount of fault-slip appears to be reduced, complicating estimates of the amount of strain accommodation by faulting.

The Wonji Fault Belt (WFB), located near the eastern rift margin and subparallel to the rift axis in the CMER, is described as a series of right-stepping, en echelon faults and dikes within the rift, which are interpreted to be <1.8 Ma (Ebinger and Casey, 2001; Woldegabriel et al., 1990; Wolfenden et al., 2004). The volcanism of the WFB is described as being bimodal in nature with silicic volcanism, sometimes forming large calderas, and basaltic volcanism forming lava flows and scoria cones (Kurz et al., 2007; Mohr, 1967; Mohr et al., 1980; Rooney et al., 2005).

Located toward the western rift margin, the SDFZ is characterized as a linear chain of nested calderas, scoria cones, and basaltic flows with associated silicic centers also oriented subparallel to the rift axis (Gasparon et al., 1993; Rooney et al., 2007; Rooney et al., 2011; Rooney et al., 2005; Woldegabriel et al., 1990).

The locations of both the WFB and SDFZ along the rift margins may be due to the existence of preexisting lithospheric weaknesses within the MER. When a low-obliquity extensional stress field is applied to these weaknesses, the result is formation of boundary faults oblique to the extension direction, followed by formation of internal faults nearly perpendicular to the direction of extension (Agostini et al., 2009; Agostini et al., 2011b; Corti, 2008). Analog modeling of stress direction oblique to preexisting lithospheric weakness suggests that the WFB in the CMER represents an intermediate stage of rift development (Agostini et al., 2011b).

We focus on the Galema range, an area of Pliocene, rift-adjacent magmatism on the Southeastern Ethiopian Plateau (Fig. 1C). Located ∼130 km southeast of Addis Ababa, it occurs to the east of the border faults that define the rift valley in the CMER (Fig. 1C). Morphologically, the Galema range displays striking similarities to magmatic features located within the rift valley of the MER. For example, the Galema range consists of a series of right stepping, en echelon mafic dikes trending at ∼10°–20°N and aligned, eroded scoria cones (Mohr, 1980; Mohr and Potter, 1976; Woldegabriel et al., 1990), similar to the WFB. The Galema range extends from the volcanic center in the north known as Gara Badda to the volcanic center in the south known as Gara Enkulo (Fig. 1C). The dikes, which constitute the range, span a distance of ∼70 km and are arranged in a dextral, en echelon formation within a ∼7-km-wide fissure system (Mohr and Potter, 1976). Maars and cinder cones exist between the two volcanic centers, along with the ∼4.5-m-thick dikes that gradually disappear at the latitude of Gara Enkulo (Mohr and Potter, 1976). The Galema range rises to an elevation ∼1000 m above the Southeastern Ethiopian Plateau and is comprised of the magmatic pile erupted from the dikes and silicic centers (Mohr and Potter, 1976). Evidence of this volcanism can be found in the volcanic edifices that appear within the Galema range. Volcanic edifices within the WFB are dominantly of two varieties: cinder cones and large silicic centers. The Galema range preserves within it volcanic necks of silicic composition, which are the exposed magmatic plumbing system of large silicic centers (Fig. 2). Moreover, the scoria cones found within the WFB are rooted by dikes; the exposed dikes evident in the Galema range are parallel to the magmatic plumbing system of these scoria cones. Field work has revealed eroded scoria cones in the southern portion of the study area (Mohr and Potter, 1976).

While the magmatic features of the Galema range and the WFB display striking morphological similarities, the direct evidence of the existence of faulting of the Galema range is less apparent than that of the WFB. This may indicate that while both areas have accommodated extensional strain, the amount of strain accommodated by faulting in the Galema range may be less profound. Indirect evidence, however, suggests that faults existed prior to the onset of magmatic intrusion. Recent work has highlighted the influence of preexisting faults and fractures within the lithosphere on the morphology of dike swarms and volcanic fields (Le Corvec et al., 2013; Mazzarini et al., 2016). These studies indicate that in areas of relatively thicker crust, alignments of dikes and monogenetic vents (i.e., scoria cones) are parallel to the trend of preexisting lithospheric fractures and faults (Mazzarini et al., 2016). Thus, the alignment of dikes and scoria cones of the Galema range, in the area of thicker crust (see section 5.5) of the Eastern Ethiopian Plateau, suggest that the ascending mafic dikes exploited preexisting lithospheric faults and fractures. As the ascending dikes are captured by the preexisting faults (Le Corvec et al., 2013), evidence of the existence of these faults may be obscured as they are cut by the intruding magma (Rubin and Pollard, 1988).

Glaciation during the middle to Late Pleistocene eroded the up to 2000-m-thick magma pile and exposed the dikes that had fed the cinder cones of the Galema range (Mohr and Potter, 1976). Taken together, it is evident that the magmatic plumbing systems of both silicic and mafic edifices are exposed within the Galema range. It is important to note that such erosion would have removed all evidence of faulting induced above propagating dikes, as mentioned earlier. The implication of such erosion would be the preservation of a magmatic plumbing system, broadly similar to that beneath the modern WFB, but devoid of clear evidence of faulting.

K-Ar dating of samples of the Galema range yielded dates of formation at ca. 2.1 to ca. 3.1 Ma (Mohr and Potter, 1976) and 1.97 ± 0.02 Ma (Kennan et al., 1990). The existence of off-rift zones of focused intrusion such as the Galema range is not well constrained within the current rifting framework but has been attributed to the possible existence of preexisting lithospheric structures (Corti, 2012; Mohr and Potter, 1976).

A field excursion in 2008 recovered 77 samples from dikes, basaltic flows, scoria cones, and silicic centers of the Galema range. These samples were cut to minimize alteration due to weathering before being further cut into ∼30 g billets. After polishing to remove saw marks and cleaning in an ultrasonic bath of deionized water, the billets were crushed in a steel jaw crusher and then powdered in a Bico ceramic disk mill. The powders were then fused with a lithium tetraborate (Li₂B₄O₇) flux into glass discs, following procedures detailed elsewhere (e.g., Rooney et al., 2012). Major-element oxide abundances and loss-on-ignition values were determined via X-ray fluorescence analysis on a Bruker S4 Pioneer instrument at Michigan State University (Supplemental Tables¹).

Analyses of clinopyroxene and plagioclase mineral phases were performed on a representative subset of samples ranging from primitive to evolved compositions, on the basis of their bulk-rock MgO content, using in situ spot analyses by laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) following protocols developed together with other studies (Rooney et al., 2016; Trestrail et al., 2016). Standard petrographic thin sections were loaded into a 15 × 15 cm two-volume HelEx sample cell, integral to a Teledyne Photon Machines Analyte G2 193 nm excimer laser ablation system at Michigan State University. Multiple spots of 110 µm diameter in individual crystals were ablated for 30 s with a laser repetition rate of 10 Hz at a fluence of 4.1 J/cm². A gas blank was collected after each sample for ∼20 s, and gas blank intensities were subtracted from ablation intensities. The ablated material was transferred to a Thermo Scientific ICAP Q quadrupole ICP-MS in a flux of 0.75 L/min of high-purity He carrier gas for analysis. The ICP-MS was tuned for highest intensities, low-oxide production rate (ThO/Th < 0.7%), and low double-charged cation formation rate (Ba²⁺/Ba <1.5%) by performing surface scans of National Institute of Standards and Technology (NIST) SRM 612 glass standard, while signal stability was subsequently tested over a 3 min surface scan of NIST SRM 612. Data were calibrated using 20 natural and synthetic geological standards (Supplemental Tables [^{footnote 1}]) mounted on the same sample tray together with the samples. These were analyzed sequentially without opening the sample chamber and interrupting automation. Data were collected using the ICP-MS in kinetic energy discrimination mode in which the ionized material circulates through a collision cell fluxed with He, reducing isobaric interferences with oxide and argide compounds, enhancing detection limits of light mass isotopes including major elements, and lowering their signal intensities, allowing for their quantification at a precision and accuracy similar to trace elements (Supplemental Tables). Data were drift corrected using surface scans of well-characterized fused rock powder standards (JB-1a and BHVO-1) analyzed approximately every hour during analysis at 10 μm/s and other ablation parameters identical to those described above. Replicated analyses of JB-1a and BHVO-1 analyzed as unknowns by surface-scan analyses achieved ±5% accuracy and day-averaged reproducibility of 5% (1σ) or better (Supplemental Tables). The raw sum of all major-element oxides with Fe expressed as FeO was calculated for each mineral analysis, in order to normalize all major- and trace-element compositions to 100 wt% oxides. Practically, 100 wt% oxides was our internal standard for all mineral analyses; such an assumption is reasonable for all stoichiometrically anhydrous minerals. If necessary, certain laser ablation signals were then filtered to remove certain sections of ablation sequences exhibiting evidence of micro-crystal inclusions and or alteration. Reproducibility of the in situ spot analyses was assessed separately by performing ten 30-second ablations on each of the basalt glass standards GSD-1G and BCR-2G for major and trace elements using analytical conditions identical to those above; on these homogeneous glass materials, reproducibility is also generally of 5% (1σ) or better (Supplemental Tables).

In hand sample, the mafic specimens collected from the Galema range can be described as massive, non-vesicular basalts with sparse, 1–5 mm–scale plagioclase crystals. A subset of samples is scoriaceous with visible plagioclase crystals at the 1–5 mm scale. As a whole, the samples show little to no evidence of alteration or magma mingling at hand-sample scale.

In thin section, a majority of the samples collected from the Galema range exhibit crystal phases dominated by plagioclase, and a lesser abundance of samples display olivine as the primary crystal phase. A subset of samples includes clinopyroxene crystals as part of the phenocryst population (Fig. 3). These crystal phases rarely display resorption or skeletal textures, and most appear to be in equilibrium with the ground mass, while none display evidence of zoning (Fig. 4) (Supplemental Tables [^{footnote 1}]). All samples are observed to contain plagioclase laths within a cryptocrystalline to aphanitic groundmass. Opaque Fe-Ti oxide crystals also occur as a fractionating phase in all samples in plane polarized light. Anhedral olivine crystals occur as part of the mineral assemblage in varying abundances in samples with >∼2.5 wt% MgO, appearing as the primary mineral phase in samples with >∼9.0 wt% MgO, though often altered to iddingsite (Figs. 3 and 4). Anhedral to subhedral clinopyroxene crystals occur as part of the phenocryst assemblage in samples with between ∼7.0 wt% and ∼8.0 wt% MgO (Figs. 3 and 4). Plagioclase crystals are present within the range of ∼1.3 to ∼9.6 MgO wt%. It is the primary mineral phase in samples from ∼1.0 to ∼4.5 MgO wt% and occurs in varying abundance with olivine between ∼4.5–9.6 wt% MgO (Fig. 3).Within the Galema range, the temporal relationship between the silicic centers and mafic dikes is not well constrained, and we therefore focus on the mafic samples in an effort to characterize their magmatic plumbing systems. Samples with >∼59 wt% SiO₂ (<∼1.4 wt% MgO), as determined by X-ray fluorescence analysis, were not considered further.

On the basis of the total alkali silica diagram (Le Maitre et al., 2005), our samples from the Galema range plot in an array from picrobasalt to rhyolite (Fig. 5). The magmatic products of the Galema range have been described previously as bimodal in silica content, composed of hawaiites (trachybasalts) and comendites (rhyolites) (Mohr, 1980). There is a noticeable lack of samples with SiO₂ content between ∼59 wt% to ∼68 wt% (Supplemental Tables [^{footnote 1}]) (Fig. 5). This lack of such sample compositions, also known as the Daly Gap, is described in many locations in the EARS (Baker, 1987; Peccerillo et al., 2003; Rooney et al., 2012; Trua et al., 1999). This gap has been interpreted as a function of differential tapping of magma chambers at depth (Peccerillo et al., 2007) and more recently as the result of a magma mush where more silicic material has evolved and removed from the magma chamber to fractionate and shallower levels (Rooney et al., 2012; Trua et al., 1999).

Variation diagrams of the major-element geochemistry of samples with >∼1.4 wt% MgO show trends consistent with fractional crystallization of olivine, plagioclase, pyroxene, and Fe-Ti oxides occurring over the course of magmatic evolution (Fig. 6). In samples with >∼9 wt% MgO content, the relatively flat trend of Fe₂O₃ and TiO₂ versus MgO is consistent with crystallization of olivine. Crystallization of pyroxene occurs in the compositional range of ∼6 to ∼9 wt% MgO as indicated by the decreasing trend of CaO/Al₂O₃ versus MgO. Below ∼5 wt% MgO, formation of Fe-Ti oxides is indicated by a break in slope of Fe₂O₃ and TiO₂ versus MgO (Fig. 6). Below ∼6 wt% MgO, plagioclase is the primary fractionating phase as indicated by the slope of the CaO/Al₂O₃ versus MgO variation diagram.

Results of in situ LA-ICP-MS analysis on the crystal phases of samples selected for crystal availability from each of the pressure regimes determined in section 5.2.1 (3028, 3029, 3061, 3068, and 3072) reveal compositions confirming the observations in thin section of the presence of plagioclase and pyroxene phases (Supplemental Tables [^{footnote 1}]). Multiple analyses of each crystal were performed, and data were averaged across each crystal because no significant core-to-rim zoning was observed (Supplemental Tables) (Fig. 4E). The size of the laser spot created an averaging effect of the crystal chemistry in each analysis, obscuring any chemical evidence of fine zoning for each crystal. After filtering the data that showed clear signs of inclusion and/or alteration, the crystal chemistry of each spot in an individual transect was analyzed for evidence of chemical zoning, and the standard deviation for each transect was calculated (Supplemental Tables) (Fig. 4E). These data were averaged to determine the average crystal chemistry. The averaging calculations allow for a more representative composition of the clinopyroxene phenocrysts in equilibrium with the melt, facilitating thermobarometric calculations.

Labradorite occurs as the most abundant plagioclase phase in samples 3029 (3.28 wt% MgO) (An_52–54), 3061 (7.94 wt% MgO) (An_63–65), and 3072 (4.86 wt% MgO) (An_51–56). Bytownite occurs as the most abundant plagioclase phase of sample 3068 (7.02 wt% MgO) (An_71–78) (Supplemental Tables [^{footnote 1}]) (Fig. 7A). The clinopyroxene phases of samples 3061 and 3028 (7.15 wt% MgO) are classified as augite with Wo_38–40, Fs_9–15 and Wo_38–42, Fs_11–16, respectively (Supplemental Tables) (Fig. 7B).

To determine whether the crystal phases are in equilibrium with the ground mass, and to choose suitable candidates for thermobarometric calculations, the in situ analyses were used in combination with the whole-rock X-ray fluorescence (XRF) analyses to calculate partition coefficient (K_D) values between crystal with melt. In sample 3061, K_D (Fe-Mg)^cpx-liq = 0.27 ± 0.02, while in sample 3028, K_D = 0.29 ± 0.01. These values are well within the K_D (Fe-Mg)^cpx-liq values of 0.29 ± 0.08 proposed by Putirka (2008), suggesting that the crystals are in equilibrium with melt.

Crystal assemblages from the analyzed samples reveal no discernable pattern with regard to sample location on the Galema range. Small groupings of samples with similar crystal assemblages are present; however, these small groups lie in close proximity (<5 km) to samples with a dissimilar assemblage. It is unclear if this is an artifact of sample availability or perhaps magmatic processes at such scale.

To first-order, dikes are emplaced perpendicular to the orientation of the least compressive stress (σ₃) (Rubin and Pollard, 1988; Rubin, 1990; Parsons and Thompson, 1991; Rubin, 1992). In rifting environments, the orientation of the least compressive stress is in the direction of the extensional strain imposed on the lithosphere. If this intruded magma is sourced from a location other than the host rock (i.e., not a result of partial melting of the crust) and is allowed to remain in the host rock through magma cooling and dike freezing, the result is an increase of the volume of the host rock (Parsons and Thompson, 1991). In this way, magmatic intrusion accommodates extensional strain. The rapid emplacement of dikes can accommodate up to millennia worth of accumulated strain within the host rock (Parsons and Thompson, 1991). The depth where magma stalls on ascent and forms magma chambers coincides with depths where lithospheric strain accommodation occurs. Locations where magma is emplaced perpendicular to the direction of least principal compressive stress (σ₃), such as areas experiencing lithospheric extension, indicate areas where the crust may accommodate extensional strain aseismically (Corti et al., 2003). Thus, an estimate of the architecture of a magmatic plumbing system and the depths at which magma accumulates can constrain where and how strain is accommodated in the lithosphere in areas of extensional tectonics.

The stalling of a magma on ascent is testable through the study of its major-element geochemistry (e.g., Peccerillo et al., 2003; Rooney, 2010; Trua et al., 1999). The depth at which a magma stalls in the lithosphere, in part, controls the phases that form from that magma, thereby influencing the composition of the phase and residual magma from which the phase fractionates (Green and Ringwood, 1967; Green, 1969; O’Hara, 1968; Thompson, 1972; Wilson, 1994). Previous studies within the region have shown that deeper dike stalling depths resulted in greater stability of clinopyroxene, in contrast to shallower conditions where olivine and plagioclase dominate the fractionating assemblages, illustrating the use of major-element geochemical analysis in the interpretation of the architecture of a magmatic plumbing system (Peccerillo et al., 2003; Rooney et al., 2007; Rooney et al., 2005; Trua et al., 1999).

To constrain the depths of magma stalling of the Galema range, and thereby establish where extensional strain is being accommodated in the lithosphere, thermodynamic modeling of the crystal fractionation sequence was performed using Excel-MELTS (Asimow and Ghiorso, 1998; Ghiorso and Sack, 1995; Gualda and Ghiorso, 2015). An initial composition for modeling was chosen as the most primitive sample from our data suite (sample 3014, MgO wt% = 11.00) (Fig. 6). Multiple models were run utilizing this initial composition with varying pressures (10 kbar–2 kbar), water contents (0%–5%), and ƒO₂ (QFM ± 1). Utilizing a constant pressure of 7 kbar, 1 wt% H₂O, and ƒO₂ value of QFM, with temperature decreasing continuously from 1350 °C, the modeled evolution of this initial composition was found to fit the observed data above ∼6.1 wt% MgO in major-element variation diagrams (Fig. 6). At 7 kbar, the MELTS model predicts minor olivine to begin fractionating at 1307 °C (10.98 wt% MgO). The dominant fractionating crystal phase of the model at this pressure is clinopyroxene, which begins fractionating at 1250 °C (9.01 wt% MgO) and continues to where the model deviates from the data at ∼6.1 wt% MgO (Fig. 6). Below ∼6.1 wt% MgO, the model with these initial conditions deviates significantly from the observed data when plotted in major-element variation diagrams. To resolve this deviation from the observed data, a second model was created, utilizing a starting composition of ∼6 wt% MgO (i.e., sample 3016, MgO wt% = 6.20) (Fig. 6). At a constant pressure of 3 kbar, 0.5 wt% H₂O, ƒO₂ at QFM, and continuous decreasing temperature from 1200 °C, the results of the model fit the observed data with a compositional range of ∼1.3 to ∼6.1 wt% MgO (Fig. 6). The 3 kbar model predicts that olivine will begin fractionating at 1181 °C (6.21 wt% MgO). At 1158 °C (5.62 wt% MgO), the MELTS model predicts the dominant fractionating crystal phase to be plagioclase, which continues fractionating through the end of the model. Minor amounts of spinel will begin fractionating at 1133 °C (4.53 wt% MgO). Minor clinopyroxene does fractionate in the model at this pressure, but plagioclase is clearly fractionating in greater proportion.

The results of the MELTS model agree with the variation diagrams and petrography of the Galema basalts. The 7 kbar model predicts the primary crystalizing phase to be olivine at 10.98 wt% MgO, while the variation diagrams show that olivine is the primary crystalizing phase above ∼9 wt% MgO (Fig. 6). This is in agreement with the petrography of samples in that compositional range (Fig. 3). The MELTS model predicts that clinopyroxene should fractionate at 9.01 wt% MgO at 7 kbar, in agreement with the slope of the CaO/Al₂O₃ versus MgO variation diagram in that compositional range (Fig. 6). A subset of samples from compositional range of ∼7 to ∼8 wt% MgO also displays clinopyroxene as a member of the phenocryst population (Fig. 3). The 3 kbar MELTS model predicts the primary fractionating crystal phase to be plagioclase at ∼5.62 wt% MgO, which is in agreement with the CaO/Al₂O₃ versus MgO variation diagram in that compositional range (Fig. 6). The petrography of samples from this compositional range also shows plagioclase as the dominant crystal phase (Fig. 3). The 3 kbar MELTS model predicts spinel fractionating at 4.53 wt% MgO, consistent with the Fe₂O₃ versus MgO variation diagram (Fig. 6) and the occurrence of oxide phases in the petrography (Fig. 3).

Spinel and plagioclase exist in the petrography of the basalts of the Galema range for compositions above ∼5.62 wt% MgO, which is not accounted for in the deeper MELTS model (Fig. 6). The presence of these phases may be the result of the shallower fractionation, subsequent to the initial stages where clinopyroxene was dominant at greater depths.

The results of these thermodynamic models are thus consistent with the existence of a polybaric fractionation system, which is centered at 25 and 11km depth (Fig. 8). The modeled pressures of fractionation, when correlated with the suite of samples analyzed, show no discernable correlation with their locations on the Galema range (Fig. 1C).

To ascertain the pressures of fractional crystallization, and to validate the Excel-MELTS model of magmatic evolution, the in situ geochemical analysis of the crystal phases of the Galema range was applied to clinopyroxene-melt thermobarometry calculations (Nimis, 1995). The samples were chosen based upon the MgO compositional range and the availability of analytically suitable clinopyroxene crystals. Liquid compositions for the selected samples were estimated using the whole-rock geochemical data for sample 3028 (60% groundmass) and 3061 (25% groundmass) following the methodology of Bédard (2014), Putirka (2008), and Putirka et al. (2003). The water content (0.1–1 wt%) for the magmas in the thermobarometry calculations were estimated based on the H₂O compositions determined in the MELTS modeling. The thermobarometric calculations determined that samples 3061 and 3028 fractionated at pressures of 6.7–8.6 kbar (24.4 km to 31.3 km depth) (Fig. 9). These calculations are within ±2.0 kbar of the pressures determined in the MELTS models (Fig. 6).

While plagioclase mineral data were collected (Supplemental Tables [^{footnote 1}]), these data were not used in thermobarometric calculations. Plagioclase-liquid thermobarometers have large errors associated with their pressure estimates, with some as high as ±3 kbar (Putirka, 2008). These estimates were not useful for our analysis of the magmatic plumbing system of the Galema range.

Geophysical observations of the lithosphere of the MER show regions of higher P-wave velocities (V_p) in the crust at 10 km beneath the surface expressions of the WFB, SDFZ, and the Galema range (Keranen et al., 2004; Rooney et al., 2014) (Fig. 10A). This higher V_p indicates the existence of higher density material, which is consistent with a considerable volume of solidified, mafic material (interpreted to be gabbroic in composition) that was intruded into the lithosphere (e.g, Keranen et al., 2004; Mahatsente et al., 1999). The differences in magnitude of V_p between the Galema range, WFB, and the SDFZ may be due to differences in the volume of higher density material stored in the lithosphere between each location. This depth at which higher V_p occurs beneath the Galema range is consistent with the results of our MELTS modeling and thermobarometric calculations. This correlation may be evidence of cooled mafic material at these depths, resulting from stalled magmatic intrusion similar to processes observed in the WFB and SDFZ (Keranen et al., 2004; Rooney et al., 2005). Residual gravity measurements are entirely consistent with the seismic velocity anomalies. Specifically, the gravity data show a distinct negative residual gravity across the entire rift. Outside the rift, there is typically no significant anomaly, with the exception of the Galema range on the eastern plateau. The anomaly beneath the Galema range has been previously interpreted to be the result of silicic flows (Mahatsente et al., 1999); however, this interpretation was not based upon field data. The magnitude of this gravity anomaly is related to the amount and density of material within the crust and does not offer a unique solution to the observations.

Recent activity in Afar indicates that extensional strain is accommodated by focused magmatic intrusion within polybaric magmatic plumbing systems. Bimodal chemistry of magmatic products within the Dabbahu magmatic segment in Afar has been attributed to polybaric fractionation resulting from multiple depths of magma storage, which is also in agreement with geophysical studies (Ayele et al., 2009; Belachew et al., 2011; Field et al., 2012; Hammond, 2014; Wright et al., 2012). Zones of focused magmatic intrusion are also found in the MER in the WFB and SDFZ. Geophysical measurements show a volume of dense, intruded material within the lithosphere at the SDFZ (Keranen et al., 2004). Geochemical analyses of the magmatic products of the SDFZ indicate that the magmas were derived from sublithospheric depths (Rooney et al., 2005). The intruded material was introduced into the crust adding volume to the lithosphere at the SDFZ, indicating that the magmatism of the SDFZ represents accommodation of strain through focused magmatic intrusion.

Similarly, the WFB in the CMER plausibly represents an area of extensional strain accommodation by focused magmatic intrusion. Intruded material beneath the WFB, which added volume to the lithosphere, has been observed in geophysical measurements (e.g., Bastow et al., 2008; Biggs et al., 2011; Keir et al., 2011; Keir et al., 2015; Keranen et al., 2004; Kim et al., 2012; Rooney et al., 2007). Seismic anisotropy indicates the structure of these intrusions to be vertical dikes aligned in the orientation perpendicular to the extensional stress field affecting the MER (Keir et al., 2005; Kendall et al., 2006). These dikes are not derived from partial melting of the lithosphere, again requiring added volume to the crust (Rooney et al., 2007; Rooney et al., 2011). The WFB thus represents an area of magmatic strain accommodation.

As noted previously, the Galema range displays a similar morphology to the Quaternary zones of focused intrusion within the MER. The morphological similarities between the Galema range and areas of focused magmatism within the rift (e.g., dikes aligned perpendicular to the σ₃ induced by the extensional stress field) suggest that the Galema range also magmatically accommodated strain, but intriguingly, exists outside the rift proper.

The polybaric fractionation system of the Galema range is also consistent with the architecture of the magma plumbing system of the WFB. Thermodynamic modeling, petrographic modeling, and thermobarometric calculations have shown that the WFB consists of a polybaric fractionation system at deep crustal levels before magma ascends to consistently shallow levels within the crust (Rooney et al., 2005; Rooney et al., 2007; Rooney et al., 2011; Trua et al., 1999). The magmas of the WFB have been interpreted, through clinopyroxene thermobarometry, as having crystallization pressures of 0.1 kbar to 5.0 kbar (Rooney et al., 2007)—a narrower and shallower range of fractionation depths than the magmas of the Galema range.

The thermobarometric calculations of the Galema range indicate greater pressures of magma stalling and fractionation in comparison to those found at the WFB. These pressures are more consistent with those determined from analyses of the SDFZ. The SDFZ displays geochemical evidence of increased fractionation of clinopyroxene over plagioclase in its magmas, suggesting deeper levels of fractionation than at the WFB (e.g., Green and Ringwood, 1967; Rooney et al., 2007; Rooney et al., 2005). Similarly, thermobarometric estimates derived from host lavas within the SDFZ indicate fractionation pressures over a range from 2.0 kbar to 5.6 kbar (Rooney et al., 2005). This increased depth of magma stalling for the less mature magmatic plumbing system of the SDFZ is more similar to the Galema range. These observations indicate a greater similarity in magmatic processes between the Galema range and the SDFZ in comparison to the WFB.

Thermal modification of the base of the lithosphere allows for magma to intrude and begin fluxing into the affected areas (Havlin et al., 2013). This thermal modification allows for subsequent dikes to intrude further before freezing, broadly modifying the lithosphere in such a way as to allow diking and magma ascent to proceed (Bastow and Keir, 2011; Bialas et al., 2010; Daniels et al., 2014; Havlin et al., 2013). Once in the lithosphere, the depth at which a magma may stall is controlled by how the buoyancy forces of the magma interact with density contrasts of the lithosphere (Menand, 2011).

Buoyancy forces affecting a dike within the lithosphere depend on the availability and amount of magma generated (Parsons and Thompson, 1991). Buoyant magma will ascend to levels within the lithosphere where density differences are equal to the magma overpressure, which is generated by the buoyancy forces (Rubin, 1990). However, in a persistent rift zone, the constant extensional stress in the lithosphere is large enough to offset the stress state that results from lithospheric density differences (Rubin, 1990). If the extensional stress field increases during dike ascension, but magma flux remains constant, the magma overpressure is essentially cancelled, and the dike is trapped at depth. Extension reduces the horizontal stress perpendicular to the rift zone (σ₃), allowing dikes to increase in width, but decreases their chance for ascent when the extensional stress field returns to normal (Rubin, 1990). If magma flux and resultant magma overpressure increased after such an event, the dike would again be allowed to ascend, but in a wider, more well developed conduit.

Magma flux into the lithosphere determines the development of a rift magmatic plumbing system. Greater available magma volume is necessary to form larger dikes, which may traverse the crust to shallow levels (Rooney et al., 2014). A reduced magma flux into the lithosphere would create a more complex magmatic plumbing system because smaller dikes would have difficulty ascending directly to the surface due to reduced thermal weakening of the lithosphere (Bialas et al., 2010; Buck, 2004; Havlin et al., 2013). The less mature magmatic plumbing system of the Galema range, in comparison to the WFB, may be due to an overall decreased magma flux into the lithosphere resulting from both the average extent and rate of lithospheric stretching and thinning. This reduced magma flux would also inhibit the ascension of any dikes trapped at depth due to a temporary increase in the local stress field.

As the lithosphere thins, the underlying asthenosphere must ascend to compensate partially for this reduced thickness, generating an increasing melt (Langmuir and Forsyth, 2007). Increased extension rate also increases the amount of adiabatic melting of the asthenosphere and can therefore increase magma flux into the lithosphere (Langmuir and Forsyth, 2007).

With progressive lithospheric thinning, the depth of melting of the asthenosphere shallows as the lithospheric cap is reduced, and adiabatic melting can occur at decreased depth (Rooney, 2010; Wang et al., 2002). Stretching of the lithosphere manifests as an overall thinning of the plate, resulting in a shallowing of both lithosphere-asthenosphere boundary (LAB) and the Moho. Active source seismic studies within the MER show a Moho depth of ∼32–36 km, but beneath the relatively less stretched Southeastern Ethiopian Plateau, upon which the Galema range lies, the Moho is at a depth of ∼42–46 km (Keranen et al., 2009; Rooney et al., 2014) (Fig. 10B). Given this increased Moho depth, we anticipate a deeper melting column beneath the Galema range (which is controlled by the depth of the LAB).

An examination of the melting depth of the magmas of the Galema range reveals depths greater than those at the WFB. Variation of SiO₂ versus FeO* for Galema range lavas (corrected for olivine fractionation using the PRIMELT3 MEGA modeling suite; Herzberg and Asimow, 2015) exhibits a distinct grouping displaying higher FeO* content than experimentally derived melts of fertile mantle peridotite at a pressure of ∼25–30 kbar (Baker and Stolper, 1994; Furman et al., 2006; Kushiro, 1996) (Fig. 11). The grouping of the Galema range lavas at high FeO* also exceeds lavas from the adjacent MER, e.g., the WFB and Akaki (Rooney et al., 2007; Rooney et al., 2014) and suggests a source melt for the Galema basalts at pressures greater than ∼30 kbar (∼110 km depth), deeper than the interpreted melting depths of the WFB and Akaki (∼73 km depth) (Fig. 11A). Compositions from the Galema range do not plot along the same trend as other areas of focused magmatic activity within the MER on an olivine-plagioclase-quartz ternary diagram. While the compositions of the other areas of focused magmatic activity (WFB, SDFZ, and Akai) plot along experimentally derived boundary curves for the 3-phase saturated liquid at 30 kbar (Hirose and Kushiro, 1993; Walker et al., 1979), Galema range compositions plot in a region that may indicate a higher melting pressure (Fig. 11B). We suggest that the source melt of the Galema range was generated beneath thicker lithosphere in comparison to melts within the MER.

Rifting initially manifests as faulting and progresses toward the development of connected half grabens (Ebinger, 2005). The classic continental rift shape of a downthrown valley floor and steep, high-angle border faults develop and progressive faulting becomes more focused within the rift valley (Corti, 2012; Corti et al., 2004; Maccaferri et al., 2014; Wolfenden et al., 2005). These faults have accommodated strain since their formation, though it is difficult to ascertain if this accommodation of strain was of a constant magnitude throughout their history.

Thinning of the lithosphere occurs in response to the extensional forces induced during continental rifting. Recent estimates indicate that within the MER, both the crust and lithosphere are thinned beneath the rift valley (Bastow and Keir, 2011; Keir et al., 2015; Keranen et al., 2009; Rooney et al., 2014). This thinning of the lithosphere results in decompression and consequent melting of the underlying asthenosphere, generating the melt necessary for continental rifting processes to occur (Ebinger, 2005; Ebinger and Casey, 2001; Ebinger and Sleep, 1998; Langmuir and Forsyth, 2007; McKenzie and Bickle, 1988). As magma began to intrude at the base of the lithosphere, perpendicular to the σ₃ generated by the regional extensional stress, the thermal gradient rose, broadly modifying the lithosphere in such a way as to allow diking and magma ascent to proceed (Bastow and Keir, 2011; Bialas et al., 2010; Daniels et al., 2014; Havlin et al., 2013).

Following the generation of melt beneath the rift margin, subsidence of the half grabens forced magma laterally along the base of the crust to outside of the border faults (Corti et al., 2003; Corti et al., 2004; Maccaferri et al., 2014). As this melt was relocated to beneath the Southeastern Ethiopian Plateau, it raised the thermal gradient and assisted in the thermomechanical modification of the lithosphere at that location (Bastow and Keir, 2011; Bialas et al., 2010; Daniels et al., 2014; Havlin et al., 2013) (Fig. 12) This modification facilitated subsequent magmatic intrusion, which ascended to form the Galema range. Additionally, the thermal modification induced by magmatic intrusion would have the effect of reducing lithospheric strength to a degree that would make strain accommodation by faulting unfavorable within the extensional stress field (Buck, 2004). Thus, activity of the nearby Asela-Sire Border Fault would have ceased or been reduced during these events. The local stress field and presence of preexisting lithospheric heterogeneities within the crust, at the time of the formation of the Galema range, dictated the pathways of ascent and determined the location and morphology of the range (Mazzarini et al., 2013; Rooney et al., 2014). Magmatic intrusion through preexisting lithospheric weaknesses will still result in extensional strain accommodation through the addition of volume to the crust if the intruded material remains within the host rock and is not sourced from the host rock itself.

After the formation of the Galema range, magmatism in this off-rift location ceased due to several possible processes related to either continued extension or reduced magma supply or both. A temporary reduction in the extension rate would lead to a reduction in adiabatic melting of the asthenosphere, limiting the amount of magma available to flux into the lithosphere. This reduced magma supply would have the effect of decreasing the magma overpressure in ascending dikes, allowing them to become trapped at depths where this pressure equals the density-controlled stress within the surrounding lithosphere. If extension and resultant magma production continued to be diminished for a long enough period, the dikes that had previously ascended into the lithosphere would freeze, effectively ending magmatic activity.

As the lithosphere continued to deform and thin along the rift axis, magmatism shifted toward more focused activity along the axis due to increased magma flux created by adiabatic melting in the newly thinned areas (Ebinger, 2005; Morton et al., 1979). Magmatism may also shift location due to a change in the topography at the base of the lithosphere (as a result of progressive thinning), which channels melt material toward the axis (Corti, 2012; Keir et al., 2015). While the shift in magmatic activity toward the rift axis may be due to a change in topography at the base of the lithosphere, we modify this model to include magmatic intrusion of the Galema range dikes and the consequent unavoidable thermomechanical modification of the lithosphere (Bastow and Keir, 2011; Bialas et al., 2010; Buck, 2004; Daniels et al., 2014; Havlin et al., 2013).

Lithospheric thinning associated with continued extension generated an increased melt fraction at shallower depths (Langmuir and Forsyth, 2007; McKenzie and Bickle, 1988). The magma generated during this more advanced state of extension intruded the lithosphere at the same broad location, which was thermomechanically modified by the previous activity of the Galema range.

The greater melt fraction increased the magma flux into the lithosphere, allowing for the formation of larger dikes. These larger dikes would be able to more easily ascend to the surface, forming more well developed magmatic plumbing system with shallower levels of fractionation, as observed at the WFB (Rooney et al., 2014) (Fig. 12). Greater melt production would allow for an increase in the amount of magmatic material erupted, increasing surface topography and simultaneously reducing the extensional stress at the bottom of the lithospheric plate through magma injection. This dynamic allows magma to ascend to shallower levels before stalling (Behn et al., 2006).

As a continental rift evolves, magmatism is influenced less by preexisting lithospheric structures and more by the lithospheric strain profile due to the thermal modification of the lithosphere (Isola et al., 2014; Robertson et al., 2016). Magmatism became the dominant form of strain accommodation in the WFB as a result of the enhancement of the magmatic plumbing system assisted by, among other processes, the earlier focused magmatism of the Galema range. Furthermore, this suggests that magmatic strain accommodation does not linearly migrate toward the rift axis during rift evolution and may in fact be more distributed and episodic.

During the early stages of continental rifting, extensional strain is accommodated through faulting and stretching of the lithosphere, resulting in the formation of the high-angle border faults that will define the rift. Following this, magmatic intrusion begins to accommodate extensional strain, which should migrate and focus along the rift axis (Wolfenden et al., 2005). Recent work has suggested that strain accommodation in the Afar region via focused magmatic intrusion did not occur until very recently and that prior to this, the only expression of magmatic activity was the emplacement of large-scale flood basalts (Stab et al., 2015). These areas of focused magmatic intrusion may represent the organization of magmatic strain accommodation organizing into axial segments (Medynski et al., 2015). The Galema range may represent an intermediate stage of continental rifting, occurring after the formation of the rift-border faults but before the organization of axial magmatic segments. This intermediate stage is characterized by pulsed accommodation of extensional strain by focused magmatic intrusion outside the rift proper and accommodation of extensional strain by faulting and stretching of the lithosphere during periods between magmatic intrusions.

The Galema range, which is located on the Southeastern Ethiopian Plateau adjacent to the rift margin, is a zone of focused dike intrusion that displays morphological similarities to zones of focused magmatic intrusion within the Main Ethiopian Rift. Our results demonstrate that the depths of stalling of the mafic magmas of the Galema range are deeper than those at the WFB, consistent with a less developed magma plumbing system for the Galema range. We find that magma generation beneath the Galema range occurs at greater depths than beneath the WFB, indicating that the lithosphere had undergone relatively less extension in comparison to the modern rift. On the basis of these observations, and existing geodynamic constraints, we propose the following model explaining the origin of off-axis zones of focused intrusion: (1) the initial development of rift-border faults controlled the migration of melt into the lithospheric mantle. The transit of the melt through the crust was controlled, in part, by the establishment of a horizontal axis of least compressive stress during graben formation and thereby facilitating rift adjacent volcanism; and (2) with continued rift evolution, the transit of melt through the lithospheric mantle occurred along the same pathways established during the earlier rift-adjacent volcanism enhancing magmatism along one rift margin. However, variations in the amount of magma production and resultant magma overpressure in ascending dikes may have allowed for the eruption of within-rift zones of focused intrusion.

We propose that the Galema range represents an intermediate stage of continental rift evolution, occurring after the formation of the rift-border faults but before the formation of axial magmatic segments. In this hybrid rifting model, strain is accommodated by oscillations between both the rift-border faults and areas of focused magmatic intrusion. However, the amount of strain partitioned between each mechanism and the timing are not yet resolved. These findings have implications for the initial development of zones of focused magmatic intrusion within rifts. In particular, this work demonstrates that it is possible to have multiple periods of extension being accommodated through focused magmatic intrusion or displacement of faults along a single margin. These observations will need to be incorporated into ongoing debates about the timing of initiation of magma-facilitated extension (Stab et al., 2015; Wolfenden et al., 2005).

This work was supported by the National Aeronautics and Space Administration Planetary Geology and Geophysics Program (grant 10-PGG10-006). We would like to thank Giacomo Corti and the anonymous reviewer for their comments, which have greatly improved this work, and Shanaka de Silva for careful editorial handling.