Deep CO2 emissions facilitated by localized shear deformation: A case study of the Karakoram fault system, western Tibet

•Along- and across-fault soil CO2 fluxes were measured in the Karakoram fault system.•Mass balance model reveals high fraction of deep carbon in the geothermal fields.•Higher CO2 fluxes are found in geothermal fields along the fault bending segment.•Localized shear deformation in the fault bend could facilitate deep CO2 emissions. Strike-slip faults play a significant role in creating deeply penetrating fractures with high permeability, thus promoting rapid migration of CO2-rich fluids to the surface. However, there are rare observations regarding how strike-slip movement could affect deep CO2 emissions. Here, we focus on the Karakoram fault system (KKFS), western Tibet, to estimate diffuse soil CO2 fluxes and to unravel potential controlling factors for CO2 emissions. Average CO2 fluxes of geothermal fields range in 22–2475 g m−2 d−1, significantly higher than the across-fault profiles (6–116 g m−2 d−1). A mass balance model based on δ13C-CO2 and CO2 concentration of soil gases reveals that deep carbon constitutes 49.1–91.5 % (average = 73.9 %) and 0.2–40.5 % (average = 25.5 %) of soil-gas carbon released from geothermal fields and across-fault profiles, respectively. Deep carbon could be produced by thermal decomposition of crustal rocks considering CO2-rich fluids with radiogenic helium isotopes. Strikingly, higher CO2 fluxes preferentially occur in geothermal fields along a bending segment of the KKFS, where localized shear deformation is prominent as documented by high slip rates over geological timescales, dense splay faults, clustering of earthquake events, and elevated strain rates. We suggest that high stress acting on the KKFS bend could enhance the deformation and fracturing of fault zone rocks, leading to production of metamorphic CO2 and efficient release of CO2-rich fluids through the highly permeable fault system. Our results could shed new light on CO2 origins and fluxes of strike-slip faults that are characterized by spatially heterogeneous strain partitioning and thus localized enhanced shear deformation.