Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Rain intensity needs to included as a parameter in the model #79

Closed
PaulMelloy opened this issue Aug 5, 2021 · 8 comments
Closed

Rain intensity needs to included as a parameter in the model #79

PaulMelloy opened this issue Aug 5, 2021 · 8 comments
Labels
documentation Improvements or additions to documentation enhancement New feature or request
Milestone

Comments

@PaulMelloy
Copy link
Collaborator

Discussed in meeting,
The amount of rainfall in a rain hour/event effects the potentially_effective_spores. However, this is not included in the model.

We could easily add this as a multiplier, we need coefficients from models to base this multiplier off.

@IhsanKhaliq could you provide references that report coefficients for rain quantity impacting spore dispersal

@PaulMelloy PaulMelloy added documentation Improvements or additions to documentation enhancement New feature or request labels Aug 5, 2021
@PaulMelloy PaulMelloy added this to the RModel_Paper milestone Aug 5, 2021
@IhsanKhaliq
Copy link
Owner

There is a rain cauchy mulitplier for rain splash dispersal reported by both Art and Coventry. I tried my best but couldn't find anything else. Would these work?

Screen Shot 2021-08-06 at 8 32 19 am

Screen Shot 2021-08-06 at 8 35 52 am

@IhsanKhaliq
Copy link
Owner

But these are related to the distance spread. We want something related to the amount of spore release?

@PaulMelloy
Copy link
Collaborator Author

Yes, something related to rainfall in millimetres influencing the number of "effective spores" spread

@IhsanKhaliq
Copy link
Owner

IhsanKhaliq commented Aug 5, 2021 via email

@PaulMelloy
Copy link
Collaborator Author

PaulMelloy commented Aug 5, 2021

On another note the model at the moment the model uses the number of new growing points produced in the last day (24h) to influence the max_interception_probability.

ascotraceR/R/one_day.R

Lines 55 to 57 in f19aec5

max_interception_probability <-
interception_probability(target_density = 5 * max(daily_vals[["paddock"]][,new_gp]),
k = spore_interception_parameter)

In the table above there is a reference that new growing points might influence max_interception_probability for the 4-5 days prior?

@PaulMelloy
Copy link
Collaborator Author

Actually now I think about this...
That would explain the 5 *
multiplier, so disregard that question

@IhsanKhaliq
Copy link
Owner

Yes, that calculation should take of that. Indirectly, max_interception_probability or deposition probability would depend on the number of susceptible growing points that are not infected yet

Screen Shot 2021-08-06 at 10 08 34 am

&

Screen Shot 2021-08-06 at 10 07 39 am

@PaulMelloy
Copy link
Collaborator Author

I am going to close this issue as it has been addressed with the rainfall_multiplier branch of this repo. More experimentation is needed to discover actual coefficients that could be used as an estimator for the volume of spores dispersed in a spread event.
As a note from eyeballing the model in the rainfall_multiplier branch against field data it seems a better fit than without a rainfall multiplier. I suggest running permutations of the model with different rainfall multiplier values against field data to find the best coefficient here

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
documentation Improvements or additions to documentation enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants