We obtained geographic and demographic data for the 2010 US Census
Bureau’s decennial Census using the package ‘tidycensus’ (Walker
2020). We then estimated the population of each watershed by
intersecting the polygons describing watershed boundaries (as
delineated at 1m resolution in Cary Institute of Ecosystem Studies,
Lagrosa, and Welty 2017) with the 2010 Census blocks using the package
‘sf’ (Pebesma 2018). We then determined the proportion of each census
block inside each watershed boundary and multiplied this by the total
population in the respective census block to estimate the population
of each block within the watershed. Finally, we summed the
proportion-adjusted populations of all Census blocks that intersected
a watershed to estimate the total population within the watershed. The
result is that the estimated population is the proportion of the block
in the watershed. All calculations were conducted in R version 3.6.2
(R Core Team 2019).
Cary Institute of Ecosystem Studies, J. Lagrosa, and C. Welty. 2017.
GIS Shapefile, Spatial boundaries and land cover summaries for eight
sub-watersheds of the Baltimore Ecosystem Study LTER ver 100.
Environmental Data Initiative.
https://doi.org/10.6073/pasta/ad0cce16ef6165913ea26b97e295f985.
Accessed 2020-04-12.
Pebesma, E., 2018. Simple Features for R: Standardized Support for
Spatial Vector Data. The R Journal 10 (1), 439-446,
https://doi.org/10.32614/RJ-2018-009.
R Core Team. 2019. R: A language and environment for statistical
computing. R Foundation for Statistical Computing, Vienna, Austria.
URL https://www.R-project.org/.
Walker, Kyle. 2020. tidycensus: Load US Census Boundary and Attribute
Data as 'tidyverse' and 'sf'-Ready Data Frames. R package version
0.9.6. https://CRAN.R-project.org/package=tidycensus.