Dear José Ramon, 

Thank you for your revised version, 

We are pleased to recommend your manuscript within the PCI for Forest and Wood Science, after you have addressed the few editing comments and typos to correct for the final version of your manuscript:

Organisation of the supplementary information, Figures and tables:

The current numbering of figures is a bit confusing. I suggest that you rename the “Appendix section (usually devoted to mathematical developments or additional raw data at the end of the main text) as a “Supplementary information” section. In this section, please add a supplementary methods section (see below) where you can include the parts of the previous texts that you removed from the original methods (original end of line 94 to line 101 and 104 to 109 which are think are important to include). 

Here in the methods, I think at least one of your Fig. A1 or A2 is needed in the main text. This will prevent having to go to supplementary to initially view the original or cleaned data with geographic map details on the case study. It can be renamed Fig. 2 since Fig 1 is the graphical method abstract. 

Then the Figure that you did not put back in the main text (A1 or A2) becomes Fig. S1,

Fig. B1, B2 and B3 become Fig. S2, S3 and S4,

Fig. C1 and C2 become Fig. S5 and S6, and Tables C1 and C2 become Tables S1 and S2,

Fig. D1, D2, D3 and D4 become Fig. S7 to S10.

Then rename the figures in the text accordingly.

The quality of the figures is rather poor (we can’t zoom in to better see the details or text within figures for example), please refer to the formatting requirements of the PCI to make sure the quality is sufficient. 

Summary: 

About the third point (sentence before last), I suggest to add “in fire and forest loss patterns” after “significant differences” based on your results.

Introduction:

Line 2: add “(SDG)” since you are using the abbreviation later after “Sustainable Development Goal”

Lines 5 to 8: can you please cite a reference after “in preserved areas”, and another one if possible after “secondary forests”.

Materials and methods

2.1 part

In Fig. 1: correct “gobal” to “global” in the top left square, replace “stats” by “statistics” the 3 times it is used, it should fit. Change “false positive” to “false positives”. In the legend of Figure 1: please add “and see the Data and code availability section for provided scripts.“ after “details”.

Lines 70: Change “Fig. A1” to “Fig. S1” or Fig. 2 (depending on which alternative you choose for or of Fig. A1 or A2 to put in the main text, see above).

Line 104: Change “Fig. A2” to “Fig. S1” or Fig. 2 (idem).

There are 2 parts of the previous text that you removed in this 2.1 part (without indicating it in tracking mode). This makes the text easier to read but as said above, please put those parts the supplementary section.

Lines 111 and 181: Replace “statistics techniques” by “statistical techniques”.

Lines 127 to 132 should be relocated at the end of part 2.1, since after the consistency check (top of Figure 1), you focus now on the MODIS dataset in this revised version. Thus the VIIRS datasets should not appear in the rest of the sections as it is implicit that you are using similar approaches within the corresponding periods for consistency checks.

Lines 134 to 140: this paragraph is not edited as requested, according to the MODIS data focus only. Suggestions are to: add “noise-free” after “MODIS”, add “the above defined” before “forest”, add “cover” after “forest”, then remove 

“stands with a canopy closure equal to or greater than 25%, which I generated from the year 2000 tree cover raster layer. Then, I”

 since it is redundant with what’s explained above and in the Fig. 1, remove also “from both datasets”.

Also please change “Last, I divided the number of points by the cell area in square kilometers, and then again by the number of years of each of the two periods of analysis, from which I obtained two data fields, one for each period of analysis, containing the average density of fire points per square kilometer per year (hereafter, fire density).” into “Last, I divided the number of points in each cell by its area (in km2) and by the number of years, which resulted in fire density (Fig. 1).” 

Lines 141 to 146: these lines can be summarized without loss of information content as follow, if you agree: “I pooled forest loss surface area representing the period 2001 to 2018, then divided it by the corresponding cell size and by 18 years to obtain the average forest loss per unit area per year.”

2.2 part: some summarizing can be suggested:

Lines 154 to 156: please replace “In particular, I focused on assessing the association between those variables considering the size of the forest clearings, using both a..” just by “with…” in the text, since the details are given afterwards on how the different size of clearings is considered.

Lines 158 to 160 until “source”: replace “I generated 18 maps of annual forest loss, one per each year of the study period, using the loss year raster as a source” with “I used the forest loss year raster to generate 18 annual maps”. 

Lines 162 to 166: replace “Afterwards…1 ha in size” with “Additionnally, annual maps of “small forest clearings” were produced with patches of less than 1 ha in size, and maps of “medium- and large-sized forest clearings” (or large clearings”) with patches larger than 1 ha (Fig. 1)”.

