This paper explores the impacts of land use change and basin topography on hydrological processes using the integrated surface-subsurface hydrological model (ISSHM), the SHUD model. The study conducted model simulations and analyses within the Shenzhen-Hong Kong basin. The experimental design is comprehensive, the arguments are credible, the charts are clear, and the writing is readable. However, I have some concerns that could make this study even more comprehensive and reliable.

1. The authors explored the relationship between slope and several hydrological processes. I did not see the process of slope calculation; therefore, I am uncertain if the slope here refers to the slope calculated from the DEM in different zones or the slope calculated based on the triangular mesh of the SHUD model. Since the resolution of the triangles is lower than that of the DEM, I believe that the slope based on triangles better reflects the impact of slope on hydrological processes. According to lines 187-189, the three zones are divided based on elevation differences, but the average slope attribute within the three zones is not displayed or discussed. Therefore, the article needs to clarify the slope calculation method.

2. Based on the SHUD model, the authors analyzed the impact of slope on hydrological processes (subsurface). Local subsurface flow is usually affected by soil characteristics, especially hydraulic conductivity, and is also influenced by slope runoff (slope + accumulating area). Therefore, the groundwater level is higher in the flat areas around the river channel, leading to a larger subsurface flow.

3. The SHUD model outputs subsurface flow, which may include the flow (Q1, Q2, Q3) in three directions of the triangle (positive outward) and the sum/net flow in the three directions (Qs = Q1 + Q2 + Q3, or eleqsub = eleqsub1 + eleqsub2 + eleqsub3). Then, the water balance of the saturated zone of this unit should be dS = Q1 + Q2 + Q3. In the long-term trend, dS should be close to zero. Here, dS is equivalent to Qs (eleqsub). Therefore, subsurface flow should be the sum of positive flows or the sum of negative flows. If the authors directly use new-flow (Qs/eleqsub) in the calculation, it no longer conforms to the meaning of subsurface flow. For example, the flow rates of the three sides of a triangle are eleqsub1 = 200, eleqsub2 = 100, and eleqsub3 = -310, then eleqsub = 200 + 100 + (-310) = -10. It indicates that the unit gains 10 units of water. In this case, the flow rate at this unit should be 300 (200 + 100) or -310, but definitely not -10. Similarly, the calculation of surface runoff may have the same concerns. I suggest that the authors clarify the reading and processing of variables in the appendix or supplementary materials.

4. Parameter calibration is a very challenging task for any ISSHMs. Although your study uses manual parameter calibration, I strongly recommend that you share more calibration details/experiences so that others can learn from your work. For example, how were sensitive parameters determined? How were parameters adjusted to gradually approach the observed results? How was it determined that the current parameter set reached "optimal" or "usable" levels?

5. Although the authors have clarified the data sources of this study, I strongly recommend that the authors share the research data of this article (including model input files, model source code used to implement this study, and calibrated parameters) in an open database (such as Zenodo). Of course, you can desensitize or not share meteorological and hydrological data that are restricted by copyright or confidentiality clauses. Sharing this data can help other researchers replicate your work and further expand the impact of your work. This is just a suggestion.

This article is a comparative analysis of two adjacent watersheds using the integrated surface-subsurface hydrological model (ISSHM) to examine the effects of topography and land use/land cover change (LULCC) on hydrological processes within the Greater Bay Area (GBA) of China. In general, the datasets used and methodology for the analyses are clear and appropriate. The results are well-organized. In the conclusion part, by mentioning the importance for hydrologic management strategies to consider the specific topography and LULCC characteristics of each watershed, the high-level significance of the study stands out.

However, if the significance of the study could be fully investigated with more in-depth thoughts about the results and conclusion, the publication would represent a more substantial contribution to scientific progress in the field.

The four scenarios using HLU and CLU seems intuitive and simple to apply, which might become a lightly novel or highlight of the methodology that could be potentially widely-applied. However, insufficient details about “all built-type land uses” and the process of “The HLU pattern involves reverting all built-type land uses in both watersheds to their pre-construction conditions” are provided. This is not addressed explicitly in the discussion or supplementary, which could become a valuable contribution to future studies.

Moderate to major revision is suggested to this manuscript.

Specific comments are listed below:

• Introduction:

1. L73-74: consider adding reference to support the statement “they are mainly based on single and spatially homogeneous watersheds”.
2. L77: may consider adding one or two sentences to describe “crucial economic zone” in China to highlight the importance of the study area.

• Methodology:

1. In section 3.4, the four scenarios using HLU and CLU seems intuitive and straightforward to apply, which might become a highlight of the methodology that could be potentially widely-applied. However, insufficient details about “all built-type land uses” and the process of “The HLU pattern involves reverting all built-type land uses in both watersheds to their pre-construction conditions” are provided. This is not addressed explicitly in the discussion or conclusion, and could be a valuable contribution in future studies.
2. L187: consider adding reference of 40 meter, at which Zone 1 and Zone 2 are divided.
3. L192-203: consider having a flowchart or a table instead of a long paragraph to demonstrate the key factors/processes, could also help the readers better understand the results.

• Results:

1. Figure 8: the text in the x and y-axis are too small to read.
2. Figure 9: the current figure is busy with many equations/texts embedded with the dots, may considering make it concise by leaving the p-values in the figure but putting the equations in the captions or texts.
3. Figure 10: similar to Figure 9

• Conclusion:

1. L373-379: the author may consider adding more specific examples or related refences about how hydrologic management or local watershed agencies could use this study to improve their methodology and strategies. Thus, the application of this publication could not only benefit not only the future studies in academia but also shed light on the practical water management or engineering world.

• Discussion (e.g., the Limitation and future work section in the current manuscript) can be substantially improved.

1. It lacks the in-depth analysis of evapotranspiration (ET), which correlates with climatic factors such as solar radiation, temperature, humidity, etc. It is worth noting that ET is different from surface runoff, subsurface flow, and infiltration. The author may consider adding sentences in the limitation/discussion section.
2. Groundwater dynamics was mentioned in other sections except the discussion part, considering address it in section 4.4.
3. L258-262: consider adding text about the statement that “topographic indices more accurately reflect hydrological responses under steady-state conditions” in discussion part.
4. Considering the study areas are the important economic zone, the author may consider relating it to other economic zone or highly-urbanized areas in other regions of the world. It may worth thinking and adding texts about how this publication could shed light on the practical water management or engineering world.

• The abstract is also suggested to revise based on the updated revision.