Estimating the Lithospheric Flexure of a Plate with Non-uniform Flexural Rigidity: A Quantitative Modelling Approach

Abstract

Lithospheric deformation is a fundamental process in plate tectonics. It is therefore critical to determine how the lithosphere responds to geological loads to better understand tectonic processes. The lithosphere can be modelled as the flexure of a thin, elastic plate over long-term (>10 5 yr) geological timescales. The partial differential equation for the flexure of an orthotropic plate is used indirectly to calculate theoretical admittance and coherence, which are then compared against the observed admittance and coherence to invert for the non-uniform flexural rigidity (or effective elastic thickness, T e ) of the plate. However, the process for accurately recovering variable lithospheric flexure remains unresolved, as the classical lithospheric model may overestimate the deflection of the plate. Here we adopt the classic lithospheric model with applied external and internal loads at the surface and Moho, respectively, and assume that the compensation material is denser than the mantle material beneath the Moho. The lithospheric flexure errors are derived mainly from the T e and Moho recovery errors in this lithospheric model. Synthetic modelling is then performed to analyse the influence of the T e and Moho errors. The analysis of synthetic modelling shows that: 1) the inverted flexure values are insensitive to the T e errors, and the T e error-induced flexure errors exhibit a rippling pattern; 2) the Moho error-induced flexure errors occur mainly in the low T e regions, and the Moho disturbance is alleviated by both an attenuation coefficient of the form and the lithospheric strength as the Moho disturbance propagates into the lithospheric flexure error; and 3) the lithospheric flexure errors are dominated by the T e and Moho errors in the high and low T e regions, respectively.