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

Adjust weights for PESTPP based on observation type #23

Open
jzwart opened this issue Sep 10, 2020 · 2 comments
Open

Adjust weights for PESTPP based on observation type #23

jzwart opened this issue Sep 10, 2020 · 2 comments

Comments

@jzwart
Copy link
Owner

jzwart commented Sep 10, 2020

Leaving the weights of the different observation types (e.g. temperature obs, flow obs) causes bias towards high value observations (or residuals) in objective function. According to PESTPP4.2.4 manual

A simple procedure that often works well is to assign observations of different types to different observation groups, and to then use a utility program such as PWTADJ1 from the PEST suite to ensure that the contribution made to the objective function by each observation group is about the same as that made by any other observation group at the commencement of the inversion process.
When calibrating a surface water model, or when including contaminant concentrations as a component of the calibration dataset of a groundwater model, it may be useful for the magnitudes of weights to reflect the magnitudes of measurements.

@jzwart
Copy link
Owner Author

jzwart commented Sep 11, 2020

Leaving flow and temperature weights all at 1 reduced RMSE for flow (6.1 RMSE compared to 6.8 for uncal), but inflated RMSE for temperature (~40 RMSE).

Currently trying out adjusting weight by log normalizing between 0 and 1. This gives more weight to low values (temperature values and low flow). See log normalizing (red) vs. linear (black) below
image
zoomed in x-axis
image

@jzwart
Copy link
Owner Author

jzwart commented Sep 18, 2020

Log normalization still biases flow. Setting temp weight to 1 and flow weights to 0.05 still biases calibration towards flow. Probably best to calibrate flow and temp separately for now.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant