Search

Teleseismic receiver-functions and Rayleigh-wave dispersion curves are jointly inverted for quantifying $S$-wave velocity profiles beneath the active volcanic zone off Mayotte. We show that the lithosphere in the east-northeast quadrant is composed of four main layers, interpreted as the volcanic edifice, the crust with underplating, the lithospheric mantle, and the asthenosphere, the latter two presenting a main low-velocity zone. The depths of the old (10–11 km) and new Moho (28–31 km) coincide with the two magma reservoirs evidenced by recent seismological and petrological methods. We propose that the main magma reservoir composed of mush with an increasing amount of liquid extends down to 54 km depth. This magma storage develops from a rheological contrast between the ductile lower and brittle upper lithospheric mantle and accounts for most of the volcanic eruption-related seismicity. Finally, the abnormally small thickness of the lithospheric mantle (33 km) is likely a result of a thermal thinning since the onset of Cenozoic magmatism.

We report in this study how ocean-bottom seismometers (OBS) can be used as passive sonars to automatically detect, localize, and track moving acoustic sources at the ocean surface. We developed single-station methods based on direction of arrival and on multi-path interference measurements capable of handling continuous erratic signals emitted by ships. Based on a Bayesian mathematical framework, we developed an azimuthal detector and a radial detector and combined them into a fully automatic tracker. We tested the developed algorithm on seismic and hydroacoustic data recorded in the Indian Ocean by an OBS deployed at 4300 m depth, 200 km west of La Réunion Island. We quantified the performances using archives of commercial-vessel trajectories in the area provided by the Automatic Identification System. Detectors demonstrate capabilities in the detection range up to 100 km from the OBS with azimuthal accuracies of a few degrees and with distance accuracies of a few hundred of meters. We expect the method to be easily transposed to any other kind of sources (such as marine mammals). [ABSTRACT FROM AUTHOR]

Author(s): James Wookey [1]; J.-Michael Kendall (corresponding author) [1]; Guilhem Barruol [2] With time, convective processes in the Earth's mantle will tend to align crystals, grains and inclusions. This mantle [...]

The lithosphere of Madagascar records a long series of tectonic processes. Structures initially inherited from the Pan-African Orogeny are overprinted by a series of extensional tectonic and magmatic events that began with the breakup of Gondwana and continued through to the present. Here, we present a Pn- tomography study in which Pn traveltimes are inverted to investigate the lateral variation of the seismic velocity and anisotropy within the uppermost mantle beneath Madagascar. Results show that the Pn velocities within the uppermost mantle vary by ±0.30 km s–1 about a mean of 8.10 km s–1. Low- Pn -velocity zones (<8.00 km s–1) are observed beneath the Cenozoic alkaline volcanic provinces in the northern and central regions. They correspond to thermally perturbed zones, where temperatures are estimated to be elevated by ∼100–300 K. Moderately low Pn velocities are found near the southern volcanic province and along an E–W belt in central Madagascar. This belt is located at the edge of a broader low S -velocity anomaly in the mantle imaged in a recent surface wave tomographic study. High- Pn -velocity zones (>8.20 km s–1) coincide with stable and less seismically active regions. The pattern of Pn anisotropy is very complex, with small-scale variations in both the amplitude and the fast-axis direction, and generally reflects the complicated tectonic history of Madagascar. Pn anisotropy and shear wave (SKS) splitting measurements show good correlations in the southern parts of Madagascar, indicating coherency in the vertical distribution of lithospheric deformation along Pan-African shear zone as well as coupling between the crust and mantle when the shear zones were active. In most other regions, discrepancies between Pn anisotropy and SKS measurements suggest that the seismic anisotropy in the uppermost mantle beneath Madagascar differs from the vertically integrated upper mantle anisotropy, implying a present-day vertical partitioning of the deformation. Pn anisotropy directions lack the coherent pattern expected for an incipient plate boundary within Madagascar proposed in some kinematic models of the region. [ABSTRACT FROM AUTHOR]

