Active Seismic Data Research Articles

Effective Active and Passive Seismics for the Characterization of Urban and Remote Areas: Four Channels for Seven Objective Functions

An efficient system for the joint acquisition and analysis of multi-component active and passive seismic data is presented. It is shown how, in spite of the limited field equipment (the system requires just a 4-channel seismograph, one 3-component and four vertical-component geophones), it is nevertheless possible to define up to seven different (but mutually related and complementary) objects used to constrain a multi-objective joint inversion capable of providing a robust subsurface shear-wave velocity (VS) profile for both geotechnical and seismic-hazard studies. The presented approach relies on acquisition techniques that require simple and straightforward field procedures useful in particular, but not solely, in the characterization of urban and remote areas where, due to logistical problems, standard acquisition procedures cannot be easily applied. Active data recorded by a single 3-component geophone are processed so to define up to five objective functions: the group-velocity spectra of the three components, the radial-to-vertical spectral ratio and the Rayleigh-wave particle motion frequency curve. Passive data are used to compute two further objects: the horizontal-to-vertical spectral ratio and the phase-velocity dispersion curve obtained via miniature array analysis of microtremors. These seven objects are jointly inverted by means of a multi-objective inversion procedure based on the Pareto criterion. Performances are assessed through a comprehensive field dataset acquired in an urban area of NW-Italy. The consistency of the overall procedure is assessed by comparing the results with the analyses accomplished by considering classical multi-channel active and passive data and methodologies (multi-component MASW, multichannel analysis of surface waves and ESAC, extended spatial auto-correlation).

Sedimentary basin investigation using receiver function: an East African Rift case study

SUMMARYWe apply receiver function (RF) methodology to map the geometry of a sedimentary basin along a ∼10-km-long profile of broadband seismometers that recorded continuously for approximately 3 months. For a subset of the stations, we apply the Neighbourhood Algorithm inversion scheme, to quantify the geometry of basin bounding fault directly beneath the stations. We compare our results with active reflection seismic data and with the lithostratigraphy from a well located along the profile. We find that the P-to-s conversions from the sediments–basement interface (SBI), recorded in RF data sets together with information on intrabasin structures, are useful for obtaining high resolution images of the basin. The depth of the SBI derived from RF inversion is consistent (within ∼0.4 km) with the estimates from active reflection seismic and the well data. This study highlights that analysis of teleseismic waveforms can retrieve relevant information on the structure of a sedimentary basin.

Shear-wave velocity estimation using a combination of ambient noise from small aperture array and small-scale active seismic measurements: a case study in the area of the natural gas fields of Northern Germany

In recent years, numerous induced seismic events have occurred in the proximity of the natural gas field in Northern Germany. To monitor the seismicity and to asses the seismic hazard potential, a local monitoring network was installed in the area. Focusing on the seismicity hazard assessment, a major challenge is the characterisation of potential site effects due to local soil characteristics. This is quantitatively performed by estimating the shear-wave velocity (V s) variation with subsurface layer thickness. Such local effects can only be covered with a coarse spatial resolution due to the limited number of monitoring stations. Profiles were determined at three test sites (Langwedel, Walle and Bomlitz) by using a combined approach of small aperture 2D array ambient noise and small-scale active 1D measurements. The high-resolution frequency-wavenumber (HRFK), spatial autocorrelation (SPAC) and multichannel analysis of surface waves (MASW) methods were applied to the recorded ambient noise and active seismic data using various array sizes supplemented by the active measurements. This jointly allowed obtaining phase velocity dispersion curves covering a wide frequency range from 2 up to 32 Hz. The inversion of the obtained dispersion curves results in average S-wave velocity profiles down to depths of 70 m, identifying thin near-surface layers of a few meters as well as thicker layers of tens of meters in greater depth. A comparison with available borehole data shows a good correlation with the layering. Additionally, to asses the impact of a seismic event at the test sites, PGV estimations for various seismic events were performed. The final results of the test surveys demonstrate that the combined approach represents a suitable tool for near-surface characterisation, which can be used to improve the seismic hazard assessment in the area of the natural gas fields in Northern Germany.

Deep Structure of Northern Apennines Subduction Orogen (Italy) as Revealed by a Joint Interpretation of Passive and Active Seismic Data

