New algorithm unclogs major bottleneck in ocean geochemical and biogeochemical modelling

Posted by mmaheigan 
· Thursday, May 16th, 2024 

Numerical models are some of the principal tools for understanding the cycling of geochemical and biogeochemical tracers in the ocean, with the latter also being important components of the Earth System Models used to project future climate change. However, in order to use these models they must first be integrated to a seasonally-repeating equilibrium with minimal drift, a computationally expensive calculation that can take months on supercomputers given the long turnover timescale – many thousands of years – of the ocean. This “spin-up” problem has long been a major bottleneck in marine geochemical and biogeochemical modelling.

In a study published last year in J. Adv. Model. Earth Sys (2023, see reference below), a new algorithm was shown to speed-up by a factor of between 10-25 the spin-up of a wide range of geochemical tracers, such as radiocarbon, protactinium/thorium and zinc. It can be applied to any model in a “black box” manner.

Now, a follow up study published recently in Sci. Adv. (2024, see reference below) extends the previous results to complex marine biogeochemical models such as those used in the Coupled Model Intercomparison Project (CMIP) that underpin IPCC reports on climate change. The algorithm can accelerate the spin-up of seasonally-forced models by over an order of magnitude, and by a factor of 5 when driven with interannually forcing as is typical in CMIP simulations.

The ability to efficiently spin-up geochemical and biogeochemical models should enable their more effective use, for example making it feasible to calibrate models against observations and performing simulations at resolutions higher than has been previously possible.

Caption: Spin-up to equilibrium of the PISCES marine biogeochemical model. PISCES is coupled to the NEMO ocean circulation model and has 24 prognostic tracers. Left: Drift in dissolved inorganic nitrate concentration (mean squared difference at all grid points between consecutive years) as a function of time. Right: Globally-integrated air-sea flux of CO2 as a function of time. The solid horizontal line is the criterion for convergence established by the Ocean Model Intercomparison Project (OMIP). The blue lines are the conventional direct integration solution and the red lines the accelerated solution using the new algorithm.

This is a joint highlight with the GEOTRACES program.

Reference:

Khatiwala, S. (2024). Efficient spin-up of Earth System Models using sequence acceleration. Science Advances, 10. Access the paper: 10.1126/sciadv.adn2839

Khatiwala, S. (2023). Fast Spin‐Up of Geochemical Tracers in Ocean Circulation and Climate Models. Journal of Advances in Modeling Earth Systems, 15. Access the paper: 10.1029/2022ms003447

Filter by Keyword

