The USGS LOADEST model uses a fitted a.sin(x)+b.cos(x) term to model seasonality in river nutrient loads. This is just a quick ReactJs / Highcharts applet to illustrate how that term behaves.