Understanding caldera degassing from a detailed investigation at Lake Rotoiti, Okataina Volcanic Centre, New Zealand

Volcanic lakes in large silicic caldera volcanoes are an important source of CO2 emissions. However, quantifying CO2 output is challenging due to the lack of observed historical CO2 flux records and the large size of the volcanic and hydrothermal systems. Twenty percent of the surface area of Okataina Caldera, Taupo Volcanic Zone, New Zealand is covered with lakes. Geothermal expressions, predominantly on the Okataina caldera margin, occur in at least six different locations, with surface expressions both on land and under water. Lake Rotoiti is located at the northwest edge of the Okataina Caldera, spans across the Tikitere and Taheke geothermal fields, and has inputs from on-land thermal springs, lake floor hydrothermal vents, and from Lake Rotorua. CO2 flux from Lake Rotoiti was assessed using the accumulation chamber method and three techniques are used to process the data for total CO2 emission of the lake: (1) sequential Gaussian simulation (sGs) method allows quantification of CO2 emission with spatial control, (2) graphical statistical approach (GSA) allows the quantification of CO2 emission from different degassing regimes, and (3) a method based on water chemistry of the lake. We find CO2 is mostly emitting at Tumoana Bay and Central Basin, and the emission rate is 271 ± 38 t d−1 of CO2 (based on sGs). We then added this estimated emission data to the existing CO2 data for individual geothermal systems distributed around the active Okataina volcanic centre and calculate that the entire caldera is emitting at least 1856 t d−1 of CO2. The total emission is dominated by the lakes, where the topography is lower. We discuss the implications of this study in terms of preferential degassing locations and the amount of degassing, in particular for CO2, in an active caldera setting.