PNNL-PREMIS / PREMIS-ghg

Soil GHG experiment
MIT License
4 stars 1 forks source link

Tree proximity analysis #14

Closed bpbond closed 5 years ago

bpbond commented 6 years ago

We expect that soil respiration will be correlated with temperature–this is almost always true, except in extremely water-limited system–and probably with moisture.

One question that might make a good poster (and later paper) would be: is soil respiration rate, or its temp/moisture sensitivity, spatially correlated with nearby trees? See for example

http://dx.doi.org/10.1007/s10533-004-5889-6 http://dx.doi.org/10.1007/s11104-006-9052-0 http://dx.doi.org/10.1016/j.soilbio.2012.04.024

Since the surface CO2 flux is derived overwhelmingly from microbes and plant roots, we might hypothesize that collars closer to more and larger trees will have different flux dynamics than those far away. This would be straightforward to test once we have a dataset documenting, for each collar, distance to nearest trees (e.g. 2 m to tree #X, 3.5 m to tree #Y, etc). And we could use only the control collars for this question.

bpbond commented 6 years ago

H1: Collars closer to more trees (tree basal area) near them will have higher fluxes. H2: Collars closer to more trees will have higher temperature sensitivity. (Mechanism for H1.) H3: This effect will be limited to leaf-on (growing) season. In dormant season, no effect. H4: This effect will be stronger in well-watered conditions (downhill, non-drought), because trees dry before soil.

bpbond commented 6 years ago

@stephpenn1 Just committed collar_proximity.csv, which is real proximity data for two GCREW collars (I measured while Licor'ing today).

I think this will be a good data processing exercise! Let's say the first goal is to produce a graph with two lines on it, one for each collar; on the x axis is distance (m) and on the y axis is cumulative number of trees. In other words, it should show the number of trees within 1 m, 2 m, 3 m, etc., with the lines constantly rising.

Can you think about the steps required for this and open a PR, e.g.:

# 1. Read proximity data
# 2. QA/error check proximity data [what are errors?]
# 3. Read tree inventory data and join by...
# (Etc.)
stephpenn1 commented 6 years ago

@bpbond, awesome about new dataset! All looks/sounds great. How creepy is it that yesterday I created a new branch also named tree-proximity-analysis!

Edit: I guess we did use those terms when we discussed on Monday.. still surprised nonetheless

bpbond commented 6 years ago

Given that we tend to text each other simultaneously, not that big a shock, honestly. 😄

stephpenn1 commented 6 years ago

Still a bit confused on what data goes on the y axis for this plot: on the x axis is distance (m) and on the y axis is cumulative number of trees ... is this to get a distribution of the trees? Just a sequence of numbers 0-10 on the y axis?

bpbond commented 6 years ago

Distance is on x, so is 1:10 basically. The y axis is the number of trees within that distance from the collar. So the curves are cumulative: if there is one trees at 1m, and the second at 2m, etc., it'll be a straight line at a 45 degree angle. Does that make sense?

bpbond commented 6 years ago

Something like this: example

stephpenn1 commented 6 years ago

tree_dist

stephpenn1 commented 6 years ago

Not a fan of the gradient legend color.. working on fixing it

bpbond commented 6 years ago

Very nice!

bpbond commented 5 years ago

Seems safe to close this.