We use seismic noise cross-correlations to obtain a 3-D tomography model of SV -wave velocities beneath the western Indian Ocean, in the depth range of the oceanic crust and uppermost mantle. The study area covers 2000 × 2000 km2 between Madagascar and the three spreading ridges of the Indian Ocean, centred on the volcanic hotspot of La Réunion. We use seismograms from 38 ocean bottom seismometers (OBSs) deployed by the RHUM-RUM project and 10 island stations on La Réunion, Madagascar, Mauritius, Rodrigues, and Tromelin. Phase cross-correlations are calculated for 1119 OBS-to-OBS, land-to-OBS, and land-to-land station pairs, and a phase-weighted stacking algorithm yields robust group velocity measurements in the period range of 3–50 s. We demonstrate that OBS correlations across large interstation distances of >2000 km are of sufficiently high quality for large-scale tomography of ocean basins. Many OBSs yielded similarly good group velocity measurements as land stations. Besides Rayleigh waves, the noise correlations contain a low-velocity wave type propagating at 0.8–1.5 km s−1 over distances exceeding 1000 km, presumably Scholte waves travelling through seafloor sediments. The 100 highest-quality group velocity curves are selected for tomographic inversion at crustal and lithospheric depths. The inversion is executed jointly with a data set of longer-period, Rayleigh-wave phase and group velocity measurements from earthquakes, which had previously yielded a 3-D model of Indian Ocean lithosphere and asthenosphere. Robust resolution tests and plausible structural findings in the upper 30 km validate the use of noise-derived OBS correlations for adding crustal structure to earthquake-derived tomography of the oceanic mantle. Relative to crustal reference model CRUST1.0, our new shear-velocity model tends to enhance both slow and fast anomalies. It reveals slow anomalies at 20 km depth beneath La Réunion, Mauritius, Rodrigues Ridge, Madagascar Rise, and beneath the Central Indian spreading ridge. These structures can clearly be associated with increased crustal thickness and/or volcanic activity. Locally thickened crust beneath La Réunion and Mauritius is probably related to magmatic underplating by the hotspot. In addition, these islands are characterized by a thickened lithosphere that may reflect the depleted, dehydrated mantle regions from which the crustal melts where sourced. Our tomography model is available as electronic supplement. [ABSTRACT FROM AUTHOR]

Over a year of seismic observations, ~5000 short duration icequakes were detected by a permanent broadband station installed at the Princess Elisabeth base, located ~180 km inland in eastern Dronning Maud Land, East-Antarctica. Icequake detection via seismic waveform pattern recognition indicates the presence of two dominating clusters of events, totalizing ~1500 icequakes. The corresponding icequake locations point towards two distinct zones of outcropping blue ice areas (BIAs) located respectively at 4 and 1 km from the seismic station, both on the leeward side of a nunatak protruding through the ice sheet. The temporal occurrence of these icequakes suggests a close genetic link with thermal contraction of ice caused by significant surface cooling controlled, in summer by variations in diurnal solar radiation and in winter by strong cooling during cold katabatic regimes. Further analysis demonstrates the dependence of these icequakes on the absolute surface temperature and on its temporal change. Besides providing information on the ice fracture mechanics and rheology, investigations of thermal icequakes may be regarded as a ground-based proxy for the monitoring of the thermal state of BIAs, and characterization of ice-sheet ablation zones. [ABSTRACT FROM AUTHOR]

We present SKS splitting measurements in the Western Indian Ocean, recorded on 20 land and 57 seafloor seismometers deployed by the RHUM-RUM experiment (Réunion Hotspot and Upper Mantle – Réunions Unterer Mantel). We discuss our splitting observations within their geodynamic settings and compare them to SKS splitting parameters predicted from an azimuthally anisotropic Rayleigh wave tomography model that includes the RHUM-RUM data. We find that anisotropic directions poorly correlate with the present-day motion of the Somali plate, which at <2.6 cm/yr may be too slow to cause strongly sheared fabric in the asthenosphere. Fast split directions ( Φ ) between La Réunion and the Central Indian Ridge (CIR) trend E–W and provide strong, first seismological evidence for near-horizontal flow in the asthenosphere that connects the Réunion mantle upwelling with the CIR, supporting a long-standing hypothesis on plume–ridge interaction. In the vicinity of the Réunion hotspot, we observe a seismic anisotropy pattern indicative of a parabolic asthenospheric flow controlled by the Réunion mantle upwelling and its consecutive asthenospheric spreading. We furthermore observe ridge-normal Φ along the CIR and ridge-parallel Φ along the Southwest Indian Ridge (SWIR), both mainly attributed to asthenospheric mantle flows. In the Mozambique Channel between East-Africa and Madagascar, we attribute E–W trending Φ to frozen lithospheric structures, recording the paleo-orientation of the spreading ridges that enabled Madagascar's separation away from Africa. Based on the synopsis of this and previous SKS splitting studies at mid-ocean ridges, we propose that ridge-normal Φ may develop at fast and intermediate spreading ridges (e.g., CIR and East Pacific Rise) and ridge-parallel Φ could be characteristic to slow spreading ridges (e.g., SWIR, Mid-Atlantic Ridge and the paleo-ridges in the Mozambique Channel). [ABSTRACT FROM AUTHOR]

