Andrew Frederiksen

Surface wave tomography applied to the north american upper mantle

Tomographic techniques that invert seismic surface waves for 3-D Earth structure differ in their definitions of data and the forward problem as well as in the parameterization of the tomographic model. However, all such techniques have in common that the tomographic inverse problem involves solving a large and mixed-determined set of linear equations. Consequently these inverse problems have multiple solutions and inherently undefinable accuracy. Smoother and rougher tomographic models are found with rougher (confined to great circle path) and smoother (finite-width) sensitivity kernels, respectively. A powerful, well-tested method of surface wave tomography (Partitioned Waveform Inversion) is based on inverting the waveforms of wave trains comprising regional S and surface waves from at least hundreds of seismograms for 3-D variations in S wave velocity. We apply this method to nearly 1400 seismograms recorded by digital broadband seismic stations in North America. The new 3-D S-velocity model, NA04, is consistent with previous findings that are based on separate, overlapping data sets. The merging of US and Canadian data sets, adding Canadian recordings of Mexican earthquakes, and combining fundamental-mode with higher-mode waveforms provides superior resolution, in particular in the US-Canada border region and the deep upper mantle.

Crustal thickness variations and isostatic disequilibrium across the North Anatolian Fault, western Turkey

We use teleseismic recordings from a dense array of seismometers straddling both strands of the North Anatolian Fault Zone to determine crustal thickness, P/S velocity ratio and sedimentary layer thickness. To do this, we implement a new grid search inversion scheme based on the use of transfer functions, removing the need for deconvolution for source normalization and therefore eliminating common problems associated with crustal-scale receiver function analysis. We achieve a good fit to the data except at several stations located in Quaternary sedimentary basins, where our two-layer crustal model is likely to be inaccurate. We find two zones of thick sedimentary material: one north of the northern fault branch, and one straddling the southern branch. The crustal thickness increases sharply north of the northern strand of the North Anatolian Fault Zone (NAFZ), where the fault nearly coincides with the trace of the Intra-Pontide Suture; the velocity ratio changes across the southern fault strand, indicating a change in basement composition. We interpret these changes to indicate that both strands of the NAFZ follow preexisting geological boundaries rather than being ideally aligned with the stress field. The thick crust north of the northern NAFZ strand is associated with low topography and so is inconsistent with simple models of isostatic equilibrium, requiring a contribution from mantle density variations, such as possible loading from underthrust Black Sea oceanic lithosphere.

Distinct crustal structure of the North American Midcontinent Rift from P wave receiver functions

Eighty-two broadband seismic stations of the Superior Province Rifting Earthscope Experiment (SPREE) collected 2.5 years of continuous seismic data in the area of the high gravity anomaly associated with the Midcontinent Rift (MCR). Over 100 high-quality teleseismic earthquakes were used for crustal P wave receiver function analysis. Our analysis reveals that the base of the sedimentary layer is shallow outside the MCR, thickens near the flanks where gravity anomalies are low, and shallows again in the MCRs center where the gravity anomalies peak. This pattern is similar to that found from local geophysical studies and is consistent with reverse faulting having accompanied the cessation of rifting at 1.1 Ga. Intermittent intracrustal boundaries imaged by our analysis might represent the bottom of the MCRs mostly buried dense volcanic layers. Outside the MCR, the Moho is strong, sharp, and relatively flat, both beneath the Archean Superior Province and the Proterozoic terranes to its south. Inside the MCR, two weaker candidate Mohos are found at depths up to 25 km apart in the rifts center. The intermediate layer between these discontinuities tapers toward the edges of the MCR. The presence of this transitional layer is remarkably consistent along the strike of the MCR, including beneath its jog in southern Minnesota, near the Belle Plaine Fault. We interpret these results as evidence for extensive underplating as a defining characteristic of the rift, which remains continuous along the Minnesota jog, where due to its orientation, it is minimally affected by the reverse faulting that characterizes the NNE striking parts of the rift.