Identification of sediment–basement structure in West Papua province, Indonesia, using gravity and magnetic data inversion as an Earth’s crust stress indicator

West Papua province in eastern Indonesia is positioned in a dynamic tectonic zone along with the collision of the Australian, Pacific, and Eurasian plates. The interaction resulted in the formation of strike-slip faults such as Koor, Sorong, Ransiki, and Yapen, that are prone to earthquakes in the region. The rocks of West Papua in the northern part are a contribution to the Pacific Ocean plate consisting of ophiolite and volcanic arcs of the archipelago, even while the rocks of the Australian plate in the southern part are dominated by quaternary and siliciclastic sedimentary. It has a wide variety of resources, including oil and gas. This study combines the interpretation of regional gravity and magnetic data derived from satellite observations to identify the subsurface structure of West Papua. This is performed since most studies were conducted on the surface and did not significantly focus on the subsurface. The composition of subsurface is determined through three-dimensional (3-D) unconstrained inversion modeling using the iterative reweighting inversion of regional gravity and magnetic anomalies as a function of density contrast and magnetic susceptibility of rocks. In depth variations, gravity inversion produces density contrast ranging from − 0.348 to 0.451 gr/cm3, whereas magnetic inversion provides rock susceptibility varying between − 0.363 and 0.223 SI. Gravity and magnetic inversions characterize the subduction of the Pacific Ocean plate in the north, extensive intrusion of igneous rocks, and low density-susceptibility contrast in the Bintuni basin as a source of oil and natural gas. The boundary between the sediment layer and the basement is believed to be 15–20 km deep, with rocks from the uplifted mantle in the north and a Silur-Devon aged Kemum formation in the south.