An oceanic two-phase plume model is developed to include bubble size distribution and bubble interactions, applied to the prediction of CO2 bubble plume and CO2 solution dynamics observed from the recent QICS field experiment in the Scottish sea at Ardmucknish Bay. Observations show bubbles form at between 2 and 12 mm in diameter, where the inclusion of the interactions within the simulations brings results of bubble plumes closer to that of the experiment. Under a given leakage flux, simulations show that the bubble size affects the maximum pCO2 dissolved in the water column, while the bubble interactions affect the vertical bubble distribution. The maximum modelled pCO2 increases from a background 360 ìatm to 400, 427 and 443 ìatm as CO2 injection rates increase from 80, 170 to 208 kg/day respectively at low tide. An increase of the leakage rate to 100% of the injection rate shows the maximum pCO2 could be 713 ìatm, approaching the mean pCO2 observed of 740 ìatm during the high leakage component of the experiment, suggesting that the flux may be greater than estimated due to the varied flux and activity across the pockmarks during the leakages. This is a publication in QICS Special Issue - International Journal of Greenhouse Gas Control, Peter Taylor et. al. Doi:10.1016/j.ijggc.2014.09.007.
