I thank the authors for their effort to improve the manuscript. Some of my concerns on the flux calculation, however, remain. The issue is that the measurement set-up is inadequate to calculate reliable fluxes, which is especially concerning during the nighttime. The authors recognize this themselves by stating in Section 3.5 that at least they are sure of the direction of the flux, but should discuss this in a more structured and defined way. Moreover, they should emphasize what this means for the interpretability of their results and also recognize this in the conclusions. As the direction of the flux has more confidence than the size of the flux, the authors could place more emphasis on the gradients rather than the flux. They can e.g. evaluate correlations such as in Figures 5 and 6 for the gradients. Please find below my remaining concerns on the methodology, and some specific comments.
Comments on methodology
Though the used alpha-factor is a good first order estimate for the correction, in reality the alpha-factor in not only height-dependent, but is also affected by canopy structure (Melman et al., 2024). The authors should indicate in Section 3.5 that the value of 0.75 is a first order of estimate, arbitrarily based on the height above the canopy. Secondly, the effect of the RSL on the AGM is different under leafy or leafless conditions (Shapkalijevski et al., 2016) (which may also play a role to explain Figure 8 and the phenomena on lines 323-326). I realize that the measurement set-up is not suited to derive an alpha factor for the two different periods, but the authors should at least discuss this in Section 3.5, and note that this causes a large uncertainty in the comparison of fluxes between the leafy and leafless periods.
The authors indicate in point 1.2 of Referee #1 that they calculate the transfer velocity at 10 minute intervals. However, choosing a short period will lead to a sampling error of large eddies, which may lead to an underestimation of the flux, especially if measured high above the surface where large eddies dominate (e.g. when measuring over forest) (Finkelstein and Sims, 2001). Usually, turbulent statistics are calculated at 30 minute intervals to prevent this error from getting too large, why did the authors choose to calculate D at 10 minute intervals? Moreover, did they apply any additional processing steps, including coordinate rotation, detrending (Gash and Culf, 1996), statistical tests and raw data screening (Vickers and Mahrt, 1997), and flux damping (Moncrieff et al., 2005)? See also Burba (2022).
I’m concerned for the quality of the nighttime fluxes, as in this case the gradients contain very mixed stability classes. The authors should make different groups in their analysis (e.g. in Figure 10) such that their daytime fluxes and nighttime fluxes can be judged individually.
Specific comments
Line 141: Please also indicate that you will also deal with other sources of uncertainty in Section 3.5.
Figures 5 and 6: I suggest to keep the correlations on the 34m level measurements, and to replace the current column with correlations at the 26m level with correlations between the gradient (dNH3) and the respective meteorological variables. The difference in relations between the 34 and 26 m levels is only limited (as mentioned in the manuscript) and should be due to flux (and thus visible in the gradient), which makes this a good alternative. This could be a good place to place some more emphasis on the found gradients (which have a higher confidence than the flux). Moreover, perhaps relations between dNH3 and temperature are less dependent on wind direction (see also point 6) and could even reveal a relation with RH.
Lines 267-273: The fact that sign of the correlation changes when incorporating data from certain wind directions aligns with comment 6 of Referee #2; Being that in this case the wind direction, NH3 concentration, and air temperature are correlated. Does that mean that the found relation is only valid under (north)easterly winds? The authors should include a discussion on multicollinearity in the manuscript. Moreover, can the authors speculate on the causality between air temperature and NH3 concentrations?
Lines 271-273: In the authors response to comment 10 of Referee #1 they also remove data from February 2nd. Did they also do this in the manuscript?
Figure 7: Please indicate which measurements where obtained on February 1st and 2nd, such that it is easily visible to support your statements on lines 271-273.
Lines 336-338: Does this mean that the authors found only a significant gradient for in total two or four of the cases?
Lines 360-361: In the previous sentence the authors refer to a difference between leafy and leafless, but in this sentence the difference seems to be about daytime or nighttime. Please rewrite for more clarity.
Figure 10: The xticks in panels (d) and (h) show values up to 1.5 kW m-2, while in Figure 4(e) and (j) maximum solar radiation is around 0.7 kW m-2. I suspect erroneous xticks. Additionally, I suspect that solar radiation may not be the sole explanatory variable here, but rather the fact that during daytime (high solar radiation) there is a higher transfer coefficient (as in Figure 8). Can the authors reflect on the causality (or multicollinearity) between solar radiation, transfer velocity and NH3 fluxes?
Section 3.5: This section should start a general statement that the calculated fluxes come with major uncertainties. I propose to start this section with: “Due to several limitations of the measurement set-up, the calculated fluxes contain major uncertainties”, followed by a systematic listing of the limitations and the induced errors, being: “(i) Instrumental set-up did not allow us to validate the roughness sublayer effect on the AGM. Therefore, we tentatively used the alpha factor of (…). This introduces an error of (…). (ii) The long measurement of ΔNH3 interval introduces two uncertainties: (ii-a) ignorance of the cross term and (ii-b) averaging over multiple stability classes. To estimate the error due to ignoring the cross-term (ii-a), we (…). This introduces an error of (…). During daytime, the effect of (ii-b) is limited because (…). However, during nighttime (…). This introduces an error of (…). Finally, (iii) the use of AGM leads to larger errors compared to REA or EC (…). This introduces an error of (…)”
The authors should end this part with an estimate of the total error for the random error and for the systematic error (as a quadratic sum to all separate errors), such that the it is easier to evaluate the reported results. The authors could even do that separately for each measurement interval (i.e. D1, D2 and N), as I suspect that during nighttime errors are much larger than during daytime. After that, they can continue with the rest of Section 3.5.
References
Burba, G. (2022). Eddy covariance method for scientific, regulatory, and commercial applications. LI-COR Biosciences.
Finkelstein, P. L., & Sims, P. F. (2001). Sampling error in eddy correlation flux measurements. Journal of Geophysical Research: Atmospheres, 106(D4), 3503-3509.
Gash, J. H. C., & Culf, A. D. (1996). Applying a linear detrend to eddy correlation data in realtime. Boundary-Layer Meteorology, 79(3), 301-306.
Melman, E. A., Rutledge-Jonker, S., Braam, M., Frumau, K. F. A., Moene, A. F., Shapkalijevski, M., ... & van Zanten, M. C. (2024). Increasing complexity in Aerodynamic Gradient flux calculations inside the roughness sublayer applied on a two-year dataset. Agricultural and Forest Meteorology, 355, 110107.
Moncrieff, J., Clement, R., Finnigan, J., & Meyers, T. (2004). Averaging, detrending, and filtering of eddy covariance time series. In Handbook of micrometeorology: A guide for surface flux measurement and analysis (pp. 7-31). Dordrecht: Springer Netherlands.
Shapkalijevski, M., Moene, A. F., Ouwersloot, H. G., Patton, E. G., & Vilà-Guerau de Arellano, J. (2016). Influence of canopy seasonal changes on turbulence parameterization within the roughness sublayer over an orchard canopy. Journal of Applied Meteorology and Climatology, 55(6), 1391-1407.
Vickers, D., & Mahrt, L. (1997). Quality control and flux sampling problems for tower and aircraft data. Journal of atmospheric and oceanic technology, 14(3), 512-526.
Webb, E. K., Pearman, G. I., & Leuning, R. (1980). Correction of flux measurements for density effects due to heat and water vapour transfer. Quarterly Journal of the Royal Meteorological Society, 106(447), 85-100. |
The manuscript entitled “Ammonia exchange flux over a tropical dry deciduous forest in the dry season in Thailand”, written by Xu et al., presents a unique study to NH3 exchange between forest and atmosphere in the tropics. It is, to my knowledge, the first to present such data and is valuable to the general NH3 community. They systematically discuss the drivers of both NH3 concentration at different levels and the flux and suggest future studies that could answer questions raised in their own study. Their methodology, however, lacks clarity. As this could potentially affect a major part of their analysis, I recommend major revision. In addition, I have some specific comments and suggestions that could further improve their study.
Comments on methodology
Specific comments:
Technical comments