Lines 176 from “I used a regular grid…”: please summarize as “This larger area than for cells used in the long-term approach was chosen to reduce the skewness of variables distributions and improve adherence to normality”.

Part 2.3

Line 197: remove “the”.

Line 198: remove “method”.

Line 204: remove “Afterward”.

Line 210: please add a reference for the Moran I statistic (suggestion: Sokal, R. R. and Oden, N. L. (1978). Spatial autocorrelation in biology. 1. Methodology. Biol J Linn Soc, 10: 199–228.).

Line 233: remove “The final stage was to produce”

Line 234: add “were produced for” after “graphics” and remove “in” before “QGIS”

Line 239: add “; see also the Data and code availability section” after 2017

Results

Line 244: replace “quite high” by “close to 1”

Line 249: replace “was good enough” by “highly correlated across time and thus appropriate”.

Line 258: add “points” after 11700.

Figure 3 legend: recall what HH and LL mean.

Line 297: replace the “Spatial modelling” subtitle with “Spatial dependence between Forest loss and Fire density”.

Lines 298 to 306: these lines are redundant with lines 199-200 and 221 to 223 in the methods, please combine lines 298 to 306 to those previous lines and remove from the results.

Lines : 309-310: please remove “The relevant statistics of the spatial error models fitted are summarized in Table 2” and add “; Table 2” after “(p<<<0.01”.

Also please refer to the formatting requirement for reporting probability for a recommended PCI manuscript format.

Line 312: add “the” before “spatial error model”.

Line 335: replace “originated” with “originating”

Lines 340 and 344: remove “see” before “Figs”

Line 357: please change “analysis” to “patterns” and remove “spatial modelling” from the subtitle.

Figure 8 and Figure 9 legends: please change “Same legend as in Figure 7” with “See Figure 7 for colour legend of hexagons”.

Line 426: add “was” before “significantly”

Discussion

Line 468: please replace “All in all” with “Overall”

Dear José Ramon, 

Three reviewers have now assessed your work. They are overall very positive and find that both the methods using remote sensing data and the patterns described are convincing, helping the understanding of links between fire due to shifting agriculture and forest loss in Dominican Republic ecosystems. They also think that your spatio-temporal approach could inspire other studies in other countries. 

They explain however that you need to address several issues before your manuscript could be recommended, and I am joining them with a few additional suggestions.

The first one concerns the methodological approach and comparison of datasets: in the long-term approach, you compare statistics and patterns from the MODIS and VIIRS datasets but use different time periods (see more details in the referees’ comments). This makes it difficult to interpret results, and the comparison should be restricted to the same periods (e.g. in Table 1). It is also proposed that after this comparison, you should focus on the longer dataset. I agree with these comments and am adding that it could also be useful, if you consider it not too redundant to the annual approach, to compare a few different multi-year periods (overlapping or not) to address your working hypothesis about an increase of association between fire and forest loss. Besides, possible benefits of one dataset compared to the other could be discussed and lead to recommending which one is more appropriate to answer the targeted questions. You may also put some of the annual approach results, using the second dataset, in supporting information.

Addressing this first issue will help address the second one that concerns mostly the form of the manuscript and the presentation of the results. I agree with the two referees who suggest restructuring part of the text (some methods being currently in the results, see details), and also reducing both methods and results, in order to clarify them (adding also sub chapters) since they are too long, and not easy to follow in some parts. I am adding one suggestion here: to please make a summarized figure of the flow of statistical analyses, to help follow the methods text, as this would greatly help to understand this part and help you summarize your text. Other comments from the referees explain specific parts that require more explanations or corrections. 

The third important point is to extend the discussion a bit more, referring to your initial questions and developing more how your approach and results could be used to monitor or predict risks and prevent forest loss threats, integrating possible alternative strategies (see comments from the third referee).

Finally, all the valuable code that you have made available on your Github website for reproducing the analyses is not put forward nor referred to in your current manuscript. I suggest that you add a "Data and code availability" paragraph at the end of the manuscript where you provide the links, explaining also how it is structured, you can also refer to it at the end of the methods.   

In your revised manuscript, please provide constructive solutions to all main and detailed comments from the referees below, as well as to the additional comments above,

Thank you for considering the PCI Forest and Wood Science for submitting your work, and I am looking forward to receiving your revised version,

Yours sincerely,

Pauline Garnier-Géré