Modeling Spatial and Temporal pCO 2 Variability at 49 ° N / 16 . 5 ° W and 56 . 5 ° N / 52 . 6 ° W in the North Atlantic Ocean

Numerical models have been employed in understanding and capturing spatiotemporal CO2 trends, inter seasonal to decadal variability, and the characterization of thermal (pCO2-T) and non-thermal (pCO2-nonT) components of surface ocean pCO2 and air-sea CO2 flux variability. We employed the MIT Ocean General Circulation Model and available data from two sub polar North Atlantic observatories at 49°N, 16.5°W and 56.5°N, 52.6°W to capture in situ pCO2 observations and deconvolute bio-physical controlling processes. The model suggests that the pCO2 cycle is marked by a summertime minimum and wintertime maximum. The physical-chemical and biological response pattern of the model is in good accordance with the observed pCO2, pCO2-T and pCO2-nonT trends. Model outputs suggest that the CO2 cycle is governed by contrasting effects of seasonal cooling and warming on one side and spring and summer net biology activities on the other side. It also predicts year-round under saturation, indicating that the region is a weak to moderate sink for CO2. Citation: Benson NU, Asuquo FE (2016) Modeling Spatial and Temporal pCO2 Variability at 49°N/16.5°W and 56.5°N/52.6°W in the North Atlantic Ocean. 7: 330. doi:10.4172/2157-7617.1000330


Introduction
Global oceans provide large uptake potential for anthropogenic atmospheric CO 2 . This is of great importance considering the oceanic stabilization effect on atmospheric CO 2 . The oceanic uptake of excess of inorganic carbon from the atmosphere has been shown to account for approximately 48% of the total anthropogenic CO 2 emissions within the past 200 years [1]. Current estimates of global oceanic CO 2 sink as well as spatial variability and trends in pCO 2 in the oceans vary significantly [2][3][4][5][6][7]. Observational studies and models suggest that in some regions of the Atlantic Ocean, there may be a weakening trend of oceanic CO 2 sinks in response to either climate change or natural variability [7][8][9][10][11]. However, a switch from natural carbon storehouses to sources is highly unlikely in the foreseeable future.
The North Atlantic is considered to represent the most intense, per unit area, ocean sink for atmospheric CO 2 in the Northern Hemisphere on interannual to decadal timescales [12]. On the global scale, the oceans are the Earth's primary sink that offsets anthropogenic CO 2 emissions and the North Atlantic is the primary gate for CO 2 entering the global ocean. However, according to the recent IPCC Global Carbon Budget 2013 assessment report, the ocean sink for anthropogenic CO 2 is estimated to be approximately equal to land sink fraction [13]. The subpolar North Atlantic is a complex and highly variable oceanic region both in space and time, and has in recent studies been confirmed as a significant sink based on automated surface ocean pCO 2 measurements [6,8,14].
Numerical ocean biogeochemical models have been and are employed in understanding and capturing physical-chemical and biological response patterns and driving mechanisms (processes) of the uptake and release of CO 2 in several ocean studies, usually in conjunction with in situ pCO 2 measurements. These coupled physicalbiogeochemical models are efficient, fast, compact and powerful scientific tools that are employed in understanding the past and present and predicting future cycles and changes in ocean carbonate chemistry system and climate, while being conscious of the ocean and climate system's response to this forcing [15] have outlined the imperativeness of accurately capturing both the pCO 2 -T (thermal) and pCO 2 -nonT (non-thermal) in order to understand the biogeochemical dynamics of the oceanic pCO 2 cycle.
The fundamental objective of this study is the characterization and evaluation of the annual/seasonal patterns of pCO 2 variability as well as a deconvolution of associated physical and biological forcing's using MITgcm biogeochemical model and in situ data from two North Atlantic time series sites. While we seek to improve our ability to quantify the North Atlantic sink as well as understand the surface ocean pCO 2 variability and governing processes, we took advantage of observed data from the porcupine abyssal plain (PAP) and the K1 central labrador sea (K1 CELAS) observatories, and attempt to investigate how well a physical-biogeochemical-ecosystem model can capture in situ observations of surface ocean pCO 2 in the North Atlantic western and eastern oceanographic locations.