AbstractThe Apennines is a well‐studied orogeny formed by the accretion of continental slivers during the subduction of the Adriatic plate, but its deep structure is still a topic of controversy. Here we illuminated the deep structure of the Northern Apennines belt by combining results from the analysis of active seismic (CROP03) and receiver function data. The result from combining these two approaches provides a new robust view of the structure of the deep crust/upper mantle, from the back‐arc region to the Adriatic subduction zone. Our analysis confirms the shallow Moho depth beneath the back‐arc region and defines the top of the downgoing plate, showing that the two plates separate at depth about 40 km closer to the trench than reported in previous reconstructions. This spatial relationship has profound implications for the geometry of the shallow subduction zone and of the mantle wedge, by the amount of crustal material consumed at trench.

High resolution, multi-2D seismic imaging of Solfatara crater (Campi Flegrei Caldera, southern Italy) from active seismic data

A multi-2D imaging of the Solfatara Crater inside the Campi Flegrei Caldera, was obtained by the joint interpretation of geophysical evidences and the new active seismic dataset acquired during the RICEN experiment (EU project MEDSUV) in 2014. We used a total of 17,894 first P-wave arrival times manually picked on pre-processed waveforms, recorded along two 1D profiles criss-crossing the inner Solfatara crater, and performed a tomographic inversion based on a multi-scale strategy and a Bayesian estimation of velocity parameters. The resulting tomographic images provide evidence for a low velocity (500–1500 m/s) water saturated deeper layer at West near the outcropping evidence of the Fangaia, contrasted by a high velocity (2000–3200 m/s) layer correlated with a consolidated tephra deposit. The transition velocity range (1500–2000 m/s) layer suggests a possible presence of a gas-rich, accumulation volume. Thanks to the mutual P-wave velocity model, we infer a detailed image for the gas migration path to the Earth surface. The gasses coming from the deep hydrothermal plume accumulate in the central and most depressed area of the Solfatara being trapped by the meteoric water saturated layer. Therefore, the gasses are transmitted through the buried fault toward the east part of the crater, where the ring faults facilitate the release as confirmed by the fumaroles. Starting from the eastern surface evidence of the gas releasing in the Bocca Grande and Bocca Nuova fumaroles, and the presence of the central deeper plume we suggest a fault situated in the central part of the crater which seems to represent the main buried conduit among them plays a key role.

Towards monitoring the englacial fracture state using virtual-reflector seismology

In seismology, coda wave interferometry (CWI) is an effective tool to monitor time-lapse changes using later arriving, multiply scattered coda waves. Typically, CWI relies on an estimate of the medium's impulse response. The latter is retrieved through simple time-averaging of receiver-receiver cross-correlations of the ambient field, that is, seismic interferometry (SI). In general, the coda is induced by heterogeneities in the Earth. Being comparatively homogeneous, however, ice bodies such as glaciers and ice sheets exhibit little scattering. In addition, the temporal stability of the time-averaged cross-correlations suffers from temporal variations in the distribution and amplitude of the passive seismic sources. Consequently, application of CWI to ice bodies is currently limited. Nevertheless, fracturing and changes in the englacial macroscopic water content alter the bulk elastic properties of ice bodies, which can be monitored with cryoseismological measurements. To overcome the current limited applicability of CWI to ice bodies, we therefore introduce virtual-reflector seismology (VRS). VRS relies on a so-called multidimensional deconvolution (MDD) process of the time-averaged crosscorrelations. The technique results in the retrieval of a medium response that includes virtual reflections from a contour of receivers enclosing the region of interest (i.e. the region to be monitored). The virtual reflections can be interpreted as artificial coda replacing the (lacking) natural scattered coda. Hence, this artificial coda might be exploited for the purpose of CWI. From an implementation point of view, VRS is similar to SI by MDD, which, as its name suggests, also relies on a multidimensional deconvolution process. SI by MDD, however, does not generate additional virtual reflections. Advantageously, both techniques mitigate spurious coda changes associated with temporal variations in the distribution and amplitude of the passive seismic sources. In this work, we apply SI by MDD and VRS to synthetic and active seismic surface-wave data. The active seismic data were acquired on Glacier de la Plaine Morte, Switzerland. We successfully retrieve virtual reflections through the application of VRS to this active seismic data. In application to both synthetic and active seismic data, we show the potential of VRS to monitor time-lapse changes. In addition, we find that SI by MDD allows for a more accurate determination of phase velocity.

