Currently, an international network of operating high-resolution microbarographs was established to record wave-induced pressure variations at the Earth’s surface. Based on these measurements, simulations are performed to analyze the characteristics of waves corresponding to the observed variations of atmospheric pressure. Such a mathematical problem involves a set of primitive nonlinear hydrodynamic equations considering lower boundary conditions in the form of pressure variations at the Earth’s surface. Selection of upward propagating acoustic-gravity waves (AGWs) generated or reflected at the Earth’s surface requires the Neumann boundary conditions involving the vertical gradients of vertical velocity at the lower boundary. To analyze the correctness of the mathematical problem, linearized equations are used for small-surface wave amplitudes excited near the ground. Using the relation for wave energy, it is proven that the solution of the boundary problem based on the nondissipative approximation is uniquely determined by the variable pressure field at the Earth’s surface. The respective dissipative problem has also a unique solution with the appropriate choice of lower boundary conditions for temperature and velocity components. To test the numerical algorithm, solutions of the linearized equations for AGW modes are used. Developed boundary conditions are implemented into the model describing acoustic-gravity wave propagation from the surface atmospheric pressure source. Atmospheric waves propagating from the observed surface pressure variations to the upper atmosphere are simulated using the obtained algorithms and the computer codes.