Open amsnyder opened 1 year ago
@galengorski
When we look at the correlation and mutual information between discharge and salt front location within each one of those river mile bins, it looks like this:
The top panel shows that in general there is a stronger correlation (more negative) when the salt front is lower in the estuary, and the correlation gets weaker as the salt front moves up. The bottom panel shows mutual information, which doesn't show much of a signal at any point in the estuary. Stronger correlation in the lower estuary make some sense as this is when discharge is likely more variable, in contrast when the salt front is further up the estuary, the discharge is low and doesn't change much.
@galengorski
This led to the question of how does past discharge affect the movement of the salt front in each river mile bin. To answer that, I did the following:
- For each day, take the mean and standard deviation of the previous n days of discharge
- This results in a running mean (and standard deviation) of the previous n days of discharge
- For each river mile bin, calculate the correlation and mutual information between the mean (and standard deviation) of the previous n days of discharge and the current salt front location
- Do this for n = 1, 7, 15, 30, 60, 120 days to determine the "time scale" of influence for each driver at each river mile bin
Here is what it looks like for correlation with mean discharge at Trenton:
I think it's interesting to look at these plots column by column. For example, the last river mile bin 82-91 show that the mean discharge from the last 60 days is a much stronger predictor than the last 15 days. Each river mile bin shows distinct time scales of influence, lower river mile bins are correlated with mean discharge of the past 7-15 days.
The plot below shows the same idea but for mutual information instead of correlation:
I think what is interesting here is that although the patterns are similar for correlation and mutual information, river mile bin 78-82 shows a very different relationship. This is an area of the river where it widens out and older water may become "trapped" during periods of low flow. This might indicate that there are non-linearities in the relationship between discharge and salt front here.
@salme146
After meeting with Galen to discuss these results I looked at the COAWST modeled depth-averaged velocities at Trenton vs. observed discharge at Trenton (USGS) and plotted them:
Given the river mile "reaches" I calculated the time scale in days for the water at Trenton to arrive at each reach and plotted the results here:
This makes sense, as the higher discharges correlate to higher flows, which means water from Trenton will arrive at the reach sooner, when flow is faster. In the summer months it takes longer for water to get to each of these locations, so the effect of discharge on the salinity in each location is lagged by a different amount.
@galengorski
That's super interesting, these figures are awesome! What do you think about converting these to a travel time distribution for each river mile bin so that way for example R6: 82-91 has distribution of times for discharge to reach that area. We could then look at when the "time scales of influence" from the heatmaps above are significantly different from the mean travel time distribution and that could clue us into areas where discharge is more and less dominant in influencing the salt front...
@salme146
If we take a look at the average time it takes water at Trenton to get to each "reach" :
Annually R1:31-58 : 6.5 days R2:58-68 : 5.7 days R3:68-70 : 5.5 days R4:70-78 : 4.8 days R5:78-82 : 4.5 days R6:82-91 : 3.7 days
August-October R1:31-58 : 13.5 days (~Spring-Neap tidal cycle) R2:58-68 : 11.7 days R3:68-70 : 11.4 days R4:70-78 : 10 days R5:78-82 : 9.2 days R6:82-91 : 7.6 days
@galengorski
Here are the same plots, but instead of looking at the mean of the previous n days of discharge at Trenton, I looked at the standard deviation of the previous n days of discharge at Trenton.
Correlation:
Mutual information:
The patterns don't look all that different from looking at the mean, which confirms that periods of low and relatively uniform flow are associated with salt front movement in the upper estuary.
@galengorski
This is a finer scaled version of the correlation heatmap, it looks like the heatmaps capture the overall patterns ok, but the finer scale shows more structure.
river mile: minimum correlations 31-58: 12 58-68: 12 68-78: 33 70-78: 87 78-82: 21 82-91: 54
@galengorski
And here is a finer scaled version for the mutual information heatmap
river mile: maximum mutual information 31-58: 66 58-68: 16 68-78: 10 70-78: 36 78-82: 82 82-91: 68
@salme146
Noes from 2/17 meeting:
Different peaks in the data (say brown line) might correspond to different dynamics and time scales of those dynamics (large discharge events, and drought periods). So it might be useful to dig a little deeper into the timescales of influence.
Sub-daily data - would have enough data to make calculations about binned flow. sub daily data would also lead towards using tidal data which is really important. daily tidal signals change on sub 6 hour time scales, spring-neap tidal time scales are 10-14 days.
@galengorski
The plots for the Schuylkill look very similar, and in the interest of not clogging this discussion with plots, I'll leave those out.
For next steps, @salme146 suggested looking at the "expected time scale of influence" might be given the average river velocity and the distance from Trenton. Comparing the expected and actual time scales of influence might give us an idea of when discharge is dominant and when other drivers are more important.
@salme146
@galengorski , I posted plots above looking at river velocity vs discharge for 2019 and then plotting travel times from Trenton to each of the reaches we picked for the analysis. I used the highest river mile to calculate the time scale. I am thinking if we are talking more about "influence" maybe I should use the middle of the reach. In any case this gives you an idea of time scales. They are shorter than I thought, but again, this is based on modeled velocities at Trenton.
@galengorski