One of the defining characteristics of hotspot volcanism is the presence of a long-lived, linear chain of age-progressive volcanoes created by the movement of the lithosphere over a stationary melting anomaly. However, the spatial distribution of volcanism at hotspots is often complex and highly variable suggesting that the relationship between magma generation and magma transport at hotspots is poorly understood. Here, I present the results of the first systematic quantitative characterization of the spatial distribution of volcanism at oceanic hotspots. In the first study I develop a novel methodology to characterize the across-strike distribution of volcanism at hotspots and apply it to a catalog of 40 oceanic hotpots. I find that only 25% (10/40) of hotspots exhibit the simple single-peak profile predicted by geodynamic models of melt generation in mantle plumes. The remaining 75% (30/40) of hotspots exhibit a dual- or multi-peak pattern. In the second study, I focus on the across-strike distribution of volcanism at the oceanic hotspots that are sourced by a deep-rooted mantle plume. 14 out of the 15 consensus plume-fed hotspots exhibit a dual-peaked across-strike profile. The spacing between these peaks display a strong negative correlation with lithospheric age, in direct contrast to models of inter-volcanic spacing controlled by elastic plate thickness. This relation suggests a different mechanism controls volcanic spacing at plume-fed hotspots. In the third chapter, I investigate variations in the average topographic profiles over time along the two longest and best-constrained oceanic hotspot tracks: Hawaii and Louisville. I find that the dual-peak across-strike profile of volcanism is a persistent feature at the Louisville hotspot over the entire length of the track examined (spanning a period of more than 65 Myr). In contrast, the dual-peak profile of volcanism at Hawaii is only evident along the most recent portion of the track (i.e., over the last 5 Myr). In total, this thesis represents a significant step foreword in the collective understanding of hotspot volcanism, and introduces a new diagnostic tool for analysis of hotspot influenced seafloor topography.

Oceanic lithosphere constitutes the upper boundary layer of the Earth’s convecting mantle. Its structure and evolution provide a vital window on the dynamics of the mantle and important clues to how the motions of Earth’s surface plates are coupled to convection in the mantle below. The three-dimensional shear-velocity structure of the upper mantle beneath the Atlantic Ocean is investigated to gain insight into processes that drive formation of oceanic lithosphere. Travel times are measured for approximately 10,000 fundamental-mode Rayleigh waves, in the period range 30-130 seconds, traversing the Atlantic basin. Paths with >30% of their length through continental upper mantle are excluded to maximize sensitivity to the oceanic upper mantle. The lateral distribution of Rayleigh wave phase velocity in the Atlantic upper mantle is explored with two approaches. One, phase velocity is allowed to vary only as a function of seafloor age. Two, a general two-dimensional parameterization is utilized in order to capture perturbations to age-dependent structure. Phase velocity shows a strong dependence on seafloor age, and removing age-dependent velocity from the 2-D maps highlights areas of anomalously low velocity, almost all of which are proximal to locations of hotspot volcanism. Depth-dependent variations in vertically-polarized shear velocity (Vsv) are determined with two sets of 3-D models: a layered model that requires constant VSV in each depth layer, and a splined model that allows VSV to vary continuously with depth. At shallow depths (~75 km) the seismic structure shows the expected dependence on seafloor age. At greater depths (~200 km) high-velocity lithosphere is found only beneath the oldest seafloor; velocity variations beneath younger seafloor may result from temperature or compositional variations within the asthenosphere. The age-dependent phase velocities are used to constrain temperature in the mantle and show that, in contrast to previous results for the Pacific, phase velocities for the Atlantic are not consistent with a half-space cooling model but are best explained by a plate-cooling model with thickness of 75 km and mantle temperature of 1400oC. Comparison with data such as basalt chemistry and seafloor elevation helps to separate thermal and compositional effects on shear velocity.

Mineral chemistry and petrography of basalts from the Kimama drill core recovered by Hotspot: Snake River Scientific Drilling Project, Idaho establish crystallization conditions of these lavas. Twenty-three basalt samples, from 20 individual lava flows were sampled from the upper 1000 m (of the 1912 m drilled) core drilled on the axis of the Snake River Plain, and represent approximately 3 m.y. of volcanism (rocks at the bottom of the hole are ~6 Ma). Rock from the upper 1000 m are typically fresh, while those lower in the core are more altered and are less likely to preserve fresh phenocrysts to analyze. Intratelluric phenocrysts (pre-eruption) are: olivine, plagioclase and Cr-spinel inclusions in olivine and plagioclase; groundmass phases (post-eruption) are: olivine, plagioclase, clinopyroxene, magnetite and ilmenite. Olivine core compositions range from Fo84-68, plagioclase cores range from An80-62, clinopyroxene ranges in composition from Wo47-34, En47-28, Fs30-15, spinel inclusions are Cr (up to 20 wt % Cr2O3) and Al-rich (up to 35 wt % Al2O3) and evolve to lower concentrations of Cr and Al and higher Fe and Ti, chromian titanomagnetite to magnetite, and ilmenite are groundmass oxide phases. Thermobarometry of Kimama core basalts indicates that the phenocryst phases crystallized at temperatures of 1155 to 1255°C at depths of 7 to 17 km, which is within or near the seismically imaged mid-crustal sill. Plagioclase hygrometry suggests that these lavas are relatively anhydrous with less than 0.4 wt % H2O. Groundmass phases crystallized at lower temperatures (<1140°C) after eruption. Oxygen fugacity inferred from Fe-Ti oxide equilibria is at or just below the QFM buffer. The origin of the basaltic rocks of the Snake River Plain has been attributed to a mantle plume or to other, shallow mantle processes. Mineral and whole rock major and trace element geochemistry of the olivine tholeiites from the Kimama core are used to distinguish between these two sources (deep or shallow mantle). Whole rock compositions were corrected for plagioclase and olivine fractionation to calculate primary liquids to estimate mantle potential temperatures. Olivine phenocrysts have the pyroxenite source characteristics of low Mn and Ca, but a peridotite source characteristic of low Ni. Thus, trace element models were used to test whether there is pyroxenite in the source of the Snake River Plain basalts, as hypothesized for Hawaii and other plume-related hotspots (e.g., Sobolev et al., 2005; Herzberg, 2011). Olivine chemistry and trace element models establish that the basalt source is a spinel peridotite, not a pyroxenite. The average mantle potential temperature obtained for these samples is 1577°C, 177°C hotter than ambient mantle, suggesting that the basaltic liquids were derived from a thermal plume. Silica activity barometry shows that melt segregation occurs between 80 and 110 km depth, which is within or very near the spinel stability field, and suggests that the lithosphere has been eroded by the plume to a maximum depth of 80 km, and recent mantle tomography suggests that it may be even thinner.