In this study, the authors conducted a comprehensive investigation of the driving factors behind compound flooding in coastal cities using a combination of hydrodynamic modeling and mathematical statistics. Their research yielded insights into the impact of rainfall and tidal levels on compound flooding, as well as the contributions of different types of floods to compound flooding. The conclusions of this study are novel and have positive implications for risk management. The chosen topic also aligns with the scope of the journal "Hydrology and Earth System Sciences." , I have several suggestions and recommendations that would help enhance the manuscript.

Specific comments

1. It is necessary to indicate whether there are tide prevention facilities in the study area and how this factor is considered in the model.
2. In the "Model construction and validation" section, a one-dimensional and two-dimensional hydrodynamic coupling model was constructed, but the analysis of flood severity primarily relied on flood volume, with limited analysis on indicators such as flood area and flood depth. This aspect needs to be supplemented.
3. The paper only mentioned the length of the original data but did not elaborate on the sample data situation for constructing the Copula function.
4. The paper is lengthy, and to be concise, compressing the research methods section would be helpful. For example, common formulas in Correlation and Copula.
5. In the section "Spatial interaction of drainage units," the analysis of the interaction forces between different drainage units is not highly relevant to the main theme of this paper.
6. In section 4.4 "Causes and prevention measures of floods in drainage units," flood prevention measures should not be discussed in the research results. It is suggested to elaborate on them in the discussion section.
7. The conclusion needs further refinement.

Technical comments

L11: The word "Currently" is repeated.

L31: Change "in this year" to a specific year.

Figure 1: Mark all drainage outlets in the figure.

Table 1: Explain why the RMSE of the edge distribution function corresponding to the optimal tidal level for 3h is not the best.

L430-431: What is the specific relationship? It needs to be clarified.

L457-459: Drainage unit 14 appears in both cases simultaneously, please verify.

L460-461: Same as above.

L508-509: This is an observation, not a conclusion.

Figure 14: The discussion of the causes of sudden changes is crucial.

L520-522: This should be in the Methods section.

In this study, the authors thoroughly explored the factors driving compound flooding in an urban area by combining hydrodynamic modeling and multivariate statistics. Their research provided insights into the analysis of the relative contribution of rainfall and tidal levels on compound flooding and the roles of different flood scenarios. The study is fairly novel and findings are beneficial for flood risk management and applications. Further, the scope and findings of the study align with the scope of the journal "Hydrology and Earth System Sciences,” and I recommend its publication but after some revisions. Below, the authors can find my suggestions and recommendations, which I believe would contribute to improving the manuscript and clarifying the methodology.

Major Comments:

1. Copulas and marginal distributions are fitted to the “tidal levels.” Is this just the astronomical tide, or do these values include wind-driven surge (in the case of Typhoon/Cyclone conditions) and river discharge components as well? The authors mention “storm surge” as a critical factor for compound flooding several times (L346, L456, L564), but it’s not clear how storm surge is taken into account when modeling “tidal levels” in this study. If the tidal range is higher and the storm surge component is relatively very small (or if only the astronomical tide is considered), fitting an extreme value model might be fundamentally incorrect since the astronomical tide is deterministic. Therefore, a clear explanation of the terminology is needed. Furthermore, are the tidal boundary conditions used at the sluice gate (as outflow boundary conditions?)? The authors mention a tidal gate, but the location is not shown in the figure.
2. When selecting the sample for joint probability analysis, the tidal level is chosen as the "highest tidal level of the corresponding day." This method implies that the calculated tidal design levels could occur within a 24-hour period from the rainfall event. In the numerical modeling process, a single event is used as a representative for the tidal levels process, but the timing of the rainfall peak and tidal peaks in generating design scenarios is not explained. Additionally, figure 8(a) shows the peak rainfall and peak tide almost coinciding, which is rare in reality (not the most probable 200-year event). The 3-hour rainfall duration shows the peak rainfall occurring approximately half an hour after the peak tide, whereas the 12-hour rainfall duration shows the peak rainfall occurring approximately 3 hours before the peak tide. How are the corresponding time lags between peak rainfall and peak tide determined for each scenario? Are they coming from a probabilistic basis? This timing significantly impacts the resultant flood volume and might be the main reason for the variations discussed around L430.
3. 41 and 42 are used to evaluate the spatial interactions between drainage units (DUs). The equations only rely on the masses and distances of two drainage units. To clarify the conclusions around L445-L450, it is essential to understand what the spatial interaction force represents. For example, the spatial interaction force between DU22 and DU6 (which are far away from each other) is 0.75. Does this relatively higher value suggest a higher level of interconnectedness or influence on each other, than even adjacent drainage units? Please explain. Also, it is not clear how the spatial interaction of the drainage units relates to the main objectives of the study.
4. It is necessary to provide additional details on the types of data used in the study. For instance, specify whether rainfall data is gridded or gauge data, and if gauge data is used, indicate the number of gauges involved. What are the locations of the tidal data and rainfall used? If so, do you assume a uniformly distributed rainfall over the entire catchment?
5. The paper lacks a proper discussion about the main assumptions and limitations of the process. For instance, using only 16 years of data to estimate a 200-year design event introduces significant uncertainty in modeling the dependence structure and the tail distributions. Since the study proposes a "universal" method, it is crucial for readers to understand these assumptions and uncertainties, especially when the method is to be applied to other study sites.

