The connection between ocean environmental parameters and transmission loss (TL) can be modeled using sound propagation models. We can analyze and improve these models using information theoretic tools such as the Fisher information matrix (FIM) and methods for model reduction, where we identify parameters that can be removed from a model without sacrificing accuracy. Critical to these methods for model reduction is the ability to evaluate derivatives of the model predictions of TL with respect to model parameters. Calculating derivatives of TL models using finite difference methods with sufficient accuracy is challenging due to the complexity of sound fields in an ocean environment; instead, we train a surrogate machine learning model to which we can apply automatic differentiation (AD). We propose a general method for model reduction in which a ML model is used as a surrogate model for physical sound propagation ocean models, which is then used to find potential simplifications of the original model. We demonstrate this method using the Pekeris waveguide model of the ocean, which calculates underwater acoustic TL in an ocean environment due to seafloor characteristics. [Work supported by Office of Naval Research]