Sources of observed data and deployments
The details of the sampling strategy, analytical methods, data quality control and calibration procedures employed for taking the diverse set of biogeochemical and physical measurements and recovered observed data at the KI CELAS and PAP observatories have been reported by [18,19] Figure 1 and Table 1). The deployments at the Porcupine Abyssal Plain time series site were designated PAP2, PAP3 and PAP4 for observations made during the second, third and fourth deployments respectively.

The physical-biogeochemical-ecosystem ocean model
Details of the MIT ocean general circulation model (MITgcm) used for the comparison of observational and model data in this study, which is regionally configured for the North Atlantic between 20°S and 81.5°N with a horizontal resolution of 0.5° latitude and 0.5° longitude have been described in previous papers and the model has been applied to the North Atlantic carbon sink by [7,[20][21][22]. However, the distinctive characteristics of the physical part of the biogeochemical model are shown in Table 2, while a schematic presentation of the ecosystem part is illustrated in Figure 2.
The ecosystem model is moderately complex taking into account the cycling of carbonates and nutrients [16,23]. It also features a pelagic ecosystem with one zooplankton class and two phytoplankton classes (diatoms and small phytoplankton) [7,15]. The atmospheric pCO 2 forcing in the model is incorporated from Mauna Loa Observations [24], and takes into account the annual and seasonal pCO 2 cycle. The dissolved and particulate carbon pools are resolved, and parameterisation of calcium carbonate cycling is included, while the sea-air exchange of CO 2 in the biogeochemical model is parameterized according to [25], with surface ocean pCO 2 determined as discussed by [23]. A positive flux is described as a net CO 2 exchange from sea to the atmosphere (outgassing) while a negative flux value implies a net CO 2 invasion from air to the sea (uptake).

Calculation of pCO 2 thermal and non-thermal components
The exchange of CO 2 between the atmosphere and ocean surface is primarily controlled by the air-sea difference in partial pressure of CO 2 (∆pCO 2 ), which provides the thermodynamical driving force and is regulated by naturally occurring physical, chemical and biological processes [26][27][28][29][30][31]. Considering the thermodynamics of the marine CO 2 system, the thermal component of surface-water pCO 2 variability is comparable in magnitude but opposite in direction to the biological effect [32]. However, it has been demonstrated that that the temperature (thermal) and biological (non-thermal) effects on seasonal pCO 2 amplitude can be derived from the total concentration of CO 2 effect, and the relative contributive impact of each processes for regulating the seasonal amplitude in pCO 2 assessed over the oceans [6]. In this synthesis, this approach was adopted in separating the thermal and non-thermal components of surface ocean pCO 2 variability.
In order to remove the temperature effect from the observed pCO 2 , the observed pCO 2 values are normalized to the annual mean temperature of seawater at the observatory, using the equation: and the effect of temperature changes on pCO 2 at a given observed temperature, T obs is computed using the equation: where T is the temperature in °C, pCO 2 is the partial pressure of carbon dioxide in µatm, and the subscripts " X " and "obs" indicate annual mean and observed values respectively, [6].

Data/model pCO 2 comparison at the NE PAP observatory
In this section, the seasonal cycles of pCO 2 , pCO 2 -T, and pCO 2 -nonT from the biogeochemical model and the PAP (PAP2 -PAP4) moored data are presented. For the most part and within specific surface mixed layer, the ocean biogeochemical model appears to be highly consistent and compares favorably with the PAP observed pCO 2 data and derived thermal/non-thermal components (Figures 3a-3c). However, the disagreement between the in situ PAP observations and model results at certain mixed layer depth may be partly attributed to the model resolution, and this obviously presents a major limitation of the MITgcm biogeochemical model.
Model results for the seasonal pCO 2 cycle relative to the observations indicate that the model captures the cycle reasonably well although with significant underestimation of the pCO 2 amplitude especially within the 35-75 m depth mixed layer during summer/early wintertime of the 2 nd (PAP2) and 3 rd (PAP3) deployments. As described below, the model fails to capture the full magnitude of the biological signal (nutrient drawdown) observed at the PAP site. However, the model captured the characteristic subpolar time series observed pCO 2 summertime minimum, wintertime maximum.
The model pCO 2 -T cycle (Figure 3b) has a pattern that clearly captures the cycle shown by the PAP data especially within the 20 -55m mixed layer, although the peak data pCO 2 -T amplitudes during late summer of the second deployment (PAP2) were higher than the model pCO 2 -T cycle. Interestingly, the model pCO 2 -T cycle at 55-75 m mixed layer depth was obviously inconsistent with the seasonal data pCO 2 -T pattern. The amplitude of the model pCO 2 -T (55-75 m) cycle was relatively higher in magnitude than the data pCO 2 -T at this depth range. In general, the model and data results reveal that the summer months were characteristically warmer with 2004 recording comparatively higher sea surface temperature values. Surface temperature values were however low at the subsurface euphotic depths especially during the winter months.
Model non-thermal forcing pattern (pCO 2 -nonT) consistently agree with the data pCO 2 -nonT cycle at 20-75 m depth mixed layer on a seasonal scale, and the amplitude broadly aligns with or comes within the range of the observed pCO 2 -nonT. However, the model pCO 2 -nonT cycle at 55-75 m mixed layer underestimates the biology signal clearly imprinted by the data pCO 2 -T cycle during the summertime of the 2 nd deployment (Figure 3c).
A combination of the seasonal patterns of pCO 2 -T and pCO 2 -nonT cycles creates the model pCO 2 cycle, which indicates a large increase in surface pCO 2 during the winter and a significant decrease during the summertime. The enhanced surface pCO 2 is due to vertical mixing of CO 2 (DIC) enriched water from the subsurface depths (mesopelagic zone) to the surface of the ocean in winter [12,18]. However, it is obvious that the model actually estimates the surface mixed layer pCO 2 levels of the PAP site, and could be considered realistic in this regard.
The modeled pCO 2 performance during the summertime corresponding to the period of the second deployment suggests a temperature-controlled cycle when actually it is the non-thermal forcing that is dominantly regulating the total pCO 2 pattern. From non-thermal component attributes of the observed/model output, it can be proposed that the pCO 2 variability of the PAP site may not be governed primarily by seasonal temperature changes only, but also by biological forcing fueled by carbon fixation and nutrient availability,    as well as deep vertical mixing of dissolved inorganic carbon enriched water from the bottom mixed layer depth to the ocean surface.
A striking observation is that the model evidently captures observed pCO 2 summertime minimum and wintertime maximum respectively except for 2003, and also accurately depicted that pCO 2 is governed by contrasting effect of temperature and biology. The model captured observed cycles of pCO 2 , thermal and biology components reasonably well within surface 20-35 m with caveats at deeper layers. More so, the model underestimated the amplitude of annual temperature cycle at 35 m depth and below, suggesting shallow depth of seasonal thermocline.

Data/model pCO 2 comparison at NW K1 central labrador sea
In this section, the pCO 2 , pCO 2 -T and pCO 2 -nonT model outputs to observed data at the K1 Labrador Sea time series are presented (Figure 4). At site K1, model pCO 2 has a rather peculiar jump during wintertime (mid-Feb. in 2004, and mid-Jan. in 2005), followed by a rather long period of flat value. It seems the model grossly underestimated both the physical and biological activities at this site during this period. The model pCO 2 is obviously out of phase with observed pCO 2 by about 3 months, peaking in November instead of February. More so, pCO 2 peak-to-peak amplitudes between the observed data and model output in early fall (October 2004) and in wintertime (January-March) were out of phase by about 60 μatm. This difference might be considered very significant, thereby indicating a moderate flaw of the model. Discrepancies might be due to model being too cold in the fall, and thus not allowing biological activity to drawdown pCO 2 as much as observed. Further, although both model and observed data show a dominance of biological activities over pCO 2 variability, caution must be made since modeled seasonal patterns were quite different from the data.
Model pCO 2 -T of K1 CELAS site looks good between November through June although it appears too cold in late summer/fall. However, it warms up late by three months (three months phase offset). K1 CELAS data pCO 2 trend may be considered to be in-phase with model pCO 2 -T pattern, although there is a significant digression during summertime. The MITgcm fairly captures the K1 CELAS wintertime pCO 2 -nonT pattern ( Figure 4). However, it partially fails in capturing the summertime biological signal observed at the K1 site within the 35-55 m depth. It appears the model pCO 2 -nonT is not warming enough in summer hence its failure to capture the observed pCO 2 -nonT drawdown and biological productivity.

Anomalies of pCO 2 , pCO 2 -T and pCO 2 -nonT
The patterns of the averaged daily anomalies of pCO 2 , pCO 2 -T and pCO 2 -nonT calculated from observed data at the PAP and K1 CELAS observatories and model outputs are presented in Figures 5 and 6 respectively. Measurements at PAP site shows a much elevated pCO 2 during summer 2003. However, the model grossly over-estimated the biological drawdown during this period, resulting in much lower pCO 2 ( Figure 5). Relatively high positive anomalies of pCO 2 (61 µatm), pCO 2 -T (58 µatm), and pCO 2 -nonT (77 µatm) were obtained based on average daily data for PAP pCO 2 , temperature and non-temperature components of the total pCO 2 considered. More so, the computed anomalies based on daily averages indicated negative deviations of 54, 58 and 84 µatm for pCO 2 , pCO 2 -T and pCO 2 -nonT respectively. For K1 CELAS site, positive anomalies of pCO 2 (35 µatm), pCO 2 -T (67 µatm), and pCO 2 -nonT (57 µatm) were obtained based on average daily data, with corresponding negative deviations of 61, 25 and 96 µatm for pCO 2 , pCO 2 -T and pCO 2 -nonT respectively ( Figure 6).

Conclusions
In this paper, the seasonal distribution of surface ocean observational data and associated thermal and non-thermal factors controlling the variability and processes of pCO 2 at the PAP (Northeastern) and K1 Central Labrador Sea (Northwestern) observatories in the North Atlantic Ocean have been investigated using a modestly highresolution ocean general circulation model. This work indicates that in the subpolar North Atlantic on a seasonal timescale the observed as well as the modeled pCO 2 distribution is marked by minimum and maximum pCO 2 for the summertime and wintertime respectively. The seasonal cycles of model pCO 2 , pCO 2 -T, and pCO 2 -nonT of the biogeochemical model compared to observational data revealed that the thermal component of surface water pCO 2 is comparable in magnitude but opposite in direction to the biological effect. For the most part and within specific depth of the mixed layer, the ocean biogeochemical model appears to reproduce seasonal pCO 2 , pCO 2 -T, and pCO 2 -nonT phases and amplitudes favorably consistent with the seasonal PAP pCO 2 , pCO 2 -T, and pCO 2 -nonT observations. However, few but significant disagreements between the in situ PAP observations and model results at certain mixed layer depths may be partly attributed to the model resolution amongst various factors. In general, model results for the seasonal pCO 2 , pCO 2 -T, and pCO 2 -nonT cycles relative to the observatory data indicates that the model captures the pCO 2 , pCO 2 -T, and pCO 2 -nonT cycles fairly well although with significant underestimation of the pCO 2 and pCO 2 -nonT amplitudes especially within the 35 -75m mixed layer depth during summer/early wintertime of the 2nd and 3rd deployments. In order words, the model fails to capture the biological signal (nutrient drawdown) observed at the PAP site.