Scatter plots demonstrating how the near-IR H2O line strength and the 10 μm silicate metric depend on the diffusivity parameter (K_zz; left), metallicity (Z; middle), and internal energy flux (T_int; right) model parameters. The size of the dots corresponds to the ${\stackrel{\sim }{\chi }}^2$ of their fits to the JWST observations. The gray rectangle corresponds to the JWST observational data and their 1σ uncertainties. The best-fit model run (star) is in 2σ agreement with the observation.