THREEDIMENSIONAL IMAGE OF CRUSTAL DENSITY MODEL: A CASE STUDY IN SOUTH KAMCHATKA
M.D. Sidorov, A.G. Nurmukhamedov
Geotechnological Scientific Research Center, Far Eastern Branch of the Russian Academy of Sciences, SeveroVostochnoe Shosse 30, PO box 56, PetropavlovskKamchatsky, 683002, Russia
Keywords: Density, model, voxel, intrusion, crust, block diagram, isosurface
Abstract
The methodology for building a 3D density model of the Earth’s crust and upper mantle for a crustal block in the South Kamchatka Peninsula is presented. The model is derived from the interpretation of the Bouguer anomaly values and data obtained by deep seismic soundings and electrical prospecting methods. 3D image creation can be viewed as a several phase procedure. At the first step, the area of interest is investigated through a network of intersecting profiles utilized for subsequent modeling depth density sections using the parameter scanning procedure. The a priori data include materials of geophysical surveys and geological, and petrophysical studies as well as results of depth estimation at singular points by the Euler deconvolution method. The next step involves formation of a 3D density matrix from the crosssection models. Coordinates in the crosssection planes are defined in the square grid nodes and converted to plane rectangular coordinates (the SK42 reference system). A file of the crosssectional database of density values is thus produced. For area shaving a sparse network of crosssections, an additional database is derived from grid files of densities mapped horizonwise. This ultimately results in the creation of the unified databases. A 3D network of cubic unit cells for threedimensional (voxel) representation of density of the geologic medium is calculated. The model of density of unit cube with edge measuring 4 km was constructed for entire investigated block of the Earth’s crust. The model was detailed by decreasing successively the unit cube edges to 2.0, 1.0, and 0.5 km. The 3D modeling results are presented as block diagrams and density isosurface maps in horizontal slices at different depths. The block diagrams show isosurfaces with densities typical of deep subsurface boundaries and structures. The isosurfaces for densities of 2.75 and 2.90 g/cm^{3 }reflect the morphology of the roof of the crystalline basement, and the lowerupper crust boundary, respectively, while the isosurfaces of 3.20 and 3.33 g/cm^{3 }represent feasible shape of the crustupper mantle contact. The lowerdensity zone ubiquitously observed in the investigated crustal block is caused by rock destruction along the deep fault. It is identified in the model with 1 km cubic unit cell side and encompasses the feeding systems of volcanoes, namely: Sopka Asacha (Asacha Volcano), Mutnovskaya Sopka (Mutnovsky Volcano), Gorelaya Sopka (Gorely Volcano), and Vilyuchinskaya Sopka (Vilyuchansky Volcano).The model with a 0.50 km cubic unit cell side shows a 2.85 g/cm^{3 }isosurface, which delineates a columnshaped highdensity region. A higher density revealed in the block is assumed to be associated with basite intrusions. 3D density models based on crosssection modeling provide additional information for geodynamic reconstructions, allowing a rough estimation of the geometry of concealed geologic features and their volumes, and can be utilized in 3D modeling to build initial models.