Many applications in seismology rely on the accurate absolute timing of seismograms. However, both seismological land stations and ocean bottom seismometers (OBSs) can be affected by clock errors, which cause the absolute timing of seismograms to deviate from a highly accurate reference time signal, usually provided by GPS satellites. Timing problems can occur in land stations when synchronization with a GPS signal is temporarily or permanently lost. This can give rise to complicated, time-dependent clock drifts relative to GPS time, due to varying environmental conditions. Seismometers at the ocean bottom cannot receive GPS satellite signals, but operate in more stable ambient conditions than land stations. The standard protocol is to synchronize an OBS with a GPS signal immediately before deployment and after recovery. The measured timing deviation, called ‘skew’, is assumed to have accumulated linearly over the deployment interval, an assumption that is plausible but usually not verifiable. In recent years, cross-correlations of ambient microseismic noise have been put to use for correcting timing errors, but have been limited to interstation distances of at most a few tens of kilometres without reducing the temporal resolution. We apply noise cross-correlations to the evaluation of clock errors in four broad-band land stations and 53 wideband and broad-band OBSs, which were installed on and around the island of La Réunion in the western Indian Ocean during the RHUM-RUM (Réunion Hotspot and Upper Mantle-Réunions Unterer Mantel) experiment. We correlate all three seismic components, plus a hydrophone channel in OBS stations. Daily cross-correlation functions are derived for intermediate distances (∼20 km) for land-to-land station pairs; stable, 10 d stacks are obtained for very large interstation distances up to >300 km for land-to-OBS and OBS-to-OBS configurations. Averaging over multiple station pairs, and up to 16 component pairs per station, improves the accuracy of the method by a factor of four compared to the single-channel approaches of prior studies. The timing accuracy of our method is estimated to be ∼20 ms standard deviation or one sample at a sampling rate of 50 Hz. In land stations, nonlinear clock drifts and clock jumps of up to 6 min are detected and successfully corrected. For 52 out of 53 OBSs, we successfully obtain drift functions over time, which validate the common assumption of linear clock drift. Skew values that were available for 29 of these OBSs are consistent with our independent estimates within their observational error bars. For 23 OBSs that lacked skew measurements, linear OBS clock drifts range between 0.2 and 8.8 ms d−1. In addition to linear drift, three OBSs are affected by clock jumps of ∼1 s, probably indicating a missing sample problem that would otherwise have gone undetected. Thus we demonstrate the routine feasibility of high-accuracy clock corrections in land and OBSs over a wide range of interstation distances. [ABSTRACT FROM AUTHOR]

As a first step to Antarctic blue whale (ABW) monitoring using passive acoustics, a method based on the stochastic matched filter (SMF) is proposed. Derived from the matched filter (MF), this filter-based denoising method enhances stochastic signals embedded in an additive colored noise by maximizing its output signal to noise ratio (SNR). These assumptions are well adapted to the passive detection of ABW calls where emitted signals are modified by the unknown impulse response of the propagation channel. A filter bank is computed and stored offline based on a priori knowledge of the signal second order statistics and simulated colored sea-noise. Then, the detection relies on online background noise and SNR estimation, realized using time-frequency analysis. The SMF output is cross-correlated with the signal's reference (SMF + MF). Its performances are assessed on an ccean bottom seismometer-recorded ground truth dataset of 845 ABW calls, where the location of the whale is known. This dataset provides great SNR variations in diverse soundscapes. The SMF + MF performances are compared to the commonly used MF and to the Z-detector (a sub-space detector for ABW calls). Mostly, the benefits of the use of the SMF + MF are revealed on low signal to noise observations: in comparison to the MF with identical detection threshold, the false alarm rate drastically decreases while the detection rate stays high. Compared to the Z-detector, it allows the extension of the detection range of ≃ 30 km in presence of ship noise with equivalent false discovery rate. [ABSTRACT FROM AUTHOR]

