Methane (CH₄) is the second most important greenhouse gases just behind carbon dioxide. Recent studies showed that half of global CH₄ emissions could come from highly variable aquatic ecosystems, including substantial emissions from lakes. However, lake CH₄ dynamics has not yet been represented in Earth system models and thus the feedback between lake methane and climate is largely unknown. In this study, we developed a lake biogeochemical model in Energy Exascale Earth System Model (E3SM) that represents CH₄ production, CH₄ oxidation and CH₄ transport through diffusion and ebullition. We used a surrogate-model based parameter calibration method to analyze the sensitive model parameters and calibrate them simultaneously. The model is validated against a global lake CH₄ dataset that is compiled from literature reviews and includes data from over 100 lakes larger than 0.1 km² across the globe. In addition, site-level model data comparison at sites with detailed CH₄ dynamics data is also performed. We also test the possibility to improve the model performance by representing CH₄ production in oxic water layers and using alternative algorithms for CH₄ oxidation. The developed model may enable realistic predictions of the global CH₄ cycle under climate change.