Last updated: 2022-01-31

Checks: 7 0

Knit directory: mapme.protectedareas/

This reproducible R Markdown analysis was created with workflowr (version 1.6.2). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.


Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

The command set.seed(20210305) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

Nice! There were no cached chunks for this analysis, so you can be confident that you successfully produced the results during this run.

Great job! Using relative paths to the files within your workflowr project makes it easier to run your code on other machines.

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version 8f3f53a. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:


Ignored files:
    Ignored:    .RData
    Ignored:    .Rhistory
    Ignored:    .Rproj.user/

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.


These are the previous versions of the repository in which changes were made to the R Markdown (analysis/rapid-assessment_centralamerica.rmd) and HTML (docs/rapid-assessment_centralamerica.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File Version Author Date Message
Rmd 8f3f53a Ohm-Np 2022-01-31 update id-68135 with latest polygon data
html 8f3f53a Ohm-Np 2022-01-31 update id-68135 with latest polygon data
Rmd b48b821 Ohm-Np 2022-01-31 replace a WDPAID pol from Panama
html b48b821 Ohm-Np 2022-01-31 replace a WDPAID pol from Panama
Rmd 0eca73c Johannes Schielein 2022-01-26 update central america analysis
html 0eca73c Johannes Schielein 2022-01-26 update central america analysis

Map and input data

The following report is a quick assessment of forest cover loss in 19 protected areas (PA) in Caentral America. It is based on publicly available data from the World Database of Protected Areas (WDPA/IUCN) which was downloaded from the Protected Planet Website. Our initial goal was to analyze 21 areas that had been provided to us by our operational departments in a PDF document in Spanisch language. The identification of areas was based on a name search by which we were able to clearly identify 19 areas of the 21. Unfortunately two areas, EL Imposible-San Benito and El Imposible-El Balsamoare not in the current database from IUCN which only has an area named El Imposible (with WDPAID 12494 and reported area 17.65 sqkm) probably the ancestor of the afforementioned two areas.

To quantify forestcover loss we utilized data from the Global Forest Watch (Hansen et al, 2013)1. forest cover loss is defined in their work “as a stand-replacement disturbance, or a change from a forest to non-forest state.”. forest cover loss can either be the result of human activities or natural factors such as forest fires or hurricanes, which are especially relevant in central America. In order to identify probable causes it is usefull to look at truecolor satellite images (see map below) and cross-check with additional data-sources such as the NOAA website that provides historical data for hurricane tracks.

It is alo important to note, that the utilized approach does not quantify total forest area in the analyzed PAs because it (currently) does not account for forest regrowth. Rather then this it compares forest areas and estimates forest cover loss in comparision to the baseline in the year 2000. In this sense the data can be utilized to e.g. assess threat levels and disturbance dynamics but it is not intended to make a complete quantification of forest stands, especially regarding regrown secondary forests.

We processed the available 19 PA polygons using mapme.forest package for the year 2000 to 2020. The map below shows the raw input data. It is an interactive map where you can toggle on and off the provided layers and zoom into the regions to inspect the analyzed areas.

Quantification of forest cover loss

The following graphs show forest cover loss in each of the 19 analyzed areas seperatly. Note that it is difficult to assess the development of all areas at once because of the high intensity of forest cover loss in one area (Reserva Biológica Indio Maíz), which makes the development of the others difficult to visualize. Therefore you might click on the legend an deactivate/activate individual areas to improve the result. A double-click allows you to isolate an individual area.

Note: For Reserva Biológica Indio Maíz the actual cause for the spike in forest cover loss was a Hurricane in 2016.

Results Download

Finally, here is a table containing desaggregated results for the forest cover area and loss area for the year 2000 - 2020. You can also download this dataset.

Note: area in hectare

Download CSV


sessionInfo()
R version 3.6.3 (2020-02-29)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.6 LTS

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.7.1
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.7.1

locale:
 [1] LC_CTYPE=C.UTF-8       LC_NUMERIC=C           LC_TIME=C.UTF-8       
 [4] LC_COLLATE=C.UTF-8     LC_MONETARY=C.UTF-8    LC_MESSAGES=C.UTF-8   
 [7] LC_PAPER=C.UTF-8       LC_NAME=C              LC_ADDRESS=C          
[10] LC_TELEPHONE=C         LC_MEASUREMENT=C.UTF-8 LC_IDENTIFICATION=C   

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] magrittr_2.0.1        rmarkdown_2.11        plotly_4.9.3         
 [4] RColorBrewer_1.1-2    htmltools_0.5.1.1     scales_1.1.1         
 [7] ggsci_2.9             leaflet.extras2_1.1.0 leaflet.extras_1.0.0 
[10] leaflet_2.0.4.1       sf_1.0-5              forcats_0.5.1        
[13] stringr_1.4.0         dplyr_1.0.7           purrr_0.3.4          
[16] readr_1.4.0           tidyr_1.1.4           tibble_3.1.6         
[19] ggplot2_3.3.4         tidyverse_1.3.1      

loaded via a namespace (and not attached):
 [1] httr_1.4.2              sass_0.4.0              viridisLite_0.4.0      
 [4] jsonlite_1.7.2          modelr_0.1.8            bslib_0.2.5.1          
 [7] assertthat_0.2.1        askpass_1.1             cellranger_1.1.0       
[10] yaml_2.2.1              pillar_1.6.4            backports_1.2.1        
[13] glue_1.6.0              digest_0.6.27           promises_1.2.0.1       
[16] rvest_1.0.0             leaflet.providers_1.9.0 colorspace_2.0-1       
[19] httpuv_1.6.1            pkgconfig_2.0.3         broom_0.7.6            
[22] haven_2.3.1             whisker_0.4             later_1.2.0            
[25] openssl_1.4.5           git2r_0.28.0            proxy_0.4-26           
[28] farver_2.1.0            generics_0.1.1          ellipsis_0.3.2         
[31] withr_2.4.2             lazyeval_0.2.2          cli_3.1.0              
[34] crayon_1.4.2            readxl_1.3.1            evaluate_0.14          
[37] fs_1.5.0                fansi_1.0.0             xml2_1.3.2             
[40] class_7.3-19            tools_3.6.3             data.table_1.13.6      
[43] hms_1.1.1               lifecycle_1.0.1         munsell_0.5.0          
[46] reprex_2.0.0            compiler_3.6.3          jquerylib_0.1.4        
[49] e1071_1.7-9             rlang_0.4.12            classInt_0.4-3         
[52] units_0.7-2             grid_3.6.3              rstudioapi_0.13        
[55] htmlwidgets_1.5.3       crosstalk_1.1.1         gtable_0.3.0           
[58] DBI_1.1.2               R6_2.5.1                lubridate_1.7.10       
[61] knitr_1.34              utf8_1.2.2              workflowr_1.6.2        
[64] rprojroot_2.0.2         KernSmooth_2.23-20      stringi_1.7.6          
[67] Rcpp_1.0.7              vctrs_0.3.8             dbplyr_2.1.1           
[70] tidyselect_1.1.1        xfun_0.24              

  1. “Hansen, M. C., P. V. Potapov, R. Moore, M. Hancher, S. A. Turubanova, A. Tyukavina, D. Thau, S. V. Stehman, S. J. Goetz, T. R. Loveland, A. Kommareddy, A. Egorov, L. Chini, C. O. Justice, and J. R. G. Townshend. 2013. “High-Resolution Global Maps of 21st-Century Forest Cover Change.” Science 342 (15 November): 850–53. Data available on-line from: http://earthenginepartners.appspot.com/science-2013-global-forest."