I'm new to calculating GPP/CR from diel oxygen curves, and have been getting negative GPPs and I am trying to figure out what this means. I know one large source of error is the reaeration coefficient. I understand there are lots of opinions about estimating reaeration coefficients- we cannot afford to do injection tests. So, I have been getting the coefficients by looking at the slope between DO deficit and DO change per hour at night based off our understanding of Heffernan and Cohen 2010/Owens 1974. The resulting coefficients are reasonable- 5 to 20 day-1. But, continuing the one-station calculations results in negative GPPs. I have noticed that, despite having a clear diel curve, our streams are almost consistently undersaturated in DO- would that be the reason I am getting negative GPPs? I am looking for advice on interpretation or different methods to try!