In the western Indian Ocean, the Réunion hot spot is one of the most active volcanoes on Earth. Temporal interactions between ridges and plumes have shaped the structure of the zone. This study investigates the mantle structure using data from the Réunion Hotspot and Upper Mantle-Réunions Unterer Mantel (RHUM-RUM) project, which significantly increased the seismic coverage of the western part of the Indian Ocean. For more than 1 year, 57 ocean bottom seismometer stations and 23 temporary land stations were deployed over this area. For each earthquake station path, the Rayleigh wave fundamental mode phase velocities were measured for the periods from 30 s to 300 s and group velocities for the period from 16 s to 250 s. A three-dimensional model of the shear velocity in the upper mantle was built in two steps. First, the path mean phase and group velocities were inverted, to obtain regionalized velocity maps for each separate period. Then, all of the phase and group velocity maps were combined and inverted at each grid point, to obtain the local S wave velocity as a function of depth, using a transdimensional inversion scheme. The three-dimensional anisotropic S wave velocity model has resolution down to 300 km in depth. The tomographic model and surface tectonics are correlated down to ∼100 km in depth. Starting at 50 km in depth, a slow velocity anomaly beneath Rodrigues Ridge and the east-west orientation of the azimuthal anisotropy show a connection between Réunion upwelling and the Central Indian Ridge. The slow velocity signature beneath La Réunion is connected at greater depths (150-300 km) with a large slow velocity zone beneath the entire Mascarene Basin. This develops along a northeast direction, following the general motion direction of the African Plate. These observations indicate nonisotropic spreading of hot plume material and dominant horizontal flow in the upper mantle beneath this area. [ABSTRACT FROM AUTHOR]

