Update SMAP L1C DA to use dynamic Tb_error information to inflate observation errstd - #190
Draft
gmao-qliu wants to merge 3 commits into
Draft
Update SMAP L1C DA to use dynamic Tb_error information to inflate observation errstd#190gmao-qliu wants to merge 3 commits into
gmao-qliu wants to merge 3 commits into
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
This PR updates read_obs_SMAP_halforbit_Tb to utilize the dynamic tb_error values provided in the SMAP L1C data to change the observation error standard deviation (errstd). The purpose is to reduce the impact of the RFI-affected Tbs in the SMAP L4 system.
Previously, errstd for SMAP L1C Tb was set at a constant 4K. By reading and applying the actual L1C Tb_error, the assimilation system now increases the error of Tb observations based on their provided Tb_error values using:
errstd = errstd * inflation_factor
where:
inflation_factor = exp(20 * (Tb_err - 1.3))
and
inflation_factor is limited to [1.0, 20.0]
Impact & Expected Behavior
This change only impacts the SMAP L1C Data Assimilation. No other sensors or general DA settings were touched.
Results will remain identical to the previous implementation when Tb_error <= 1.3K.
Assimilation results will diverge from the baseline in instances where Tb_error > 1.3K.
Testing Done
Verified that the code compiles and runs without structural errors.
Confirmed Tb_errors are predominantly < 1.3K in normal circumstances.
Confirmed expected numerical deviations in assimilation for cases where Tb_error > 1.3K (e.g., due to low-level RFI).