Receiver function imaging in strongly laterally heterogeneous crust: Synthetic modeling of BOLIVAR data

We resolve a large (~20 km) discrepancy in Moho depth determined from PdS receiver functions (RFs) and from active source seismic profiling in the complex Caribbean-South American plate boundary zone in eastern Venezuela. As part of the BOLIVAR experiment 20 broadband stations were deployed along an active source profile to record teleseisms. Using the extremely heterogeneous crustal model obtained from active source data, we generated 2D finite-difference elastic wave synthetics and from them calculated receiver functions and CCP stacks. We compare the observations with synthetic sections that have been spatially sampled at 0.25 km to 40 km. The densely sampled synthetics show that several events in the field data that were originally interpreted as the Moho are multiple reflections within sedimentary basins. Where the Moho has the steepest dip under the plate boundary the CCP stacks fail to image the Moho well, regardless of the density of spatial sampling. A suitable spatial sampling criterion for clearly imaging the lower crust and Moho is to overlap Fresnel zones by 50% at Moho depth, which for the 1 Hz receiver functions examined here, requires an instrument spacing of 15–20 km, with the actual field data density ranging from 20 km to 100 km.