Ocean acidification (OA) generally receives far less consideration than other climate stressors and related hazards, such as global warming and extreme weather events. Canada is uniquely vulnerable to OA given its extensive coastal oceans, the oceanographic processes in its three basins, accelerated warming and sea-ice melt, and extensive coastal communities and maritime economic sectors. Canada’s coastline is also home to extensive and diverse First Nations peoples with distinct histories, rights, title, laws, governance and whose traditions and cultures are extrinsically linked to the sea. However, there are currently very limited pathways to support OA action, mitigation, and/or adaptation in Canada, particularly at the policy level. Here, we present a first synthesis of the current state of OA knowledge across Canada's Pacific, Arctic, and Atlantic regions, including monitoring, modelling, biological responses, socioeconomic and policy perspectives, and examples of existing OA actions and efforts at local and provincial levels. We also suggest a step-wise pathway for actions to enhance the coordinated filling of OA knowledge gaps and integration of OA knowledge into decision-making frameworks. The goals of these recommendations are to improve our ability to respond to OA in Canada, and minimize risks to coastal marine environments and ecosystems, vulnerable sectors, and communities.

Coastal waters are increasingly susceptible to low-oxygen conditions where circulation or nutrient inputs have been altered by anthropogenic activities, leading to increased occurrence of coastal hypoxia worldwide. Low-oxygen conditions also occur naturally in coastal settings, typically in bottom water or landward of sills where circulation is weak. Here we identify and characterize a persistent, low-oxygen layer in the mid-water column of a remote fjord of the northeast Pacific Ocean that uncharacteristically lacks an entrance sill. This feature is sustained by weak reoxygenation processes constrained by fjord and watershed morphology and geography, and expands seaward into a neighboring channel during summer and autumn. Marine CO2 system conditions characteristic of adverse habitat for local calcifiers also occurred within the low-oxygen layer, demonstrating that such features are multi-stressor habitats more typical of deeper bottom waters or eutrophied settings. Seasonal movement of this layer overlapped habitat margins for some sensitive local species, indicating risk from interannual variability or response to ongoing climate change. These results provide an improved understanding of fjord habitat that can be applied to other fjords where similar characterization of morphological features and water properties can be evaluated.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesise datasets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (E-FOS) are based on energy and cement production data. Emissions from land-use change (E-LUC) are estimated by bookkeeping models based on land-use data. The global atmospheric CO2 growth rate (G(ATM)) is computed from changes in concentration measured at surface stations. The global net uptake of CO2 by the ocean (S-OCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO(2)-products. The global net uptake of CO2 by the land (S-LAND) is estimated with dynamic global vegetation models. Additional lines of evidence are provided by atmospheric inversions, atmospheric oxygen measurements, ocean interior observation-based estimates, and Earth System Models. This year, we introduced corrections on the E-LUC, S-OCEAN and S-LAND estimates. The sum of all sources and sinks results in the carbon budget imbalance (B-IM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as +/- 1 sigma. For the year 2024, E-FOS increased by 1.1 % relative to 2023, with fossil emissions at 10.3 +/- 0.5 GtC yr(-1) (including the cement carbonation sink, 0.2 GtC yr(-1)), E-LUC was 1.3 +/- 0.7 GtC yr(-1), for total anthropogenic CO2 emissions of 11.6 +/- 0.9 GtC yr(-1) (42.4 +/- 3.2 GtCO(2) yr(-1)). Also, for 2024, G(ATM) was 7.9 +/- 0.2 GtC yr(-1) (3.73 +/- 0.1 ppm yr(-1)), 2.2 GtC above the 2023 growth rate. S-OCEAN was 3.4 +/- 0.4 GtC yr(-1) and S-LAND was 1.9 +/- 1.1 GtC yr(-1), leaving a large negative B-IM (-1.7 GtC yr(-1)), suggesting that the total sink or G(ATM) is strongly overestimated in 2024. The global atmospheric CO2 concentration averaged over 2024 reached 422.8 +/- 0.1 ppm. Preliminary data for 2025 suggest an increase in E-FOS relative to 2024 of +1.0 % (0.2 % to 1.7 %) globally, and atmospheric CO2 concentration increasing by 2.1 ppm reaching 425.6 ppm, 53 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959-2024, with a near-zero overall budget imbalance, although discrepancies of up to around 1 GtC yr(-1) persist for the representation of annual to decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows: (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the mean ocean sink.

Given the United Nations’ ambitious goal—endorsed by more than 100 nations, including Canada—of protecting 30% of the world’s marine ecosystems by 2030 (UN, 2023), the need to establish, expand, and track the effectiveness of protected areas is becoming more pressing each year. Effective planning for and monitoring of marine protected areas (MPAs) rely on the availability of quality baseline ecological and oceanographic information. Many of the sites in the Canadian Pacific Ocean best suited for protection, due to the confluence of ecological, cultural, and political significance, are remote. This poses two challenges: (1) there is little baseline information to guide planning, and (2) they are challenging to monitor once established.

Large changes in marine CO2 chemistry manifest in areas with weakly-buffered seawater where ocean acidification (OA) acts in concert with natural CO2 additions. These settings can exhibit periods of extreme OA in the form of multiple co-occurring stressors, including calcite undersaturation and low pH. Such conditions were observed in the northern Strait of Georgia, on the northeast Pacific coast, where extreme OA spanned a 3-year period. Here, we utilized an 8-year, highly-resolved record of seawater CO2 partial pressure and total dissolved inorganic carbon to decompose the drivers of this extreme OA. We find that variability in storm season intensity shaped the extent of conservative mixing and biogeochemical drivers such that manifests of extreme OA arise in this setting. Extreme OA manifested during years with weak storm seasons due to direct and indirect biogeochemical factors and the reduced impact of conservative mixing. This sensitivity to the storm season intensity highlights how vulnerable the northern Strait of Georgia is to subtle changes in environmental forcing and provides some predictive capacity for OA conditions over the coming year. These results illustrate that OA is not a “slow burn” process within weakly-buffered settings, but rather invokes periods of intensification with poorly understood biological implications.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC) are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The global net uptake of CO2 by the ocean (SOCEAN, called the ocean sink) is estimated with global ocean biogeochemistry models and observation-based fCO2 products (fCO2 is the fugacity of CO2). The global net uptake of CO2 by the land (SLAND, called the land sink) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The sum of all sources and sinks results in the carbon budget imbalance (BIM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2023, EFOS increased by 1.3 % relative to 2022, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (10.3 ± 0.5 GtC yr−1 when the cement carbonation sink is not included), and ELUC was 1.0 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1 ± 0.9 GtC yr−1 (40.6 ± 3.2 GtCO2 yr−1). Also, for 2023, GATM was 5.9 ± 0.2 GtC yr−1 (2.79 ± 0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 2.3 ± 1.0 GtC yr−1, with a near-zero BIM (−0.02 GtC yr−1). The global atmospheric CO2 concentration averaged over 2023 reached 419.31 ± 0.1 ppm. Preliminary data for 2024 suggest an increase in EFOS relative to 2023 of +0.8 % (−0.2 % to 1.7 %) globally and an atmospheric CO2 concentration increase by 2.87 ppm, reaching 422.45 ppm, 52 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2023, with a near-zero overall budget imbalance, although discrepancies of up to around 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the mean ocean sink. This living-data update documents changes in methods and datasets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2024 (Friedlingstein et al., 2024).

High mortality rates of cultured Pacific oysters ( Crassostrea gigas ) during the summer months have regularly occurred on oyster farms in British Columbia, Canada over the last 10 years, but little is known about the microbial and environmental conditions that contribute to such mortality events. The objective of the study was to determine correlative factors associated with the onset of a summer mortality event in oysters (mean ± SD shell height: 14.2 ± 0.5 mm) grown in suspended culture at four stocking densities (150, 300, 450, 600 oysters tray −1 ) from May 11 to September 17, 2018. Variables examined included both biotic (oyster growth, mortality, reproductive development, and microbiome (approximately every week); Vibrio and harmful algal species abundance) and abiotic (temperature, salinity, turbidity, dissolved oxygen, p CO 2 , pH, and aragonite saturation) ones. Both the absolute abundance of V. aestuarianus and the relative abundance of Vibrio spp. increased with observed oyster mortality and declining health. Mortality was highest on August 12 and associated with a prior period of elevated temperatures (i.e., increasing temperatures from early July to early August) and increased oyster growth/reproductive development. At that time, systemic mixed microbial infections and necrotic gill tissue in histological cross sections were observed in 19% of oysters that appeared healthy macroscopically. Cumulative percent mortalities per tray ranged from 34 to 75%, the highest-density trays having significantly less mortality and smaller shell width, shell length, and gonad length than lower-density trays. This study demonstrates the significant impact of summer mortality on Pacific oysters and highlights the biotic (host growth, reproductive development, and microbiome composition as well as Vibrio spp. abundance) and abiotic (water temperature) factors associated with the observed mortality in this region.

As elsewhere in the global ocean, the Gulf of Alaska is experiencing the rapid onset of ocean acidification (OA) driven by oceanic absorption of anthropogenic emissions of carbon dioxide from the atmosphere. In support of OA research and monitoring, we present here a data product of marine inorganic carbon chemistry parameters measured from seawater samples taken during biannual cruises between 2008 and 2017 in the northern Gulf of Alaska. Samples were collected each May and September over the 10 year period using a conductivity, temperature, depth (CTD) profiler coupled with a Niskin bottle rosette at stations including a long-term hydrographic survey transect known as the Gulf of Alaska (GAK) Line. This dataset includes discrete seawater measurements such as dissolved inorganic carbon and total alkalinity, which allows the calculation of other marine carbon parameters, including carbonate mineral saturation states, carbon dioxide (CO2), and pH. Cumulative daily Bakun upwelling indices illustrate the pattern of downwelling in the northern Gulf of Alaska, with a period of relaxation spanning between the May and September cruises. The observed time and space variability impart challenges for disentangling the OA signal despite this dataset spanning a decade. However, this data product greatly enhances our understanding of seasonal and interannual variability in the marine inorganic carbon system parameters. The product can also aid in the ground truthing of biogeochemical models, refining estimates of sea–air CO2 exchange, and determining appropriate CO2 parameter ranges for experiments targeting potentially vulnerable species. Data are available at https://doi.org/10.25921/x9sg-9b08 (Monacci et al., 2023).

Attribution of the ocean acidification (OA) signal in estuarine carbonate system observations is necessary for quantifying the impacts of global anthropogenic CO2 emissions on water quality, and informing managers of the efficacy of potential mitigation options. We present an analysis of observational data to characterize dynamics and drivers of seasonal carbonate system variability in two seagrass habitats of Puget Sound, WA, USA, and estimate how carbon accumulations due to anthropogenic CO2 emissions (C-anth) interact with these drivers of carbonate chemistry to determine seasonally resolved rates of acidification in these habitats. Three independent simulations of C-anth accumulation from 1765 to 2100 were run using two previously published methods and one novel method for C-anth estimation. Our results revealed persistent seasonal differences in the magnitude of carbonate system responses to anthropogenic CO2 emissions caused by seasonal metabolic changes to the buffering capacity of estuarine waters. The seasonal variability of pH(T) and pCO(2) is increased (while that of Omega(aragonite) is decreased) and acidification rates are accelerated when compared with open-ocean estimates, highlighting how feedbacks between local metabolism and C-anth can control the susceptibility of estuarine habitats to OA impacts. The changes in seasonal variability can shorten the timeline to exceedance of established physiological thresholds for endemic organisms and existing Washington State water quality criteria for pH. We highlight how C-anth estimation uncertainties manifest in shallow coastal waters and limit our ability to predict impacts to coastal organisms and ecosystems from anthropogenic CO2 emissions.

Abstract Arctic outflow winds bring cold air from the continent to the coastline through mountain passes. Using observational data and a 2‐D model, we show that a February 2019 outflow event caused the upper 100 m in Bute Inlet, British Columbia (within the traditional territory of the Homalco Nation) to cool up to 1.9°C and gain up to 4.1 mLL−1 of oxygen. The cold, oxygenated water persisted for almost 1 year within the 1,023–1,023.5 kgm−3 isopycnal range (∼50–150 m). Atmospheric (from 1929 to 2022) and oceanographic (from 1951 to 2022) data showed a statistically significant relationship between continental air temperature at Tatlayoko Lake and temperature and oxygen in Bute Inlet. This local mechanism that counters some effects of climate change could create a biological refugium as surrounding waters warm and lose oxygen at a faster rate. The number of outflow events decreased from 1951 to 2018, and increased since.

In recent years, mortality of cultured Pacific oysters (Crassostrea gigas) during July and August has regularly exceeded 50% in many oyster farms in Baynes Sound, British Columbia (BC), Canada. However, little is known about the bacterial diversity and environmental/biological conditions that may be contributing to such mortality events in the region. A study of Pacific oysters in intertidal culture trays was conducted during the summer of 2017 in Baynes Sound to monitor the occurrence of oyster mortality, while identifying potentially pathogenic bacteria present in the oysters, examining oyster reproductive development, and recording various environmental parameters (temperature, salinity, dissolved oxygen, chlorophyll a, carbon chemistry, phytoplankton species/abundance, zooplankton species/abundance) at three farm sites in Buckley Bay, Fanny Bay, and Ship's Point. Continuous mortality was observed throughout July and August, cumulative mortalities ranging from 16.5 to 35.1% at the three sites. Mortality over time was examined with a binomial generalized linear mixed-effects regression model, which demonstrated a correlation with date, tides, temperature, gonad somatic index, and gonad developmental phase (Rc2 = 0.598). Average cumulative mortality per module was highly correlated with average oyster gonad length and the proportion of female oysters in a multiple linear-regression model (p = 0.002, R2 = 0.824). A total of 158 bacterial isolates were cultured and identified based on sequencing of the recA gene. Among those isolates, Vibrio aestuarianus and V. harveyi, with 27 and three isolates, respectively, are well-documented pathogens of Pacific oysters in other regions. This study is the first to concurrently evaluate a broad array of factors potentially associated with Pacific oyster summer mortality in BC. We provide evidence that summer mortality of Pacific oysters in Baynes Sound is occurring with an etiology similar to summer mortalities of adult Pacific oysters reported in other countries.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (E-FOS) are based on energy statistics and cement production data, while emissions from land-use change (E-LUC), mainly deforestation, are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (G(ATM)) is computed from the annual changes in concentration. The ocean CO2 sink (S-OCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO(2) products. The terrestrial CO2 sink (S-LAND) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The resulting carbon budget imbalance (B-IM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as +/- 1 sigma. For the year 2022, E-FOS increased by 0.9% relative to 2021, with fossil emissions at 9.9 +/- 0.5 GtC yr(-1) (10.2 +/- 0.5 GtC yr(-1) when the cement carbonation sink is not included), and E-LUC was 1.2 +/- 0.7 GtC yr(-1), for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1 +/- 0.8 GtC yr(-1) (40.7 +/- 3.2 GtCO(2) yr(-1)). Also, for 2022, G(ATM) was 4.6 +/- 0.2 GtC yr(-1) (2.18 +/- 0.1 ppm yr(-1); ppm denotes parts per million), S-OCEAN was 2.8 +/- 0.4 GtC yr(-1), and S-LAND was 3.8 +/- 0.8 GtC yr(-1), with a B-IM of 0.1 GtC yr(-1) (i.e. total estimated sources marginally too low or sinks marginally too high). The global atmospheric CO2 concentration averaged over 2022 reached 417.1 +/- 0.1 ppm. Preliminary data for 2023 suggest an increase in E-FOS relative to 2022 of +/- 1:1% (0.0% to 2.1 %) globally and atmospheric CO2 concentration reaching 419.3 ppm, 51% above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959-2022, with a near-zero overall budget imbalance, although discrepancies of up to around 1 Gt Cyr(-1) persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living-data update documents changes in methods and data sets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2023 (Friedlingstein et al., 2023).

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1 (40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with a BIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (E-FOS) are based on energy statistics and cement production data, while emissions from land-use change (E-LUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (G(ATM)) is computed from the annual changes in concentration. The ocean CO2 sink (S-OCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (S-LAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (B-IM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as +/- 1 sigma. For the first time, an approach is shown to reconcile the difference in our E-LUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, E-FOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 +/- 0.5 GtC yr(-1) (9.3 +/- 0.5 GtC yr(-1) when the cement carbonation sink is included), and E-LUC was 0.9 +/- 0.7 GtC yr(-1), for a total anthropogenic CO2 emission of 10.2 +/- 0.8 GtC yr(-1) (37.4 +/- 2.9 GtCO(2)). Also, for 2020, G(ATM) was 5.0 +/- 0.2 GtC yr-1 (2.4 +/- 0.1 ppm yr(-1)), S-OCEAN was 3.0 +/- 0.4 GtC yr(-1), and S-LAND was 2.9 +/- 1 GtC yr(-1), with a B-IM of -0.8 GtC yr(-1). The global atmospheric CO2 concentration averaged over 2020 reached 412.45 +/- 0.1 ppm. Preliminary data for 2021 suggest a rebound in E-FOS relative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959-2020, but discrepancies of up to 1 GtC yr(-1) persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; Le Quere et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at (Friedlingstein et al., 2021).

Information on marine CO2 system variability has been limited along the northeast Pacific Inside Passage despite the region's rich biodiversity, abundant fisheries, and developing aquaculture industry. Beginning in 2017, the Alaska Marine Highway System M/V Columbia has served as a platform for surface underway data collection while conducting twice weekly ∼1600 km transits between Bellingham, Washington, and Skagway, Alaska. Marine CO2 system patterns were evaluated using measurements made over a 2-year period, which revealed the seasonal cycle as the dominant mode of temporal variability. The amplitude of this signal varied spatially and was modulated by the relative influences of tidal mixing, net community production, and the magnitude and character of freshwater input. Surface water pHT (total hydrogen ion scale) and aragonite saturation state (Ωarag) were determined using carbon dioxide partial pressure (pCO2) data with alkalinity derived from a regional salinity-based relationship, which was evaluated using intervals of discrete seawater samples and underway pH measurements. High-pCO2, low-pHT, and corrosive Ωarag conditions (Ωarag<1) were seen during winter and within persistent tidal mixing zones, and corrosive Ωarag values were also seen in areas that receive significant glacial melt in summer. Biophysical drivers are shown to dominate pCO2 variability over most of the Inside Passage except in areas highly impacted by glacial melt. pHT and Ωarag extremes were also characterized based on degrees of variability and severity, and regional differences were evident. Computations of the time of detection identified tidal mixing zones as strategic observing sites with relatively short time spans required to capture secular trends in seawater pCO2 equivalent to the contemporary rise in atmospheric CO2. Finally, estimates of anthropogenic CO2 showed notable spatiotemporal variability. Changes in total hydrogen ion content ([H+]T), pHT, and Ωarag over the industrial era and to an atmospheric pCO2 level consistent with a 1.5 ∘C warmer climate were theoretically evaluated. These calculations revealed greater absolute changes in [H+]T and pHT in winter as opposed to larger Ωarag change in summer. The contemporary acidification signal everywhere along the Inside Passage exceeded the global average, with specific areas, namely Johnstone Strait and the Salish Sea, standing out as potential bellwethers for the emergence of biological ocean acidification (OA) impacts. Nearly half of the contemporary acidification signal is expected over the coming 15 years, with an atmospheric CO2 trajectory that continues to be shaped by fossil–fuel development.

Assessing the status of ocean acidification across ocean and coastal waters requires standardized procedures at all levels of data collection, dissemination, and analysis. Standardized procedures for assuring quality and accessibility of ocean carbonate chemistry data are largely established, but a common set of best practices for ocean acidification trend analysis is needed to enable global time series comparisons, establish accurate records of change, and communicate the current status of ocean acidification within and outside the scientific community. Here we expand upon several published trend analysis techniques and package them into a set of best practices for assessing trends of ocean acidification time series. These best practices are best suited for time series capable of characterizing seasonal variability, typically those with sub-seasonal (ideally monthly or more frequent) data collection. Given ocean carbonate chemistry time series tend to be sparse and discontinuous, additional research is necessary to further advance these best practices to better address uncharacterized variability that can result from data discontinuities. This package of best practices and the associated open-source software for computing and reporting trends is aimed at helping expand the community of practice in ocean acidification trend analysis. A broad community of practice testing these and new techniques across different data sets will result in improvements and expansion of these best practices in the future.