Real-time microseismic overburden surveillance at the Grane PRM field offshore Norway

Real-time microseismic detection at offshore hydrocarbon fields is on its way to becoming a standard monitoring tool. Recently, increased focus on injection and overburden surveillance for an improved health, safety and environment (HSE) and cost saving has led to this development. Several hydrocarbon fields are already equipped with permanent reservoir monitoring (PRM) systems with seismic sensors permanently installed at the seafloor (Caldwell et al., 2015), and similar installations are planned or under consideration for some other offshore fields. PRM systems are in principle designed for acquiring active time-lapse seismic data 1-2 times per years, and as such, they are not used during most of their lifetime. But apart from active seismic, PRM systems can also be used for recording passive seismic data. With appropriate processing and analysis methods, such as microseismic event detection, the continuous stream of passive data can be converted into useful real-time subsurface information. This results in improved HSE, and therefore a more valuable PRM system. In 2014, a mini-PRM system with 172 multi-component sensors was installed at the Oseberg field, offshore Norway. The mini-PRM system indeed demonstrated the feasibility of injection and overburden surveillance using real-time passive seismic (Bussat et al., 2016). Despite high installation costs, a cost/benefit evaluation shows net benefits for the Oseberg system owing to better control on waste injection rates. The next step, and topic of the current paper, is to scale up and transfer the learnings from the small Oseberg system to the large Grane PRM system, so as to enable processing and analysis of large amounts of passive data and allow for real-time injection monitoring at Grane. Compared to the Oseberg case, the microseismic data processing is distributed and optimized to be able to use many more sensors and monitor larger subsurface volumes with increased resolution. Moreover, the detection method is improved from an originally strict semblance-based method to also include signal-to-noise (S/N) analysis of the semblance-weighted stack. Crucial noise filtering is integrated into the real-time processing flow, enhancing the sensitivity of the system and ensuring the best possible detection/localization at any time. Key personnel receive an alert immediately after an event detection. This makes the monitoring fully integrated into the follow-up of the daily operation.

Retrieval of Body-Wave Reflections Using Ambient Noise Interferometry Using a Small-Scale Experiment

We report the retrieval of body-wave reflections from noise records using a small-scale experiment over a mature oil field. The reflections are obtained by cross-correlation and stacking of the data. We used the stacked correlograms to create virtual source-to-receiver common shot gathers and are able to obtain body-wave reflections. Surface waves that obliterate the body-waves in our noise correlations were attenuated following a standard procedure from active source seismics. Further different strategies were employed to cross-correlate and stack the data: classical geometrical normalized cross-correlation (CCGN), phase cross-correlation (PCC), linear stacking**** and phase weighted stacking (PWS). PCC and PWS are based on the instantaneous phase coherence of analytic signals. The four approaches are independent and reveal the reflections; nevertheless, the combination of PWS and CCGN provided the best results. Our analysis is based on 2145 cross-correlations of 600 s data segments. We also compare the resulted virtual shot gathers with an active 2D seismic line near the passive experiment. It is shown that our ambient noise analysis reproduces reflections which are present in the active seismic data.

Seismic monitoring leveraging existing telecom infrastructure at the SDASA: Active, passive, and ambient-noise analysis

We analyze active and passive seismic data recorded by the Stanford distributed acoustic sensing array (SDASA) located in conduits under the Stanford University campus. For the active data we used low-energy sources (betsy gun and sledge hammer) and recorded data using both the DAS array and 98 three-component nodes deployed along a 2D line. The joint analysis of shot profiles extracted from the two data sets shows that some surface waves and refracted events are consistently recorded by the DAS array. In areas where geophone coupling was suboptimal because of surface obstructions, DAS recordings are more coherent. In contrast, surface waves are more reliably recorded by the geophones than the DAS array. Because of the noisy environment and weak sources, neither data set shows clear reflections. We demonstrate the repeatability of DAS recordings of local earthquakes by comparing two weak events (magnitude 0.95 and 1.34) with epicenters 100 m apart that occurred only one minute from each other. Analyzing another local, and slightly stronger, earthquake (magnitude 2.0) we show how the kinematics of both the P-arrival and S-arrival can be measured from the DAS data. Interferometric analysis of passive data shows that reliable virtual-source responses can be extracted from the DAS data. We observe Rayleigh waves when correlating aligned receivers, and Love waves when correlating receivers belonging to segments of the array parallel to each other. Dispersion analysis of the virtual sources shows the expected decrease in surface-wave velocity with increasing frequency.

