Numerical modeling of the conductive heat transfer in western Anatolia
Using the available data, the crustal temperatures were modeled along two profiles in the western Anatolia. Modeling was achieved by solving a two-dimensional steady-state heat conduction equation in a heterogeneous medium by successive overrelaxation method of finite differences. The Bouguer gravity data were used to construct the crustal structures along the profiles. The average crustal thickness was taken as 30 km. An exponential decrease with depth in heat production was assumed. The temperatures in the grabens and regions under them are obtained as relatively higher than those in the surrounding regions. The thickness of the seismogenic zone and depth to the Curie isotherm in the region were estimated as 10 and 15 km respectively. Also, the temperature at the bottom of the crust was calculated as 1075-1100ºC approximately. Sedimentary fills in the grabens affect the temperature distribution in the upper crust.