Abstract
Here, we investigate the seasonal variability in the dissolved inorganic carbon (DIC) cycle in the Northwest Pacific using a high-resolution biogeochemical and carbon model coupled with an operational ocean model. Results show that the contribution to DIC from air–sea CO2 exchange is generally offset by vertical mixing at the surface at all latitudes, with some seasonal variation. Biological processes in subarctic regions are evident at the surface, whereas in the subtropical region these processes take place within the euphotic layer and then DIC consumption deepens southward with latitude. Such latitudinal differences in biological processes lead to marked horizontal and vertical contrasts in the distribution of DIC, with modulation by horizontal and vertical advection–diffusion processes.
Similar content being viewed by others
1 Introduction
The atmospheric partial pressure of CO2 (pCO2) has been increasing at a rate of ~ 1.8 ppm by volume (ppmv) per year in recent decades as a result of human activities such as fossil-fuel burning, deforestation, and cement production (Takahashi et al. 2009; IPCC 2013). In the pre-industrial era, the ocean was generally a net source of CO2 emissions to the atmosphere because of the mineralization of land-derived organic matter in addition to that produced by in situ production, and CaCO3 precipitation (Mackenzie et al. 2004). Rising atmospheric CO2 concentrations caused by fossil-fuel combustion and land-use changes (Mackenzie et al. 2004; Bauer et al. 2013; IPCC 2013; IGBP, IOC, SCOR 2013) reversed the direction of the air–sea CO2 flux, leading the global ocean to become a net sink of anthropogenic CO2. The present thickness of the upper thermocline, where large amounts of anthropogenic CO2 emissions are stored, is estimated to be of the order of a few hundred meters (Mackenzie et al. 2004). The oceanic coastal zone changed from being a source to a sink during the industrial era (Mackenzie et al. 2004; Bauer et al. 2013). Several estimates of CO2 sinks and sources in ocean provinces (Cai et al. 2006) and/or spatially explicit typology (Laurelle et al. 2010) showed that marginal seas in the tropics are sources of CO2, whereas those in temperate regions and at high latitudes act as sinks (Cai et al. 2006; Laurelle et al. 2010).
Data-based estimates of variability and trends in ocean CO2 uptake are limited by the short record of observations. Although high-quality measurements of CO2 in surface waters and air commenced in the early 1960s, the amount of available information is still limited (Wanninkhof et al. 2013). The principal observational approaches for estimating sea–air fluxes of CO2 are to measure ΔpCO2 from ships (Takahashi et al. 2009; Nakaoka et al. 2013) and moorings (Sutton et al. 2017), and apply a parameterization using a function of wind speed (Wanninkhof et al. 2013). Other approaches rely on simulations made by ocean biochemistry general circulation models (OBGCMs) with parameterization of biogeochemical processes and total dissolved inorganic carbon (DIC) measurements in the ocean interior, and/or atmospheric data. However, gaps remain in the understanding of ocean CO2 uptake, especially the spatiotemporal variability of the seasonal inorganic/organic carbon cycle, because CO2 concentrations and other related oceanic variables are difficult to observe simultaneously, frequently, and widely. The seasonal variability in pCO2 shows differences at a local scale (Takahashi et al. 2009; Sutton et al. 2017). Model estimates of temporal trend detection (Keller et al. 2014; Lovenduski et al. 2015) show the influence of both decadal/interannual and seasonal variabilities and suggest that the time of emergence of a trend signal is basically around 10 years in the surface but the tropical area needs more time.
The detailed processes that generate variation in the DIC of the ocean interior are still uncertain. Several studies have proposed possible mechanisms for the oceanic annual carbon cycle (Palmer and Totterdell 2001; Takahashi et al. 2002; Xiu and Chai 2013). Palmer and Totterdell (2001) discussed physical and biological mechanisms that contribute to the global annual mean carbon cycle using an ecosystem model without considering the contribution of air–sea CO2 exchange. They reported that the effects of vertical mixing were largely offset by biological processes in the latitudinal range of 25−60° N over ocean surfaces and that the effects of advection were mostly offset by biological processes at latitudes of < 20° N. Takahashi et al. (2009) focused on the relative importance of temperature and biological effects to the global seasonal cycle of air–sea CO2 exchange by evaluating monthly climatological maps of air–sea CO2 flux and pCO2. Xiu and Chai (2013) investigated the seasonal and decadal variability of the upper-ocean carbon cycle in the North Pacific using a physical–biogeochemical model. Their results showed that the seasonal variability in pCO2 and CO2 flux in the North Pacific followed the change in sea surface temperature closely, with high and low values in summer and winter, respectively, and that surface pCO2 variations at the modeled sites correspond to well-known observational monitoring points controlled primarily by anthropogenic CO2 and modulated by decadal variations.
The Japan Coastal Ocean Predictability Experiment (JCOPE; http://www.jamstec.go.jp/jcope/) is an operational eddy-resolving physical ocean model for the Northwest Pacific, the Japan Sea, the Okhotsk Sea, and the East China Sea (Miyazawa et al. 2009, 2014). Ishizu et al. (2019) recently developed a biogeochemical and carbon model coupled with the JCOPE (JCOPE_EC). This model generally reproduces the observed seasonal variability of chlorophyll-a (Chl-a), dissolved inorganic nitrogen (DIN), phosphorus (DIP), and inorganic carbon (DIC), and total alkalinity (ALK), but includes monthly climatological damping for DIN, DIP, DIC, and ALK (Ishizu et al. 2019). The damping forcibly constrains the calculated biological parameters (DIN, DIP, DIC, and ALK) around the monthly climatological values with a timescale of 30 days (Ishizu et al. 2019), meaning that those authors were unable to discuss the mechanism of the inorganic carbon cycle. In this study, we therefore performed a simulation using JCOPE_EC without any climatological damping and examined the physical and biological mechanisms represented by the model dynamics, focusing on the critical roles of horizontal and vertical advection–diffusion processes in generating seasonal variation in DIC.
We present the results from JCOPE_EC (without climatological damping) and discuss the mechanisms responsible for the simulated seasonal inorganic carbon cycle for the Northwest Pacific. Details of the model configuration are given in Section 2, model accuracy is described in Section 3, the processes involved in the inorganic carbon cycle in the model are discussed in Section 4, and the conclusions of the study are provided in Section 5.
2 Model and data
2.1 Model configuration
The JCOPE_EC (Ishizu et al. 2019) is an off-line tracer model driven by physical processes simulated by an operational eddy-resolving ocean general circulation model (JCOPE2M; Miyazawa et al. 2017) based on the Princeton Ocean Model with a generalized sigma coordinate (Mellor 2001). The model is a three-dimensional high-resolution regional model covering the Northwest Pacific (108–180° E, 10.5−62° N) with a horizontal resolution of 1/12° (4.4–9.1 km) and 46 vertical active levels. The model structure is the same as that described by Ishizu et al. (2019), but our model differs in that the governing equations for DIN, DIP, DIC, and ALK (equations 4, 5, 17, and 18 of Ishizu et al. 2019) remove climatological damping.
We determined the biogeochemical model parameters using multi-optimized operations (Ishizu et al. 2019) separately for subarctic and subtropical regions (Table 1). These biogeochemical parameters are the maximum growth rate from photosynthesis (Vmax), the phytoplankton mortality rate at 0 °C (MP), the phytoplankton respiration rate (R), the maximum grazing rate of zooplankton (Gz), the zooplankton mortality rate (M)z, the decomposition rate (VPN), and the optimum light intensity for phytoplankton (Iopt). In addition to the parameters described in Ishizu et al. (2019), we introduced latitudinal changes in Vmax, MP, R, Mz, VPN, and Iopt according to the results of several sensitivity experiments as follows:
where Vmax, MP, R, Mz, VPN, and Iopt change latitudinally from 0.28 to 0.97 day−1, 0.054 to 0.155 (mmol Nm−1) m−3 day−1, 0.0011 to 0.00256 day−1, 0.044 to 0.12 (°C)−1, 0.0954 to 0.47 day−1, and from 20 to 100 W m−2, respectively; Vmax0, R0, MP0, Mz0, VPN0, Iopt0, Vmax1, R1, MP1, Mz1, VPN1, and Iopt1 are the tunable parameters (Table 1); Latbnd.vmax, Latbnd, and Latslp are the coefficients representing the values at latitudinal boundaries and the latitudinal slopes for these parameters, respectively (Table 1 and Eqs. 1, 3, 4, 5, and 6).
The first version of JCOPE_EC had a model bias resulting in a large decrease in ALK in summer due to anomalously high CaCO3 production with a large increase in Chl-a during summer in subarctic regions (Ishizu et al. 2019). To suppress this large decrease in ALK in this model, we set the CaCO3 to non-photosynthetic POC production ratio to a much smaller value (0.00035) in the version of the model used here (Table 1). This improvement allows the model to well represent the seasonal variability of ALK in our target region (Section 3).
The model was driven by forcing from daily oceanic (JCOPE2M) and six-hourly atmospheric (NCEP/NCAR) reanalysis data for a 1-year period (2015). The initial concentrations of phytoplankton were set to 0.1 and 0.0 mmol N m−3 for depths above and below 150 m, respectively. The initial zooplankton concentrations were set to 10% of the phytoplankton concentration. The initial detritus concentration was set to 0.0 mmol N m−3. The variables DIN, DIP, and DIC were initialized using the climatology for January, and ALK was initialized using the annual climatology (Ishizu et al. 2019).
2.2 Model validation
To validate phytoplankton concentrations in the model, MODIS-Aqua Ocean Color Data for Chl-a from 2015 were used, as downloaded from the website of the Physical Oceanography Distributed Active Archive Center (PODAAC). To validate the model results for DIN, DIP, and DIC, we used monthly climatological DIN, DIP, and DIC data (World Ocean Atlas 2013 (WOA13); Yasunaka et al. 2013) and Japan Meteorological Agency (JMA) observational data for 2015, as in Ishizu et al. (2019). There are no applicable monthly climatological datasets for ALK in our target region (Goyet et al. 2000; Key et al. 2004; Takatani et al. 2014, Takahashi et al. 2014). We therefore used ALK observational data obtained by JMA in 2015 for comparison with model results (Section 3).
3 Results
3.1 Accuracy of modeled Chl-a, DIN, DIC, and ALK
Model results presented here are slightly less accurate than those of the model with the climatological conditions described by Ishizu et al. (2019), except for ALK. However, the model used here still simulates observed conditions well, capturing the basic seasonal variations of DIN, DIP, DIC, and ALK in comparison with the climatological and in situ data (Figs. 1, 2, 3, 4, and 5). Although the satellite data in the subarctic region exhibit distinctive double peaks in spring and autumn (Fig. 5(a)), the concentration of Chl-a given by the model shows only a single peak in summer in the subarctic region (Figs. 1 and 5(a, d)). Ishizu et al. (2019) suggested that this single peak in summer is caused by a lack of iron restriction (Tsuda et al. 2003) of the model. The seasonal variabilities in DIN and DIC in the model outputs are represented well (Fig. 2(g, h) and the 50° N line in Fig. 5(e)) compared to the DIN and DIC climatology (Figs. 2(c, d) and 5(c−e)). It is difficult to judge whether the distribution of ALK is accurate, but modeled ALK concentrations are generally constant throughout the year (Fig. 4(d−g)) and are consistent with JMA observations (Fig. 4(a−c)).
The reproducibility of the model is poorer than that of Ishizu et al. (2019) for DIN and DIC (Table 2), but the reproducibility of ALK is much improved, especially in the subarctic region, the Kuroshio Extension, and the Japan Sea. The correlation between observed ALK and simulated ALK is low and negative (R = − 0.24) in the Japan Sea because the simulated ALK values there have near-uniform values with depth (not shown).
3.2 Ocean acidification indices pHin situ, pH25, and Ωarg
The ocean acidification indices pHin situ, pH25, and aragonite saturation (Ωarg) were calculated from model results for temperature, salinity, DIC, and ALK. The pHin situ values change throughout the year (Figs. 6 and 7). The summer pHin situ values in the subarctic region (150–175° E, 50−60° N) are slightly higher (8.05–8.10) than those shown in Fig. 9 of Ishizu et al. (2019); 7.95–8.00; Figs. 6(a−d) and 7a). Summer pHin situ values for regions north of 35° N decreased by 0.05 from winter to summer (Figs. 6(a−d) and 7a). The general seasonal variability is consistent with the previous model version used by Ishizu et al. (2019).
Horizontal distributions of pH25 exhibit consistent seasonal variability across all latitudes (Figs. 6(e−h) and 7b) but with different amplitudes. The pH25 values at higher latitudes are lower than at lower latitudes, with an increase of 0.15 in summer. Their amplitudes gradually decrease southward. The pH25 values in the subtropical region south of 20° N are generally constant throughout the year.
The seasonal pattern of Ωarg (Fig. 6(i–l)) is similar to that of pH25 (Fig. 6(e–h)). The Ωarg values in the same area of pH25 north of 35° N become 0.5–1.0 larger in summer and in autumn. The summer increase in Ωarg is less evident in the south.
Correlation coefficients for modeled pHin situ, pH25, and Ωarg are higher than those reported by Ishizu et al. (2019) because the accuracy of the modeled ALK concentrations is much improved compared to them.
4 Discussion
4.1 Processes affecting the seasonal inorganic carbon cycle
To identify processes affecting the inorganic carbon cycle, we examined the physical and biological mechanisms underlying the seasonal inorganic carbon cycle in the Northwest Pacific using the model results. We separately evaluate each process included in the governing equation as follows:
where the subscripts A, xy_dif, z_dif, Bio, and air–sea represent the time derivatives of DIC induced by advection, horizontal diffusion (i.e., horizontal mixing), vertical mixing, biological processes, and air–sea CO2 exchange (positive values indicate a transfer of CO2 from air to sea), respectively. Note that we refer to these as “DIC variation terms” here, and the DIC variations induced by air–sea CO2 exchange are included only for the top (surface) level and not for the levels below it. The total DIC variation term on the left side of Eq. (7) is a summation of all the terms on the right-hand side of the equation.
Air–sea CO2 exchange shows negative values (CO2 emission to the atmosphere) in winter north of 35° N (Fig. 8(a)) and positive values (CO2 absorption from the atmosphere) south of 35° N, and vice versa in summer for these regions. The subarctic region releases CO2 to the atmosphere in winter north of 35° N (Fig. 8(a)), and absorbs CO2 from the atmosphere in summer north of 35° N (Fig. 8(c)). In contrast, the subtropical region south of 35° N intensely absorbs CO2 south of the Kuroshio Extension in winter (Fig. 8(a)). A weak release of CO2 from the ocean to the atmosphere occurs south of 40° N in summer. The areas of strongest release are located in the northeast of the Kuril Islands in the Oyashio region and in the Okhotsk Sea in autumn–winter (Fig. 8(a, d)). (Further details on air–sea CO2 exchange represented in the model are described in the Appendix.)
Horizontal distributions and monthly mean balances of the DIC variation terms in the surface layer induced by all subprocesses (Figs. 8 and 9(a−d)) indicate that the air–sea CO2 exchange shown in Fig. 8(a−d) is generally offset by vertical mixing. Biological processes, however, made a subordinate contribution to the total DIC variation (Figs. 8 and 9(c, d, g, h, j, k, m)). Negative biological process values indicate consumption of DIC through photosynthesis. Biological process contributions vary spatially and with depth (Figs. 8 and 9; Appendix Figs. 16–18). At the surface level, the DIC consumption induced by biological processes is high in the subarctic region around 50° N throughout the year (Figs. 8(m−p) and 9(a−d)). The peak consumption of DIC moves deeper beneath the surface southward (Fig. 9(d, e–g, j, m); Appendix Figs. 17 and 18). The highest DIC consumption occurs at 50–100 m depth at 30–40° N (Fig. 9(g, j, k); Appendix Figs. 16(m–p) and 17(m–p)) and at 200 m depth in the subtropical region south of 30° N (Fig. 9(m); Appendix Fig. 18(m–p)). The contributions from biological processes at the surface and at 200 m depth are of opposite sign (Fig. 8(m–p); Appendix Fig. 18(m–p)). Instead of the contributions of the vertical mixing or biological processes below the surface, we see that advection processes are relatively contributed to the total DIC variation (Fig. 9(i–k, n, m)).
The relative contribution of biological processes compared with the other terms at 165° E is calculated as
as shown in Fig. 10. The highest DIC consumption occurs above 100 m depth north of 40° N throughout the year, gradually deepening in the range 30–40° N, and spreading vertically at 100–350 m depth in the subtropical region at latitudes south of 30° N. These patterns of DIC consumption/production may be caused by the latitudinal difference of the Chl-a maximum depth (Ishizu et al. 2019; Sauzede et al. 2015). Vertical distributions of Chl-a along 165° E from Ishizu et al. (2019) and global ocean Chl-a data from Sauzede et al. (2015) indicate that the Chl-a maximum deepens southward; the Chl-a maximum is located in the surface layer in subarctic regions and at ~ 150 m depth in subtropical regions. The highest Chl-a consumption in the subtropical region from our results (Fig. 10) occurs at greater depths (200–500 m), but the magnitude is relatively low and makes a negligible contribution to the total variation in DIC (Fig. 9(n−p)).
4.2 Mechanisms of the seasonal inorganic carbon cycle
The relative contributions of terms in the governing equation of DIC (Figs. 8 and 9; Appendix Figs. 16–18) suggest possible mechanisms for the seasonal carbon cycle, as follows. CO2 (DIC) is absorbed from the atmosphere to the ocean during winter south of 40° N. The corresponding absorbed volume of DIC is conveyed from the surface to the subsurface by vertical mixing (Figs. 8(i, j) and 9(a−c)). In summer and autumn, CO2 (DIC) is released to the atmosphere by air–sea interactions, and vertical mixing generally offsets carbon emissions (Figs. 8(c, d) and 9(a−c)). In the subarctic region, this interaction is opposite that in the subtropical region (Figs. 8(a–k) and 9(a−d)), with a negative contribution from biological processes (i.e., a sink) throughout the year (Figs. 8(m−p) and 9(a−d)). Below 50 m depth, the advection term becomes more prominent in the subtropical region and in the Kuroshio Extension region, connecting with total DIC time variation (Fig. 9(i–k, n, m)). Working against the small negative contribution of biological processes below 200 m depth in our model (Fig. 9(n−p)) is advection, which can produce an increase or decrease in DIC in the subtropical region (Fig. 9(n)).
The modeled carbon cycle also includes latitudinal differences in total DIC time variation, as shown in Figs. 12 and 13(a). Differences in the contribution from biological processes are related to the total DIC time variation both horizontally and vertically (Figs. 12 and 13(a)). A distinctive area of DIC decrease exists near 0–200 m depth, and deepens southward from the subarctic region to the subtropical region (Figs. 11(a−d) and 12b). A comparison between the annual mean DIC time variation and the Chl-a maximum depth (Ishizu et al. 2019; Fig. 12) shows some similarity above 200 m depth. Advection dominates the trend in annual mean DIC time variation in the subtropical region south of 30° N below 200 m depth (Fig. 9(n)). These results suggest that the positive and negative patterns of annual mean DIC time variation are caused by ocean currents (Fig. 9(n)).
The density range of this distinctive zone of contrast is comparable with that of the North Pacific Ocean Central Mode Water (CMW; σθ = 26.0–26.5 kg m−3; Oka and Suga 2005; Oka and Qui 2012; Fig. 12b). Those waters are formed by winter ventilation around thermocline fronts, including the Kuroshio Extension front, the Kurhoshio Bifurcation front, and the subarctic front, and are then spread by advection (Oka and Suga 2005; Oka and Qui 2012). We therefore suggest that the positive and negative contrast in the temporal DIC variations depends on uptake in the ventilation areas north of 40° N and could be transported by advection (Fig. 12). A sensitivity experiment was performed in which the maximum depth of Chl-a was decreased by adjusting the biological parameters related to photosynthesis (Iopt), to check whether advection affects the contrast in DIC time variation below the surface layer. Features of the DIC time variation south of 25° N were almost identical in the base and sensitivity experiments (not shown), supporting the dominant role of advection in the subtropical region.
4.3 Comparison with previous work on the marginal seas
The balance of the monthly mean DIC time variation terms in the East China and Japan seas (Fig. 13(a, b)) shows that air–sea CO2 exchange is generally offset by vertical mixing. These balances differ from the simulation results for the Yellow and East China Seas of Luo et al. (2015), which indicates that the contributions of vertical mixing, advection, and biological processes to the inorganic carbon cycle largely offset each other on the continental shelves and vary seasonally. The difference between the results of Luo et al. (2015) and those presented here may be due to missing processes in our model that may be required to represent local variability in these regions (e.g., the DIC input from river discharge), although the reproducibilities for the Japan and East China seas are high, with high correlation coefficients for several variables, including Chl-a, DIN, DIC, pH25, and Ωarg (Table 2).
The reproducibility for the Okhotsk Sea could not be evaluated because of a lack of observational data. Variations in DIC in Okhotsk Sea simulated by our model (Fig. 13(c)) are similar to those for the subarctic region of the Northwest Pacific, where the contributions of air–sea CO2 exchange and vertical mixing are generally offset by each other (Fig. 13(c)) and biological processes make a subordinate contribution.
The relatively small contribution of biological processes in the DIC cycle in marginal seas can be explained by the latitude-dependent functions of biological parameters (Section 2), as we adjusted the biological parameters to focus mainly on the Northwest Pacific in the subtropical region and on subarctic regions (Ishizu et al. 2019; Eqs. 1–6). Tittensor et al. (2010) showed that biological communities in the marginal seas are dissimilar to those of the Pacific at the same latitude, suggesting that there is a limitation for the ability of the parameters optimizing method applied in the present model (Section 2.1).
4.4 Comparison with a study based on observation data
Yasunaka et al. (2013) discussed the inorganic carbon cycle using surface DIC climatological data, omitting physical processes such as advection and vertical mixing. Our results agree with those of Yasunaka et al. (2013) in terms of the importance of biological processes from spring to summer, especially in the subarctic region. However, in our model, biological processes are a subordinate contributor to variation in DIC compared with both vertical mixing and air–sea exchange throughout the year. The subsurface vertical structure of the biological contribution varies with latitude.
Average daily net community production (NCP) values were estimated by Yasunaka et al. (2013) from March to July as > 14 mmol m−2 C day−1 in the Kuroshio Extension region (140–170° E, 30–40° N) and 2–6 mmol C m−2 day−1 in the subarctic region. The corresponding values of NCP in our model can be calculated as \( {\int}_0^{\mathrm{MLD}}\frac{\partial \left[\mathrm{DIC}\right]}{\partial t} dz \), where MLD indicates the mixed layer depth, which is defined as the depth at which the density is 0.125 kg m−3 greater than the density at the surface (Ohishi et al. 2019; Suga et al. 2004; Ohno et al. 2004). The modeled NCP (Fig. 14) shows positive values in the subtropical region of less than 2 mmol C m−2 day−1 and negative values in the subarctic region of less than − 4 mmol C m−2 day−1, respectively (Fig. 14), which are smaller than their respective estimated values (Yasunaka et al. 2013). In addition, our model predicts that the distinctive positive/negative (maximum/minimum) values spread sparsely and weakly only at the boundary between the Kuroshio Extension and the subarctic region. Such differences between our model and that of Yasunaka et al. (2013) can be explained by the horizontal and vertical advection/diffusion processes represented in our model, which should be considered for better understanding of the inorganic carbon cycle in our target regions.
5 Conclusions
To understand the seasonal inorganic carbon cycle in the Northwest Pacific, we performed simulations using a biogeochemical and carbon model coupled with an operational ocean model. The reproducibility of the model was sufficient to evaluate processes related to the seasonal inorganic carbon cycle. We found that contributions to the inorganic carbon cycle from air–sea CO2 exchange generally offset those of vertical mixing at the surface in the Northwest Pacific. Biological processes are a subordinate contributor to variation in DIC and show a latitudinal dependence in the euphotic layer. Advection actively contributes to variation in DIC below the layer where biological processes contribute.
A schematic representation of the main inorganic carbon cycle in the Northwest Pacific is shown in Fig. 15. DIC is absorbed from the atmosphere during winter south of 40° N and is released to the atmosphere north of 40° N (Fig. 15a). The DIC introduced in the subtropical region is conveyed from the surface to the subsurface by vertical mixing, and the DIC released in the subarctic region is compensated for by vertical mixing from subsurface layers (Fig. 15). In summer, the opposite pattern occurs because vertical mixing is weaker and the corresponding influence of mixing on variation in DIC decreases (Fig. 15b). Unlike the contribution from seasonal air–sea CO2 exchange, the negative contribution induced by biological processes (i.e., a sink) below the surface occurs throughout the year (Fig. 15) and becomes stronger to the subarctic region. Below the surface, the contribution of advection is prominent in the subtropical and Kuroshio Extension regions. This contribution of advection connects with the increase/decrease in DIC in the subtropical region (Fig. 15).
Ocean circulation (Tally et al. 1993; Suga and Hanawa 1995a; Yasuda et al. 1997; Qui and Chen 2006) transports DIC-rich and DIC-poor water and redistributes them over time. The horizontal pattern of long-term DIC trends varies with depth. The impacts of long-term trends in ocean acidification indices (pH and aragonite saturation) also show horizontal and vertical dependencies. Future studies will use modeling experiments to evaluate decadal variations in the NPZD and carbon cycle. Results of these experiments are expected to improve our understanding of the variability in carbon, biological processes, and ocean acidification.
References
Bauer JF, Cai WJ, Raymond PA, Bianchi TS, Hopkinson CH, Pierre A, Regnier P (2013) The changing carbon cycle of the coastal ocean. Nature 504:61–70. https://doi.org/10.1038/nature12857
Cai W-J, Dai HM, Wang CY (2006) Air–sea exchange of carbon dioxide in ocean margins: a province-based synthesis. Geophysical Res Letter 33:L12603. https://doi.org/10.1186/1472-6785-8-15
Goyet C, Healy R, Ryan J (2000) Global distribution of total inorganic carbon and total alkalinity below the deepest winter mixed layer depths, ORNL/CDIAC-127, NDP-076. In: Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory. U.S. Department of Energy, Oak Ridge, TN, USA, pp 1–40. https://doi.org/10.3334/CDIAC/otg.ndp076
IGBP, IOC, SCOR (2013) Ocean acidification summary for policymakers – third symposium on the Ocean in a High-CO2 World International Geosphere-Bioshere Programme. Stockholm, Sweden
IPCC, Summery for Policymakers. In: Climate Change 2013: The physical science basis. Contribution of working group I to the fifth assessment report of the Intergovernmental Panel on Climate Change [Stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J.Boschung, A. Nauels, Y. Xia, V. Bex and P.M. Midgley (eds.)]. Cambridge University press, United Kingdom and New York, NY, USA. 2013, pp.1–29
Ishizu M, Miyazawa Y, Tsunoda T, Guo X (2019) Development of a biogeochemical and carbon model related to ocean acidification indices with an operational ocean model product in the northwestern Pacific. Sustainability 11:2677. https://doi.org/10.3390/su11092677
Keller, M. K.; Joos, F.; Raible, C. C.: Time of emergence of trends in ocean biogeochemistry, Biogeosciences, 11, 3647–3659, 2014, www.biogeosciences.net/11/36472014/, doi.https://doi.org/10.5194/bg-11-3647-2014
Key MR, Kozyer A, Sabine LC, Lee K, Wanninkhof R, Bullister LJ, Feely AR, Millero JF, Mordy C, Peng TH (2004) A global ocean carbon climatology: results from global data analysis project (GLODAP). Glob Biogeochem Cycles 18:GB4031. https://doi.org/10.1029/2004GB002247
Laruelle GG, Durr HH, Slomp PC, Borges VA (2010) Evaluation of sinks and sources of CO2 in the global coastal ocean using a spatiallyexplicit typology of estuaries and continental shelves. Geophys Res Lett 37:L15607. https://doi.org/10.1029/2010GL043691
Li PQ, Franks JSP, Landry RM, Goericke R, Taylor GA (2010) Modeling phytoplankton growth rates and chlorophyll to carbon ratios in California coastal and pelagic ecosystems. J Geophys Res 115:G04003. https://doi.org/10.1029/2009JG001111
Lovenduski NS, Long CM, Lindsay K (2015) Natural variability in the surface ocean carbonate ion concentration. Biogeosciences 12:6321–6335. https://doi.org/10.5194/bg-12-6321-2015
Luo X, Wei H, Liu Z, Zhao L (2015) Seasonal variability of air–sea CO2 fluxes in the yellow and East China seas: a case study of continental shelf sea carbon cycle model. Cont Shelf Res 107:69–78
Mackenzie FT, Lerman A, Anderson AJ (2004) Past and present of sediment and carbon biogeochemical cycling models. Biogeosciences 1:11–32
Mellor GL (2001) An equation of state for numerical models of oceans and estuaries. Journal of Atmospheric and Oceanic technology 199(8):609–611
Miyazawa Y, Zhang R, Guo X, Tamura H, Ambe D, Lee JS, Yoshinari H, Setou T (2009) Water mass variability in the western North Pacific detected in a 15-year eddy resolving ocean reanalysis. J Oceanogr 65:737–756
Miyazawa Y, Yamashita N, Taniyasu S, Yamazaki E, Guo X, Varlamov MS, Miyama T (2014) Ocean dispersion simulation of perfluoroalkyl substances in the Western North Pacific associated with the great East Japan earthquake of 2011. J Oceanogr 70:535–547
Miyazawa Y, Varlamov MS, Miyama T, Guo X, Hihara T, Kiyomatsu K, Kachi M, Kurihara Y, Murakami H (2017) Assimilation of high-resolution sea surface temperature data into an operational nowcast/forecast system around Japan using a multi-scale three-dimensional variational scheme. Ocean Dyn 67:713–728
Nakaoka S, Telszewski M, Nojiri Y, Yasunaka S, Miyazaki C, Mukai M, Usui N (2013) Estimation temporal and spatial variation of ocean surface pCO2 in the North Pacific using a self-organizing map neural network technique. Biogeosciences 10:6093–6106. https://doi.org/10.5194/bg-10-6093-2013
Ohishi S, Aiki H, Tozuka T, Cronin FM (2019) Frontolysis by surface heat flux in the eastern Japan Sea: importance of mixed layer depth. J Oceanogr 75:283–297
Ohno Y, Kobayashi T, Iwasaka N, Suga T (2004) The mixed layer depth in the North Pacific as detected by Argo floats. Geophysical Research Letter 63:125–134
Oka E, Qui B (2012) Progress of North Pacific mode water research in the past decade. J Oceanogr 68:5–20. https://doi.org/10.1007/s10872-011-0032-5
Oka E, Suga T (2005) Differential formation and circulation of North Pacific central mode water. J Phys Oceanogr 35:1997–2011
Orr CJ, Epitalon MJ (2015) Improved routines to model the ocean carbonate system: mocsy 2.0. Geosci. Model Dev. 8:485–499. https://doi.org/10.5194/gmd-8-485-2015
Palmer JR, Totterdell IJ (2001) Production and export in a global ocean ecosystem model. Deep-Sea Research I 48:1169–1198
Qui B, Chen S (2006) Decadal variability in the formation of the North Pacific subtropical mode water; oceanic versus atmospheric control. Joural of Physical Oceanography 36:1365–1380
Sauzede, R.; Lavigne, H.; Claustre, H.; Schmechtig, C.; Ortenzio, D. F.; Guinet, C.; Pesant, S. Vertical distribution of chlorophyll-a concentration and phytoplankton community composition from in situ flurorescence profiles: a first database for the global ocean. Earth Syst Sci Data, 7, 2015, 261–273., doi: 10.5194/essd–7-261–2015
Stelmakh L, Gorbunova IT (2018) Carbon-to-chlorophyll-a ratio in the phytoplankton of the Black Sea surface layer: variability and regulatory factors. Ecologica Montenegria 17:60–73
Suga T, Motoki K, Aoki Y, MacDonald A (2004) The North Pacific climatology of winter mixed layer and mode waters. J Phys Oceanogr 34:3–22
Sutton, J. A.; Wanninkhof, R.; Sabine, L. C., Feely, A. R.; Cronin, F. M; Weller, A. R.; Variability and trends in surface seawater pCO2 and CO2 flux in the Pacific Ocean. Geophysical Research Letter, 2017, 44, 5627–5636, doi:https://doi.org/10.1002/2017GL073814
Takahashi T, Sutherland CS, Sweeney C, Poisson A, Metzl N, Tilbrook B, Bates N, Wanninkhof R, Feely AR, Sabin C, Olafsson J, Nojiri Y (2002) Global Sea–air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects. Deep-Sea Research II 49:1601–1622
Takahashi T, Sutherland CS, Wanninkhof R, Sweeney C, Feely AR, Chipman WD, Hales B, Friedrich G, Chaves F, Sbaine C, Watson A, Bakker CED, Schuster U, Metzl N, Yoshika-Wainoue H, Ishii M, Midorikawa T, Nojiri Y, Kortzinger A, Steinhoff T, Hoppema M, Olafsson J, Arnarson ST, Tilbrook B, Johannessen T, Olsen A, Bellerby R, Wong CS, Delille B, Bates NR, de Baar JWH (2009) Climatological mean and decadal change in surface ocean pCO2 and net sea–air CO2 flux over the global oceans. Deep-Sea Research II 56:554–577
Takahashi T, Sutherland CS, Chipman WD, Goddard GJ, Cheng H, Newberger T, Sweeney C, Munro RD (2014) Climatological distributions of pH, pCO2, total CO2, alkalinity, and CaCO3 saturation in the global surface ocean, and temporal changes at selected locations. Mar Chem 164:95–125
Takatani Y, Enyo K, Iida, Kojima A, Sasano D, Kosugi N, Midorikawa, Suzuki T, Ishii M (2014) Relationships between total alkalinity in surface water and seas surface dynamic height in the Pacific Ocean. Journal of Geophysical Research Oceans 119:2806–2814. https://doi.org/10.1002/2013JC009739
Tittensor D, Mora C, Jetz W, Lotze H, Ricard D, Vanden Berghe E, Worm B (2010) Global patterns and predictors of marine biodiversity across taxa. Nature 466(7310):1098–1101. https://doi.org/10.1038/nature09329
Tsuda, A.; Takeda, S.; Saito, H.; Nishioka, J.; Nojiri, Y.; Kudo, I.; Kiyosawa, H.; Imai, K.; Ono, T.; Shimamoto, A.; Tsumune, D.; Yoshimura, T.; Aono, T.; Hinuma, A.; Kinugasa, M.; Suzuki, K.; T. A mesoscale iron enrichment in the western subarctic Pacific induces large centric diatom bloom. Science, 2003, 958–961. doi: https://doi.org/10.1126/science.1082000
Wanninkhof R, Park HG, Takahashi T, Sweeney C, Feely R, Nojiri Y, Gruber N, Doney CS, McKinely AG, Lenton A, Le Quere C, Heinze C, Schwinger J, Garven H, Khatiwala S (2013) Global Ocean carbon uptake: magnitude, variability and trends. Biogeosciences 10:1983–2000. https://doi.org/10.5194/bg-10.1983-2013
Weiss R (1974) The provisions of social relationships. In: Rubin Z (ed) Doing unto others. Prentice Hall, Englewood Cliffs, pp 17–26
While J, Totterdell I, Martin M (2012) Assimilation of pCO2 data into a global coupled physical-biogeochemical ocean model. Journal of Geophysical Research 117:C03037. https://doi.org/10.1029/2010JC006815
Xiu P, Chai F (2013) Variability of oceanic carbon cycle in the North Pacific from seasonal to decadal scales. J. Geophys. Res Oceans 119:5270–5288. https://doi.org/10.1002/2013JC009505
Yasunaka S, Nojiri Y, Nakaoka S, Ono T, Mukai H, Usui N (2013) Monthly maps of sea surface dissolved inorganic carbon in the North Pacific: basin-wide distribution and seasonal variation. J Geophys Res Oceans 118:3843–3850. https://doi.org/10.1002/jgrc.20279,2013
Yasunaka S, Nojiri Y, Nakaoka S, Ono T, Whitney F, Telszewski M (2014) Mapping of sea surface nutrients in the North Pacific: basinwide distribution and seasonal to interannual variability. J. Geophys. Res. Oceans 119:7756–7771. https://doi.org/10.1002/2014JC010318
Acknowledgments
We acknowledge the continuing support from the Sasakawa Peace Foundation of the Ocean Policy Research Institute (OPRI-SPF). This study is a part of the Japan Coastal Ocean Predictability Experiment (JCOPE) sponsored by the Japan Agency for Marine–Earth Science and Technology (JAMSTEC). Comments from two anonymous reviewers helped to improve earlier versions of the manuscript. World Ocean Atlas 2001 (WOA01) and 2013 (WOA13) data were downloaded from websites of the US National Oceanographic Data Center (NODC): https://www.nodc.noaa.gov/OC5/WOA01/pr_woa01.html and https://www.nodc.noa.gov/OC5/woa13/, respectively. Monthly climatological data for DIN, DIP, and DIC from Yasunaka et al. (2013) were downloaded from http://soop.jp/index.html. In situ observational data from the Japan Meteorological Agency (JMA) were downloaded from https://www.data.jma.go.jp/kaiyou/db/vessel_obs/data-report/html/ship/ship.php. Partial pressure data for carbon dioxide from the Ayasato observatory were downloaded from https://ds.data.jma.go.jp/ghg/kanshi/obs/co2_monthave_ryo.html. Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua Ocean Color Data for 2015 were downloaded from the Physical Oceanography Distributed Active Archive Center (PODAAC) ftp site: ftp://podaac-ftp.jpl.nasa.gov/allData/modis/L3/aqua/4um/v2014.0/4km/daily/2015/.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
We calculated air–sea CO2 fluxes (Fig. 8(a–d)) as \( {\rho}_w{V}_p{k}_{0C}\left({P}_{{\mathrm{CO}}_2}-{P}_{{\mathrm{CO}}_2}^{\mathrm{atm}}\right) \) (Kantha, 2004; Ishizu et al. 2019), where ρw is the density of seawater (kg m−3); k0C is the solubility of CO2 in seawater (mol kg−1 atm−1; Weiss 1974); Vp is the piston velocity, which depends on wind speed (U10) and the Schmidt number (Sc) and is expressed in units of m s-1 (\( {V}_{\mathrm{p}}=8.61\times {10}^{-7}\ {U}_{10}^2\ {\left({S}_c/660\right)}^{-1/2} \); Wanninkhof, 1992); and \( {P}_{{\mathrm{CO}}_2} \) and \( {P}_{{\mathrm{CO}}_2}^{\mathrm{atm}} \) are the partial pressures of CO2 in surface waters and the atmosphere (μatm), respectively.
Compared with the global air–sea CO2 fluxes evaluated on the basis of monthly climatological data (Takahashi et al. 2002, 2009; Yasunaka et al. 2013), our model provides a finer spatiotemporal resolution of variability, with both strong and weak contrasts in air–sea CO2 fluxes in the target region (not shown). The seasonal changes in air–sea CO2 flux simulated by our model (Fig. 8(a−d); note that the values in Fig. 8(a−d) do not represent the actual air–sea CO2 flux, but the contribution of air–sea CO2 exchange to the total DIC balance) are similar to the seasonal climatology (Takahashi et al. 2002; Yasunaka et al. 2013), although the annual mean air–sea CO2 flux (Fig. 19) indicates a smaller sink in the Kuroshio Extension region between 30° N and 40° N, and the transition between sink and source areas is less distinct than that in the climatology (Takahashi et al. 2002; Yasunaka et al. 2013). For example, the climatology (Yasunaka et al. 2013) shows an extensive sink area between 30° N and 40° N of > 6 mmol m−2 C day−1, but our model predicts a mixture of sink and source areas in the range from − 5 to 5 mmol m−2 C day−1. One reason for the relatively weak sink areas predicted by our model may be uncertainty in the calculation of pCO2, as revealed by the MATLAB program CO2sys.m (Orr and Epitalon 2015; Mackenzie et al. 2004). The calculation of pCO2 performed here gives higher pCO2 values in the ocean compared with the pCO2 climatology (Nakaoka et al. 2013). While et al. (2012) also referred to large estimation in the calculation of pCO2.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ishizu, M., Miyazawa, Y., Tsunoda, T. et al. Seasonal variability in the inorganic ocean carbon cycle in the Northwest Pacific evaluated using a biogeochemical and carbon model coupled with an operational ocean model. Climatic Change 162, 877–902 (2020). https://doi.org/10.1007/s10584-020-02779-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10584-020-02779-2