Future CH4 budgets as modelled by a fully coupled Earth system model using prescribed GHG concentrations vs. interactive CH4 sources and sinks

We have used the NASA Goddard Institute for Space Studies (GISS) Earth system model GISS-E2.1 to study the future budgets and trends of global and regional CH4 under different emission scenarios. GISS-E2.1 is one of the few ESMs that can be driven by anthropogenic CH4 emissions, as well as interactive natural sources such as wetlands, and can simulate the tropospheric CH4 chemistry. In frame of the recent short-lived climate forcers (SLCFs) assessment by the Arctic Monitoring and Assessment Programme (AMAP), we used the GISS-E2.1 model with prescribed long-lived greenhouse gas (GHG) concentrations. In the present study, we have supplemented these simulations using the interactive CH4 sources and sinks in order to quantify the model performance and the sensitivity to CH4 sources and sinks. We have used the Current Legislation (CLE) and the Maximum Feasible Reduction (MFR) emission scenarios from the Eclipse V6b emission database to simulate the future chemical composition and climate impacts from 2015 to 2050. We have also simulated 1995-2014 in order to evaluate the model performance following the AMAP-SLCF protocol. The prescribed GHG version underestimates the Global Atmospheric Watch (GAW) surface CH4 observations during the period between 1995 and 2023 by 1% -8.4%-2.0%, with a correlation (r) of 0.71 -0.41 0.99. The largest underestimations are over the continental emission regions such as North America, Europe, and Asia, while biases are smallest over oceans. On the other hand, the simulation with interactive sources and sinks underestimates the GAW observations more than the prescribed simulation, by 18.5% -25% -10.4%, with a lower r of 0.36 -0.82 0.93. Opposite to the prescribed simulation, the biases are largest over oceans and smaller over the continents, however they are still larger over land than the prescribed simulation. The interactive simulation, with large sources virtually over land and strong sink over oceans, has a land/ocean ratio larger than 1 while the prescribed simulation has this ratio equal to 1 as it distributes the global prescribed CH4 concentration equally in longitude over a given latitude. This clearly shows that the interactive sources and sinks should be represented in models in order to realistically simulate the chemical composition and the oxidative capacity of the atmosphere. As expected, the MFR scenario simulates lower global surface CH4 concentrations and burdens compared to the CLE scenario, however in both cases, global surface CH4 and burden continue to increase through 2050 compared to present day. In the CLE scenario, increases are largest over the equatorial belt, in particular over India and East China, while the MFR scenario shows increases over the whole Southern Hemisphere, however much smaller compared to CLE. Finally, the interactive simulation shows that the chemical CH4 sink increases in the CLE scenario, while it slightly decreases in the MFR, leading to a larger CH4 lifetime in the MFR scenario compared to in the CLE scenario.