The manuscript presents LM4-SHARC v1.0, a novel parameterization scheme that integrates catchment-scale soil, groundwater, and river interactions into the GFDL land model. This study addresses a critical limitation in current Earth System Models (ESMs) by improving hydrologic predictions and enhancing the representation of energy-water flux dynamics. The results, validated using observational data from Providence Creek, are compelling. I support publication with minor revisions. Below are my detailed comments:

General Comments

The manuscript's most significant issue lies in the methods section, which lacks clarity despite being central to the paper. From my understanding, the model uses baseflow observations and Equation (19) to derive parameters a and b, which are then used in Equations (20) and (21) to determine other parameters. During the simulation, Equation (1) calculates ql, and Equation (6) calculates Q. However, the purpose of the other equations remains unclear. The authors should explicitly outline the role of each equation, the model's calibration process, and the input-output structure. Additionally, the distinction between analytical and numerical solutions is difficult to grasp and warrants better explanation.

Furthermore, the section on energy transport methods could be shortened, as the hydrological components are more critical to the study.

Regarding model evaluation, while improvements are demonstrated for the Providence Creek catchment, they could stem from calibration against local streamflow data. Without such data in other regions, how can the authors ensure similar improvements? The discussion on future global tuning using remote sensing data is insufficient. Since this scheme is intended for global implementation, the manuscript should address how effective parameters can be derived at a global scale at this time. Expanding this section would significantly strengthen the manuscript.

Specific Comments

1. L86: "The properties of the groundwater were…"
   
   If soil and bedrock types are stable, why should groundwater parameters vary over time? The reasoning for developing time-varying parameters is not sufficiently convincing.
2. L111: "Then, the catchment-scale hydrologic structure…"
   
   What is meant by "inter-tile connection"? Has this been fully implemented in the land model?
3. L127: "Sierra National Forest…"
   
   What criteria led to the selection of this catchment?
4. L190-205:
   
   This paragraph is unclear. The relationship between Equations (4)(5)(6) and Equations (2)(19), as well as their roles within the model, should be explicitly explained.
5. L402: "The k was set to 1…"
   
   Why was k=1? Even in a headwater catchment, different types of hillslopes could exist.
6. L541: "Baseflow observation…"
   
   The baseflow observation method derived from Szilagyi and Parlange (1998) should include a discussion of its accuracy. Since baseflow is a minor yet essential component in model calibration, its uncertainty warrants consideration.
7. L576: "The evaluation was performed for four years…"
   
   Why wasn’t the model evaluated against streamflow/baseflow observations?

Conclusion

This manuscript is a valuable contribution to Earth system modeling. With improved clarity in the methods section, expanded discussion on global parameterization, and more detailed explanations for specific concerns, it would be well-suited for publication.

This paper presents a new parameterization scheme that cooperate the groundwater component to the ESM model, which is a good starting point of considering the effect of groundwater on the energy component in the ESM model. 

In general, the manuscript is written in very good and comprehensive english and it is easy to read and follow. However, I find the manuscript is too long and with some redundant information that reduces the readability. I do agree with the first reviewer that the part describing the energy transport makes the manuscript discontinued, which might be better to put to as a supplementary. I believe with some more efforts for revision, this manuscript could reach the standard of a publication. 

Here follows some specific comments. 

1. Figure 1 appears to be difficult to understand. especially when it appears, there was no proper explanation of what is a "height band" until I read section 2.4. I suggest to have a good explain of  the following terms at the begining: what is soil columns, tiles, height band and characterisitc hillslope, repectively and what are their potential relationship.  and also the meaning of the black dots and what does the size of the dots mean in each soil layer is not clear. 

2. line 183 states rho is r_j ^j , which does not make sense.  and if rho is density, then could you check the unit consisensy of this equation?

3. line 196 eq 6 seems is the spatial derivative of the "stream" discharge instead of the time derivative of the "steam" discharge?  195 has typo "steam discharge" as well. 

4. Figure 3: the analytical area is not explained anywhere (or I missed it), if you present it there, it is better to give some explanations. 

5. Line492 "Following the extraction criteria(section 3.3.2)....  there is section 3.3.2

6. Figure 7 (a)(b)(c): the first subfigure with legend "VWC", which is not explained anywhere. according to the caption it is the soil moisture, then should it SMC? (according to line 581). And what is the meaning of fprec in the third figure? how to interpret it?

7. Line 636-648 + Figure 9 : It is not clear why reducing spin-up time matters. it is an interesting discovery, but I don't find it important.  Do you have any important reason to keep this discussion in the main body of the manuscript? Otherwise consider move it as a supplementary mateiral. 

8. Section 4.4 is very important since it finally hits the point that the modeled LAI has difference after involving the groundwater to the model, which lead to the potential improvement of the modeled transpiration. It would be better if there are more statistics given to show how much improvement between the lm4 and lm4-sharc. 

9. I see all the catchments the authors selected are headwater catchment. so before discussing applying the sharc scheme to global scale, could the author elaborate more on how to apply the scheme to downstream catchments where the drainage networks are more complicated?  I think the accuracy of the analytical solution of those catchments should be more challenging.