Integrated tomographic methods for seismic imaging and monitoring of volcanic caldera structures and geothermal areas

In this paper we present innovative methodologies for seismic monitoring of volcanic structures in space and time (4D) which can possibly evolve toward an unrest stage. They are based on repeated phase and amplitude measurements done on active and/or passive seismic data including shots, vibrations, earthquakes and ambient noise in order to characterize the structure of the volcano and track its evolution through time. The characterization of the medium properties is performed through the reconstruction of an image of the elastic and anelastic properties of the propagation medium crossed by seismic waves. This study focuses on the application of specific tomographic inversion methods to obtain high quality tomographic images. The resolution of the tomographic models is influenced by the number and spatial distribution of data. The expected resolution thus guides the setup of, for example, active seismic surveys. To recognize and monitor changes in the properties of the propagation medium without performing an active survey we identify a fast proxy based on the time evolution of the Vp/Vs ratio. The advantages and limitations of the methods are discussed through synthetic tests, resolution analysis and case studies in volcanic areas such as the Campi Flegrei (southern Italy) and The Geysers geothermal area (California).

The composition and structure of the deep crust of the Capricorn Orogen

ABSTRACTUnderstanding how the Australian continent came together requires an understanding of structure in all levels of the lithosphere. Deep seismic reflection profiles across several Proterozoic orogens have revealed entirely buried tectonic elements, termed seismic provinces. Although undoubtedly important, the nature of these seismic provinces is typically not well characterised. The Capricorn Orogen is one such region, where the upper crust is relatively well known from geological and geophysical studies, but much of the deep crust is buried beneath Proterozoic basins. Here we combine geophysical datasets, including active and passive source seismic data and gravity data, to image the density, seismic velocity and compositional structure of the deep crust of the Capricorn Orogen. Crustal structure interpreted from deep seismic reflection studies is re-scaled using velocity information from receiver function studies. This modified geometry is used to construct a density model that satisfies Bouguer gravity data. Finally, after correcting for temperature and pressure dependencies, the velocity and density information is used to generate a compositional model of the orogen. This model indicates a varied structure with at least four distinct blocks between the Yilgarn and Pilbara cratons, bounded by major shear zones. We suggest that this variation is linked to multiple accretion events during the amalgamation of the West Australian Craton.

Spatiotemporal changes of seismic attenuation caused by injected CO2 at the Frio‐II pilot site, Dayton, TX, USA

AbstractA continuous active source seismic monitoring data set was collected with crosswell geometry during CO2 injection at the Frio‐II brine pilot, near Liberty, TX. Previous studies have shown that spatiotemporal changes in the P wave first arrival time reveal the movement of the injected CO2 plume in the storage zone. To further constrain the CO2 saturation, particularly at higher saturation levels, we investigate spatial‐temporal changes in the seismic attenuation of the first arrivals. The attenuation changes over the injection period are estimated by the amount of the centroid frequency shift computed by local time‐frequency analysis. We observe that (1) at receivers above the injection zone seismic attenuation does not change in a physical trend; (2) at receivers in the injection zone attenuation sharply increases following injection and peaks at specific points varying with distributed receivers, which is consistent with observations from time delays of first arrivals; then, (3) attenuation decreases over the injection time. The attenuation change exhibits a bell‐shaped pattern during CO2 injection. Under Frio‐II field reservoir conditions, White's patchy saturation model can quantitatively explain both the P wave velocity and attenuation response observed. We have combined the velocity and attenuation change data in a crossplot format that is useful for model‐data comparison and determining patch size. Our analysis suggests that spatial‐temporal attenuation change is not only an indicator of the movement and saturation of CO2 plumes, even at large saturations, but also can quantitatively constrain CO2 plume saturation when used jointly with seismic velocity.

Crust and upper-mantle structure of Wanganui Basin and southern Hikurangi margin, North Island, New Zealand as revealed by active source seismic data

Crust and upper-mantle structure of Wanganui Basin and southern Hikurangi margin, North Island, New Zealand as revealed by active source seismic data

The shallow magma chamber of Stromboli Volcano (Italy)

