Nuclear warfare could have devastating impacts on millions of people, yet it has been suggested that regional or global nuclear conflict may be possible in the future (Toon et al., 2019). In addition to the calamitous impacts of nuclear conflict on a local level, research conducted with a range of climate models finds a global cooling in response to various conflict scenarios (Coupe et al., 2019; Malone et al., 1985; Mills et al., 2014; Pausata et al., 2016; Robock et al., 2007; Turco et al., 1983). This global cooling is driven by fires started by the nuclear weapons. These fires inject smoke into the upper troposphere, where rapid lofting can spread the sunlight-absorbing soot particles into the stratosphere (Turco et al., 1983). Recent research implies that even a small nuclear conflict may have impacts on the global climate system, affecting the state and circulation of the atmosphere (Robock et al., 2007), increasing the sea ice extent in both hemispheres (Mills et al., 2014), and reducing plant productivity and crop yields in regions far from the conflict location (Özdogan et al., ˘ 2013; Toon et al., 2019; Xia & Robock, 2013). While less studied, the potential impacts of nuclear conflict on the ocean are many. Numerous physical, chemical, and biological processes in the ocean are temperature dependent, and sunlight is a critical ingredient for photosynthesizing phytoplankton at the base of the marine food web. Using a climate model with an interactive ocean, Mills et al. (2014) evaluated the ocean physical response to a potential India/Pakistan nuclear war that lofts 5 Tg of black carbon particles into the stratosphere; they find a 0.8◦ C decrease in globally averaged sea surface temperature, with smaller temperature reductions at depth. Recently Toon et al. (2019) used an Earth system model that includes a representation for phytoplankton to evaluate the ocean biological response to nuclear conflict; they report a 5–15% decrease in phytoplankton productivity under a range of conflict scenarios. Such findings prompt further investigation into how nuclear conflict and the resulting global cooling may alter the chemical state of the ocean. Perturbations in the ocean's carbonate chemistry are of particular interest, owing to their importance for ocean acidification. Ocean acidification is an ongoing, large-scale environmental problem driven by fossil fuel emissions of carbon dioxide (CO2). Cumulatively since the preindustrial era, the ocean has absorbed 41% of the carbon emitted by human industrial activities (McKinley et al., 2017). While this ocean absorption of carbon has partially mitigated anthropogenic global warming, it has fundamentally altered the carbonate chemistry of the ocean, increasing the concentration of hydrogen ions ([H+]) while decreasing the concentration of carbonate ions ([CO2− 3 ]). Observations collected at time series sites across the global ocean find statistically significant reductions in the potential hydrogen (pH = −log([H+])) and the saturation state of the calcium carbonate mineral aragonite (Ωarag, which is proportional to [CO2− 3 ]) over the past few decades (Bates et al., 2014). These changes are a direct consequence of the ocean absorption of anthropogenic carbon; carbonate chemistry dictates that the excess carbon will react with water and CO2− 3 to decrease ocean pH and Ω (Feely et al., 2004). Both of these changes may have negative consequences for marine organisms, in particular for those that precipitate calcium carbonate shells (e.g., coccolithophores, pteropods, foraminifera, corals, molluscs, and echinoderms), as the precipitation is hindered by low pH, and because decreases in Ω favor shell dissolution (Doney et al., 2009). To date, there have been no studies of the effects of nuclear conflict on ocean acidification, though past modeling studies on the ocean's response to volcanic forcing and to proposed geoengineering schemes have intimated that ocean carbonate chemistry is highly sensitive to these types of external forcings. Using a fully coupled carbon-climate model, Frölicher et al. (2011) find that volcanic-induced cooling following the 1991 Mt. Pinatubo eruption led to immediate increases in the flux of carbon from atmosphere to ocean and consequently, increases in the total dissolved inorganic carbon (DIC) concentration in the surface ocean. Eddebbar et al. (2019) demonstrate that air-to-sea CO2 fluxes are significantly enhanced following the eruptions of Agung, El Chichón, and Pinatubo in a large ensemble of simulations with an Earth system model. Matthews et al. (2009) conduct solar radiation management climate engineering simulations with an intermediate complexity model of the coupled climate-carbon system; they find changes in ocean pH and Ωarag as a result of the anomalous cooling. Similarly, Lauvset et al. (2017) indicate that radiation management geoengineering leads to changes in North Atlantic pH in a fully coupled Earth system model, but they do not explore changes in Ωarag. While these studies are suggestive of the carbonate chemistry response to nuclear conflict, the external forcing perturbations are of a different magnitude and duration than those imposed by nuclear conflict. Further, it is difficult to mechanistically understand the ocean carbonate chemistry response to such external forcing perturbations in fully coupled models, where the terrestrial response to forcing additionally influences the atmospheric CO2 concentration. Here, we use a state-of-the art Earth system model to simulate the ocean carbonate chemistry response to a range of nuclear conflict scenarios. We decouple the ocean carbon cycle from that of the terrestrial carbon cycle via a direct prescription of the atmospheric CO2 boundary condition used for air-sea CO2 flux, that is, changes in the terrestrial biosphere have no influence on the atmospheric CO2 that the ocean sees. As we will demonstrate, we find large perturbations in ocean pH and Ωarag as a result of nuclear conflict. These perturbations have relatively long duration (order of 10 years) and are driven by decreases in temperature and subsequent increases in the ocean carbon inventory. 2. Methods We analyse output generated by the Community Earth System Model (CESM) version 1.3, a state-of-the-art coupled climate model consisting of atmosphere, ocean, land, and sea ice components (Hurrell et al., 2013). The atmosphere component of CESM in our simulations is the Whole Atmosphere Community Climate Model (WACCM; Marsh et al., 2013) with nominal 2◦ resolution, 66 vertical levels, and a model top at ∼145 km; it uses the Rapid Radiative Transfer Model for GCMs (RRTMG; Iacono et al., 2000) for the radiative transfer. The Community Aerosol and Radiation Model for Atmospheres (Bardeen et al., 2008) is coupled with WACCM to simulate the injection, lofting, advection, and removal of soot aerosols in the troposphere and stratosphere, and their subsequent impact on climate (Coupe et al., 2019; Toon et al., 2019). The ocean component of CESM is the Parallel Ocean Program version 2 (Danabasoglu et al., 2012) with nominal 1◦ resolution and 60 vertical levels. The biogeochemical ocean component of CESM is the Biogeochemical Elemental Cycling model that represents the lower trophic levels of the marine ecosystem, full carbonate system thermodynamics, air-sea CO2 fluxes, and a dynamic iron cycle (Doney et al., 2006; Moore et al., 2004, 2013; Moore & Braucher, 2008; Long et al., 2013; Lindsay et al., 2014). LOVENDUSKI ET AL. 2 of 9 Geophysical Research Letters 10.1029/2019GL086246 The ocean in the coupled CESM simulation is initialized from rest with World Ocean Circulation (WOCE) temperature and salinity (Gouretski & Koltermann, 2004). Biogeochemical tracers are initialized to observationally based climatologies where possible (Lauvset et al., 2016); where these were not available (such as dissolved iron and phytoplankton biomass), the model is initialized with fields interpolated from an existing CESM simulation. The new, fully coupled simulation was spun up for 4 years to an approximate steady state with a constant atmospheric CO2 mixing ratio of 370 ppm, representative of the mixing ratio in the year 2000. Due to the relatively short spin-up period, the globally integrated air-sea CO2 flux is not in steady state (drifting at a rate of 0.14 Pg C year−2) when the perturbation forcing is applied. We therefore present our results as anomalies from the drifting control integrations. Three control simulations of 20-year duration are generated using round-off level differences in atmospheric initial conditions. As each of these control simulations has different phasing of internal variability (e.g., El Niño-Southern Oscillation), we use the standard deviation across this ensemble to identify statistically significant perturbations due to nuclear conflict. We report on the anomalies generated from four simulations of nuclear conflict with varying amounts of soot injection: three India/Pakistan conflict scenarios that inject 5, 27, and 47 Tg of soot, respectively, and one US/Russia conflict scenario that injects 150 Tg of soot. The initial soot injection amounts are generated from plausible scenarios for nuclear conflict following advice from a number of military and policy experts; the reader is referred to Toon et al. (2019) for further details on scenario development. In each case, we prescribe that the conflict begins on 15 May of the 5th year of the first control simulation, and we integrate the model for a 15-year period following the injection. We assume that the smoke generated by mass fires from nuclear conflict is injected into the upper troposphere above the target sites (in the U. S./Russia case, smoke is spread evenly over the two nations), as in Toon et al. (2019). WACCM lofts much of this smoke higher into the stratosphere via solar heating of black carbon aerosols in the smoke, where the black carbon aerosols persist for about a decade. The resulting annual mean, post-conflict (May to the following April) anomalies in aerosol optical depth are shown in Figure 1a. These optical depth changes result in a 10–40% reduction in incoming solar energy (Toon et al., 2019). While we discuss the anomalies generated from all four of these conflict simulations, we describe two in greater detail throughout this manuscript: the U. S./Russia case, as it is the largest climate perturbation overall, and the India/Pakistan 47-Tg case, as it is the largest climate perturbation generated by a regional nuclear conflict. Ocean biogeochemistry in the version of CESM used for our simulations has been extensively validated in the literature (Brady et al., 2019; Freeman et al., 2018; Harrison et al., 2018; Krumhardt et al., 2017; Lindsay et al., 2014; Lovenduski et al., 2015, 2016; Long et al., 2013, 2016; Moore et al., 2013; McKinley et al., 2016; Negrete-García et al., 2019). Of particular note for our study, the simulated surface ocean carbonate ion concentration from a long, preindustrial control simulation of CESM compares favorably with reconstructed observations, albeit with lower interannual variance than has been measured at subtropical time series sites (Lovenduski et al., 2015). In Figure S1 in the supporting information, we illustrate the comparison between observationally based estimates of surface ocean pH and Ωarag (from GLODAPv2; Lauvset et al., 2016) and the CESM control ensemble mean. In this comparison, we note that the observational estimates have been extensively interpolated and are intended to represent year 2002 carbonate chemistry parameters, whereas CESM has been integrated under an atmospheric CO2 mixing ratio that corresponds to year 2000 forcing. We find high correspondence between the spatial patterns of modeled and observed pH and Ωarag, giving us confidence that CESM is capable of representing the mean state of these two variables. 3. Results Globally averaged surface ocean pH increases in response to each of the nuclear conflicts, where the magnitude of the pH anomaly scales with the amount of soot injected (Figure 1b). In each case, the pH anomaly exceeds the interannual standard deviation of pH in the control ensemble mean (gray shading in Figure 1b). We observe the largest increases in surface ocean pH in response to the U. S./Russia 150-Tg case; here the globally averaged surface ocean pH anomaly exceeds 0.05, corresponding to a ∼10% decrease in the global mean hydrogen ion concentration. Under each scenario, the pH anomaly peaks 2–4 years after the conflict and persists for ∼10 years. With the exception of the high-latitude oceans, the pH increase following the nuclear conflict is pervasive across the surface ocean (Figures 2a– 2c). In the 47-Tg India/Pakistan scenario, we observe local pH anomalies exceeding 0.06 units on average in years 2–5 post conflict (Figure 2c); the anomalies are largest in the North Atlantic, North Pacific, and Equatorial Pacific. These large, abrupt changes in surface ocean pH may have important consequences for calcifying organisms, as shell precipitation can be affected by the ambient hydrogen ion concentration in seawater (Kroeker et al., 2013). Since the beginning of the industrial revolution, global ocean pH has dropped by an estimated 0.1 units (Ciais & Sabine, 2013). The anomalies in pH generated by our simulations exceed 50% of this historical change and occur over a much shorter time period. Whether and how organisms respond to the initial and rapid alleviation of low pH, followed by an immediate return to the current pH state in the global ocean, is as yet unknown (see, e.g., Haigh et al., 2015). In contrast to our results for pH, we observe decreases in surface ocean Ωarag following nuclear conflict (Figure 1c), which should tend to inhibit the maintenance of shells and skeletons in calcified organisms. While minimal changes in Ωarag are simulated for the 5-Tg India/Pakistan case, the other three cases produce large decreases in saturation state, on the order of 0.1 to 0.3 units (Figure 1c). In each of these three cases, the anomalies exceed the interannual standard deviation of Ωarag in the control ensemble mean (gray shading in Figure 1c). The peak response in these three cases occurs 3–5 years post conflict, a year or so later than the pH response. While for pH the globally averaged anomaly is negligibly small, 10-years post conflict; anomalies in globally averaged Ωarag persist beyond our 15-year simulation time frame for all conflict scenarios. The decreases in aragonite saturation state span the tropics and subtropics, with the exception of the central and eastern Equatorial Pacific region (Figures 2d– 2f). Local decreases in saturation state exceed 0.5 units in the western North Atlantic and western North Pacific under the 47-Tg India/Pakistan scenario (Figure 2f). Importantly, the simulated decreases in saturation state are highly pronounced in regions that host diverse coral reef ecosystems (for instance, the western and southwestern Pacific and the Caribbean), and like pH, the changes in saturation state occur fairly rapidly. Projections from climate models suggest that coral reef ecosystems across the world will experience aragonite saturation state declines from their preindustrial value of 3.5 to 3.0 by the end of the century (Ricke et al., 2013); alarmingly, our simulations project similar Ωarag declines over a 3- to 5-year period, which then persist for years after the initial forcing dissipates. The opposite-signed anomalies in pH and Ωarag induced by nuclear conflict seem puzzling at first, as for "typical" anthropogenic ocean acidification scenarios, both of these variables simultaneously decrease. Why would nuclear conflict cause opposing responses in pH and saturation state? To understand these opposing responses, we need to consider the carbonate chemistry system in seawater and its sensitivity to changing temperature. Gaseous CO2 reacts with seawater to form carbonic acid (H2CO3), which then dissociates to form H+ and bicarbonate (HCO− 3 ). The hydrogen ion then reacts with CO2− 3 to form additional HCO− 3 , CO2 + H2O− ↽−−−−−−⇀−H2CO3. (1) H2CO3− ↽−−−−−−⇀−H+ + HCO− 3 . (2) H+ + CO2− 3 − ↽−−−−−−⇀−HCO− 3 . (3) The equilibrium constants for these reactions (typically expressed as K0, K1, and K2, respectively; Sarmiento & Gruber, 2006) are sensitive to changes in temperature, for example, the cooling induced by nuclear conflict. We need to also consider the dissolution reaction for mineral calcium carbonate (CaCO3) in seawater, CaCO3(s)− ↽−−−−−−⇀−Ca2+ sat + CO2− 3,sat, (4) where [Ca2+]sat and [CO2− 3 ]sat are the concentrations of dissolved calcium and carbonate in equilibrium with mineral CaCO3, and the solubility product (Ksp) for this reaction is also sensitive to temperature (Sarmiento & Gruber, 2006). Further, the saturation state for a calcium carbonate mineral in seawater (here: aragonite), can be expressed as Ωarag = [Ca2+][CO2− 3 ] Ksp , (5) where both [CO2− 3 ] and Ksp are affected by changes in temperature (Ca2+ is highly abundant in seawater, and thus changes in temperature do not affect its concentration enough to matter for CaCO3 dissolution; Emerson & Hedges, 2008; Sarmiento & Gruber, 2006). Thus, we can decompose the anomalies in pH and Ωarag into the component driven by temperature-induced changes in the carbonate chemistry equilibrium constants (K0, K1, K2, and Ksp) and the component driven by all other changes to the carbonate chemistry system, such as changes in the DIC concentration, the alkalinity, or the salinity. We approximate the temperature sensitivity of the equilibrium constants using a program developed for CO2 system calculations (CO2SYS; van Heuven et al., 2011) via finite difference approximation. The component driven by all other changes to the carbonate system is computed as the residual of the other two terms. The pH response to nuclear conflict is the sum of two opposing drivers: an increase in pH driven by a decrease in sea surface temperature that alters the carbonate chemistry equilibrium constants and a decrease in pH driven by an increase in the DIC concentration of the upper ocean. Figure 1b illustrates the temporal evolution of the components of the global pH anomalies from the India/Pakistan 47-Tg simulation driven by changes in the equilibrium constants versus all other changes in the carbonate chemistry system. The equilibrium constant-driven pH anomaly is positive, peaking 2–3 years after the conflict, whereas the “other” component of the pH anomaly is negative, peaking 3–5 years after the conflict. The resulting total pH anomaly is positive, indicating that it is more strongly influenced by changes in the equilibrium constants than other changes. In the India/Pakistan 47-Tg case, globally averaged temperature reaches a minimum 2 to 3-years post conflict; the model initially produces 3.5◦C–4◦C anomalies at the surface that rewarm toward pre-conflict values for the duration of the simulation (Figure 3a). In contrast, surface ocean salinity-normalized DIC anomalies peak 3 to 5-years post conflict (Figure 3b), mainly as a result of the enhanced solubility of CO2 in colder seawater. While decreasing biological export production also contributes to increased DIC in the surface ocean, this signal is small relative to the change driven by enhanced air-to-sea CO2 flux (e.g., Figure S2). The delay in DIC relative to temperature anomalies is a result of the long (order months to years) timescale for CO2 to fully equilibrate with the surface mixed layer (Emerson & Hedges, 2008). The cold, high DIC surface anomalies slowly propagate into the global ocean thermocline; we observe 1◦ C and 10 mmol m−3 anomalies in temperature and DIC, respectively, at a depth of 300 m that persist beyond the length of our simulation (Figure 3). As there are no significant anomalies in global mean alkalinity or salinity post conflict (not shown), we conclude that the DIC perturbation drives the “other” component of the pH anomalies. We find similar behavior for these components in the other conflict scenarios (not shown). The negative Ωarag anomalies post conflict are driven by a combination of lower temperatures and higher DIC concentrations. Colder surface temperatures tend to increase Ksp, while higher surface DIC concentrations tend to decrease [CO2− 3 ], resulting in lower Ωarag values post conflict. Figure 1c illustrates that the DIC (other) component dominates the total Ωarag anomaly for the India/Pakistan 47-Tg simulation. As for pH, the equilibrium constant component peaks earlier than the other component; this is due to the timing of the temperature and DIC perturbations (Figure 3). The spatial patterns of the post-conflict surface pH and Ωarag anomalies in the India/Pakistan 47-Tg scenario (Figures 2c and 2f) result from perturbations in local surface ocean temperature and DIC (Figure S3). Negative temperature anomalies and positive DIC anomalies are pervasive in the tropics and extratropics, with the exception of the eastern Equatorial Pacific, where a large and long-lasting El Niño-like event develops following the conflict (Coupe, et al., manuscript in review). This strong reduction in the equatorial trade winds greatly weakens upwelling in the cold tongue region, producing near-zero surface temperature anomalies and a reduction in vertical DIC supply here (Figure S3). In the Southern Ocean, temperature and DIC are not much affected by the nuclear conflict, likely a result of enhanced upwelling of warm water from the subsurface (Harrison, et al., manuscript in preparation). Taken together, the aforementioned changes in temperature and DIC lead to increases in pH and decreases in Ωarag over most of the ocean surface (Figure S4). The changes in surface ocean pH that we simulate for nuclear conflict resemble the simulated response of pH to volcanic eruptions, but are an order of magnitude larger. Figure S5 illustrates the anomaly in surface ocean pH in the first year following the eruptions of Agung, El Chichón, and Mt. Pinatubo, as estimated by the CESM Large Ensemble (Kay et al., 2015), which uses the same physical and biogeochemical ocean components as in our nuclear conflict simulations. The ensemble mean isolates the evolution of the Earth system under historical external forcing, including the aerosol loading following volcanic eruptions (Eddebbar et al., 2019), and averages across the various representations of internal variability (Deser et al., 2012; we note that ensembles are not necessary for the nuclear conflict scenarios since the much larger magnitude of forcing provides a higher signal-to-noise ratio). The anomaly in the ensemble mean shown here thus cleanly captures the response of surface ocean pH to volcanic eruptions. Here we show the anomaly in preindustrial pH (pH anomalies in equilibrium with preindustrial atmospheric CO2, which is computed simultaneously with contemporary pH at model run time), as the contemporary pH anomalies include also the response to increasing atmospheric CO2 from one year to the next. The similarity in the spatial patterns of volcanically induced pH anomalies and those produced under nuclear conflict is striking (cf. Figures S5 and 2c), suggesting that volcanic forcing produces similar temperature, DIC, and thus pH anomalies (including the El Niño-like response to volcanic forcing in the eastern Equatorial Pacific, described in Eddebbar et al., 2019). However, the eruption-driven pH anomaly is both smaller (an order of magnitude) and of shorter duration (∼2 years) than in the India/Pakistan 47-Tg simulation. Unfortunately, a similar analysis of volcanic Ωarag anomalies in the CESM Large Ensemble was not possible as preindustrial [CO2− 3 ] was not saved to disk. 4. Conclusions and Discussion We report on the surface ocean pH and Ωarag anomalies generated from four simulations of nuclear conflict using the CESM with full ocean carbonate system thermodynamics. Globally averaged surface ocean pH increases in response to each conflict, with the largest increases in the North Atlantic, North Pacific, and Equatorial Pacific Ocean. The pH anomalies persist for 10 years post conflict and are primarily driven by changes in the carbonate chemistry equilibrium constants as a result of decreases in sea surface temperature. In contrast, CESM simulates globally averaged decreases in surface ocean Ωarag in response to nuclear conflict, with the largest decreases in the tropics and subtropics. The Ωarag anomalies persist beyond the length of our 15-year simulations and are driven by a combination of changes in the carbonate chemistry equilibrium constants and the solubility-driven increases in DIC. We further demonstrate that the surface pH anomalies induced by nuclear conflict resemble those induced by volcanic eruptions in the same modeling system. The simulated changes in global and regional pH and Ωarag as a result of nuclear conflict are large and abrupt. In the most extreme forcing scenario (U. S./Russia 150 Tg), over a period of ∼5 years, global surface ocean pH increases by 0.06 units, and Ωarag decreases by 0.3 units. To put these numbers into perspective, this simulated rate of change of pH is 10 times larger than the rate of change we have observed over the past two decades as a result of ocean acidification (−0.0018 year−1; Lauvset et al., 2015). Worryingly, surface ocean Ωarag decreases more than six times faster than has been observed in the open ocean over the past three decades (−0.0095 year−1 at the Bermuda Atlantic time series; Bates et al., 2014). While the cooling associated with nuclear conflict rapidly and briefly alleviates the decline in pH associated with ocean acidification, the increase in solubility causes the ocean to absorb ∼11 Pg of excess carbon in a 10-year period, leading to a rapid drop in Ωarag. Whether and how calcifying organisms might respond to such rapid and opposing changes in pH and Ωarag is as yet unknown. In order to measure organism response to ocean acidification, a majority of laboratory studies perform CO2 bubbling perturbation experiments, which simultaneously decrease the pH and Ωarag in the surrounding seawater solution (Pörtner et al., 2014). This simultaneous change in two carbonate chemistry parameters challenges our ability to isolate the organism response to changes in pH or changes in Ωarag alone. A recent laboratory sensitivity study of marine bivalve larvae used chemical manipulation experiments to decouple these two parameters; they found that larval shell development and growth were negatively impacted by decreasing Ω and unaffected by changes in pH (Waldbusser et al., 2014). If these sensitivities are sustained in other organisms, we might conclude that calcifying organisms would be severely affected by nuclear conflict. Our findings shed light on the ocean biogeochemical response to other forms of extreme external forcing, such as volcanic eruptions (Eddebbar et al., 2019; Frölicher et al., 2011) and solar radiation management climate engineering (Lauvset et al., 2017; Matthews et al., 2009). They may further inform the study and understanding of the role of ocean acidification in marine extinction following the Chicxulub impact event (Henehan et al., 2019). Importantly, our results suggest that even a regional nuclear conflict can have an impact on global ocean acidification, adding to the list of the many, far-reaching consequences of nuclear conflict for global society.