www.waters.gu.se
An uncertainty framework for marine indicators based on monitoring data
Jacob Carstensen
Department of Bioscience Aarhus University HELCOM workshop, Uppsala 18 May 2015
- Variations in time, space and
methodology
An uncertainty framework for marine indicators based on monitoring - - PowerPoint PPT Presentation
HELCOM workshop, Uppsala 18 May 2015 An uncertainty framework for marine indicators based on monitoring data - Variations in time, space and methodology Jacob Carstensen Department of Bioscience Aarhus University www.waters.gu.se
www.waters.gu.se
Jacob Carstensen
Department of Bioscience Aarhus University HELCOM workshop, Uppsala 18 May 2015
methodology
www.waters.gu.se
Decision-making is inherently uncertain!
www.waters.gu.se
Why bother about uncertainty?
CIS guideline #7
WFD classification and uncertainty
www.waters.gu.se
How can we determine the confidence in status classification?
www.waters.gu.se
How can we determine the confidence in status classification?
www.waters.gu.se
The simple and ”convenient” approach
2 4 6 8 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
Spreadsheet solution Mean = 1.61 µg L-1 St.dev. = 1.56 µg L-1 Mean = 1.52 µg L-1 St.dev. = 1.35 µg L-1
www.waters.gu.se
The simple and ”convenient” approach
2 4 6 8 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
Spreadsheet solution Mean = 1.61 µg L-1 St.dev. = 1.56 µg L-1 Mean = 1.52 µg L-1 St.dev. = 1.35 µg L-1 Mean = 1.56 µg L-1 St.dev. = 1.46 µg L-1
www.waters.gu.se
The simple and ”convenient” approach
2 4 6 8 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
Spreadsheet solution Mean = 1.61 µg L-1 St.dev. = 1.56 µg L-1 Mean = 1.52 µg L-1 St.dev. = 1.35 µg L-1 Mean = 1.56 µg L-1 St.dev. = 1.46 µg L-1
www.waters.gu.se
www.waters.gu.se
The uncertainties to be considered
Fixed Random
www.waters.gu.se
70% 10% 20%
Including fixed factors to reduce random variation
www.waters.gu.se
70% 10% 20% 4.5% 0.5% 95%
Including fixed factors to reduce random variation
www.waters.gu.se
Reducing uncertainty by including seasonal variation
108% and 113%
0.1 1 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
www.waters.gu.se
Reducing uncertainty by including seasonal variation
108% and 113%
0.5 1 1.5 2 2.5 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
(season): 101% and 99%
0.1 1 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
www.waters.gu.se
Reducing uncertainty by including interannual variation
108% and 113%
0.5 1 1.5 2 2.5 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
(season) 101% and 99%
0.1 1 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
0.5 1 1.5 2 2007 2008 2009 2010 2011 2012 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
(season & year): 93% and 95%
www.waters.gu.se
Splitting the residual variation between random components
0.1 1 10 2007 2008 2009 2010 2011 2012 2013 Chlorophyll a (µg L-1) Station 3 m depth Station 7 m depth
(season & year): 93% and 95% Component VAR SE
Station 0% Station×month 0% Station×year 0% year×month 0.1472 0.3837 47% Residual 0.2730 0.5225 69%
www.waters.gu.se
Component VAR SE
Station 0.1212 0.3481 42% Station×month 0.0249 0.1578 17% Station×year 0.0321 0.1791 20% year 0.0160 0.1265 13% year×month 0.1255 0.3543 43% Residual 0.4113 0.6414 90% Total variation 0.7310 0.8550 135%
Chlorophyll a uncertainty components
To reduce uncertainty it is important to focus on the largest sources
www.waters.gu.se
www.waters.gu.se
Another example: Eelgrass shoot density
Area Site # stations # years # divers # observations North Sound Vitsandsbrygga 1 2 2 12 Höganäs 2 12 2 144 Central Sound Landskrona 2 17 4 288 Bjärred 4 16 3 288 Lomma 3 13 3 168 Limhamn 2 4 2 120 South Sound Bunkeflo 3 11 3 204 Klagshamn 5 17 4 450 Bredgrund 3 11 2 167 Ö. Haken 1 1 1 6 Lilla Hammar 1 1 1 6 South coast Segelskär 1 1 1 6 Fredshög 1 11 3 60
www.waters.gu.se
Area Site V[GRADIENT] V[YEAR] V[PERSON ] V[G×Y] V[PATCHINESS] North Vitsandsbrygga
Sound Höganäs 0.2119 0.0164
0.1761 Central Landskrona 0.0262 0.0082 0.0697 0.0082 0.2275 Sound Bjärred 0.0520 0.0292 0.0718 0.0111 0.1370 Lomma 0.0777 0.0067 0.0123
Limhamn
0.0167 0.0121 0.4577 South Bunkeflo 0.5013
0.0508 0.1132 Sound Klagshamn 0.1142 0.0368 0.0137
Bredgrund 0.0734 0.0388 0.1986 0.0117 0.0391 Ö. Haken
Lilla Hammar
South coast Segelskär
Fredshög
0.0917
Model 1: Analysis of individual sites
Eelgrass shoot density
www.waters.gu.se
Model 2: all sites combined
Component VAR SE
YEAR 0.0098 0.0990 10% STATION 0.1273 0.3570 43% YEAR×STATION 0.0133 0.1153 12% PERSON 0.0599 0.2447 28% PATCHINESS 0.2151 0.4638 59% Total variation 0.4254 0.6522 92%
500 1000 1500 2000 Shoot density (m-2) A)
200 400 600 800 1000 1200 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 Shoot density (m-2) B)
www.waters.gu.se
Component VAR SE
YEAR 0.0063 0.0794 8% STATION 0.0181 0.1345 14% YEAR×STATION 0.0222 0.1490 16% PERSON 0.0768 0.2771 32% PATCHINESS 0.0944 0.3072 36% Total variation 0.2178 0.4667 59%
Model 3: Explaining the gradient by depth
500 1000 1500 2000 Shoot density (m-2) A)
200 400 600 800 1000 1200 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 Shoot density (m-2) B)
www.waters.gu.se
Estimation of indicator variance – crossed design
Estimating indicator variance – crossed design
Number of years Number of stations Number of replicates
www.waters.gu.se
Uncertainties at the indicator level propagate through the integrated assessment
www.waters.gu.se
Uncertainties at the indicator level propagate through the integrated assessment
www.waters.gu.se
Uncertainties at the indicator level propagate through the integrated assessment
www.waters.gu.se
Conclusions
uncertainty that must be considered
precision requires a large dataset with a structure that allows for identification of these components
components can be reduced by including fixed, explanatory factors in the model
integrated assessment and this allows for estimating the confidence in classification and the risk of misclassification
www.waters.gu.se
So, how do I calculate uncertainties in my spreadsheet?
www.waters.gu.se
So, how do I calculate uncertainties in my spreadsheet?
www.waters.gu.se
Uncertainty
www.waters.gu.se
Uncertainty Uncertainty
www.waters.gu.se
Uncertainty Uncertainty