BGC-val: a model- and grid-independent Python toolkit to evaluate marine biogeochemical models
Lee de Mora1, Andrew Yool2, Julien Palmieri2, Alistair Sellar3, Till Kuhlbrodt4,Ekaterina Popova2, Colin Jones5, and J. Icarus Allen1
- 1Plymouth Marine Laboratory, Prospect Place, The Hoe, Plymouth, PL1 3DH, UK
- 2National Oceanography Centre, University of Southampton Waterfront Campus, European Way, Southampton SO14 3ZH, UK
- 3Met Office Hadley Centre, Exeter, EX1 3PB, UK
- 4NCAS, Department of Meteorology, University of Reading, Reading, RG6 6AH, UK
- 5NCAS, School of Earth and Environment, University of Leeds, Leeds, LS2 9JT, UK
The biogeochemical evaluation toolkit, BGC-val, is a model- and grid-independent Python toolkit that has been built to evaluate marine biogeochemical models using a simple interface. Here, we present the ideas that motivated the development of the BGC-val software framework, introduce the code structure, and show some applications of the toolkit using model results from the Fifth Climate Model Intercomparison Project (CMIP5). A brief outline of how to access and install the repository is presented in Appendix A, but the specific details on how to use the toolkit are kept in the code repository.
The key ideas that directed the toolkit design were model and grid independence, front-loading analysis functions and regional masking, interruptibility, and ease of use. We present each of these goals, why they were important, and what we did to address them. We also present an outline of the code structure of the toolkit illustrated with example plots produced by the toolkit.
After describing BGC-val, we use the toolkit to investigate the performance of the marine physical and biogeochemical quantities of the CMIP5 models and highlight some predictions about the future state of the marine ecosystem under a business-as-usual CO2concentration scenario (RCP8.5).
Big in the benthos: Future change of seafloor community biomass in a global, body size‐resolved model
Andrew Yool, Adrian P. Martin, Thomas R. Anderson, Brian J. Bett, Daniel O. B. Jones, Henry A. Ruhl
National Oceanography Centre, University of Southampton Waterfront Campus, Southampton, UK
Deep‐water benthic communities in the ocean are almost wholly dependent on near‐surface pelagic ecosystems for their supply of energy and material resources. Primary production in sunlit surface waters is channelled through complex food webs that extensively recycle organic material, but lose a fraction as particulate organic carbon (POC) that sinks into the ocean interior. This exported production is further rarefied by microbial breakdown in the abyssal ocean, but a residual ultimately drives diverse assemblages of seafloor heterotrophs. Advances have led to an understanding of the importance of size (body mass) in structuring these communities. Here we force a size‐resolved benthic biomass model, BORIS, using seafloor POC flux from a coupled ocean‐biogeochemistry model, NEMO‐MEDUSA, to investigate global patterns in benthic biomass. BORIS resolves 16 size classes of metazoans, successively doubling in mass from approximately 1 μg to 28 mg. Simulations find a wide range of seasonal responses to differing patterns of POC forcing, with both a decline in seasonal variability, and an increase in peak lag times with increasing body size. However, the dominant factor for modelled benthic communities is the integrated magnitude of POC reaching the seafloor rather than its seasonal pattern. Scenarios of POC forcing under climate change and ocean acidification are then applied to investigate how benthic communities may change under different future conditions. Against a backdrop of falling surface primary production (−6.1%), and driven by changes in pelagic remineralization with depth, results show that while benthic communities in shallow seas generally show higher biomass in a warmed world (+3.2%), deep‐sea communities experience a substantial decline (−32%) under a high greenhouse gas emissions scenario. Our results underscore the importance for benthic ecology of reducing uncertainty in the magnitude and seasonality of seafloor POC fluxes, as well as the importance of studying a broader range of seafloor environments for future model development.
How to cite
Yool A, Martin A, Anderson T, Bett B, Jones D, Ruhl H, (2017). Big in the benthos: Future change of seafloor community biomass in a global, body size-resolved model. Global change biology, 23 (9), pp. 3554-3566, 2017, https://doi.org/10.1111/gcb.13680
Improved Aerosol Processes and Effective Radiative Forcing in HadGEM3 and UKESM1
J. P. Mulcahy1, C. Jones2, A. Sellar1, B. Johnson1, I. A. Boutle1, A. Jones1, T. Andrews1, S. T. Rumbold3, J. Mollard4, N. Bellouin4, C. E. Johnson1, K. D. Williams1, D. P. Grosvenor2, and D. T. McCoy5
- 1Met Office, Exeter, UK,
- 2National Centre for Atmospheric Science, University of Leeds, Leeds, UK
- 3National Centre for Atmospheric Science, University of Reading, Reading, UK
- 4Department of Meteorology, University of Reading, Reading, UK
- 5Institute of Climate and Atmospheric Sciences, University of Leeds, Leeds, UK
Aerosol processes and, in particular, aerosol‐cloud interactions cut across the traditional physical‐Earth system boundary of coupled Earth system models and remain one of the key uncertainties in estimating anthropogenic radiative forcing of climate. Here we calculate the historical aerosol effective radiative forcing (ERF) in the HadGEM3‐GA7 climate model in order to assess the suitability of this model for inclusion in the UK Earth system model, UKESM1. The aerosol ERF, calculated for the year 2000 relative to 1850, is large and negative in the standard GA7 model leading to an unrealistic negative total anthropogenic forcing over the twentieth century. We show how underlying assumptions and missing processes in both the physical model and aerosol parameterizations lead to this large aerosol ERF. A number of model improvements are investigated to assess their impact on the aerosol ERF. These include an improved representation of cloud droplet spectral dispersion, updates to the aerosol activation scheme, and black carbon optical properties. One of the largest contributors to the aerosol forcing uncertainty is insufficient knowledge of the preindustrial aerosol climate. We evaluate the contribution of uncertainties in the natural marine emissions of dimethyl sulfide and organic aerosol to the ERF. The combination of model improvements derived from these studies weakens the aerosol ERF by up to 50% of the original value and leads to a total anthropogenic historical forcing more in line with assessed values.
How to cite
The Low‐Resolution Version of HadGEM3 GC3.1: Development and Evaluation for Global Climate
Till Kuhlbrodt1, Colin G. Jones2, Alistair Sellar3, Dave Storkey3, Ed Blockley3, Marc Stringer1, Richard Hill3, Tim Graham3, Jeff Ridley3, Adam Blaker4, Daley Calvert3, Dan Copsey3, Richard Ellis5, Helene Hewitt3, Patrick Hyder3, Sarah Ineson3, Jane Mulcahy3, Antony Siahaan4,6, and Jeremy Walton3
- 1National Centre for Atmospheric Science, Department of Meteorology, University of Reading, Reading, UK,
- 2National Centre for Atmospheric Science, School of Earth and Environment, University of Leeds, Leeds, UK,
- 3Met Office Hadley Centre, Exeter, UK,
- 4National Oceanography Centre, Southampton, UK,
- 5Centre for Ecology and Hydrology, Wallingford, UK,
- 6Now at the British Antarctic Survey, Cambridge, UK
A new climate model, HadGEM3 N96ORCA1, is presented that is part of the GC3.1 configuration of HadGEM3. N96ORCA1 has a horizontal resolution of ~135 km in the atmosphere and 1° in the ocean and requires an order of magnitude less computing power than its medium‐resolution counterpart, N216ORCA025, while retaining a high degree of performance traceability. Scientific performance is compared to both observations and the N216ORCA025 model. N96ORCA1 reproduces observed climate mean and variability almost as well as N216ORCA025. Patterns of biases are similar across the two models. In the northwest Atlantic, N96ORCA1 shows a cold surface bias of up to 6 K, typical of ocean models of this resolution. The strength of the Atlantic meridional overturning circulation (16 to 17 Sv) matches observations. In the Southern Ocean, a warm surface bias (up to 2 K) is smaller than in N216ORCA025 and linked to improved ocean circulation. Model El Niño/Southern Oscillation and Atlantic Multidecadal Variability are close to observations. Both the cold bias in the Northern Hemisphere (N96ORCA1) and the warm bias in the Southern Hemisphere (N216ORCA025) develop in the first few decades of the simulations. As in many comparable climate models, simulated interhemispheric gradients of top‐of‐atmosphere radiation are larger than observations suggest, with contributions from both hemispheres. HadGEM3 GC3.1 N96ORCA1 constitutes the physical core of the UK Earth System Model (UKESM1) and will be used extensively in the Coupled Model Intercomparison Project 6 (CMIP6), both as part of the UK Earth System Model and as a stand‐alone coupled climate model.
How to cite
New insight from CryoSat-2 sea ice thickness for sea ice modelling
David Schröder1, Danny L. Feltham1, Michel Tsamados2, Andy Ridout2, and Rachel Tilling3
- 1Centre for Polar Observation and Modelling, Department of Meteorology, University of Reading, Reading, RG6 6BB, UK
- 2Centre for Polar Observation and Modelling, Department of Earth Sciences, University College London, London, WC1E 6BT, UK
- 3Centre for Polar Observation and Modelling, School of Earth and Environment, University of Leeds, Leeds, LS2 9JT, UK
Estimates of Arctic sea ice thickness have been available from the CryoSat-2 (CS2) radar altimetry mission during ice growth seasons since 2010. We derive the sub-grid-scale ice thickness distribution (ITD) with respect to five ice thickness categories used in a sea ice component (Community Ice CodE, CICE) of climate simulations. This allows us to initialize the ITD in stand-alone simulations with CICE and to verify the simulated cycle of ice thickness. We find that a default CICE simulation strongly underestimates ice thickness, despite reproducing the inter-annual variability of summer sea ice extent. We can identify the underestimation of winter ice growth as being responsible and show that increasing the ice conductive flux for lower temperatures (bubbly brine scheme) and accounting for the loss of drifting snow results in the simulated sea ice growth being more realistic. Sensitivity studies provide insight into the impact of initial and atmospheric conditions and, thus, on the role of positive and negative feedback processes. During summer, atmospheric conditions are responsible for 50 % of September sea ice thickness variability through the positive sea ice and melt pond albedo feedback. However, atmospheric winter conditions have little impact on winter ice growth due to the dominating negative conductive feedback process: the thinner the ice and snow in autumn, the stronger the ice growth in winter. We conclude that the fate of Arctic summer sea ice is largely controlled by atmospheric conditions during the melting season rather than by winter temperature. Our optimal model configuration does not only improve the simulated sea ice thickness, but also summer sea ice concentration, melt pond fraction, and length of the melt season. It is the first time CS2 sea ice thickness data have been applied successfully to improve sea ice model physics.
How to cite
Representation of fire, land-use change and vegetation dynamics in the Joint UK Land Environment Simulator vn4.9 (JULES)
Chantelle Burton1,2, Richard Betts1,2, Manoel Cardoso3, Ted R. Feldpausch2, Anna Harper2,Chris D. Jones1, Douglas I. Kelley4, Eddy Robertson1, and Andy Wiltshire1
- 1Met Office Hadley Centre, Exeter, EX1 3PB, UK
- 2College of Life and Environmental Science, University of Exeter, Exeter, EX4 4SB, UK
- 3Brazilian Institute for Space Research (INPE), Earth System Science Center (CCST), São José dos Campos, Brazil
- 4Centre for Ecology and Hydrology, Wallingford, OX10 8BB, UK
Disturbance of vegetation is a critical component of land cover, but is generally poorly constrained in land surface and carbon cycle models. In particular, land-use change and fire can be treated as large-scale disturbances without full representation of their underlying complexities and interactions. Here we describe developments to the land surface model JULES (Joint UK Land Environment Simulator) to represent land-use change and fire as distinct processes which interact with simulated vegetation dynamics. We couple the fire model INFERNO (INteractive Fire and Emission algoRithm for Natural envirOnments) to dynamic vegetation within JULES and use the HYDE (History Database of the Global Environment) land cover dataset to analyse the impact of land-use change on the simulation of present day vegetation. We evaluate the inclusion of land use and fire disturbance against standard benchmarks. Using the Manhattan metric, results show improved simulation of vegetation cover across all observed datasets. Overall, disturbance improves the simulation of vegetation cover by 35 % compared to vegetation continuous field (VCF) observations from MODIS and 13 % compared to the Climate Change Initiative (CCI) from the ESA. Biases in grass extent are reduced from −66 % to 13 %. Total woody cover improves by 55 % compared to VCF and 20 % compared to CCI from a reduction in forest extent in the tropics, although simulated tree cover is now too sparse in some areas. Explicitly modelling fire and land use generally decreases tree and shrub cover and increases grasses. The results show that the disturbances provide important contributions to the realistic modelling of vegetation on a global scale, although in some areas fire and land use together result in too much disturbance. This work provides a substantial contribution towards representing the full complexity and interactions between land-use change and fire that could be used in Earth system models.
How to cite
Studying the impact of biomass burning aerosol radiative and climate effects on the Amazon rainforest productivity with an Earth system model
Florent F. Malavelle1, Jim M. Haywood1,2, Lina M. Mercado3,4, Gerd A. Folberth2,Nicolas Bellouin5, Stephen Sitch3, and Paulo Artaxo6
- 1CEMPS, University of Exeter, Exeter, EX4 4QE, UK
- 2UK Met Office Hadley Centre, Exeter, EX1 3PB, UK
- 3CLES, University of Exeter, Exeter, EX4 4RJ, UK
- 4Centre for Ecology and Hydrology, Wallingford, OX10 8BB, UK
- 5Department of Meteorology, University of Reading, Reading, RG6 6BB, UK
- 6Department of Applied Physics, Institute of Physics, University of São Paulo, São Paulo, Brazil
Diffuse light conditions can increase the efficiency of photosynthesis and carbon uptake by vegetation canopies. The diffuse fraction of photosynthetically active radiation (PAR) can be affected by either a change in the atmospheric aerosol burden and/or a change in cloudiness. During the dry season, a hotspot of biomass burning on the edges of the Amazon rainforest emits a complex mixture of aerosols and their precursors and climate-active trace gases (e.g. CO2, CH4, NOx). This creates potential for significant interactions between chemistry, aerosol, cloud, radiation and the biosphere across the Amazon region. The combined effects of biomass burning on the terrestrial carbon cycle for the present day are potentially large, yet poorly quantified. Here, we quantify such effects using the Met Office Hadley Centre Earth system model HadGEM2-ES, which provides a fully coupled framework with interactive aerosol, radiative transfer, dynamic vegetation, atmospheric chemistry and biogenic volatile organic compound emission components. Results show that for present day, defined as year 2000 climate, the overall net impact of biomass burning aerosols is to increase net primary productivity (NPP) by +80 to +105 TgC yr−1, or 1.9 % to 2.7 %, over the central Amazon Basin on annual mean. For the first time we show that this enhancement is the net result of multiple competing effects: an increase in diffuse light which stimulates photosynthetic activity in the shaded part of the canopy (+65 to +110 TgC yr−1), a reduction in the total amount of radiation (−52 to −105 TgC yr−1) which reduces photosynthesis and feedback from climate adjustments in response to the aerosol forcing which increases the efficiency of biochemical processes (+67 to +100 TgC yr−1). These results illustrate that despite a modest direct aerosol effect (the sum of the first two counteracting mechanisms), the overall net impact of biomass burning aerosols on vegetation is sizeable when indirect climate feedbacks are considered. We demonstrate that capturing the net impact of aerosols on vegetation should be assessed considering the system-wide behaviour.
How to cite
The Cost of Reducing the North Atlantic Ocean Biological Carbon Pump
- 1Plymouth Marine Laboratory, Plymouth, UK
- 2National Oceanography Centre, University of Southampton, Waterfront Campus, Southampton, UK
- 3Plymouth Marine Laboratory and National Centre for Earth Observation, Plymouth, UK
To predict the impacts of climate change it is essential to understand how anthropogenic change alters the balance between atmosphere, ocean, and terrestrial reservoirs of carbon. It has been estimated that natural atmospheric concentrations of CO2 are almost 200 ppm lower than they would be without the transport of organic material produced in the surface ocean to depth, an ecosystem service driven by mechanisms collectively referred to as the biological carbon pump. Here we quantify potential reductions in carbon sequestration fluxes in the North Atlantic Ocean through the biological carbon pump over the twenty-first century, using two independent biogeochemical models, driven by low and high IPCC AR5 carbon emission scenarios. The carbon flux at 1000 m (the depth at which it is assumed that carbon is sequestered) in the North Atlantic was estimated to decline between 27 and 43% by the end of the century, depending on the biogeochemical model and the emission scenario considered. In monetary terms, the value of this loss in carbon sequestration service in the North Atlantic was estimated to range between US$170–US$3000 billion in abatement (mitigation) costs and US$23–US$401 billion in social (adaptation) costs, over the twenty-first century. Our results challenge the frequent assumption that coastal habitats store more significant amounts of carbon and are under greater threat. We highlight the largely unrecognized economic importance of the natural, blue carbon sequestration service provided by the open ocean, which is predicted to undergo significant anthropogenic-driven change.
How to cite
Vegetation distribution and terrestrial carbon cycle in a carbon cycle configuration of JULES4.6 with new plant functional types
Anna B. Harper1, Andrew J. Wiltshire2, Peter M. Cox1, Pierre Friedlingstein1, Chris D. Jones2,Lina M. Mercado3,4, Stephen Sitch3, Karina Williams2, and Carolina Duran-Rojas1
- 1College of Engineering, Mathematics, and Physical Sciences, University of Exeter, Exeter EX4 4QF, UK
- 2Met Office Hadley Centre, Fitzroy Road, Exeter EX1 3PB, UK
- 3College of Life and Environmental Sciences, University of Exeter, Exeter EX4 4PS, UK
- 4Centre for Ecology and Hydrology, Wallingford OX10 8BB, UK
Dynamic global vegetation models (DGVMs) are used for studying historical and future changes to vegetation and the terrestrial carbon cycle. JULES (the Joint UK Land Environment Simulator) represents the land surface in the Hadley Centre climate models and in the UK Earth System Model. Recently the number of plant functional types (PFTs) in JULES was expanded from five to nine to better represent functional diversity in global ecosystems. Here we introduce a more mechanistic representation of vegetation dynamics in TRIFFID, the dynamic vegetation component of JULES, which allows for any number of PFTs to compete based solely on their height; therefore, the previous hardwired dominance hierarchy is removed.
With the new set of nine PFTs, JULES is able to more accurately reproduce global vegetation distribution compared to the former five PFT version. Improvements include the coverage of trees within tropical and boreal forests and a reduction in shrubs, the latter of which dominated at high latitudes. We show that JULES is able to realistically represent several aspects of the global carbon (C) cycle. The simulated gross primary productivity (GPP) is within the range of observations, but simulated net primary productivity (NPP) is slightly too high. GPP in JULES from 1982 to 2011 is 133 Pg C yr−1, compared to observation-based estimates (over the same time period) between 123 ± 8 and 150–175 Pg C yr−1. NPP from 2000 to 2013 is 72 Pg C yr−1, compared to satellite-derived NPP of 55 Pg C yr−1 over the same period and independent estimates of 56.2 ± 14.3 Pg C yr−1. The simulated carbon stored in vegetation is 542 Pg C, compared to an observation-based range of 400–600 Pg C. Soil carbon is much lower (1422 Pg C) than estimates from measurements (> 2400 Pg C), with large underestimations of soil carbon in the tropical and boreal forests.
We also examined some aspects of the historical terrestrial carbon sink as simulated by JULES. Between the 1900s and 2000s, increased atmospheric carbon dioxide levels enhanced vegetation productivity and litter inputs into the soils, while land use change removed vegetation and reduced soil carbon. The result is a simulated increase in soil carbon of 57 Pg C but a decrease in vegetation carbon of 98 Pg C. The total simulated loss of soil and vegetation carbon due to land use change is 138 Pg C from 1900 to 2009, compared to a recent observationally constrained estimate of 155 ± 50 Pg C from 1901 to 2012. The simulated land carbon sink is 2.0 ± 1.0 Pg C yr−1 from 2000 to 2009, in close agreement with estimates from the IPCC and Global Carbon Project.
How to cite
Technical note: A simple theoretical model framework to describe plant stomatal “sluggishness” in response to elevated ozone concentrations
Chris Huntingford1, Rebecca J. Oliver1, Lina M. Mercado2,1, and Stephen Sitch2
- 1Centre for Ecology and Hydrology, Benson Lane, Wallingford, Oxfordshire, OX10 8BB, UK
- 2College of Life and Environmental Sciences, University of Exeter, Amory Building, Rennes Drive, Exeter, EX4 4RJ, UK
Elevated levels of tropospheric ozone, O3, cause damage to terrestrial vegetation, affecting leaf stomatal functioning and reducing photosynthesis. Climatic impacts under future raised atmospheric greenhouse gas (GHG) concentrations will also impact on the net primary productivity (NPP) of vegetation, which might for instance alter viability of some crops. Together, ozone damage and climate change may adjust the current ability of terrestrial vegetation to offset a significant fraction of carbon dioxide (CO2) emissions. Climate impacts on the land surface are well studied, but arguably large-scale modelling of raised surface level O3 effects is less advanced. To date most models representing ozone damage use either O3 concentration or, more recently, flux-uptake-related reduction of stomatal opening, estimating suppressed land–atmosphere water and CO2 fluxes. However there is evidence that, for some species, O3 damage can also cause an inertial “sluggishness” of stomatal response to changing surface meteorological conditions. In some circumstances (e.g. droughts), this loss of stomata control can cause them to be more open than without ozone interference. To both aid model development and provide empiricists with a system on to which measurements can be mapped, we present a parameter-sparse framework specifically designed to capture sluggishness. This contains a single time-delay parameter , characterizing the timescale for stomata to catch up with the level of opening they would have without damage. The larger the value of this parameter, the more sluggish the modelled stomatal response. Through variation of , we find it is possible to have qualitatively similar responses to factorial experiments with and without raised O3, when comparing to reported measurement time series presented in the literature. This low-parameter approach lends itself to the inclusion of ozone-induced inertial effects being incorporated in the terrestrial vegetation component of Earth system models (ESMs).
How to cite
Modelling marine DOC degradation time scales
Luca Polimene1,∗, Richard B. Rivkin2, Ya-Wei Luo3, Eun Young Kwon4, Marion Gehlen5, M. Angelica Peña6, Nannan Wang3, Yantao Liang7, Hermanni Kaartokallio8 and Nianzhi Jiao3
- 1Plymouth Marine Laboratory, Prospect Place, UK
- 2Department of Ocean Sciences, Memorial University of Newfoundland, Canada
- 3State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, China
- 4Center for Climate Physics, Institute for Basic Science, South Korea
- 5Laboratoire des Sciences du Climat et de l’Environnement, Institut Pierre Simon Laplace, CEA-CNRS-UVSQ, France
- 6Institute of Ocean Sciences, Fisheries and Oceans Canada, Canada
- 7Key Laboratory of Biofuels, Shandong Provincial Key Laboratory of Energy Genetics, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, China
- 8Finnish Environment Institute, Marine Research Centre, Finland
Marine dissolved organic carbon (DOC) is formed of a large number of highly diverse molecules. Depending on the environmental conditions, a fraction of these molecules may become progressively resistant to bacterial degradation and accumulate in the ocean for extended time scales. This long-lived DOC (the so-called recalcitrant DOC, RDOC) is thought to play an important role in the global carbon cycle by sequestering carbon into the ocean interior and potentially affecting the climate. Despite this, RDOC formation is underrepresented in climate models. Here we propose a model formulation describing DOC recalcitrance through two state variables: one representing the bulk DOC concentration and the other representing its degradability (κ) which varies depending on the balance between the production of ‘new’ DOC (assumed to be easily degradable) and bacterial DOC utilization assumed to leave behind more recalcitrant DOC. We propose this formulation as a means to include RDOC dynamics into climate model simulations….
How to cite
Large but decreasing effect of ozone on the European carbon sink
Rebecca J. Oliver1, Lina M. Mercado1,2, Stephen Sitch2, David Simpson3,4, Belinda E. Medlyn5,Yan-Shih Lin5, and Gerd A. Folberth6
- 1Centre for Ecology and Hydrology, Benson Lane, Wallingford, OX10 8BB, UK
- 2College of Life and Environmental Sciences, University of Exeter, EX4 4RJ, Exeter, UK
- 3EMEP MSC-W Norwegian Meteorological Institute, PB 43, NO-0313, Oslo, Norway
- 4Dept. Space, Earth & Environment, Chalmers University of Technology, Gothenburg, SE-41296 Sweden
- 5Hawkesbury Institute for the Environment, Western Sydney University, Locked Bag 1797, Penrith NSW 2751 Australia
- 6Met Office Hadley Centre, Exeter, UK
The capacity of the terrestrial biosphere to sequester carbon and mitigate climate change is governed by the ability of vegetation to remove emissions of CO2 through photosynthesis. Tropospheric O3, a globally abundant and potent greenhouse gas, is, however, known to damage plants, causing reductions in primary productivity. Despite emission control policies across Europe, background concentrations of tropospheric O3 have risen significantly over the last decades due to hemispheric-scale increases in O3 and its precursors. Therefore, plants are exposed to increasing background concentrations, at levels currently causing chronic damage. Studying the impact of O3 on European vegetation at the regional scale is important for gaining greater understanding of the impact of O3 on the land carbon sink at large spatial scales. In this work we take a regional approach and update the JULES land surface model using new measurements specifically for European vegetation. Given the importance of stomatal conductance in determining the flux of O3 into plants, we implement an alternative stomatal closure parameterisation and account for diurnal variations in O3 concentration in our simulations. We conduct our analysis specifically for the European region to quantify the impact of the interactive effects of tropospheric O3 and CO2 on gross primary productivity (GPP) and land carbon storage across Europe. A factorial set of model experiments showed that tropospheric O3 can suppress terrestrial carbon uptake across Europe over the period 1901 to 2050. By 2050, simulated GPP was reduced by 4 to 9% due to plant O3 damage and land carbon storage was reduced by 3 to 7%. The combined physiological effects of elevated future CO2 (acting to reduce stomatal opening) and reductions in O3 concentrations resulted in reduced O3 damage in the future. This alleviation of O3 damage by CO2-induced stomatal closure was around 1 to 2% for both land carbon and GPP, depending on plant sensitivity to O3. Reduced land carbon storage resulted from diminished soil carbon stocks consistent with the reduction in GPP. Regional variations are identified with larger impacts shown for temperate Europe (GPP reduced by 10 to 20%) compared to boreal regions (GPP reduced by 2 to 8%). These results highlight that O3 damage needs to be considered when predicting GPP and land carbon, and that the effects of O3 on plant physiology need to be considered in regional land carbon cycle assessments.