AbstractIn this work, we integrate artificial and natural seismic sources data to obtain high‐resolution images of the shallow inner structure of Stromboli Volcano. Overall, we used a total of 21,953 P readings from an active seismic experiment and an additional 2731 P and 992 S readings deriving from 269 local events. The well‐defined Vp, Vs, and Vp/Vs tomograms have highlighted the following: (i) the region where magma cumulates at shallow depths (2–4 km below sea level (bsl)), forming an elongated NE‐SW high‐velocity body (Vp ≥ 6.0 km/s and Vs ≥ 3.5 km/s), with a very fast velocity core (6.5 ≤ Vp < 7.0 km/s) of ~2 km3; (ii) the presence of some near‐vertical pipe‐like structures, characterized by relatively high P velocities values, mainly linked to past activity (e.g., Strombolicchio); and (iii) a near‐vertical pipe‐like volume with high Vp/Vs (1.78 ÷ 1.85), located beneath to the craters (down to ~1.0 km bsl), overlying a deeper region (1.0 to 3.0 km bsl) with low Vp/Vs (1.64 ÷ 1.69), interpreted as the actual and preferential pathway of magma toward the surface. Our results demonstrate the importance of combining passive and active seismic data to improve, in a tomographic inversion, the resolution of the volcanic structures and to discover where magma may be stored.

Source characteristics and geological implications of the January 2016 induced earthquake swarm near Crooked Lake, Alberta

On 2016 January 12, an intraplate earthquake with an initial reported local magnitude (ML) of 4.8 shook the town of Fox Creek, Alberta. While there were no reported damages, this earthquake was widely felt by the local residents and suspected to be induced by the nearby hydraulic-fracturing (HF) operations. In this study, we determine the earthquake source parameters using moment tensor inversions, and then detect and locate the associated swarm using a waveform cross-correlation based method. The broad-band seismic recordings from regional arrays suggest a moment magnitude (M) 4.1 for this event, which is the largest in Alberta in the past decade. Similar to other recent M ∼ 3 earthquakes near Fox Creek, the 2016 January 12 earthquake exhibits a dominant strike-slip (strike = 184°) mechanism with limited non-double-couple components (∼22 per cent). This resolved focal mechanism, which is also supported by forward modelling and P-wave first motion analysis, indicates an NE–SW oriented compressional axis consistent with the maximum compressive horizontal stress orientations delineated from borehole breakouts. Further detection analysis on industry-contributed recordings unveils 1108 smaller events within 3 km radius of the epicentre of the main event, showing a close spatial-temporal relation to a nearby HF well. The majority of the detected events are located above the basement, comparable to the injection depth (3.5 km) on the Duvernay shale Formation. The spatial distribution of this earthquake cluster further suggests that (1) the source of the sequence is an N–S-striking fault system and (2) these earthquakes were induced by an HF well close to but different from the well that triggered a previous (January 2015) earthquake swarm. Reactivation of pre-existing, N−S oriented faults analogous to the Pine Creek fault zone, which was reported by earlier studies of active source seismic and aeromagnetic data, are likely responsible for the occurrence of the January 2016 earthquake swarm and other recent events in the Crooked Lake area.

Sedimentary basin exploration with receiver functions: Seismic structure and anisotropy of the Dublin Basin (Ireland)

Teleseismic receiver functions (RFs) were used to investigate the seismic structure of the southern margin of the Dublin Basin, a potential geothermal site. Through an inversion-based approach, the elastic properties and seismic anisotropy of sedimentary basin units were examined, using data from a linear array of closely spaced seismic stations. Our results were compared with sonic logs and lithostratigraphies from two nearby boreholes, NGE1 and NGE2 and colocated active seismic data. Including a high-frequency RF (up to 8 Hz) allowed us to compute S-wave velocity models with a vertical resolution [Formula: see text]. The results indicated the presence of a subvertical lateral discontinuity in [Formula: see text], in correspondence with the main basin-bounding fault (Blackrock-Newcastle Fault [BNF]). North of this discontinuity, a shallow low-velocity layer thickens (from 0.7 to 1.0 km thick) toward the inner basin, in agreement with the geometry of the shallowest reflector found by active seismics. A good correlation was also found between the sonic log at NGE1 and our velocity model. Station DB02 showed an increase in [Formula: see text] at a depth of approximately 0.7 km and a decrease in [Formula: see text] at approximately 1.4 km in depth. Two velocity jumps with matching polarities were also observed in the NGE1 sonic log at the contact between the Upper and Lower Calp formations (positive jump, 688 m deep), and between a calcarenite and a sandstone layers (negative jump, 1337 m deep). Moreover, the main velocity contrasts in our model agree with the major lithostratigraphic boundaries inferred from borehole-drilled samples. Two juxtaposed anisotropic layers are identified close to the BNF. Directions of the slow axis of anisotropy are consistent with the borehole structural data. From these observations, the presence of aligned open cracks within the sandstones, possibly fluid-filled, was inferred up to a depth of 2.3 km in the vicinity of the BNF.

