Search

Ocean swell interacting with Antarctic ice shelves produces sustained (approximately, 2×106 cycles per year) gravity-elastic perturbations with deformation amplitudes near the ice front as large as tens to hundreds of nanostrain. This process is the most energetically excited during the austral summer, when sea ice-induced swell attenuation is at a minimum. A 2014–2017 deployment of broadband seismographs on the Ross Ice shelf, which included three stations sited, approximately, 2 km from the ice front, reveals prolific swell-associated triggering of discrete near-ice-front (magnitude≲0) seismic subevents, for which we identify three generic types. During some strong swell episodes, subevent timing becomes sufficiently phase-locked with swell excitation, to create prominent harmonic features in spectra calculated across sufficiently lengthy time windows via a Dirac comb effect, for which we articulate a theoretical development for randomized interevent times. These events are observable at near-front stations, have dominant frequency content between 0.5 and 20 Hz, and, in many cases, show highly repetitive waveforms. Matched filtering detection and analysis shows that events occur at a low-background rate during all swell states, but become particularly strongly excited during large amplitude swell at rates of up to many thousands per day. The superimposed elastic energy from swell-triggered sources illuminates the shelf interior as extensional (elastic plate) Lamb waves that are observable more than 100 km from the ice edge. Seismic swarms show threshold excitation and hysteresis with respect to rising and falling swell excitation. This behavior is consistent with repeated seismogenic fracture excitation and growth within a near-ice-front damage zone, encompassing fracture features seen in satellite imagery. A much smaller population of distinctly larger near-front seismic events, previously noted to be weakly associated with extended periods of swell perturbation, likely indicate calving or other larger-scale ice failures near the shelf front.

Recent efforts to characterize small (Mw

Observations of teleseismic earthquakes using broadband seismometers on the Ross Ice Shelf (RIS) must contend with environmental and structural processes that do not exist for land-sited seismometers. Important considerations are: (1) a broadband, multi-mode ambient wavefield excited by ocean gravity wave interactions with the ice shelf; (2) body wave reverberations produced by seismic impedance contrasts at the ice/water and water/seafloor interfaces and (3) decoupling of the solid Earth horizontal wavefield by the sub-shelf water column. We analyze seasonal and geographic variations in signal-to-noise ratios for teleseismic P-wave (0.5–2.0 s), S-wave (10–15 s) and surface wave (13–25 s) arrivals relative to the RIS noise field. We use ice and water layer reverberations generated by teleseismic P-waves to accurately estimate the sub-station thicknesses of these layers. We present observations consistent with the theoretically predicted transition of the water column from compressible to incompressible mechanics, relevant for vertically incident solid Earth waves with periods longer than 3 s. Finally, we observe symmetric-mode Lamb waves generated by teleseismic S-waves incident on the grounding zones. Despite their complexity, we conclude that teleseismic coda can be utilized for passive imaging of sub-shelf Earth structure, although longer deployments relative to conventional land-sited seismometers will be necessary to acquire adequate data.

In weathered bedrock aquifers, groundwater is stored in pores and fractures that open as rocks are exhumed and minerals interact with meteoric fluids. Little is known about this storage because geochemical and geophysical observations are limited to pits, boreholes, or outcrops or to inferences based on indirect measurements between these sites. We trained a rock physics model to borehole observations in a well-constrained ridge and valley landscape and then interpreted spatial variations in seismic refraction velocities. We discovered that P-wave velocities track where a porosity-generating reaction initiates in shale in three boreholes across the landscape. Specifically, velocities of 2.7 ± 0.2 km/s correspond with growth of porosity from dissolution of chlorite, the most reactive of the abundant minerals in the shale. In addition, sonic velocities are consistent with the presence of gas bubbles beneath the water table under valley and ridge. We attribute this gas largely to CO2 produced by 1) microbial respiration in soils as meteoric waters recharge into the subsurface and 2) the coupled carbonate dissolution and pyrite oxidation at depth in the rock. Bubbles may nucleate below the water table because waters depressurize as they flow from ridge to valley and because pores have dilated as the deep rock has been exhumed by erosion. Many of these observations are likely to also describe the weathering and flow path patterns in other headwater landscapes. Such combined geophysical and geochemical observations will help constrain models predicting flow, storage, and reaction of groundwater in bedrock systems.

The Ross Ice Shelf (RIS) is host to a broadband, multimode seismic wavefield that is excited in response to atmospheric, oceanic and solid Earth source processes. A 34-station broadband seismographic network installed on the RIS from late 2014 through early 2017 produced continuous vibrational observations of Earth's largest ice shelf at both floating and grounded locations. We characterize temporal and spatial variations in broadband ambient wavefield power, with a focus on period bands associated with primary (10–20 s) and secondary (5–10 s) microseism signals, and an oceanic source process near the ice front (0.4–4.0 s). Horizontal component signals on floating stations overwhelmingly reflect oceanic excitations year-round due to near-complete isolation from solid Earth shear waves. The spectrum at all periods is shown to be strongly modulated by the concentration of sea ice near the ice shelf front. Contiguous and extensive sea ice damps ocean wave coupling sufficiently so that wintertime background levels can approach or surpass those of land-sited stations in Antarctica.

