A two-dimensional computational ELF technique has been applied to interpret variations of peak frequencies of Schumann resosnances (SR) observed at the MIT experimental site (West Greenwich, Rhode Island). The spatial-temporal dynamics of global thunderstorm activity has been simulated on the basis of general meteorological data. It is shown that the proposed models provide a reasonable qualitative agreement between computed and observed variations for SR I to IV. Some inverse task diagrams has been constructed as an instrument for distinguishing between day-to-day thunderstorm scenarios.