  Nonlinear acoustic evolution is often discussed in the context of
wave-steepening that leads to shock formation, and is of special interest in
applications where the shock continues to strengthen due to a narrowing of its
channel or the stratification of the medium. Accurate scalings govern low
amplitude waves and strong shocks, but connecting these phases, or describing
waves that are nonlinear from the outset, generally requires simulation. We
address this problem using the fact that waves within a plane-parallel,
isentropic and gravitationally stratified atmosphere are described by exact
simple-wave solutions, thanks to the conservation of Riemann invariants in a
freely falling reference frame. Our solutions enable us to discriminate waves
that reflect from those that form shocks, and to capture wave and shock
evolution using an ordinary differential equation. For several relevant values
of the adiabatic index $\gamma$ the solutions are explicit; furthermore,
nonlinear wave reflection from a free surface can be described analytically for
$\gamma=3$. Comparison to hydrodynamic simulations shows that our analytic
shock approximation is accurate up to moderate ($\sim$ few--15) Mach numbers,
where the accuracy increases with the adiabatic index. Our solutions also imply
that an initially subsonic pulse is unable to unbind mass from the atmosphere
without significantly increasing its entropy.
