Identifying processes within the Earth System that have modulated atmospheric pCO₂ during each glacial cycle of the late Pleistocene stands as one of the grand challenges in climate science. The growing array of surface ocean pH estimates from the boron isotope proxy across the last glacial termination may reveal regions of the ocean that influenced the timing and magnitude of pCO₂ rise. Here we present two new boron isotope records from the subtropical‐subpolar transition zone of the Southwest Pacific that span the last 20 kyr, as well as new radiocarbon data from the same cores. The new data suggest this region was a source of carbon to the atmosphere rather than a moderate sink as it is today. Significantly higher outgassing is observed between ~16.5 and 14 kyr BP, associated with increasing δ¹³C and [CO₃]²⁻ at depth, suggesting loss of carbon from the intermediate ocean to the atmosphere. We use these new boron isotope records together with existing records to build a composite pH/pCO₂ curve for the surface oceans. The pH disequilibrium/CO₂ outgassing was widespread throughout the last deglaciation, likely explained by upwelling of CO₂ from the deep/intermediate ocean. During the Holocene, a smaller outgassing peak is observed at a time of relatively stable atmospheric CO₂, which may be explained by regrowth of the terrestrial biosphere countering ocean CO₂ release. Our stack is likely biased toward upwelling/CO₂ source regions. Nevertheless, the composite pCO₂ curve provides robust evidence that various parts of the ocean were releasing CO₂ to the atmosphere over the last 25 kyr.

Seafloor pockmarks of varying size occur over an area of 50,000 km² on the Chatham Rise, Canterbury Shelf and Inner Bounty Trough, New Zealand. The pockmarks are concentrated above the flat‐subducted Hikurangi Plateau. Echosounder data identify recurrent episodes of pockmark formation at ~100,000‐year frequency coinciding with Pleistocene glacial terminations. Here we show that there are structural conduits beneath the larger pockmarks through which fluids flowed upward toward the seafloor. Large negative Δ¹⁴C excursions are documented in marine sediments deposited next to these subseafloor conduits and pockmarks at the last glacial termination. Modern pore waters contain no methane, and there is no negative δ¹³C excursion at the glacial termination that would be indicative of methane or mantle‐derived carbon at the time the Δ¹⁴C excursion and pockmarks were produced. An ocean general circulation model equipped with isotope tracers is unable to simulate these large Δ¹⁴C excursions on the Chatham Rise by transport of hydrothermal carbon released from the East Pacific Rise as previous studies suggested. Here we attribute the Δ¹⁴C anomalies and pockmarks to release of ¹⁴C‐dead CO₂ and carbon‐rich fluids from subsurface reservoirs, the most likely being dissociated Mesozoic carbonates that subducted beneath the Rise during the Late Cretaceous. Because of the large number of pockmarks and duration of the Δ¹⁴C anomaly, the pockmarks may collectively represent an important source of ¹⁴C‐dead carbon to the ocean during glacial terminations.

A transect of seafloor heat probe measurements on the Hikurangi Margin shows a significant increase of thermal gradients upslope of the updip limit of gas hydrate stability at the seafloor. We interpret these anomalously high thermal gradients as evidence for a fluid pulse leading to advective heat flux, while endothermic cooling from gas hydrate dissociation depresses temperatures in the hydrate stability field. Previous studies predict a seamount on the subducting Pacific Plate to cause significant overpressure beneath our study area, which may be the source of the fluid pulse. Double‐bottom simulating reflections are present in our study area and likely caused by uplift based on gas hydrate phase boundary considerations, although we cannot exclude a thermogenic origin. We suggest that uplift may be associated with the leading edge of the subducting seamount. Our results provide further evidence for the transient nature of fluid expulsion in subduction zones.

Moderate elevated vertical methane (CH₄) flux is associated with sediment accretion and raised fluid expulsion at the Hikurangi subduction margin, located along the northeast coast of New Zealand. This focused CH₄ flux contributes to the cycling of inorganic and organic carbon in solid phase sediment and pore water. Along a 7 km offshore transect across the Porangahau Ridge, vertical CH₄ flux rates range from 11.4 mmol·m⁻²·a⁻¹ off the ridge to 82.6 mmol·m⁻²·a⁻¹ at the ridge base. Stable carbon isotope ratios (δ¹³C) in pore water and sediment were variable across the ridge suggesting close proximity of heterogeneous carbon sources. Methane stable carbon isotope ratios ranging from −107.9‰ to −60.5‰ and a C1:C2 of 3000 indicate a microbial, or biogenic, source. Near ridge, average δ¹³C for pore water and sediment inorganic carbon were ¹³C-depleted (−28.7‰ and −7.9‰, respectively) relative to all core subsamples (−19.9‰ and −2.4‰, respectively) suggesting localized anaerobic CH₄ oxidation and precipitation of authigenic carbonates. Through the transect there was low contribution from anaerobic oxidation of CH₄ to organic carbon pools; for all cores δ¹³C values of pore water dissolved organic carbon and sediment organic carbon averaged −24.4‰ and −22.1‰, respectively. Anaerobic oxidation of CH₄ contributed to pore water and sediment organic carbon near the ridge as evidenced by carbon isotope values as low as to −42.8‰ and −24.7‰, respectively. Carbon concentration and isotope analyses distinguished contributions from CH₄ and phytodetrital carbon sources across the ridge and show a low methane contribution to organic carbon.