<div><div>The upper ocean box model is forced by observed atmospheric pCO2 and temperature. It calculates the pCO2 and air-sea CO2 flux. <span style="font-size:13px">The model runs 1960-2017 (Dec 2017) . If linear atmospheric boundary condition is selected, this is applied only 1980-2017. The output that is the focus of the plotting is 1980-2017. </span>For complete documentation, including the equation being solved, see section 2.5 of McKinley et al. (2020). </div>
<div>
The zip file McKinley2020_boxmodel.zip contains code and documentation for the Upper Ocean Box Model along with information about how to make modifications to do sensitivity tests. See "Data Files" for access to McKinley2020_boxmodel.zip.<br />
</div>
<div>====================================================</div>
<div>Files included in the zip file McKinley2020_boxmodel.zip:</div>
<div> </div>
<div>2 main files</div>
<div>Ocean_Atm_Carbon_Box_McKinley2020.m</div>
<div>plot_boxmodel_McKinley2020.m</div>
<div> </div>
<div>3 helper routines</div>
<div>remove_fixed_trend.m</div>
<div>make_annual.m</div>
<div>calc_pCO2_McKinley2020.m </div>
<div> </div>
<div>1 data file for atmospheric pCO2</div>
<div>adjusted_atmpco2_fromML_v3.mat</div>
<div> </div>
<div>"README for Upper ocean diagnostic box model.rtf" (McKinley et al. 2020). Text reproduced below with additional information about data sources.</div>
<div> </div>
<div>====================================================</div>
<div>
<p>In Matlab2017a.</p>
<p>Coding and documentation by Galen A. McKinley, Columbia University / Lamont Doherty Earth Observatory. <a href="mailto:mckinley@ldeo.columbia.edu">mckinley@ldeo.columbia.edu</a></p>
</div>
<div>DESCRIPTION from "README for Upper ocean diagnostic box model.rtf"</div>
<div> </div>
<div>In Ocean_Atm_Carbon_Box_McKinley2020.m</div>
<div>Atmospheric pCO2 boundary conditions are loaded, and the linear forcing, as a change from change from Jan 1980 to Dec 2017 is calculated. </div>
<div>The magnitude and duration of El Chichon and Mt Pinatubo volcano temperature anomalies are set and plotted. </div>
<div>The model equation is integrated, including carbon chemistry at each timestep. </div>
<div> </div>
<div>In plot_boxmodel_McKinley2020.m, set key input variables:</div>
<div> On lines 28-33, set model temperature, overturning, depth of box, piston velocity (kw), deep abiotic DIC concentration, Alkalinity. </div>
<div> On line 37 set real or linear atmospheric pCO2 (do_linear_pco2 = 0, 1, respectively).</div>
<div> On line 38, set volcano forcing of SST or not (vary_temp = 1, 0, respectively). </div>
<div> </div>
<div>In plot_boxmodel_McKinley2020.m, set years to plot and to save on lines 44-45.</div>
<div> </div>
<div>In plot_boxmodel_McKinley2020.m, line 50, set to save out *mat files (tosave=1), into an existing directory. </div>
<div> </div>
<div>In plot_boxmodel_McKinley2020.m, plot results starting on line 90. Save results starting on line 164.</div>
<div> </div>
<div>====================================================</div>
<div> </div>
<div>PARAMETER INFORMATION</div>
<div> </div>
<div>List variables hard coded in Ocean_Atm_Carbon_Box_McKinley2020.m</div>
<div>% A % ocean Area</div>
<div>% rho % density (kg/m3), default is 1027</div>
<div>% timestep % choose 1/12 for monthly</div>
<div> </div>
<div>Input variables for each run of Ocean_Atm_Carbon_Box_McKinley2020.m</div>
<div>do_linear_pco2atm=0; % if = 1, will replace NOAA ESRL global pCO2atm with linear change from 1980-2018 </div>
<div>vary_temp=1; % if = 1, volcano imposed, mean T set by Temp; details in Ocean_Atm_Carbon_Box_McKinley2020.m</div>
<div> % if = 0, TempC throughout</div>
<div>% TempC % Temperature, set for box (C)</div>
<div>% nu_in % rate of ocean overturning (m3/yr)</div>
<div>% dz % Box depth (m)</div>
<div>% kw_in % piston velocity, cm/hr input</div>
<div>% deepDIC % deep ocean DIC </div>
<div>% ALK % Alkalinity, global constant, mmol/m3</div>
<div> </div>
<div>---</div>
<div>pCO2_atm (uatm) is read in from file adjusted_atmpco2_fromML_v4.mat (Mean adjusted in years 1959-1979 to merge from Mauna Loa only prior to 1980 to global after. Water vapor correction has been applied using global SST to find uatm.) Mauna Loa xCO2 data came from <a href="https://www.esrl.noaa.gov/gmd/ccgg/trends/data.html ">https://www.esrl.noaa.gov/gmd/ccgg/trends/data.html </a> from 1959-1979, and then ESRL global marine boundary layer <a href="https://www.esrl.noaa.gov/gmd/ccgg/mbl/data.php ">https://www.esrl.noaa.gov/gmd/ccgg/mbl/data.php </a> . From annual mean, we interpolate back to monthly to create a smoothly evolving timeseries that is deseasonalized. The slight adjutment is to bring the mean of Mauna Loa and the global MBL into agreement at the 1979-1980 transition.
<div> </div>
<div><span style="font-size:13px">There is additional detail on the data processing to derive pCO2 from NOAA xCO2 in supplementary of McKinley, G. A., Fay, A. R., Eddebbar, Y. A., Gloege, L. and Lovenduski, N. S.: External Forcing Explains Recent Decadal Variability of the Ocean Carbon Sink, AGU Advances, 1(2), 1, doi:10.1029/2019AV000149, 2020. </span>
<p>adjusted_atmpco2_fromML_v3.mat contains variables:</p></div>
<p> adjusted_atmpco2_1959_2018 60x1 480 double <br />
atm_pco2_wvcorr_global_AW 468x1 3744 double <br />
atm_pco2_wvcorr_global_AW_annual 39x1 312 double <br />
time_year 40x1 320 double</p>
</div>
<div> </div>
<div>Output from Ocean_Atm_Carbon_Box_McKinley2020.m</div>
<div>[pco2_ocean,pco2atm,trend_pco2atm,flux,time,dt]</div>
<div> </div>
<div>pco2_ocean % Global mean Ocean pCO2 (uatm)</div>
<div>pco2atm % Global mean atmospheric pCO2 (uatm)</div>
<div>trend_pco2atm % atmospheric pCO2 trend for 1980 to 2017 (uatm/yr)</div>
<div>flux % air-sea CO2 flux (PgC/yr)</div>
<div>time % time of simulation </div>
<div>dt % timestep (fraction of year)</div>
<div> </div>
<div>If plot_boxmodel_McKinley2020.m, these variables can be set there. Plotting of output is automated. </div>
<div>% TempC % Temperature, set for box (C)</div>
<div>% NU_IN % rate of ocean overturning (m3/yr)</div>
<div>% dz % Box depth (m)</div>
<div>% kw_in % piston velocity, cm/hr input</div>
<div>% deepDIC % deep ocean DIC </div>
<div>% ALK % Alkalinity, global constant, mmol/m3<br />
</div></div>
Upper Ocean Box Model
840371
Upper Ocean Box Model
2021-02-09T17:09:01-05:00
2021-02-09T17:09:01-05:00
2021-05-10T17:30:43-04:00
urn:bcodmo:dataset:840371
Upper Ocean Box Model which solves for the time change of Dissolved Inorganic Carbon (DIC) in single upper ocean box
The box model solves for the time change of Dissolved Inorganic Carbon (DIC) in single upper ocean box.
The upper ocean box model is forced by observed atmospheric pCO2 and temperature. It calculates the pCO2 and air-sea CO2 flux.
false
McKinley, G. A. (2021) Upper Ocean Box Model which solves for the time change of Dissolved Inorganic Carbon (DIC) in single upper ocean box. Biological and Chemical Oceanography Data Management Office (BCO-DMO). (Version 1) Version Date 2021-02-16 [if applicable, indicate subset used]. doi:10.26008/1912/bco-dmo.840371.1 [access date]
true
1
10.26008/1912/bco-dmo.840371.1
false
2021-02-16
Datapackage.json
Frictionless Data Package
https://www.bco-dmo.org/dataset/840371/datapackage.json
application/vnd.datapackage+json
PDF
https://www.bco-dmo.org/dataset/840371/Dataset_description.pdf
application/pdf
ISO 19115-2 (NOAA Profile)
https://www.bco-dmo.org/dataset/840371/iso
application/xml
http://www.isotc211.org/2005/gmd-noaa
Dublin Core
https://www.bco-dmo.org/dataset/840371/dublin-core
application/xml
http://purl.org/dc/elements/1.1/
840371
http://lod.bco-dmo.org/id/dataset/840371
1960 - 2017
1960
1960
2017
2017
OSPREY
http://www.opengis.net/def/crs/OGC/1.3/CRS84
<http://www.opengis.net/def/crs/OGC/1.3/CRS84> POLYGON((-180 -90, 180 -90, 180 -90, 180 90, -180 90, -180 -90))
-180
-90
180
90
0.000000000000
0.000000000000