A mathematical model of lithosphere–atmosphere coupling for seismic events