research_blog/2022/12/13/TRAPS-update

[utterances-bot]

TRAPS Update | Weekly Research Update

This week after unsuccessfully attempting to shift Birch Bay WWTP seawards, I decided to remove the point source entirely. After doing so, LiveOcean was able to run for two days. I also calculated some simple summary statistics (means and standard deviations) of the climatology values for rivers and point sources. These mean values were then used to replace the climatology of rivers with weird biogeochemistry data. More details below. Shifting Birch Bay WWTP Seawards I decided to shift Birch Bay WWTP seawards by two grid cells. Previously, shifting this souce seawards by one grid cell was not enough to eliminate all instabilities. I was hoping that shifting Birch Bay WWTP seawards by two grid cells would be sufficient to prev

https://ajleeson.github.io/research_blog/2022/12/13/TRAPS-update.html

[parkermac]

I'm definitely not going kayaking anywhere near Birch Bay.

[parkermac]

Aurora, this is an excellent write up. I am still perplexed by Birch Bay, but I think for now there is no harm in just leaving it out.

I really like the summary statistics and climatologies. Does your list of rivers include the ones already in LiveOcean? I would be interested to know the annual flow rate of the sum of all the rivers that are new. How does the sum of new "tiny" Puget Sound rivers compare to the annual average of about 1000 m3/s?

Also I would be interested to know the sum of all the N-loading from the WWTP's (annual average) to compare with the loadings currently in LiveOcean.

With all of these queries I am trying to get a sense of how big things are. How much loading and river flow were we missing by not having the TRAPS? This is the kind of information you want when designing the model forcing: is it big enough to make a difference that warrants the difficulty of getting the information?

Final question: how soon can I start making my own TRAPS forcing to include in some year-long parameter tests I am doing? And over what years would the forcing have inter-annual variation, as opposed to being filled with the climatologies you created?

[ajleeson]

Hi Parker, thanks for the feedback. I like the idea of comparing the magnitude of the new inputs to the current inputs.

The list of rivers no longer includes the ones already in LiveOcean, so the summary statistics you see here are only the new rivers. I'll note that the rivers don't all discharge into Puget Sound (e.g. it includes the Willamette down in Oregon).

The annual average flow of the sum of all tiny rivers is 5,822 m3/s. The largest river is the Willamette. If I look at only the rivers discharging into the Salish Sea (Puget Sound + Straits), then the annual average flow for the sum of those rivers is 4,892 m3/s.

The annual N-loading of the sum of all WWTPs is 64,444 kg/day. I'll need to do a bit more work to get the comparison to LiveOcean's current loadings. Something to look forward to.

As for adding TRAPS forcing to your own runs, I had shared an older version of the TRAPS code with Jilian which works for her Hood Canal hydrodynamics model. However, I haven't updated it to include any recent updates (such as fixing the nans in the biogeochemistry). I should be able to clean up the new code and push it to Github by the end of the week. I'll send you an email when it's done. Right now it is set up to only run with the climatologies and not with the individual years.

[parkermac]

Aurora,

Thanks for the notes.

This one seems surprising to me: "The annual average flow of the sum of all tiny rivers is 5,822 m3/s. The largest river is the Willamette. If I look at only the rivers discharging into the Salish Sea (Puget Sound + Straits), then the annual average flow for the sum of those rivers is 4,892 m3/s."

The reason this is surprising to me is that the annual average of all the rivers coming into the Salish Sea in LiveOcean (the largest being the Fraser) is about 6000 m3/s. Is your 4,892 m3/s really excluding the LiveOcean rivers? I'm just double checking because that would mean we ignored almost 50% of the freshwater forcing.

Parker

On Wed, Dec 21, 2022 at 11:02 AM ajleeson @.***> wrote:

Hi Parker, thanks for the feedback. I like the idea of comparing the magnitude of the new inputs to the current inputs.

The list of rivers no longer includes the ones already in LiveOcean, so the summary statistics you see here are only the new rivers. I'll note that the rivers don't all discharge into Puget Sound (e.g. it includes the Willamette down in Oregon).

The annual average flow of the sum of all tiny rivers is 5,822 m3/s. The largest river is the Willamette. If I look at only the rivers discharging into the Salish Sea (Puget Sound + Straits), then the annual average flow for the sum of those rivers is 4,892 m3/s.

The annual N-loading of the sum of all WWTPs is 64,444 kg/day. I'll need to do a bit more work to get the comparison to LiveOcean's current loadings. Something to look forward to.

As for adding TRAPS forcing to your own runs, I had shared an older version of the TRAPS code with Jilian which works for her Hood Canal hydrodynamics model. However, I haven't updated it to include any recent updates (such as fixing the nans in the biogeochemistry). I should be able to clean up the new code and push it to Github by the end of the week. I'll send you an email when it's done. Right now it is set up to only run with the climatologies and not with the individual years.

— Reply to this email directly, view it on GitHub https://github.com/ajleeson/research_blog/issues/20#issuecomment-1361916655, or unsubscribe https://github.com/notifications/unsubscribe-auth/AC5WRX7GOAYXVO25A3IHXXDWONH5BANCNFSM6AAAAAAS6WIOG4 . You are receiving this because you commented.Message ID: @.***>

--

Parker MacCready

Leo Maddox Endowed Professor in Oceanography

Research Professor, School of Oceanography I live and work on the land and waters of the Squaxin Island Tribe. Email: @.** URL: faculty.washington.edu/pmacc LiveOcean Daily Forecasts http://faculty.washington.edu/pmacc/LO/LiveOcean.html Cell: (360) 359-1936, Office: (206) 685-9588 pronouns: he/him*

[ajleeson]

Hi Parker,

Yes, I did double check now. This is excluding the LiveOcean rivers. I will note that the largest new rivers coming into the Salish Sea are mostly Canadian rivers. Below I've included the flowrates for the largest 12 individual rivers that I'm adding.

Without any new Canadian rivers, the average annual flowrate being added into the Sound is 189 m3/s.

[parkermac]

Aurora,

Thanks for checking on this. I wonder how many of these enter the Salish Sea and how many go right into the Pacific.

I'm leaving on vacation tomorrow morning for a week. Have a great holiday!

Parker

On Wed, Dec 21, 2022 at 2:05 PM ajleeson @.***> wrote:

Hi Parker,

Yes, I did double check now. This is excluding the LiveOcean rivers. I will note that the largest new rivers coming into the Salish Sea are mostly Canadian rivers. Below I've included the flowrates for the largest 12 individual rivers that I'm adding.

[image: image] https://user-images.githubusercontent.com/15829099/209009713-61a8aa03-9a5a-408b-86fe-be733aacae7a.png

Without any new Canadian rivers, the average annual flowrate being added into the Sound is 189 m3/s.

— Reply to this email directly, view it on GitHub https://github.com/ajleeson/research_blog/issues/20#issuecomment-1362164847, or unsubscribe https://github.com/notifications/unsubscribe-auth/AC5WRX2IR7MXL5JHGWVMPMLWON5L5ANCNFSM6AAAAAAS6WIOG4 . You are receiving this because you commented.Message ID: @.***>

--

Parker MacCready

Leo Maddox Endowed Professor in Oceanography

Research Professor, School of Oceanography I live and work on the land and waters of the Squaxin Island Tribe. Email: @.** URL: faculty.washington.edu/pmacc LiveOcean Daily Forecasts http://faculty.washington.edu/pmacc/LO/LiveOcean.html Cell: (360) 359-1936, Office: (206) 685-9588 pronouns: he/him*

[ajleeson]

Happy holidays!