Monitoring the temperature-dependent elastic and anelastic properties in isotropic polycrystalline ice using resonant ultrasound spectroscopy

Abstract. The elastic and anelastic properties of ice are of interest in the study of the dynamics of sea ice, glaciers, and ice sheets. Resonant ultrasound spectroscopy allows quantitative estimates of these properties and aids calibration of active and passive seismic data gathered in the field. The elastic properties and anelastic quality factor Q in laboratory-manufactured polycrystalline isotropic ice cores decrease (reversibly) with increasing temperature, but compressional-wave speed and attenuation prove most sensitive to temperature, indicative of pre-melting of the ice. This method of resonant ultrasound spectroscopy can be deployed in the field, for those situations where shipping samples is difficult (e.g. remote locations), or where the properties of ice change rapidly after extraction (e.g. in the case of sea ice).

The shallow elastic structure of the lunar crust: New insights from seismic wavefield gradient analysis

AbstractEnigmatic lunar seismograms recorded during the Apollo 17 mission in 1972 have so far precluded the identification of shear‐wave arrivals and hence the construction of a comprehensive elastic model of the shallow lunar subsurface. Here, for the first time, we extract shear‐wave information from the Apollo active seismic data using a novel waveform analysis technique based on spatial seismic wavefield gradients. The star‐like recording geometry of the active seismic experiment lends itself surprisingly well to compute spatial wavefield gradients and rotational ground motion as a function of time. These observables, which are new to seismic exploration in general, allowed us to identify shear waves in the complex lunar seismograms, and to derive a new model of seismic compressional and shear‐wave velocities in the shallow lunar crust, critical to understand its lithology and constitution, and its impact on other geophysical investigations of the Moon's deep interior.

PARTOS - Passive and Active Ray TOmography Software: description and preliminary analysis using TOMO-ETNA experiment’s dataset

<p>In this manuscript we present the new friendly seismic tomography software based on joint inversion of active and passive seismic sources called PARTOS (Passive Active Ray TOmography Software). This code has been developed on the base of two well-known widely used tomographic algorithms (LOTOS and ATOM-3D), providing a robust set of algorithms. The dataset used to set and test the program has been provided by TOMO-ETNA experiment. TOMO-ETNA database is a large, high-quality dataset that includes active and passive seismic sources recorded during a period of 4 months in 2014. We performed a series of synthetic tests in order to estimate the resolution and robustness of the solutions. Real data inversion has been carried out using 3 different subsets: i) active data; ii) passive data; and iii) joint dataset. Active database is composed by a total of 16,950 air-gun shots during 1 month and passive database includes 452 local and regional earthquakes recorded during 4 months. This large dataset provides a high ray density within the study region. The combination of active and passive seismic data, together with the high quality of the database, permits to obtain a new tomographic approach of the region under study never done before. An additional user-guide of PARTOS software is provided in order to facilitate the implementation for new users.</p>

Acquisition of a Unique Onshore/Offshore Geophysical and Geochemical Dataset in the Northern Malawi (Nyasa) Rift

The Study of Extension and maGmatism in Malawi aNd Tanzania (SEGMeNT) project acquired a comprehensive suite of geophysical and geochemical datasets across the northern Malawi (Nyasa) rift in the East Africa rift system. Onshore/offshore active and passive seismic data, long‐period and wideband magnetotelluric data, continuous Global Positioning System data, and geochemical samples were acquired between 2012 and 2016. This combination of data is intended to elucidate the sedimentary, crustal, and upper‐mantle architecture of the rift, patterns of active deformation, and the origin and age of rift‐related magmatism. A unique component of our program was the acquisition of seismic data in Lake Malawi, including seismic reflection, onshore/offshore wide‐angle seismic reflection/refraction, and broadband seismic data from lake‐bottom seismometers, a towed streamer, and a large towed air‐gun source.