Crustal structure beneath broad-band seismic stations in the Mediterranean region

We have analysed receiver functions to derive simple models for crustal structure below 12 broad-band seismological stations from the MIDSEA project and 5 permanent broad-band stations in the Mediterranean region including northern Africa. To determine an accurate Moho depth we have reduced the trade-off between crustal velocities and discontinuity depth using a new grid search method, which is an extension of recently published methods to determine crustal thickness. In this method the best fitting synthetic receiver function, containing both the direct conversion and the reverberated phases, is identified on a model grid of varying Moho depth and varying Poisson's ratio. The values we found for Moho depth range from around 20 km for intra-oceanic islands and extended continental margins to near 45 km in regions where the Eurasian and African continents have collided. More detailed waveform modelling shows that all receiver functions can be well fit using a 2- or 3-layer model containing a sedimentary layer and/or a mid-crustal discontinuity. On comparing our results with Moho maps inferred from interpolated reflection and refraction data, we find that for some regions the agreement between our receiver function analysis and existing Moho maps is very good, while for other regions our observations deviate from the interpolated map values and extend beyond the geographic bounds of these maps.