We present two independent, automated methods for estimating the absolute horizontal misorientation of seismic sensors. We apply both methods to 44 free-fall ocean-bottom seismometers (OBSs) of the RHUM-RUM experiment (http://www.rhum-rum.net/). The techniques measure the 3-D directions of particle motion of (1) P-waves and (2) Rayleigh waves of earthquake recordings. For P-waves, we used a principal component analysis to determine the directions of particle motions (polarizations) in multiple frequency passbands. We correct for polarization deviations due to seismic anisotropy and dipping discontinuities using a simple fit equation, which yields significantly more accurate OBS orientations. For Rayleigh waves, we evaluated the degree of elliptical polarization in the vertical plane in the time and frequency domain. The results obtained for the RHUM-RUM OBS stations differed, on average, by 3.1° and 3.7° between the methods, using circular mean and median statistics, which is within the methods' estimate uncertainties. Using P-waves, we obtained orientation estimates for 31 ocean-bottom seismometers with an average uncertainty (95 per cent confidence interval) of 11° per station. For 7 of these OBS, data coverage was sufficient to correct polarization measurements for underlying seismic anisotropy and dipping discontinuities, improving their average orientation uncertainty from 11° to 6° per station. Using Rayleigh waves, we obtained misorientation estimates for 40 OBS, with an average uncertainty of 16° per station. The good agreement of results obtained using the two methods indicates that they should also be useful for detecting misorientations of terrestrial seismic stations. [ABSTRACT FROM AUTHOR]

We deployed five broadband three-components seismic stations in the Iles Eparses in the south-west Indian Ocean and on Mayotte Island, between April 2011 and January 2014. These small and remote oceanic islands suffer the effects of strong ocean swells that affect their coastal environments but most islands are not instrumented by wave gauges to characterize the swells. However, wave action on the coast causes high levels of ground vibrations in the solid earth, so-called microseismic noise. We use this link between the solid earth and ocean wave activity to quantify the swells locally. Spectral analyses of the continuous seismic data show clear peaks in the 0.05–0.10 Hz frequency band (periods between 10 and 20 s), corresponding to the ocean wave periods of the local swells. We analyze an example of austral swell occurring in August 2013 and a cyclonic event (Felleng) that developed in January 2013, and quantify the ground motion at each station induced by these events. In both cases, we find a linear polarization in the horizontal plane with microseismic amplitude directly correlated to the swell height (as predicted by the global swell model WaveWatchIII), and a direction of polarization close to the predicted swell propagation direction. Although this analysis has not been performed in real time, it demonstrates that terrestrial seismic stations can be efficiently used as wave gauges, and are particularly well suited for quantifying extreme swell events. This approach may therefore provide useful and cheaper alternatives to wave buoys for monitoring swells and the related environmental processes such as beach erosion or coral reef damages. [ABSTRACT FROM AUTHOR]

The Piton de la Fournaise basaltic volcano, on La Réunion Island in the western Indian Ocean, is one of the most active volcanoes in the world. This volcano is classically considered as the surface expression of an upwelling mantle plume and its activity is continuously monitored, providing detailed information on its superficial dynamics and on the edifice structure. Deeper crustal and upper mantle structure under La Réunion Island is surprisingly poorly constrained, motivating this study. We used receiver function techniques to determine a shear wave velocity profile through the crust and uppermost mantle beneath La Réunion, but also at other seismic stations located on the hotspot track, to investigate the plume and lithosphere interaction and its evolution through time. Receiver functions (RFs) were computed at permanent broadband seismic stations from the GEOSCOPE network (on La Réunion and Rodrigues), at IRIS stations MRIV and DGAR installed on Mauritius and Diego Garcia islands, and at the GEOFON stations KAAM and HMDM on the Maldives. We performed non-linear inversions of RFs through modelling of P-to-S conversions at various crustal and upper mantle interfaces. Joint inversion of RF and surface wave dispersion data suggests a much deeper Mohorovičić discontinuity (Moho) beneath Mauritius (∼21 km) compared to La Réunion (∼12 km). A magmatic underplated body may be present under La Réunion as a thin layer (≤3 km thick), as suggested by a previous seismic refraction study, and as a much thicker layer beneath other stations located on the hotspot track, suggesting that underplating is an important process resulting from the plume-lithosphere interaction. We find evidence for a strikingly lowvelocity layer starting at about 33 km depth beneath La Réunion that we interpret as a zone of partial melt beneath the active volcano. We finally observe low velocities below 70 km beneath La Réunion and below 50 km beneath Mauritius that could represent the base of the oceanic lithosphere. [ABSTRACT FROM AUTHOR]

The Terre Adélie and George V Land (East Antarctica) represent key areas for understanding tectonic relationships between terranes forming the Neoarchean-Palaeoproterozoic Terre Adélie Craton (TAC) and the neighbouring lithospheric blocks, together with the nature of its boundary. This region that represents the eastern border of the TAC is limited on its eastern side by the Mertz shear zone (MSZ) separating more recent Palaeozoic units from the craton. The MSZ, that recorded dextral strike-slip movement at 1.7 and 1.5 Ga, is likely correlated with the Kalinjala or Coorong shear zone in South Australia, east of the Gawler Craton and may therefore represent a frozen lithospheric-scale structure. In order to investigate the lithospheric structure of the TAC and the MSZ, we deployed from 2009 October to 2011 October four temporary seismic stations, which sampled the various lithospheric units of the TAC and of the neighbouring Palaeozoic block, together with the MSZ. We used receiver function method to deduce Moho depths and seismic anisotropy technique to infer the upper mantle deformation. Results from receiver functions analysis reveal Moho at 40-44 km depth beneath the TAC, at 36 km under the MSZ and at 28 km beneath the eastern Palaeozoic domain. The MSZ therefore delimits two crustal blocks of different thicknesses with a vertical offset of the Moho of 12 km. Seismic anisotropy deduced from SKS splitting at stations on the TAC shows fast polarisation directions (Φ) trending E-W, that is, parallel to the continental margin, and delay times (δt) ranging from 0.8 to 1.6 s. These results are similar to the splitting parameters observed at the permanent GEOSCOPE Dumont D'Urville station (DRV: Φ 95°N, δt 1.1 s) located in the Palaeoproterozoic domain of TAC. On the MSZ, the small number of good quality measurements limits the investigation of the deep signature of the shear zone. However, the station in the Palaeozoic domain shows Φ trending N60°E, which is significantly different to the Φ trending measurements from stations on the TAC, suggesting that the MSZ may also represent a major frontier between the Neoarchean-Palaeoproterozoic and Palaeozoic terranes. [ABSTRACT FROM AUTHOR]

Author(s): James Wookey (corresponding author) [1]; J.-Michael Kendall [1]; Guilhem Barruol [2] Wookey et al . reply Saul and Vinnik's comment raises an interesting point concerning the influence of shear-coupled [...]

SUMMARY We performed forward modelling of seismic anisotropy beneath a migrating strike-slip plate boundary to: (1) test if such a geodynamic context might explain teleseismic shear wave splitting data in the vicinity of the central part of the San Andreas Fault and (2) constrain the power of such data to unravel vertical and lateral variations in deformation patterns in the upper mantle. The modelling involves five steps: (1) thermomechanical modelling, using a finite-element code, of the deformation field, (2) viscoplastic self-consistent modelling of the resulting olivine and pyroxene crystal preferred orientations, (3) calculation of the elasticity tensors for different domains of the finite elements (FEs) model, (4) forward modelling of seismic wave propagation through the model using ray theory, finite-frequency theory and a full-wave approach and (5) performing splitting measurements on the synthetic seismograms. SKS splitting data in central California are best fitted by a model with a hotter geotherm within 60 km of the plate boundary accounting for the opening of an asthenospheric window due to the northward migration of the Mendocino Triple Junction. The westward motion of the plate boundary cannot however explain the rotation of fast polarization directions east of the San Andreas Fault in central California. Comparison between modelled and measured individual shear wave splitting also implies that the homogeneity of the two-layer models accounting for the observations in the vicinity of the San Andreas Fault indicates a sharp transition between lithospheric and asthenospheric deformations beneath this plate boundary. The ability of different synthetic approaches to localize horizontally and vertically the plate boundary-related deformation differs significantly. Splitting in ray theory synthetics closely follow variations in olivine crystallographic preferred orientations in the model. In contrast, splitting analysis on full-wave synthetics, which should be more representative of actual long-period SKS waves, results in smooth apparent lateral variations of anisotropy; the exact location and width of the plate boundary may only be retrieved by comparing fast-polarization profiles obtained using a multichannel analysis on waves with different periods. [ABSTRACT FROM AUTHOR]

Upper-mantle flow beneath the South Pacific is investigated by analysing shear wave splitting parameters at eight permanent long-period and broad-band seismic stations and 10 broad-band stations deployed in French Polynesia from 2001 to 2005 in the framework of the Polynesian Lithosphere and Upper Mantle Experiment (PLUME). Despite the small number of events and the rather poor backazimuthal coverage due to the geographical distribution of the natural seismicity, upper-mantle seismic anisotropy has been detected at all stations except at Tahiti where two permanent stations with 15 yr of data show an apparent isotropy. The median value of fast polarization azimuths (N67.5°W) is parallel to the present Pacific absolute plate motion direction in French Polynesia (APM: N67°W). This suggests that the observed SKS fast polarization directions result mainly from olivine crystal preferred orientations produced by deformation in the sublithospheric mantle due to viscous entrainment by the moving Pacific Plate and preserved in the lithosphere as the plate cools. However, analysis of individual measurements highlights variations of splitting parameters with event backazimuth that imply an actual upper-mantle structure more complex than a single anisotropic layer with horizontal fast axis. A forward approach shows that a two-layer structure of anisotropy beneath French Polynesia better explains the splitting observations than a single anisotropic layer. Second-order variations in the measurements may also indicate the presence of small-scale lateral heterogeneities. The influence of plumes or fracture zones within the studied area does not appear to dominate the large-scale anisotropy pattern but may explain these second-order splitting variations across the network. [ABSTRACT FROM AUTHOR]

We present a regional surface waveform tomography of the Pacific upper mantle, obtained using an automated multimode surface waveform inversion technique on fundamental and higher mode Rayleigh waves, to constrain the structure down to ∼400 km depth. We have improved on previous implementations of this technique by robustly accounting for the effects of uncertainties in earthquake source parameters in the tomographic inversion. We have furthermore improved path coverage in the South Pacific region by including Rayleigh wave observations from the French Polynesian Pacific Lithosphere and Upper Mantle Experiment deployment. This improvement has led to imaging of vertical low-velocity structures associated with hotspots within the South Pacific Super-Swell region. We have produced an age-dependent average cross-section for the Pacific Ocean lithosphere and found that the increase in with age is broadly compatible with a half-space cooling model of oceanic lithosphere formation. We cannot confirm evidence for a Pacific-wide reheating event. Our synthetic tests show that detailed interpretation of average trends across the Pacific Ocean may be misleading unless lateral resolution and amplitude recovery are uniform across the region, a condition that is difficult to achieve in such a large oceanic basin with current seismic stations. [ABSTRACT FROM AUTHOR]