Numerical Prediction of the Behavior of CO2 Bubbles Leaked from Seafloor and Their Convection and Diffusion near Southeastern Coast of Korea

Among various carbon capture and storage technologies to mitigate global warming and ocean acidification due to greenhouse gases, ocean geological storage is considered the most feasible for Korea due to insufficient inland space to store CO2. However, the risk of CO2 leakage and the behavior and environmental effects of the leaked CO2 need to be assessed for its successful implementation. Therefore, the behavior of CO2 bubbles/droplets dissolving into the surrounding seawater and the diffusion of dissolved CO2 by ocean flows should be accurately predicted. However, finding corresponding research has been difficult in Korea. Herein, the behavior and convection-diffusion of CO2 that was assumed to have leaked from the seafloor near the southeastern coast of Korea were numerically predicted using a multi-scale ocean model for the first time. In the simulation region, one of the pilot projects of CO2 ocean geological storage had started but has been temporarily halted. In the ocean model, hydrostatic approximation and the Eulerian–Lagrangian two-phase model were applied for meso- and small-scale regions, respectively. Parameters for the simulations were the leakage rate and the initial diameter of CO2. Results revealed that all leaked and rising CO2 bubbles were dissolved into the seawater before reaching the free surface; further, the change in the partial pressure of CO2 did not exceed 500 ppm during 30 days of leakage for all cases.

CO2 leakage can cause loss of benthic biodiversity in submarine sands

One of the options to mitigate atmospheric CO₂ increase is CO₂ Capture and Storage in sub-seabed geological formations. Since predicting long-term storage security is difficult, different CO₂ leakage scenarios and impacts on marine ecosystems require evaluation. Submarine CO₂ vents may serve as natural analogues and allow studying the effects of CO₂ leakage in a holistic approach. At the study site east of Basiluzzo Islet off Panarea Island (Italy), gas emissions (90-99% CO₂) occur at moderate flows (80-120 L m^-2 h^-1). We investigated the effects of acidified porewater conditions (pH_T range: 5.5-7.7) on the diversity of benthic bacteria and invertebrates by sampling natural sediments in three subsequent years and by performing a transplantation experiment with a duration of one year, respectively. Both multiple years and one year of exposure to acidified porewater conditions reduced the number of benthic bacterial operational taxonomic units and invertebrate species diversity by 30-80%. Reduced biodiversity at the vent sites increased the temporal variability in bacterial and nematode community biomass, abundance and composition. While the release from CO₂ exposure resulted in a full recovery of nematode species diversity within one year, bacterial diversity remained affected. Overall our findings showed that seawater acidification, induced by seafloor CO₂ emissions, was responsible for loss of diversity across different size-classes of benthic organisms, which reduced community stability with potential relapses on ecosystem resilience.

Differing responses of the estuarine bivalve Limecola balthica to lowered water pH caused by potential CO2 leaks from a sub-seabed storage site in the Baltic Sea: An experimental study

Sub-Seabed CCS is regarded as a key technology for the reduction of CO₂ emissions, but little is known about the mechanisms through which leakages from storage sites impact benthic species. In this study, the biological responses of the infaunal bivalve Limecola balthica to CO₂-induced seawater acidification (pH7.7, 7.0, and 6.3) were quantified in 56-day mesocosm experiments. Increased water acidity caused changes in behavioral and physiological traits, but even the most acidic conditions did not prove to be fatal. In response to hypercapnia, the bivalves approached the sediment surface and increased respiration rates. Lower seawater pH reduced shell weight and growth, while it simultaneously increased soft tissue weight; this places L. balthica in a somewhat unique position among marine invertebrates.

Small-scale modelling of the physiochemical impacts of CO2 leaked from sub-seabed reservoirs or pipelines within the North Sea and surrounding waters

A two-fluid, small scale numerical ocean model was developed to simulate plume dynamics and increases in water acidity due to leakages of CO2 from potential sub-seabed reservoirs erupting, or pipeline breaching into the North Sea. The location of a leak of such magnitude is unpredictable; therefore, multiple scenarios are modelled with the physiochemical impact measured in terms of the movement and dissolution of the leaked CO2. A correlation for the drag coefficient of bubbles/droplets free rising in seawater is presented and a sub-model to predict the initial bubble/droplet size forming on the seafloor is proposed. With the case studies investigated, the leaked bubbles/droplets fully dissolve before reaching the water surface, where the solution will be dispersed into the larger scale ocean waters. The tools developed can be extended to various locations to model the sudden eruption, which is vital in determining the fate of the CO2 within the local waters.

Multi-scale modeling of CO2 dispersion leaked from seafloor off the Japanese coast

A numerical simulation was conducted to predict the change of pCO(2) in the ocean caused by CO(2) leaked from an underground aquifer, in which CO(2) is purposefully stored. The target space of the present model was the ocean above the seafloor. The behavior of CO(2) bubbles, their dissolution, and the advection-diffusion of dissolved CO(2) were numerically simulated. Here, two cases for the leakage rate were studied: an extreme case, 94,600 t/y, which assumed that a large fault accidentally connects the CO(2) reservoir and the seafloor; and a reasonable case, 3800 t/y, based on the seepage rate of an existing EOR site. In the extreme case, the calculated increase in DeltapCO(2) experienced by floating organisms was less than 300 ppm, while that for immobile organisms directly over the fault surface periodically exceeded 1000 ppm, if momentarily. In the reasonable case, the calculated DeltapCO(2) and pH were within the range of natural fluctuation.