Abstract. An improved divergence method has been developed to estimate annual methane (CH4) emissions from TROPOspheric Monitoring Instrument (TROPOMI) observations. It has been applied to the period of 2018 to 2021 over the Middle East, where the orography is complicated, and the mean mixing ratio of methane (XCH4) might be affected by albedos or aerosols over some locations. To adapt to extreme changes of terrain over mountains or coasts, winds are used with their divergent part removed. A temporal filter is introduced to identify highly variable emissions and further exclude fake sources caused by retrieval artifacts. We compare our results to widely used bottom-up anthropogenic emission inventories: Emissions Database for Global Atmospheric Research (EDGAR), Community Emissions Data System (CEDS) and Global Fuel Exploitation Inventory (GFEI) over several regions representing various types of sources. The NOX emissions from EDGAR and Daily Emissions Constrained by Satellite Observations (DECSO), and the industrial heat sources identified by Visible Infrared Imaging Radiometer Suite (VIIRS) are further used to better understand our resulting methane emissions. Our results indicate possibly large underestimations of methane emissions in metropolises like Tehran (up to 50 %) and Isfahan (up to 70 %) in Iran. The derived annual methane emissions from oil/gas production near the Caspian Sea in Turkmenistan are comparable to GEFI but more than two times higher than EDGAR and CEDS in 2019. Large discrepancies of distribution of methane sources in Riyadh and its surrounding areas are found between EDGAR, CEDS, GFEI and our emissions. The methane emission from oil/gas production in the east to Riyadh seems to be largely overestimated by EDGAR and CEDS, while our estimates, and also GFEI and DECSO NOX indicate much lower emissions from industry activities. On the other hand, regions like Iran, Iraq, and Oman are dominated by sources from oil and gas exploitation that probably includes more irregular releases of methane, with the result that our estimates, that include only invariable sources, are lower than the bottom-up emission inventories.