We present a shear wave model of the West Antarctic upper mantle to ∼200 km depth with enhanced regional resolution from the 2016-2018 UK Antarctic Seismic Network. The model is constructed from the combination of fundamental mode Rayleigh wave phase velocities extracted from ambient noise (periods 8-25 s) and earthquake data by two-plane wave analysis (periods 20-143 s). We seek to (i) image and interpret structures against the tectonic evolution of West Antarctica, and (ii) extract information from the seismic model that can serve as boundary conditions in ice sheet and glacial isostatic adjustment modelling efforts. The distribution of low velocity anomalies in the uppermost mantle suggests that recent tectonism in the West Antarctic Rift System (WARS) is mainly concentrated beneath the rift margins and largely confined to the uppermost mantle (

New broadband seismic data from Botswana and South Africa have been combined with existing data from the region to develop improved P and S wave velocity models for investigating the upper mantle structure of southern Africa. Higher craton‐like velocities are imaged beneath the Rehoboth Province and parts of the northern Okwa Terrane and the Magondi Belt, indicating that the northern edge of the greater Kalahari Craton lithosphere lies along the northern boundary of these terranes. Lower off‐craton velocities are imaged beneath the Damara‐Ghanzi‐Chobe Belt, and may in part reflect thinning of the lithosphere beneath the incipient Okavango Rift. Lower velocities are also imaged to the north and northwest of the Bushveld Complex beneath parts of the Okwa Terrane, Magondi Belt, and Limpopo Belt, indicating that cratonic upper mantle in some areas beneath these terranes may have been modified by the 2.05‐Ga Bushveld and/or 1.1‐Ga Umkondo magmatic events.

We present an upper mantle P-wave velocity model for the Ross Sea Embayment (RSE) region of West Antarctica, constructed by inverting relative P-wave travel-times from 1881 teleseismic earthquakes recorded by two temporary broadband seismograph deployments on the Ross Ice Shelf, as well as by regional ice- and rock-sited seismic stations surrounding the RSE. Faster upper mantle P-wave velocities ( ∼ + 1 % ) characterize the eastern part of the RSE, indicating that the lithosphere in this part of the RSE may not have been reheated by mid-to-late Cenozoic rifting that affected other parts of the Late Cretaceous West Antarctic Rift System. Slower upper mantle velocities ( ∼ − 1 % ) characterize the western part of the RSE over a ∼500 km-wide region, extending from the central RSE to the Transantarctic Mountains (TAM). Within this region, the model shows two areas of even slower velocities ( ∼ − 1.5 % ) centered beneath Mt. Erebus and Mt. Melbourne along the TAM front. We attribute the broader region of slow velocities mainly to reheating of the lithospheric mantle by Paleogene rifting, while the slower velocities beneath the areas of recent volcanism may reflect a Neogene-present phase of rifting and/or plume activity associated with the formation of the Terror Rift. Beneath the Ford Ranges and King Edward VII Peninsula in western Marie Byrd Land, the P-wave model shows lateral variability in upper mantle velocities of ±0.5% over distances of a few hundred km. The heterogeneity in upper mantle velocities imaged beneath the RSE and western Marie Byrd Land, assuming no significant variation in mantle composition, indicates variations in upper mantle temperatures of at least 100 °C. These temperature variations could lead to differences in surface heat flow of ∼±10 mW/m2 and mantle viscosity of 102 Pa s regionally across the study area, possibly influencing the stability of the West Antarctic Ice Sheet by affecting basal ice conditions and glacial isostatic adjustment.

We present a 3D shear-wave velocity model of the southern African upper mantle developed using 30–200 s period Rayleigh waves recorded on regional seismic networks spanning the subcontinent. The model shows high velocities (∼4.7–4.8 km/s) at depths of 50–250 km beneath the Archean nucleus and several surrounding Paleoproterozoic and Mesoproterozoic terranes, placing the margin of the greater Kalahari Craton along the southern boundary of the Damara Belt and the eastern boundaries of the Gariep and Namaqua-Natal belts. At depths ≥250 km, there is little difference in velocities beneath the craton and off-craton regions, suggesting that the cratonic lithosphere extends to depths of about 200–250 km. Upper mantle velocities beneath uplifted areas of southern Africa are higher than the global average and significantly higher than beneath eastern Africa, indicating there that is little thermal modification of the upper mantle present today beneath the Southern African Plateau.

We report new P and S wave velocity models of the upper mantle beneath southern Africa using data recorded on seismic stations spanning the entire subcontinent. Beneath most of the Damara Belt, including the Okavango Rift, our models show lower than average velocities (−0.8% Vp; −1.2% Vs) with an abrupt increase in velocities along the terrane's southern margin. We attribute the lower than average velocities to thinner lithosphere (~130 km thick) compared to thicker lithosphere (~200 km thick) immediately to the south under the Kalahari Craton. Beneath the Etendeka Flood Basalt Province, higher than average velocities (0.25% Vp; 0.75% Vs) indicate thicker and/or compositionally distinct lithosphere compared to other parts of the Damara Belt. In the Rehoboth Province, higher than average velocities (0.3% Vp; 0.5% Vs) suggest the presence of a microcraton, as do higher than average velocities (1.0% Vp; 1.5% Vs) under the Southern Irumide Belt. Lower than average velocities (−0.4% Vp; −0.7% Vs) beneath the Bushveld Complex and parts of the Mgondi and Okwa terranes are consistent with previous studies, which attributed them to compositionally modified lithosphere resulting from Precambrian magmatic events. There is little evidence for thermally modified upper mantle beneath any of these terranes which could provide a source of uplift for the Southern African Plateau. In contrast, beneath parts of the Irumide Belt in southern and central Zambia and the Mozambique Belt in central Mozambique, deep‐seated low velocity anomalies (−0.7% Vp; −0.8% Vs) can be attributed to upper mantle extensions of the African superplume structure.

Ice shelves play an important role in buttressing land ice from reaching the sea, thus restraining the rate of grounded ice loss. Long-period gravity-wave impacts excite vibrations in ice shelves that can expand pre-existing fractures and trigger iceberg calving. To investigate the spatial amplitude variability and propagation characteristics of these vibrations, a 34-station broadband seismic array was deployed on the Ross Ice Shelf (RIS) from November 2014 to November 2016. Two types of ice-shelf plate waves were identified with beamforming: flexural-gravity waves and extensional Lamb waves. Below 20 mHz, flexural-gravity waves dominate coherent signals across the array and propagate landward from the ice front at close to shallow-water gravity-wave speeds (~70 m s−1). In the 20–100 mHz band, extensional Lamb waves dominate and propagate at phase speeds ~3 km s−1. Flexural-gravity and extensional Lamb waves were also observed by a 5-station broadband seismic array deployed on the Pine Island Glacier (PIG) ice shelf from January 2012 to December 2013, with flexural wave energy, also detected at the PIG in the 20–100 mHz band. Considering the ubiquitous presence of storm activity in the Southern Ocean and the similar observations at both the RIS and the PIG ice shelves, it is likely that most, if not all, West Antarctic ice shelves are subjected to similar gravity-wave excitation.

Although the East African rift system formed in cratonic lithosphere above a large‐scale mantle upwelling, some sectors have voluminous magmatism, while others have isolated, small‐volume eruptive centers. We conduct teleseismic shear wave splitting analyses on data from 5 lake‐bottom seismometers and 67 land stations in the Tanganyika‐Rukwa‐Malawi rift zone, including the Rungwe Volcanic Province (RVP), and from 5 seismometers in the Kivu rift and Virunga Volcanic Province, to evaluate rift‐perpendicular strain, rift‐parallel melt intrusion, and regional flow models for seismic anisotropy patterns beneath the largely amagmatic Western rift. Observations from 684 SKS and 305 SKKS phases reveal consistent patterns. Within the Malawi rift south of the RVP, fast splitting directions are oriented northeast with average delays of ~1 s. Directions rotate to N‐S and NNW north of the volcanic province within the reactivated Mesozoic Rukwa and southern Tanganyika rifts. Delay times are largest (~1.25 s) within the Virunga Volcanic Province. Our work combined with earlier studies shows that SKS‐splitting is rift parallel within Western rift magmatic provinces, with a larger percentage of null measurements than in amagmatic areas. The spatial variations in direction and amount of splitting from our results and those of earlier Western rift studies suggest that mantle flow is deflected by the deeply rooted cratons. The resulting flow complexity, and likely stagnation beneath the Rungwe province, may explain the ca. 17 Myr of localized magmatism in the weakly stretched RVP, and it argues against interpretations of a uniform anisotropic layer caused by large‐scale asthenospheric flow or passive rifting.

Buried deep beneath the Antarctic polar ice sheet, the geological structure and tectonic activity of East Antarctica have long remained unknown. The apparent lack of tectonic seismicity was thought to be anomalous relative to other continental interiors and has been attributed to a lack of intraplate stress due to the surrounding spreading ridges and low absolute plate velocity or to the weight of ice sheets increasing the normal stress. Here we report 27 intraplate tectonic earthquakes detected by the AGAP/GAMSEIS seismic array during 2009 in the interior of East Antarctica, which represents locally recorded seismicity in the region. The earthquakes are primarily extensional and located at shallow to mid-crustal depths beneath sedimentary basins aligned linearly adjacent to the Gamburtsev Subglacial Mountains. The basins may be part of an ancient continental rift system, which provides a zone of pre-existing tectonic weakness that focuses the seismicity. These events, when combined with events in published catalogues of Antarctic seismicity, indicate levels of seismicity in East Antarctica of the same order of magnitude as that of other stable cratons, such as the Canadian Shield.

The responses of the Ross Ice Shelf (RIS) to the September 16, 2015 8.3 (Mw) Chilean earthquake tsunami (> 75 s period) and to oceanic infragravity (IG) waves (50–300 s period) were recorded by a broadband seismic array deployed on the RIS from November 2014 to November 2016. Here we show that tsunami and IG-generated signals within the RIS propagate at gravity-wave speeds (∼70 m/s) as water-ice coupled flexural-gravity waves. IG band signals show measureable attenuation away from the shelf front. The response of the RIS to Chilean tsunami arrivals is compared with modeled tsunami forcing to assess ice shelf flexural-gravity wave excitation by very long period (VLP; > 300 s) gravity waves. Displacements across the RIS are affected by gravity-wave incident direction, bathymetry under and north of the shelf, and water layer and ice shelf thicknesses. Horizontal displacements are typically about 10 times larger than vertical displacements, producing dynamical extensional motions that may facilitate expansion of existing fractures. VLP excitation is continuously observed throughout the year, with horizontal displacements highest during the austral winter with amplitudes exceeding 20 cm. Because VLP flexural-gravity waves exhibit no discernable attenuation, this energy must propagate to the grounding zone. Both IG and VLP band flexural-gravity waves excite mechanical perturbations of the RIS that likely promote tabular iceberg calving, consequently affecting ice shelf evolution. Understanding these ocean-excited mechanical interactions is important to determine their effect on ice shelf stability to reduce uncertainty in the magnitude and rate of global sea level rise.

Rayleigh wave group and phase velocity measurements obtained from ambient noise and earthquake data at 51 broadband stations were used to construct the first 3-D crustal and upper mantle shear wave velocity model of Botswana. The model shows low crustal velocities associated with the Passarge and Nosop sedimentary basins, whereas the Kaapvaal, Zimbabwe, Maltahohe, and Congo Cratons are recognized by high mantle velocities. The lowest upper mantle shear wave velocity, beneath northeastern Botswana, is associated with the southwestern branch of the East African Rift System. This low-velocity mantle anomaly appears to be linked to the crust of the Okavango Rift Zone and the location of the 3 April 2017 Mw 6.5 earthquake in central Botswana. We suggest that fluids or melt at the base of the crust from the southward continuation of the East African Rift Zone triggered the intraplate earthquake in an extensional tectonic setting.

International audience; The crust and upper mantle of the Madagascar continental fragment remained largely unexplored until a series of recent broadband seismic experiments. An island-wide deployment of broadband seismic instruments has allowed the first study of phase velocity variations, derived from surface waves, across the entire island. Late Cenozoic alkaline intraplate volcanism has occurred in three separate regions of Madagascar (north, central and southwest), with the north and central volcanism active until

Seismicity patterns and focal mechanisms in southeastern Tanzania, determined from data recorded on temporary and permanent AfricaArray seismic stations, have been used to investigate the propagation direction of the Eastern branch of the East African Rift System southward from the Northern Tanzania Divergence Zone (NTDZ). Within the NTDZ, the rift zone is defined by three segments, the Eyasi segment to the west, the Manyara segment in the middle, and the Pangani segment to the east. Results show that most of the seismicity (~ 75%) extends to the south of the Manyara segment along the eastern margin of the Tanzania Craton, and at ~ 6–7° S latitude trends to the SE along the northern boundary of the Ruvuma microplate, connecting with a N–S zone of seismicity offshore southern Tanzania and Mozambique. A lesser amount of seismicity (~ 25%) is found extending from the SE corner of the Tanzania Craton at ~ 6–7° S latitude southwards towards Lake Nyasa. This finding supports a model of rift propagation via the Manyara segment to the southeast of the Tanzania Craton along the northern boundary of the Ruvuma microplate. However, given the limited duration of the seismic recordings used in this study, the possibility of another zone of extension developing to the south towards Lake Nyasa (Malawi) cannot be ruled out. Focal mechanisms along the boundary between the Victoria and the Ruvuma microplates and offshore southeastern Tanzania show a combination of normal and strike slip faulting indicating mainly extension with some sinistral motion, consistent with the mapped geologic faults and a clockwise rotation of the Ruvuma microplate.

The seismic velocity structure of Antarctica is important, both as a constraint on the tectonic history of the continent and for understanding solid Earth interactions with the ice sheet. We use Rayleigh wave array analysis methods applied to teleseismic data from recent temporary broadband seismograph deployments to image the upper mantle structure of central and West Antarctica. Phase velocity maps are determined using a two-plane-wave tomography method, and are inverted for shear velocity using a Monte-Carlo approach to estimate three-dimensional velocity structure. Results illuminate the structural dichotomy between the East Antarctic Craton and West Antarctica, with West Antarctica showing thinner crust and slower upper mantle velocity. West Antarctica is characterized by a 70-100 km thick lithosphere, underlain by a low velocity zone to depths of at least 200 km. The slowest anomalies are beneath Ross Island and the Marie Byrd Land dome, and are interpreted as upper mantle thermal anomalies possibly due to mantle plumes. The central Transantarctic Mountains are marked by an uppermost mantle slow velocity anomaly, suggesting that the topography is thermally supported. The presence of thin, higher velocity lithosphere to depths of about 70 km beneath the West Antarctic Rift System limits estimates of the regionally averaged heat flow to less than 90 mW/m2. The Ellsworth-Whitmore block is underlain by mantle with velocities that are intermediate between those of the West Antarctic Rift System and the East Antarctic Craton. We interpret this province as Precambrian continental lithosphere that has been altered by Phanerozoic tectonic and magmatic activity.

New P- and S-wave velocity models of the upper mantle from 100 to 400 km depth beneath the central portions of West Antarctica, obtained by inverting relative travel-times from teleseismic earthquakes recorded on Polar Earth Observing Network (POLENET/ANET) and UK Antarctic Network (UKANET) seismic stations between 2007 and 2017, reveal a heterogeneous upper mantle. A low velocity anomaly (−1.0% Vp; −2.0% Vs) imaged beneath Marie Byrd Land is attributed to thermally perturbed upper mantle of possible plume origin, and a low velocity anomaly imaged beneath the Pine Island Glacier and the mouth of Thwaites Glacier is interpreted as a rift-related thermal structure that may include warm mantle flowing from Marie Byrd Land. High velocity anomalies (≤0.8% Vp; 1.5% Vs) imaged in the central portion of the West Antarctic Rift System indicate the presence of lithosphere unmodified by tectonic activity since the Late Cretaceous formation of the rift system. Within the region of high velocities, localized low velocity anomalies beneath parts of the Bentley Subglacial Trench are suggestive of focused Cenozoic rifting. The models also show variable velocity structure beneath the Haag-Ellsworth Whitmore crustal block and low velocities beneath the Thurston Island-Eights Coast crustal block. The heterogenous upper mantle structure of central West Antarctica indicates that upper mantle temperatures could vary by 100 K or more over distances of less than 100 km, which may add complexity to solid earth-ice interactions and influence basal ice sheet conditions.

Continuous seismic observations across the Ross Ice Shelf reveal ubiquitous ambient resonances at frequencies >5 Hz. These firn‐trapped surface wave signals arise through wind and snow bedform interactions coupled with very low velocity structures. Progressive and long‐term spectral changes are associated with surface snow redistribution by wind and with a January 2016 regional melt event. Modeling demonstrates high spectral sensitivity to near‐surface (top several meters) elastic parameters. We propose that spectral peak changes arise from surface snow redistribution in wind events and to velocity drops reflecting snow lattice weakening near 0°C for the melt event. Percolation‐related refrozen layers and layer thinning may also contribute to long‐term spectral changes after the melt event. Single‐station observations are inverted for elastic structure for multiple stations across the ice shelf. High‐frequency ambient noise seismology presents opportunities for continuous assessment of near‐surface ice shelf or other firn environments.

Author Posting. © The Author(s), 2016. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 205 (2016): 785-795, doi:10.1093/gji/ggw036.

Broadband seismic stations were deployed across the Ross Ice Shelf (RIS) in November 2014 to study ocean gravity wave-induced vibrations. Initial data from three stations 100 km from the RIS front and within 10 km of each other show both dispersed infragravity (IG) wave and ocean swell-generated signals resulting from waves that originate in the North Pacific. Spectral levels from 0.001 to 10 Hz have the highest accelerations in the IG band (0.0025–0.03 Hz). Polarization analyses indicate complex frequency-dependent particle motions, with energy in several frequency bands having distinctly different propagation characteristics. The dominant IG band signals exhibit predominantly horizontal propagation from the north. Particle motion analyses indicate retrograde elliptical particle motions in the IG band, consistent with these signals propagating as Rayleigh-Lamb (flexural) waves in the ice shelf/water cavity system that are excited by ocean wave interactions nearer the shelf front.

SUMMARY New estimates of crustal thickness, Poisson’s ratio and crustal shear wave velocity have been obtained for 39 stations in Angola, Botswana, the Democratic Republic of Congo, Malawi, Mozambique, Namibia, Rwanda, Tanzania and Zambia by modellingP-wave receiver functions using the H–κ stacking method and jointly inverting the receiver functions with Rayleigh-wave phase and group velocities. These estimates, combined with similar results from previous studies, have been examined for secular trends in Precambrian crustal structure within the southern African subcontinent. In both Archean and Proterozoic terranes we find similar Moho depths [38–39 ± 3k mSD (standard deviation)], crustal Poisson’s ratio (0.26 ± 0.01 SD), mean crustal shear wave velocity (3.7 ± 0.1 km s −1 SD), and amounts of heterogeneity in the thickness of the mafic lower crust, as defined by shear wave velocities ≥4.0 km s −1 . In addition, the amount of variability in these crustal parameters is similar within each individual age grouping as between age groupings. Thus, the results provide little evidence for secular variation in Precambrian crustal structure, including between Meso- and Neoarchean crust. This finding suggests that (1) continental crustal has been generated by similar processes since the Mesoarchean or (2) plate tectonic processes have reworked and modified the crust through time, erasing variations in structure resulting from crustal genesis.

The Cameroon Volcanic Line (CVL) is a 1800 km long volcanic chain, extending SW-NE from the Gulf of Guinea into Central Africa, that lacks the typical age progression exhibited by hot spot-related volcanic tracks. This study investigates the upper mantle seismic structure beneath the CVL and surrounding regions to constrain the origin of volcanic lines that are poorly described by the classic plume model. Rayleigh wave phase velocities are measured at periods from 20 to 182 s following the two-plane wave methodology, using data from the Cameroon Seismic Experiment, which consists of 32 broadband stations deployed between 2005 and 2007. These phase velocities are then inverted to build a model of shear wave velocity structure in the upper mantle beneath the CVL. Results show that phase velocities beneath the CVL are reduced at all periods, with average velocities beneath the CVL deviating more than –2% from the regional average and +4% beneath the Congo Craton. This distinction is observed for all periods but is less pronounced for the longest periods measured. Inversion for shear wave velocity structure indicates a tabular low velocity anomaly directly beneath the CVL at depths of 50 to at least 200 km and a sharp vertical boundary with faster velocities beneath the Congo Craton. These observations demonstrate widespread infiltration or erosion of the continental lithosphere beneath the CVL, most likely caused by mantle upwelling associated with edge-flow convection driven by the Congo Craton or by lithospheric instabilities that develop due to the nearby edge of the African continent.

A quick rebound for Antarctic crust Earth's crust deforms under the load of glaciers and ice sheets. When these masses are removed, the crust rebounds at a time scale determined by the viscosity of the upper mantle. Using GPS, Barletta et al. found that the viscosity of the mantle under the West Antarctic Ice Sheet is much lower than expected. This means that as ice is lost, the crust rebounds much faster than previously expected. Although estimates of total ice loss have to be revised upward, the surprising finding indicates that the ice sheet may stabilize against catastrophic collapse. Science , this issue p. 1335

Located at the southernmost sector of the Western Branch of the East African Rift System, the Malawi Rift exemplifies an active, magma-poor, weakly extended continental rift. To investigate the controls on rifting, we image crustal and uppermost mantle structure beneath the region using ambient-noise and teleseismic Rayleigh-wave phase velocities between 9 and 100 s period. Our study includes six lake-bottom seismometers located in Lake Malawi (Nyasa), the first time seismometers have been deployed in any of the African rift lakes. Noise levels in the lake are lower than that of shallow oceanic environments and allow successful application of compliance corrections and instrument orientation determination. Resulting phase-velocity maps reveal slow velocities primarily confined to Lake Malawi at short periods (T 25 s), a prominent low-velocity anomaly exists beneath the Rungwe Volcanic Province at the northern terminus of the rift basin. Estimates of phase-velocity sensitivity indicates these low velocities occur within the lithospheric mantle and potentially uppermost asthenosphere, suggesting that mantle processes may control the association of volcanic centres and the localization of magmatism. Beneath the main portion of the Malawi Rift, a modest reduction in velocity is also observed at periods sensitive to the crust and upper mantle, but these velocities are much higher than those observed beneath Rungwe.

Online Material: Table of station parameters; figures of mean acceleration power differences, interpolated noise maps. Seismographic coverage of Antarctica prior to 2007 consisted overwhelmingly of a handful of long running and sporadically deployed transient stations, many of which were principally collocated with scientific research stations. Despite very cold temperatures, sunless winters, challenging logistics, and extreme storms, recent developments in polar instrumentation driven by new scientific objectives have opened up the entirety of Antarctica to year‐round and continuous seismological observation (e.g., Nyblade et al. , 2012). Motivations for these recent studies include improved understanding of seismogenic, volcanic, tectonic and glaciological processes, heat flow, dynamic glaciological/ocean interactions, and mantle viscosity. Such studies contribute generally to improvements in understanding the geophysical, geological, and glaciological history of the continent and how these processes interact with the past and present state of the glaciological and climate system (e.g., Winberry et al. , 2009; Hansen et al. , 2010; West et al. , 2010; Winberry et al. , 2011; Heeszel et al. , 2013; Lough et al. , 2013; Chaput et al. , 2014; Accardo et al. , 2014), including processes relevant to glacial isostatic adjustment and sea level rise (Intergovernmental Panel on Climate Change [IPCC] Report, 2007). In addition, microseisms arising from ocean wave activity contain useful climate proxy information on the state and variability of the relatively poorly sensed southern oceans (Aster et al. , 2008; Stutzmann et al. , 2009; Aster et al. , 2010), and such observations are sensitive to sea ice concentration and areal coverage in the polar regions (Grob et al. , 2011; Tsai and McNamara, 2011; Koch et al. , 2013). This characterization of the seismic noise environment of Antarctica, documentation of instrument performance, and comparisons of installation conditions (e.g., ice vaults vs. rock sites) is intended to facilitate optimization …

Large earthquakes can trigger seismicity in remote regions. Analysis of seismic data from Antarctica reveals ice quakes coincident with passing seismic waves from the 2010 Chile earthquake, suggesting that the ice sheet is sensitive to such triggers.

[1] P-to-S receiver functions (PRFs) from the Polar Earth Observing Network (POLENET) GPS and seismic leg of POLENET spanning West Antarctica and the Transantarctic Mountains deployment of seismographic stations provide new estimates of crustal thickness across West Antarctica, including the West Antarctic Rift System (WARS), Marie Byrd Land (MBL) dome, and the Transantarctic Mountains (TAM) margin. We show that complications arising from ice sheet multiples can be effectively managed and further information concerning low-velocity subglacial sediment thickness may be determined, via top-down utilization of synthetic receiver function models. We combine shallow structure constraints with the response of deeper layers using a regularized Markov chain Monte Carlo methodology to constrain bulk crustal properties. Crustal thickness estimates range from 17.0±4 km at Fishtail Point in the western WARS to 45±5 km at Lonewolf Nunataks in the TAM. Symmetric regions of crustal thinning observed in a transect deployment across the West Antarctic Ice Sheet correlate with deep subice basins, consistent with pure shear crustal necking under past localized extension. Subglacial sediment deposit thicknesses generally correlate with trough/dome expectations, with the thickest inferred subice low-velocity sediment estimated as ∼0.4 km within the Bentley Subglacial Trench. Inverted PRFs from this study and other published crustal estimates are combined with ambient noise surface wave constraints to generate a crustal thickness map for West Antarctica south of 75°S. Observations are consistent with isostatic crustal compensation across the central WARS but indicate significant mantle compensation across the TAM, Ellsworth Block, MBL dome, and eastern and western sectors of thinnest WARS crust, consistent with low density and likely dynamic, low-viscosity high-temperature mantle.

Numerous volcanoes exist in Marie Byrd Land, a highland region of West Antarctica. High heat flow through the crust in this region may influence the stability of the West Antarctic Ice Sheet 1‐4 . Volcanic activity progressed from north to south in the Executive Committee mountain range between the Miocene and Holocene epochs, but there has been no evidence for recent magmatic activity 5‐7 . Here we use a recently deployed seismic network to show that in 2010 and 2011, two swarms of seismic activity occurred at 25‐40 km depth beneath subglacial topographic and magnetic highs, located 55 km south of the youngest subaerial volcano in the Executive Committee Range. We interpret the swarm events as deep long-period earthquakes based on their unusual frequency content. Such earthquakes occur beneath active volcanoes, are caused by deep magmatic activity and, in some cases, precede eruptions 8‐11 . We also use radar profiles to identify a prominent ash layer in the ice overlying the seismic swarm. Located at 1,400 m depth, the ash layer is about 8,000 years old and was probably sourced from the nearby Mount Waesche volcano. Together, these observations provide strong evidence for ongoing magmatic activity and demonstrate that volcanism continues to migrate southwards along the Executive Committee Range. Eruptions at this site are unlikely to penetrate the 1.2 to 2-km-thick overlying ice, but would generate large volumes of melt water that could significantly affect ice stream flow.

Region‐specific empirically based ground‐truth (EBGT) criteria used to estimate the epicentral‐location accuracy of seismic events have been developed for the Main Ethiopian Rift and the Tibetan plateau. Explosions recorded during the Ethiopia–Afar Geoscientific Lithospheric Experiment (EAGLE), the International Deep Profiling of Tibet, and the Himalaya (INDEPTH III) experiment provided the necessary GT0 reference events. In each case, the local crustal structure is well known and handpicked arrival times were available, facilitating the establishment of the location accuracy criteria through the stochastic forward modeling of arrival times for epicentral locations. In the vicinity of the Main Ethiopian Rift, a seismic event is required to be recorded on at least 8 stations within the local ![Graphic][1] crossover distance and to yield a network‐quality metric of less than 0.43 in order to be classified as EBGT595% (GT5 with 95% confidence). These criteria were subsequently used to identify 10 new GT5 events with magnitudes greater than 2.1 recorded on the Ethiopian Broadband Seismic Experiment (EBSE) network and 24 events with magnitudes greater than 2.4 recorded on the EAGLE broadband network. The criteria for the Tibetan plateau are similar to the Ethiopia criteria, yet slightly less restrictive as the network‐quality metric needs to be less than 0.45. Twenty‐seven seismic events with magnitudes greater than 2.5 recorded on the INDEPTH III network were identified as GT5 based on the derived criteria. When considered in conjunction with criteria developed previously for the Kaapvaal craton in southern Africa, it is apparent that increasing restrictions on the network‐quality metric mirror increases in the complexity of geologic structure from craton to plateau to rift. [1]: /embed/inline-graphic-1.gif

SUMMARY An expanded model of the 3-D shear wave velocity structure of the uppermost mantle beneath easternAfricahasbeendevelopedusingearthquakesrecordedbytheAfricaArrayEastAfrican SeismicExperimentinconjunctionwithdatafrompermanentstationsandpreviouslydeployed temporary stations. The combined data set comprises 331 earthquakes recorded on a total of 95 seismic stations spanning Kenya, Uganda, Tanzania, Zambia and Malawi. In this study, data from 149 earthquakes were used to determine fundamental-mode Rayleigh wave phase velocities at periods ranging from 20 to 182s using the two-plane wave method, and then combined with the similarly processed published measurements and inverted for a 3-D shear wave velocity model of the uppermost mantle. New features in the model include (1) a low-velocity region in western Zambia, (2) a high-velocity region in eastern Zambia, (3) a low-velocity region in eastern Tanzania and (4) low-velocity regions beneath the Lake Malawi rift. When considered in conjunction with mapped seismicity, these results support a secondarywesternriftbranchstrikingsouthwestwardsfromLakeTanganyika,likelyexploiting therelatively weak lithosphere ofthesouthernKibaran Beltbetween theBangweulu Blockand the Congo Craton. We estimate a lithospheric thickness of ∼150–200km for the substantial fast shear wave anomaly imaged in eastern Zambia, which may be a southward subsurface extension of the Bangweulu Block. The low-velocity region in eastern Tanzania suggests that the eastern rift branch trends southeastwards offshore eastern Tanzania coincident with the purported location of the northern margin of the proposed Ruvuma microplate. Pronounced velocity lows along the Lake Malawi rift are found beneath the northern and southern ends of the lake, but not beneath the central portion of the lake.

[1] New estimates of seismic anisotropy from shear wave splitting measurements in eastern Africa reveal a pattern of seismic anisotropy dominated by a NE alignment of fast polarization directions with local changes around the thick Archean lithosphere of the Tanzania craton. The overall pattern is consistent with mantle flow from the African superplume but not with absolute plate motion, a plume head, or fossil anisotropy in the lithosphere. In combination with tomographic images of the African superplume, this finding suggests that plateau uplift, volcanism, and continental breakup along the Afro-Arabian rift system is strongly influenced by flow from the lower mantle and indicates a connection between lower mantle processes and the tectonic deformation of the Earth's surface.

Onshore gravity measurements were conducted over the inland traverse routes from Syowa Station (69.0S, 39.6E; SYO) to Dome-F (77.4°S, 39.6°E) by the Japanese Antarctic Research Expeditions (1992, 1997, and 1998). The crustal density structure between Lutzow-Holm Bay (LHB) and the inland plateau was inferred from Bouguer gravity anomalies along these routes by assuming an initial model of the P-wave structure derived from seismic refraction surveys in the LHB. A decrease in the Bouguer anomalies down to − 200 mGal toward the inland plateau indicate that the crust thickens from 38–40 km at LHB to 48–50 km beneath Dome-F. During the International Polar Year in 2007–2008, the GAmburtsev Mountain SEISmic experiment (GAMSEIS) deployed many broadband stations over the large highland on the ice sheet from the crest of the Gambursev Subglacial Mountains (GSM) to the vicinity of Dome-F. S-wave receiver functions and Rayleigh phase velocities determined using the GAMSEIS data indicate that the cratonic crust surrounding the GSM is 40–47 km thick. These thickness estimates agree with those beneath Dome-F from the JARE's gravity surveys and are also consistent with average Pre-Cambrian terrains. Beneath the GSM, the crustal thickness increases to 55–58 km and has been interpreted as providing isostatic compensation for the high mountain elevations. Accordingly, a long-distance crustal model extending over 3000 km from LHB to Dome-F and to the GSM was compiled, providing significant information for future studies of the tectonic evolution of the Antarctic craton.

Subglacial seismicity reveals information about glacier behaviour. Analysis of repeated seismic events beneath an Antarctic outlet glacier is consistent with sliding of debris-laden ice over a bedrock asperity with an event frequency that is modulated by the ocean tides.

SUMMARY The Cameroon Volcanic Line (CVL) is an 1800-km-long line of Cenozoic volcanoes that does not show a chronological progression consistent with hotspot-related volcanism. We investigate seismic anisotropy to determine the upper-mantle lattice preferred orientation and constrain the mantle flow pattern using a temporary array of 32 broad-band seismographs deployed throughout Cameroon between 2005 and 2007 along with two additional permanent seismographs in adjacent countries. We determine the fast direction and lag time beneath each station by stacking SKS and SKKS splitting measurements from multiple events. The results indicate four regions with different splitting parameters. The Congo Craton in southern Cameroon and the Garoua rift region in northeast Cameroon have northeast–southwest-oriented fast directions and split times of about 1 s. Between the Congo Craton and the CVL, in central Cameroon, the fast directions are variable and have small splitting times of 0.3 s or less. Along the CVL, where previous studies show a strong slow velocity anomaly in the mantle, the fast direction is oriented approximately north–south, with splitting times of about 0.7 s. We interpret measurements from southern Cameroon and northeast Cameroon as indications of lattice-preferred orientation frozen into the Congo Craton and subcontinental lithosphere related to relict plate motion and deformation. The distinct pattern of splitting along the CVL suggests the existence of small-scale convection in the asthenosphere related to the formation of the CVL, perhaps driven by the adjacent cold edge of the Congo Craton.

SUMMARY The shear wave velocity structure of the upper mantle beneath the East African plateau has been investigated using teleseismic surface waves recorded on new broadband seismic stations deployed in Uganda and Tanzania, as well as on previously deployed stations in Tanzania and Kenya. Rayleigh wave phase velocities at periods between 20 and 182 s, measured with a two-plane wave method, have been used to create phase velocity maps, and dispersion curves extracted from the maps have been inverted to obtain a quasi-3-D shear wave velocity model of the upper mantle. We find that phase velocities beneath the Tanzania Craton and areas directly north and west of the craton are faster, at all periods, than those beneath the Western and Eastern branches of the East African Rift System. At periods 400 K. Consequently, it is likely that the velocity reduction reflects a combination of thermal and compositional changes, and also possibly the presence of partial melt. The width and thickness of the low velocity anomaly is greater than typically expected for a plume head and is more easily attributed to an upward continuation of the lower mantle African superplume structure into the upper mantle.

SUMMARY The thickness of the mantle transition zone beneath Cameroon has been mapped using data from the 2005–2007 Cameroon Broadband Seismic Experiment to evaluate models for the originof the Cameroon Volcanic Line (CVL).Some 2200 receiver functions have been stacked using a 3-D velocity model, revealing Ps conversions from the mantle transition zone discontinuities at depths of ∼410 and 660 km. Results yield a nearly uniform transition zone thickness (251 ± 10 km) that is similar to the global average, implying that any thermal anomalies in the upper mantle beneath the CVL do not extend as deep as the transition zone. This finding, when combined with regional P and S velocity models of the mantle, supports an explanation attributing the origin of the CVL to upper mantle processes such as an edge flow convection cell in the upper mantle along the northwestern side of the Congo Craton lithosphere.

We present new one-dimensional SH-wave velocity models of the upper mantle beneath the Kalahari craton in southern Africa obtained from waveform inversion of regional seismograms from an Mw = 5.9 earthquake located near Lake Tanganyika recorded on broadband seismic stations deployed during the 1997–1999 Southern African Seismic Experiment. The velocity in the lithosphere beneath the Kalahari craton is similar to that of other shields, and there is little evidence for a significant low velocity zone beneath the lithosphere. The lower part of the lithosphere, from 110 to 220 km depth, is slightly slower than beneath other shields, possibly due to higher temperatures or a decrease in Mg number (Mg#). If the slower velocities are caused by a thermal anomaly, then slightly less than half of the unusually high elevation of the Kalahari craton can be explained by shallow buoyancy from a hot lithosphere. However, a decrease in the Mg# of the lower lithosphere would increase the density and counteract the buoyancy effect of the higher temperatures. We obtain a thickness of 250 ± 30 km for the mantle transition zone, which is similar to the global average, but the velocity gradient between the 410 and 660 km discontinuities is less steep than in global models, such as PREM and IASP91. We also obtain velocity jumps of between 0.16 ± 0.1 and 0.21 ± 0.1 km/s across the 410 km discontinuity. Our results suggest that there may be a thermal or chemical anomaly in the mantle transition zone, or alternatively that the shear wave velocity structure of the transition zone in global reference models needs to be refined. Overall, our seismic models provide little support for an upper mantle source of buoyancy for the unusually high elevation of the Kalahari craton, and hence the southern African portion of the African Superswell.

SUMMARY Broad-band seismic data from the southern African seismic experiment and the AfricaArray network are used to investigate the seismic velocity structure of the upper mantle beneath southern Africa, and in particular beneath the Kaapvaal Craton. A two-plane approximation method that includes a finite frequency sensitivity kernel is employed to measure Rayleigh wave phase velocities, which are inverted to obtain a quasi-3-D shear wave velocity model of the upper mantle. We find phase velocities for the Kaapvaal Craton and surrounding mobile belts that are comparable to those reported by previous studies, and we find little evidence for variation from east to west across the Namaqua–Natal Belt, a region not well imaged in previous studies. A high-velocity upper-mantle lid is found beneath the Kaapvaal Craton and most of southern Africa. For the Kaapvaal Craton, the thickness of the lid (∼150–200 km) is consistent with the lid thicknesses reported in many previous studies. The cratonic lid is underlain by a ∼100-km thick low-velocity zone with a 3.9 per cent maximum velocity reduction. By comparing the velocity model to those published for other Archean cratons, we find few differences, and therefore conclude that there is little evidence in the shear wave velocity structure of the mantle to indicate that the southern African plateau is supported by an upper-mantle thermal anomaly.

The Gamburtsev Subglacial Mountains (GSM), located in central East Antarctica, are one of the most enigmatic tectonic features on Earth. Buried beneath several kilometers of ice, the mountains are characterized by peak elevations reaching ~ 3000 m above sea level. In this study, new data from the Gamburtsev Antarctic Mountains Seismic Experiment (GAMSEIS) are presented, which substantially improve constraints on the crustal and upper mantle structure in this region. S-wave receiver functions and Rayleigh wave phase velocities are used to analyze data from the GAMSEIS deployment and to improve estimates of crustal thickness beneath the East Antarctic craton and the GSM. Our results indicate that the cratonic crust surrounding the GSM is ~ 40–45 km thick. Beneath the GSM, the crust thickens to ~ 55–58 km and provides isostatic support for the high mountain elevations. It has been suggested that thicker crust beneath the GSM may reflect magmatic underplating associated with a mantle plume. However, considering our results with those from other previous and ongoing studies, we instead favor models in which the GSM are an old continental feature associated with either Proterozoic or Paleozoic tectonic events.

The Cameroon Volcanic Line (CVL) is a major geologic feature that cuts across Cameroon from the south west to the north east. It is a unique volcanic lineament which has both an oceanic and a continental sector and consists of a chain of Tertiary to Recent, generally alkaline volcanoes stretching from the Atlantic island of Pagalu to the interior of the African continent. The oceanic sector includes the islands of Bioko (formerly Fernando Po) and Sao Tome and Principe while the continental sector includes the Etinde, Cameroon, Manengouba, Bamboutos, Oku and Mandara mountains, as well as the Adamawa and Biu Plateaus. In addition to the CVL, three other major tectonic features characterize the region: the Benue Trough located northwest of the CVL, the Central African Shear Zone (CASZ), trending N70 degrees E, roughly parallel to the CVL, and the Congo Craton in southern Cameroon. The origin of the CVL is still the subject of considerable debate, with both plume and non-plume models invoked by many authors (e.g., Deruelle et al., 2007; Ngako et al, 2006; Ritsema and Allen, 2003; Burke, 2001; Ebinger and Sleep, 1998; Lee et al, 1994; Dorbath et al., 1986; Fairhead and Binks, 1991; King and Ritsema,more » 2000; Reusch et al., 2010). Crustal structure beneath Cameroon has been investigated previously using active (Stuart et al, 1985) and passive (Dorbath et al., 1986; Tabod, 1991; Tabod et al, 1992; Plomerova et al, 1993) source seismic data, revealing a crust about 33 km thick at the south-western end of the continental portion of the CVL (Tabod, 1991) and the Adamawa Plateau, and thinner crust (23 km thick) beneath the Garoua Rift in the north (Stuart et al, 1985) (Figure 1). Estimates of crustal thickness obtained using gravity data show similar variations between the Garoua rift, Adamawa Plateau, and southern part of the CVL (Poudjom et al., 1995; Nnange et al., 2000). In this study, we investigate further crustal structure beneath the CVL and the adjacent regions in Cameroon using 1-D shear wave velocity models obtained from the joint inversion of Rayleigh wave group velocities and P-receiver functions for 32 broadband seismic stations. From the 1-D shear wave velocity models, we obtain new insights into the composition and structure of the crust and upper mantle across Cameroon. After briefly reviewing the geological framework of Cameroon, we describe the data and the joint inversion method, and then interpret variations in crustal structure found beneath Cameroon in terms of the tectonic history of the region.« less