Minor Comments

1. L11 “currently” repeats.
2. In the abstract, the authors mention the Kendall return period for the combined event of rainfall and tidal levels, greater than the “Or” return period and less than the “And” return period. The Kendall return period is typically expected to fall between the "Or" and "And" scenarios for combined events. It's not a significant finding of the study to be mentioned in the abstract.
3. L32, and L36, could you provide references?
4. Figure 2: what represents the purple arrow from correlation analysis to “effects of social factors on flood”?
5. L61, in underdeveloped areas, fluvial flooding may indeed be more prevalent due to the natural terrain and lack of infrastructure to mitigate such events. However, it is not accurate to say that "only" fluvial flooding exists in these areas. Could you explain or provide references?
6. Provide references for the “Pilgrim & Cordery rainfall model”, “Maximum Possible Weighting Function” and “Pilgrim & Cordery rainfall model”
7. Figure 6 seems incorrect. The theoretical distributions (GEV, Gamma) should be smooth unless you don’t change the distribution parameters.
8. What types of information does Figure 7 (b) provide? What are the insights we can gain into joint distribution probabilities? The color scale is also not given.
9. Although the peak tidal levels could be negative (depending on the datum used), the gamma distribution is selected to fit the peak tidal levels for some rainfall durations. The gamma distribution is lower bounded by zero thus making it impossible to sample negative realizations from the fitted distribution (if necessary).
10. The extreme sample selection is not well explained. Does it consider annual maxima? How to ensure the events are independent? What is the sample size?
11. Figure 10: The study area is in the legend (red). But cannot be seen in the figure.
12. Figure 12: Do these different colors represent any information? If not use a single color.
13. It's not clear how this observation is made “Despite differences in rainfall peak timing, rainfall, and tide peak timings leading to variations in calculated flooding volumes, the impact degree index is not affected by these differences”?
14. It's not clear how this conclusion is made “This study also highlights that under the same recurrence interval, rainfall events with larger peak timings are more destructive than those with earlier peak timings”? Have you run many events for the same return period by only changing the peak timing?
15. The authors suggest a set of equations (from eq. 30 to eq. 40), but poorly explained. Consider elaborating more about the governing process of equations. Additionally, it's mentioned that “calculate the range of rainfall and tidal level design values for different durations from 2-yr RP to 200-yr RP using the following formula” (L286,), but the formula doesn’t give the range of rainfall, and since it divides the range by 200yr Rainfall value. So, ∆X_t will be a dimensionless parameter that is related to the amount of variability (range) of rainfall. Same for Tides. Same in L304, ∆Vx_t and ∆Vy_t do not quantify the variation of flooding volume. Accordingly, check Dx_t and Dy_t.
16. The paragraph from L543 to 558 is more about discussing the results. Consider moving it into the discussion section.
17. The paper is lengthy. Consider moving some of the extra results to the supplementary materials. For example, Tables 1 and 2. Commonly used equations, such as those for correlation and distributions, can be omitted from the text. Instead, it is sufficient to cite the relevant references for these well-known equations.
18. Check the way of citing publications in the manuscript according to the HESS guideline (e.g., L40, L42, L50).