“…Moving forward, there is an opportunity to improve representation of storage and transmission of soil water in land models, while maintaining computational efficiency in mind, by (1) explicitly representing variably saturated flow, as is possible in the mixed form of Richards' equation [Celia et al, 1990;Maxwell and Miller, 2005;Kumar et al, 2009], to improve simulations of shallow groundwater dynamics; (2) explicitly represent airflow (vapor diffusion) through the soil to improve simulations of bare soil evaporation [Parker et al, 1987;Painter, 2011;Zeng et al, 2011;Smits et al, 2012]; and (3) explicitly represent macropore flow to simulate the nonuniform wetting of the soil matrix and the heterogeneity of flow paths at larger spatial scales Germann, 1981, 1982;SimunEk et al, 2003;Weiler, 2005;McDonnell et al, 2007;Maxwell and Kollet, 2008a;Nimmo, 2010;Yu et al, 2014]. Some integrated models have incorporated these processes-for example, the coupling of ParFlow and CLM represents threedimensional variably saturated flow and has now been applied at continental scales [Maxwell et al, 2015], CLM now parameterizes the diffusion of water vapor through a dry surface layer to improve simulations of bare soil evaporation [Swenson and Lawrence, 2014], and the LM3 model has a simple representation of macropores .…”