Title: Upscaling Methane Flux From Plot Level to Eddy Covariance Tower Domains in Five Alaskan Tundra Ecosystems

Spatial heterogeneity in methane (CH ₄ ) flux requires a reliable upscaling approach to reach accurate regional CH ₄ budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH ₄ flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH ₄ flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH ₄ flux in early growing seasons. The simulated monthly CH ₄ flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m ⁻²  s ⁻¹ and 0.782 using the DF algorithm, but 0.007 μmol m ⁻²  s ⁻¹ and 0.758 using HF and 0.007 μmol m ⁻²  s ⁻¹ and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH ₄ flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH ₄ flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH ₄ flux for all five tower domains despite relatively weak differences in simulated CH ₄ flux among three footprint algorithms. The CLM-Microbe model can simulate CH ₄ flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH ₄ flux in terrestrial ecosystems.