abundance acidification additionality advection africa air-sea air-sea interactions algae alkalinity allometry ammonium AMO AMOC anoxic Antarctic Antarctica anthro impacts anthropogenic carbon anthropogenic impacts appendicularia aquaculture aquatic continuum aragonite saturation arctic Argo argon arsenic artificial seawater Atlantic atmospheric CO2 atmospheric nitrogen deposition authigenic carbonates autonomous platforms bacteria bathypelagic BATS BCG Argo benthic bgc argo bio-go-ship bio-optical bioavailability biogeochemical cycles biogeochemical models biogeochemistry Biological Essential Ocean Variables biological pump biophysics bloom blue carbon bottom water boundary layer buffer capacity C14 CaCO3 calcification calcite carbon carbon-climate feedback carbon-sulfur coupling carbonate carbonate system carbon budget carbon cycle carbon dioxide carbon export carbon fluxes carbon sequestration carbon storage Caribbean CCA CCS changing marine chemistry changing marine ecosystems changing marine environments changing ocean chemistry chemical oceanographic data chemical speciation chemoautotroph chesapeake bay chl a chlorophyll circulation CO2 coastal and estuarine coastal darkening coastal ocean cobalt Coccolithophores commercial community composition competition conservation cooling effect copepod copepods coral reefs CTD currents cyclone daily cycles data data access data assimilation database data management data product Data standards DCM dead zone decadal trends decomposers decomposition deep convection deep ocean deep sea coral denitrification deoxygenation depth diatoms DIC diel migration diffusion dimethylsulfide dinoflagellate dinoflagellates discrete measurements distribution DOC DOM domoic acid DOP dust DVM ecology economics ecosystem management ecosystems eddy Education EEZ Ekman transport emissions ENSO enzyme equatorial current equatorial regions ESM estuarine and coastal carbon fluxes estuary euphotic zone eutrophication evolution export export fluxes export production extreme events faecal pellets fecal pellets filter feeders filtration rates fire fish Fish carbon fisheries fishing floats fluid dynamics fluorescence food webs forage fish forams freshening freshwater frontal zone functional role future oceans gelatinous zooplankton geochemistry geoengineering geologic time GEOTRACES glaciers gliders global carbon budget global ocean global warming go-ship grazing greenhouse gas greenhouse gases Greenland ground truthing groundwater Gulf of Maine Gulf of Mexico Gulf Stream gyre harmful algal bloom high latitude human food human impact human well-being hurricane hydrogen hydrothermal hypoxia ice age ice cores ice cover industrial onset inland waters in situ inverse circulation ions iron iron fertilization iron limitation isotopes jellies katabatic winds kelvin waves krill kuroshio lab vs field land-ocean continuum larvaceans lateral transport LGM lidar ligands light light attenuation lipids low nutrient machine learning mangroves marine carbon cycle marine heatwave marine particles marine snowfall marshes mCDR mechanisms Mediterranean meltwater mesopelagic mesoscale mesoscale processes metagenome metals methane methods microbes microlayer microorganisms microplankton microscale microzooplankton midwater mitigation mixed layer mixed layers mixing mixotrophs mixotrophy model modeling model validation mode water molecular diffusion MPT MRV multi-decade n2o NAAMES NCP nearshore net community production net primary productivity new ocean state new technology Niskin bottle nitrate nitrogen nitrogen cycle nitrogen fixation nitrous oxide north atlantic north pacific North Sea nuclear war nutricline nutrient budget nutrient cycles nutrient cycling nutrient limitation nutrients OA observations ocean-atmosphere ocean acidification ocean acidification data ocean alkalinity enhancement ocean carbon storage and uptake ocean carbon uptake and storage ocean color ocean modeling ocean observatories ocean warming ODZ oligotrophic omics OMZ open ocean optics organic particles oscillation outwelling overturning circulation oxygen pacific paleoceanography PAR parameter optimization parasite particle flux particles partnerships pCO2 PDO peat pelagic PETM pH phenology phosphate phosphorus photosynthesis physical processes physiology phytoplankton PIC piezophilic piezotolerant plankton POC polar polar regions policy pollutants precipitation predation predator-prey prediction pressure primary productivity Prochlorococcus productivity prokaryotes proteins pteropods pycnocline radioisotopes remineralization remote sensing repeat hydrography residence time resource management respiration resuspension rivers rocky shore Rossby waves Ross Sea ROV salinity salt marsh satellite scale seafloor seagrass sea ice sea level rise seasonal seasonality seasonal patterns seasonal trends sea spray seawater collection seaweed secchi sediments sensors sequestration shelf ocean shelf system shells ship-based observations shorelines siderophore silica silicate silicon cycle sinking sinking particles size SOCCOM soil carbon southern ocean south pacific spatial covariations speciation SST state estimation stoichiometry subduction submesoscale subpolar subtropical sulfate surf surface surface ocean Synechococcus technology teleconnections temperate temperature temporal covariations thermocline thermodynamics thermohaline thorium tidal time-series time of emergence titration top predators total alkalinity trace elements trace metals trait-based transfer efficiency transient features trawling Tris trophic transfer tropical turbulence twilight zone upper ocean upper water column upwelling US CLIVAR validation velocity gradient ventilation vertical flux vertical migration vertical transport warming water clarity water mass water quality waves weathering western boundary currents wetlands winter mixing zooplankton

Copyright © 2025 - OCB Project Office, Woods Hole Oceanographic Institution, 266 Woods Hole Rd, MS #25, Woods Hole, MA 02543 USA Phone: 508-289-2838  •  Fax: 508-457-2193  •  Email: ocb_news@us-ocb.org

link to nsflink to noaalink to WHOI

Funding for the Ocean Carbon & Biogeochemistry Project Office is provided by the National Science Foundation (NSF) and the National Aeronautics and Space Administration (NASA). The OCB Project Office is housed at the Woods Hole Oceanographic Institution.