Stochastic Modeling of the Thermal Structure to Decipher the Lithospheric Thickness: Application to Dharwar Craton

A wide range of speculation regarding lithospheric thermal thickness in a region calls for quantifying the uncertainties associated with model parameters. In this paper, we present an attempt to decipher it by solving the differential equations governing the heat transfer through the lithosphere. Adequately quantifying the temperature field is very important, since it is highly dependent on the controlling thermal parameters. Thermal conductivity is modeled as a random parameter with a known mean value and correlation structure. We compute the temperatures in the crust and lithospheric mantle, along with their error bounds, using an analytical solution to the steady-state heat equation and the Adomian method of decomposition. We apply these solutions in the Archean Dharwar craton, which is divided into western Dharwar craton (WDC) and eastern Dharwar craton (EDC). We obtain a Moho temperature estimate of 304 +/- 58 degrees C in WDC, whilst in EDC it is 375 +/- 74 degrees C for a 20% coefficient of variability in thermal conductivity. The range of lithospheric thickness obtained for WDC is 150-280km, and 110-160km for the EDC. These values are in agreement with other geophysical and geological findings in the region. The thermal lithosphere obtained from this study is based on the mantle solidus and the surface wave seismological studies considering radial anisotropy reveal that the lithosphere is characterized by higher shear wave velocities, and decrease at the lithosphere-asthenosphere boundary. Average shear velocity model for two different cratonic blocks (EDC and WDC) is used to show the LAB thickness for the eastern and western Dharwar craton. The two independent studies show a correlation in the lithospheric thickness.