From 48f97dfdc2407ddff7319a3f935a6cd0085213de Mon Sep 17 00:00:00 2001 From: shannarucker <42379382+shannarucker@users.noreply.github.com> Date: Thu, 6 Feb 2025 23:18:57 +0000 Subject: [PATCH] =?UTF-8?q?Deploying=20to=20gh-pages=20from=20@=20LimnoTec?= =?UTF-8?q?h/SCI@7d0431578845b94746e2f6c7aa41f942ad38d2ca=20=F0=9F=9A=80?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- articles/figures/rsa_join.png | Bin 23608 -> 21701 bytes articles/pre-processing.html | 20 ++++++---- pkgdown.yml | 2 +- reference/assess_dumpsites.html | 50 ++++++++++-------------- reference/assess_trash.html | 59 ++++------------------------- reference/assess_wq_nutrients.html | 12 ++++-- search.json | 2 +- 7 files changed, 49 insertions(+), 96 deletions(-) diff --git a/articles/figures/rsa_join.png b/articles/figures/rsa_join.png index a25d3c499307ae5cbe376dd467fe0d3d8522bf86..b72da14b4eb51aa2e95f3fdbad64f71843ee6a55 100644 GIT binary patch literal 21701 zcmcG$by!qi+crFOH&POk(x52KAR&lIgCLU93X;;@Al>QEASj`9hk%3(h;&OMIdp$( z@OR(u{T#>hJb%33_j!zN8DQ_V*Iw(o&g(qSYr<5OWkF5l@J1B%#?eAIG?R+6nfEVjG15X1Sd&-)=V>4 zS!XCY<-~1ERGfk~z=QS1YeUz<8|Exp(wyxX*%xp54+ZBLaZ^v4ALBk?DQU+Rc=GC5 zPMDlj>iKzhs)*xCKT)tTmD*fcQ$c1O+-uv>e6wnQ!Iy3EvcLa{uh2+~&5Re>BF4MxVg4We(1vu-x63WR<@AS|*^HBQIL= zit2dIr=|Ku_6)tlK4d%Bj!Q+Wv|69}F99aNYWWrpZD6mDYtJmksjpu`>CP z9|be$khSs&)kbfLGv!S6$&2>xp3!6bKFLhfq=JQ?WzB`+%~uWQe(>f|=ev=Q75i%| zMPt9x%gQ=wJb9|z^%V|%UtG-D0{LnDHM_VN7Y;q=j*E*=bcXFRLWQH;Y+93AJZKx5 zm(Q0uyqHP4llqUWn8D!3&rcSktCSvv9)%@+b*(Ujs9{BCbiVgG<9Y)Pto5a9xc1j( zI2DXCn^_s=NQwWvV=_C$cp8)v#Z#iLc0gcjxi>)wnS0vB>f)=$V9V#Oxj2pSr)=d_ ztjEf{>7VnI;{|S~A1x7XSOcflL=T%oPUeN*ir5WI+09+{x<{z4@@>)L89az@op7V6@v+tDUp2B@)3?w8o(_oUMU%a3m zEi;UM@b!h&vFC)grshfi{fHItG&oo(;b+lL6!}Ba()Xt3t|;&*dD>$;LeMaYPWnAF zPHD?Nq=qp_I+Un{7!fj*_hio}NB5bf<_lEVALGzcAkId>HB-g}$)3ZH3=53HgG6hM3i$-JMy4K%kdY`_OvD?fsUPu*chcr$b$m6KBx36FxXE zuarl+GnJ$>TZk=ouQ#EQnp*2UXQ7RTMK*0kbBkjw&$=8fluZm43ww`{H+ixAcE_U! z%fPx}NdOZpNB;y7|3slzJQjPe#7~#vg?~*)+=Y4GH`?{B-ppbBhM_rDMDbztU@?!F z7{h;9ya*R?RXR6ljty-}6?hmxOhVH0%PgwN z;~102Hlj;lcyH0H&$HAW#|h(^u&A|QxIW$huBWzVh!duX48|HyeZYn$q2=k zGy3`!mE*U<9@IMuGUF#@KI;BGJlDS>BVikZr|Xa1#J0PLG=Alxq)4EBtPPAFa8IXi zTBUdw?d!$wfl05~g#YGZ+i{$%h)7P-XZ?2ANoNgT&(w@xA9Bo*uuhy%tH?oC#IR6G zS+^o5a{X9wY0xawi)R07T?M+Fp*jVE$SYRxpc*d_eIyvsmbjI*4+fbuiC8^ddr@|( zGM|sypugL}7squ8V~R1FLQCMD9X4pbofQqU>mnM(%C#QWsr(jEcD?O=?J&eX)$R;_ zHe6A-w1@AmzUy`TQ!~zFdS8k2 zCcI}swZ0o@4%AHJM|6>ed$To*UZtUJScSV6Sh)v7*cKN?Ud2QQxfa!oYl85!q}352 zO))Ufb+F(`ir5L4#qO_s#wNKgv7b-oYma*=kwG#v?-#1e$!RB_&pM^JxB~@W7oYd3 zA}oe{_h`dAxe(PmIY-*Vjonv%weTcCzcRVKxWfu5+!_lk7!fpEp$$gy3lbE$j z4}P629N;)#c_;rYp7kiQ=}kf5tr5zkZt<|ExJIG_=$Lr9+-u z*1<3~U7=|3F^F{%i_=#uEZ;C_&vr@ay_Vi$f{CeL;P=!sK{RPbopU3$tp2oh-J9Nx z3=Lg#Z~xNAlBJ)V)6^+WLNXM&lokYjoT*-zcBAJxj7scC7y?JagMIEmil=l6|`CHdCf&CJ1pBbQ6kQw;mg3g#@m!xbSlxNE=M zd#iD;>$dCltlc#k1qG4P8JIUXUZ7SEWh z9xSxn)d1fG6FWO%kK^sf>guEj1cINrHOdaw3L@Fw)%i|tc{wp#yuwrvSbJ}Yy-$tg zLOtzRT%o%!_mS7#oN1t(P?3!Ehh4B(?m~P&Y83uXoz^42Sh_pj;KIkNSDUpZEj#8* zO!$C&m9eZAg>ZzAE#$W)yPQFy($pq%2DJToP-7$Kc)K_+WV{&B)ks9rnG{F#1ture zrsC-@uc1L^;I@2g`S90AjlF(>I(ByK3?h)XZre8c<>j%6pRMIWw|@Ni;hLR??z6su zu)^{MN%@Y11l{@Z&e&U~oHuP&GGY<@lT*fyuZ5JMD~DM2yYucTr-}oz@z2XWA&|_> zo!qbcUN(+{p7Uf-Oe`$qZ8P)nsmw$~L(i~nTMj8!EX;pT?XR3#?|=6`v&s;0#ANU~ zh;=3->D1}=)4kJ#-uHF0NMbO1=9%mRwl;*6@RzTb`+h1X!rm9cD?id#HC~JOT!|lz zj7r<`-w_hho|y?zF6=Ogmgw4*a*D(;4){D0>hzV~N)gxCx(b8a0!)Nag${VdXixj1 z*o5$(2X+jpeu^-u(7AXjR9!4iE+h=m>M81FyRPYcC1qf6kB0VY$1n7iOb?dSsB0i8 zJHtsJQQ>eHdXs_MmP6g{)Z7JTM@PrW0+COs|5&~C4p1pl3tm|hK$9xZyQ(X8tX%SdCrw1>r{kk2|9INrUyoD+tb`G zQ5)A@pnO!R%f;zO%Ju0vSMQh2qR-O|M`yQ~!g;~<1JP}vR3%vSyFaj?3noTJkLvU= zx6h(x-gYN)X_x;oyuw)KYo^*=h_&ROKUg2`eZv&qMEH^`V%?mI1t-qbS7avixz^ck ze*}3FcTihXm-PLz?yhG$lpexf?Tmp0Sg-xs>HZ`f(#Dy#{*0R!xKH2JPRcFt{XCC{ zBX9^xyw{IM5*aPZ@ZjOY?8b}7u~0YjruS(9qYJQ;)sj5wG-81CWV$nCthI}6P1 z2p)wjKC*p8dNz+;K_%E@TR3 z^}HsCRgz(EZ!hJZBU||U_Z6pL~;{8erxHB8dnB1n~E$Z%Z54U&vR*MRApGwew!m)BE;uv^(m}-#$}EG;=RS! zxQj}B2dsr?uvK7b{a}jKPKII9wR*0LkwSLIGcH)3uEJ;Mu?Y!zjchsNB<(0K8|Zqt zikBbK6w4omn~qP}Akx57q(bXw!pATAEhp^dZzp=z$mbtHW^Eo6%68WVO2rk@A%t9S zxkI7IP@`R{r#AU}anZO3uFYFRkpHy=Pi~Q_&@7ByG()g*%1e=#-rZ3&;z*Zv4Z_lK z|CCg31XbC%q&-P+M?$cDoh!a1|*!oePKdvI~1_|ht6wuBefXnIq2d*js| zXYwQC_a3w=Fq=Mns0s1VQvOl^bMxsP+U_Vi8NCcJjB}f0#bjRNKck=7PkL1NEe2>0 z3$g|>L<01$5qvr=-c2gh!N|~w{}jr;3&$#JMBom5kSO`Iw=3+pg73VAfMaaFgz+?ByRq6{bC;7bAw3iOvE5ug%=?L`XSPoo2DxL9wZz zkooo?yzv$vNpRNp-Si&#!QW(Iu>N-k$GKqQ?M3Ol5%D1V2Y^(?HlIv{pXNEgCU0Ru z`TMv3YCzUP$6{_Yv|XOTD-Z;N45$|kZ!#XJEF|ON;*1{9)cyUv2Om87U!UATX03wx zP8x8`!D;=Q*ZxUw|9>2iDHF$r8Qk)aW%I=uuIO3F(D^d!fsZ`~Bu4lR=QEf7Z!V+- zi_KGomM$n%P%z8vq}3nz_Fz6D(RIgjn8S26k9@RgokZ!NP#5Z>GTz0O%;ddX+u5i> zhGqiM%%!hcIWE7K{-OUEjhOj$6^oLJs73L(7Xrr&JY(Cl6I|Amz zXOOz*+H$S<%lSh`j_r)!j)a}qthhae&vrxTVQN_IPBq8L214s%e@N)}?+PB0IA$C| z^pKI(jwOM^XTr9hynKv)_ms3^jIhi2xx;gbM?&8x`B+p!p8N4IVn8`ahLyyI{*)=WMTCZSm0rK+8k}P-+s9rS#?v0!O@{BzH`$+b(G4d! zcj#*jo(CO!2yN~T*=JEWo3>?OK8Zu=OG-_&v0CAVc-aYSiJ7J*@+I1x-Ns9UW z4Lfw{mQEa-t>P1Y_j4E8N1taxrQuzbtU9sWgR72q;3rE4Zb4Yq6Yh5~rzSsEyyH<)XwF@Ga zYogCmP^!(sO-C%2x42;O8;-xyW22nM{Ua7V@INDqD$6f4)YNJesjt1<7dV}p&hKqE z+2DF#om@hvK=9DxO+K!IKPUKT^B#erJCwHkHB>|5-G`QZr3S7D92{tw{GT001x%_q zlA%)(rvI1ie+R~2_8Td>cS6`)@+UoRYJ!%9pIqBm~!An<)Va}sd7({Z+kMgU!GNWT+N}G<` z{t(>8$7_|20x5IO?CfNuOe&mRrBuyWl{ZRxsO3g?-?NA^qqyZ!&}ZD!myE1;Dn{Ji&*0+k6hEKFQU4*Lmo-_G%kqDp5*AuYSs-TKel=odBLK*F>&k53pYb ztl&L6G^2g)!zVA3=<-q;HO{1t>i&<2PowUd)pr@EHyPn6<7y?kK0_5vuObBQ=4wRc zp3!#mp#TWR1gNpJ=__4HwCkFP|}izd6;^0En~=kr^^dcnANaP&*mq z-)-TNmooRx0(xUD9@ILUldi)&s4{*DbN?;wyVK;=DgAM2W2~VIY!mlyu+qp<_ zSuu_rIJY|X@@vY&M0xNeEmh{)(hEsQ^zWL_-Mb}uwzB2A+yV_6ytuw@@o8f4xpr?+ zxjns{vO{Cu{=lQ)C&ezC2`JNK{UWgacy>JfxtK}45o2HNSMlrmzhu`@LUs|a$RP7S zZ&A{SN_I}X?>Ol${C&?P+@(2S0i^oSE>l_T#-jzdi_%ptajQnMFc+netAX{wb>Cq_ z!Vdel*zz9mg}uO+yNkAP1P;9sGP&nGIn2D>Rnd$m{O&E4*Xj~Rp@NH&3_^c+SErrF zbee8Gu#&t*i*SaP!V=NU_Z%#pZ@E07a_89bHaO_U+}=y|5qMW^P6@9 zwsD8^djj)6qA9Ej4Q|pC#fb9!AJ-`W-d5P8wiPNt;qEM(Mc?iR>_D@tkuc+1N~HjU z^tfuGb_8aKkKp|>N$ir5@UE`N@8mr?YBS)A4N7rvhn9%c%A09^`%Ddop_luy%q?yT zJn#y?P1`(WO~5|Thu8Q+qG7?yro&kgV=_-;TstAUkux2>=U<`Kv-pym;@S)I=$y__ zQMhq6C0QFIeHZqQ5zN_M*wd4k2{kdJq&~*m0>)($S!amq+-iFE3f&={Ce)(O_&Uma zPdh~4q^hfn)Xsj|B95u2jUUl16YW~-rfd)1dNUQ<+<{v%g8PPAwPkL7I0U&E_`(-D zLfcub5W^1yz6@f8AZBc2jeFFQC6+ODlEuV5{;bJ|PXm=cxb?v2{GY>kjA`fAA2Do# zQS(JtU@4_ues^ByjPbLG_Z@EzsblUpzN3sD}Sl13E22qjE6?|OxgiD@Lnb_eByNs$ny zSabNSFH7Z@Uw|}u(DsTP8u5gtD4~z~@p+6U+uUv0IDDOYaXg!sOX?w=n5K2dF)RZ8 zFC3`9+3s21u~JkG)g@nE?NAr-Wr3wew{~mCvj~iCAps4hVdYaW_C<<4rLOkQ+iu`7 z$LdCyBf(;*DLw=E7P$sVvc5Ilj#WDHy^EA#`a?H~(%pO4J$EaAfbYwO=k7HAp@ix< zBOg%8tDCLz-yg!=Pl$Q)`Kj>1Yw(Za5TJol*j@)(N_{ja6D_3ubr9vPoKu4-Vtf=C zxY~G-&x11A?J}LDzH|QLS_S5|en2+W!XEqzS{C=coTd@J=YT5$k zRbYcEU0L#2$dg+q(0{<;)5O-LIb{yUZG0{@zIO}x{jzI{FB4j%6%`#5en#l6#1vkS zEYrMGXw}G3y$ZM4+t!$58hjZ~_QjJp^sM%Seb%Ys3pazsm(|!8v+io~i;@Pa1zPM1 zdzbycP%~kxR2vnK6|H=KKH+m0)!-k)9#?T62g+qR3q!PEv|M|kUp+$u$RVKGXWoR2M9ma%L~1Wygq4=^!j8Z{n$MrVva=g3 zoPiGx#)K{nIalhV|8_VG)t7ui?A~P1xLeRSb5!8UxPbLaJNoRM|4Ccvht{b%PpOHO z?Cr$~$px$i+2}r2J=^9x2kuQSjY5XRI5>QvLi!(?XDd9MFIyX69<xuS`PeU|_ZTF{_wTG0iVURlxt@QV5{D_kOBqf)3nnjv3R*69)}YVfynATJ z78IeLA2#63{qU#flb1W&)E%pH8jhuH$BEa+lp8*hVAFYnl9nu8(1tWKGgBIyQB+K< zvp^~RTZKut=UEk<&2b;`Q06J*jsn{!=#s3@6q{2S{$a8YrvqjE^|xx+BO%c}b#o?s zOf*thiS`xeF?ZavXP$)b(k_K|^f2OCQ1B>>>rz(fY}=T~vu@u_=t18Z-|y2XubF6q zU82{k`98*5GKjN!_!L)kL5LCZz1w4UX2{m_)-Y|}=x2P76$}Uq$8EkHC8^WW8c(4> z(H4yNs;VwJsT{)7P1GNT_b-+WMON7fKOYS^qmvc2cpcL`Sl*_6_e0Uh*!wx-)r{BG zG0B<$L3PX+UoK7jN?BO>sZmX@+gO9=smSX5FXqpVDL|Ni|-{YzO1PuC{} zBja@6|)(1UsUt;m_!j}`DAn!hbX1KYhCo89(fRJ65q|H2Ji8Ih6YwE?GMuJ zs(AavnYI^arVlkaNyrKI@lt>ouU&wee>$Gg0#se@t�b?;+T+U7RL!9<9NqJel7r*W6V9LhU zUS5F3K?3(gUlDnG{PZ*>_*Ch;#2F&RT&{9Q$%`-Ybh%iS7!z(Ci*T<5BbVjfP)%ks z>_k0|2-dvLa3I7y_E!Nx>W7ElDk)iP$_Z_G4Rwi<`rs9r)cFgHOP4Y1gH+5Z?YL(; zBaM_C#YaCW7F?_Wt?}XH#MMUQ$ZoJ>{}TT394B+3HemhVkwlo~&h|kG z;+F=>Vrye(0T9s_`-^vQ{u2B3D}Skof-Ws&#$cj}tllKG|I41B_wn9joPQF}|H}hO z2oa-oQE!RA(zmNZ<0X^TwpSJG(Z4aVKyL1I>$J2(m%QsX*Xkv0gA`%cLWMn_RSpIC zB5QjjrCX&7U=RvG0M&HU!Dg5WQ>h^Tv4ttzMkdF#F7m|P?0)qJnx_d;p4AgHs++ut zCpM1ox2}d$uWugXy0`M|SoGr~G<6CSjUV==UYkB`+;rwh@yEDzU_*m8)XbJtjQA@0Gc6(|GO z0%T50g>uSGIq>7qd&=KvPMCsN9DB2w- z&+_|Xm#2`RIfpxZfmvKycz$iTS~&rLnFGv(RsnY86I1*jkW38W%Y$BeceJ>{4ei%^ z*+O@E4Bn&BArx1y&erH);GC%8YA5_mR$1!2KIB_6>JtEq0$#iI=kg27sow?-#sTU_ zf1%6-%QQye>p0V{db*%TTe@cYq0Y1l_f7Ouk53YxSWE!<#nN9vMn)$4xp!1Ll70Hc zY%I^Z(}J9ZiZmqPW{SJ#3G`64XOH!qE8aC0+8}wlKQSU(YBs}xG}@x#)Lz=X*-ynki#`DkS@9bi6&|dN`tu6tTn+-tB8iKrh5npE|ii}6mp1eqm`w&f3#!i4MTnt9$Ro!onDAHcc; z)e$ja0kcZ?CEn6qg~Pf-7`#J)ws!{#r7r}Q1s@+t2S>+>SDmxXmgS;!w+aFyJ zB~SR|Zxz=M+z&&MtYFSWL_|RD*Fb#&M!&JK5fdLD0P?{WulL)g!{^iPe?F%`y+k4N zh~>T%5yyfrUS|jEG1GG<$oWJrrkA!&S`+k74?5%vUA#i94jD!6Sk^-zn|B8ZcDTb# zbtB{al@Hj-Xg`iOYFzn`z9qY<5uPC7(T8o5q;QcB339A9l)c=Zq;8Q z$-xL|5L!Itw?LL~OtG=C+2}DVE9R`OrlzKr4QCUTR%A)sFZ_2rSH&!CZLuKYC-VWq z4olcb?O}fvzIZrh7}910fn2Ps_((fA+(+#1S6B_R{jj{sxeY}^hPN5zZ-gSEiJOFn zw!CLw26S@yNnBp8?bvp#Xxi-Xk4I;H6u@wNIKh7b0m zJ3s_j#km7eBc7Nmr1jH_uRt&W!}Jr|u0?}+pK;K;uQ5(L^zk4!{n6GGQY{wW`$ksl zv`z|SVc$ZSHx1Y$U(J8~BwW2)ydZNGhta!gOXYR1(4JR)Uf+e${f+92LIb{^VFk?F z$ZG*i)-Kk9}oWiD6d#MxICh`PJ?H^cW9KI&AgWpI_zZnfF<3b#!*c z1#@Xr!@~GKtxG-C5n#x4<}#QHyn0A| z0%;4}A4GPhps!sre}3jGx-KVb$E2i$&Rt&}yPj{=*e`u2@jP2E6bACKuCDIxO1f=M zT^%);x2mctW5X>bppiu0e@v3{c~PR@2Z3X=ucR1CVa>1kYCvtSFJeUZ>7j5VyVEWJ zBH{U?+7aQg-d1t+j^)G-t+=^V$OG~EWmj}aZ>+4cB zHay5(@JY7{w`R_}!I@oB4R%`l0lohGn$i%LB099_=zJOpI?N@J$_xXAG>(5-Zgi ze!pN_sxJ!!gRQ3rC!1MKS_jJm=v$P#J$GRz+}xQD-#j00L*s4$d2VLfKT)6Ur;@}C zCsZL+dSk1ozv^^fVg7w1PLbpMmjw4up8LN(blrjmmS=O8Q4vH9N)j;x{eUnRmiC#G zh~jl@NxRfe-QXl3PeQfFrIVr=A3S zF_lDJsSg||Yoes_Hv?WG8QEn_)^1#m<4o5B$QGKgN}$E4$kA5Rch5~ra4r(%N`j%R*0#haQ;P@vh00BA`v@dt0tR# z(fy;fICY~TL`P0OZ87Ov9gjMoLx|}1`JP=W!=NB*J*dc`fX+K*5dpHMhmW=6J~#Xe z&~K(lAQ9|(zqbJnrvO5?02+3>xJewW1>-aHWyD`H9z(Kv*bh2hQO0gaZvp_8H>d}| z!hX!>`OO!|=y|B)nh%GwKxCcLd1 zxz7IwWsTs8G~mV&Yzvl;yh>cySUCS-j37uIdF1w%Wlde73iAJGv{PPD+*M#z=!=mA z5xEUn(sSSPnpF?Efm3TGc60G2Hex?rJ_WqMy+%x`DJ3g_WQ=5E{~>|GYd2s0aaEZ@zBfjK+*$=G zRR0OI2$4YR4aNV)9r*@Qwbl8#XPga_F7LU`wMf+lSLt^L;wgiL@m4yQM_A&pkQL%f z?UMH!s&UUr<827lDUF57*0(Xk3$-cz*ozYsn|N3e+%q&u{U0oJm~Tl!VwPS6Qb@jg z(?*7a>af_E2q1r!0w;u0lGG&t^EA0~PYo7wVbJA*r1ln`x?4*&@BJ{85?XeK_>#Ne zZKv-aKXJGAX_c@2Z<(iM9LN|#QQ%hk83jfqHhwX%v*PmArtb8wT@FSX*(N_;lW zioeWvO0fgY=Rs|HDDreJBh6r$R$7f=On~sgSB&0$8e0X9!UCarkND%Q&RDiPF4-hq zL3NLYNe*z`Xy>l-z0c2m1X0X+++Da969hYNw(iZWQqAvE+#n|y29m$~A2=fk zTEk4Y4nEU^mKxes<(A^D;uQPyno1=om~|xMAr{3{W7UjUX06%kbht;9K6$n|K(hJr5F2HjI>l|%8Klc>WEk2&9#$h)kTk@2+fU9}Wp+I+#^nqE~^opI) zAC~kgP<{*o>>CTIl12p}akWR#C7 z-HKX1Xv8)QvP)L|S58UZBKn-q{TOxh+qYHK!U>=K0WE$7(sCn`r1x&KdlcN#DE0Qz zslIzr5-s8w9_cgZ6DnD^F?5!$f|H+>sHqR)9A`gNuQ+vck#E`a>S0X|0>91D69}Fp zV3aCnq*c8Hy&HGL4o-uKH}?g9NWGrvEw=y$KTAZ<_wTBxPaGY;l8WC@L-G$Vs_0}W zRbHpd-`d|@0CcO#>WbzsI>KYcn$`x@k%gFRX>A^(ts5Y<1zkY=LlULcwyXtNwv{53 zothFpCh*-i+(6wR#L5Vx2FTw&=j9X)?=}3T!Su#AGSb)@!d%{W0+X7Cn+|pedJ9j{COn-j`N3c7B8FrudBx^Zv!b6; zIn;5E+!|7!Ge}PrxUCTiRskRtZ}6AdVIri5*zN6`Octm`3jXo43?br9dH>-QVR#-E zmddU&@m{7xDJIXFn=_{H6m8~MI_}FJV(7=Cy2j_<=3nU5&%0U6iQFe6Gu9anNC>rn zqVChyc;tIN>i~(L-plQo0Q2PW5di3U^@_Z+P&V%VL<_&l%c?T@c7J ziEt!;qowY9wz{MqzQC_KKP%%AVFar=qvrvXqUsk;CkNnJ=;4WbNRNwBTk4%}L#HLY z5a1#VJd^b}K^`H)74GvdU&`||@FNT6nLs%_1~QI%Ghjny%2B&{Pwm-1wRGg0A4t!> zeoN^{T6E}5rK{9fTPqA3V`gTCmUS!E>cz0ov1Fn@x-);tALM(CNs8LbrQq+2j^pQy zAohK2dviXex4g-UQx7<^6^P?lPQ~?jDfO6n*PKMb%jMDTW-cM#_Vx2FGK8bcS}WD| z`+VRF1V+?=_-RI|IF=_d^Q64k$#^j=FcF*up#YLQo<&pCNS5ZU_~aKwVbmXhwheOQiS?=8Ec({E!9a344 z-&Vs>o#)7zh=n~eW=g3zz___j41rNq65hW*`4o5Y-t+&Ku3E(Xw_|@Wr$tZaplOWNu8|$*=@Bade_(#n+CDg@pzS4CDW!PD{$^$R%oz;LndqLL|)KIkc9^i&mAJf|7JmLV<1IVL`_ytSq5B zIG@YkeNP9}-{I%k)EJxvXhHR_-!it3;+iP1rdOAz4dF=Dbf(h(hg~3)uQ__)Ue)q6 z(0{3W?7ecO2BU%~klzLUgMyFyGmh8`lru5MzyRdrF4{xfgJt4afcG;ZBI2Y!kUcO% zY^k4-LyG0~g2`39ie>r@R(CkabTL-;(5$0QVR3tlzY>Y{ftpC=)a)CbIyglPbMO)& zor6*AVu+}xtWyY-Q~|vE+mCQFAp>zzQgO^K1g`%mR+Bf)q1Xf5p(^H(sg`akMAaWsQZH(YG3Y4 z-dj#yMb?}wEG$NI&0ff2cQ+k4lR*=el&&t#tmnQ`Mg0`p{{6o+*1^?ZSwd+pins3qlu6q}j~URmdhu@t5dLu5`}Ky`(rUh? zTGYA(Wk8!R0Crp3P$m&5;7Du5rKGU^s0GaDw;Oi4qJ1t*XeGUA@hEu%draG}&kV2e zK#472-ruU?bEc7(pU)^M$>4Lf53BGzbyy>CcUzGy1oltDIAuqddo@VvPPg3Hj4o+s zJ!F6IbxN13c6Hd9|E;EnI^r!jb@Fs}!b{z=3G{(i^)Snise>`JhI}^Tn zFE}^sKRCrcmNXtDQ#|%BN%oJp8Rdfin{PpWekj73XjD{GtujOsfq=eS1itZZ4Hb!? zcYj6~HYdJX+Sqi0gZ)A~#zp=gzO|j5vz@7bJXQjlf@3Uzo0>oRO6j|9`AFH$zKADS7s5k2}vvLih2L%n8>m!yfgUg-2MQCe@shBjEC^`c1B_4)Q=rKs7}>Boxed0atA9|sjJ zt*n_DYkhsa=VU{Dy($oE+uEdH&X8JBhg-x539;wo6gl9bMHzM`DnTZ zN?-yIou5vl0}Wiy6Unob6HRAD^+NpPv`l!Bw|{u|f8jKHiF%d7-pjP35#22P!9im4 zA88oLeL)Kg#z@W`Np%7m%pNXr+-ms%6~e36=nKfuX^)+z+Wjc8g*BJ?-b60Dx+RYw z%f~#P=%$i>Z=HSvu^tqkEKu#}0s7;JP!HqRuYNlg4R2FcYcBC7%()hsWDhqqcT}Fm z&=M!y)BG>gsB18adWBoz1xD!%#ph!I-JW8rAD$1E={rZQT|%*snkl)Z{Td>eW$_4I zY`gZ22d;?0E%Grl%i!DrqV3b&@87=;h!F^qr# z;)}Yuo=o{%lyH?Tv>=dl z;*oSZ|C>>6eEsGKkQAj0-)neJJ+^t(5duc1y2|6@dT^9Atbtp>Af-{bW6pO2^yz0s zX+57aey>>KGtTn73Lc_#tnc2{a#IBJZoc++x$(DiQd3j2b8_NRi_bT=grJ|&)A2I$i;G8;@p(ly zH$y?h!e-q4Ve3tXY#L4CGmHXVI!H#HSbF|FU-P4~nJ8)ji#y7N8s4GGh4GN@l5&ND z*$W+xx3GFw&w0D_{*zw~_b7N2_obrf5icXKh#)NqoE(oX9UTqv zP~oigiezE%?7!sgFB7LSr#T-JZ!(AIVl*jEb@$QlCZ^3x1b9-jDwebAVQs2mBJy;z z>b>j~jh^InO8eRHkA8G`{>`?4)nMg#SER_X=TU|PsMwWVGJVUWzAV?u?=2#`OZ6L! z?x&W!>=;y5SLZY}(cM+~@NItnP87XF_}Q6z4D3y}mah9!OGzxZ1LrTaYj?z9Hv-sOzICet_MG8+UGJZ#E4bw0z(8;M-C4u_oH z=Tmc8EiDWa6BFMD2D%uoPvjdDn$PWanvUCFRm@)=&vkQ^HTu8!YQA6rhvEe= zs?}A6^s(LvRP=Cp7aZNw<p_>BXnFzIA#MY-P3sjcpiq;i-1Usmd&jX% znrzp$?66Q$*8DBm&^rv!!q|*PHB%K9aQbQcZn~RxA}j_^VvT}lZJ^1;KUKiu_kh^W zH_+U()Q}3HBY0|)$W;eIm!4J~oc8McEj&jV? z*M;#c&r9UYG8xKu%cceudToA|88n51nEk3Zfg|Pq!LTdQ$5a`?)nK=R)Pl*|CDQLbco4y*T`G<32MRv0_4Mm(yT5-A zbs)Q!sf-08X36%L{WeA+Xe2m(Q%%2kBf{o{u`=%JOM79TN?tS%N|+HH)8H=)%}$r; zw_`A(BSTVuh6tJHvDp_$*09p_h3vN*BO#lC-ZZ&QyyS?GTRGo5l@nD1MTEcVqbrGI zFaY+JCtS%&sBm3bG$gPaVi9zB7?B!XUevtpB^??FHoa8guvD!r?}`oLP_Rp|LOtPA^ zKydS2N`Z&O4n@1-2|Tmcnk~3yIRRKVzxx&+YAoBsh#(VDdFEM!O?fk%ycJ3Ky`&KF zRmRGIUjKyGY&17-|BzLeC?&Z8fdA@s{20nt)dc2S`B2PQawMcLIixb(>lsXtpr zqe`1~DMz`v4|H=}?QaWke) zgH>akYWyk*Z(s2PD8DsHLZa;&k0##~ef-!~=6%Su9P+w-0a)-cddqE6r5ZrN5L)U; zMgcYO*m`&|cYW@I#Qgxlt)>Xtq*W^#Rd4A`eTbs_^V2hJ`Mx|oP@55mVk5`7tB0x$ z_um-pzvX(1t%m@f(J3B^zLclpD2f6ETD>|S7K217w>E>1s|coeZdt}F8yl)M&$+ps zX=&Z^m~-@nOKo!{mx-@xneAUsdrrzzH<-2@JK_rp0-(vhNXTI5xrcGIIkC9U<+?Gi zrmGtcd{8hi|F0%<&4Q*phPVM6g`Pp>#^;mHc{}CpudCX1F@Np2VA>oMdksqzj?ot z$M<@6u>m*DkJ@&y39-EoOAOmySLGBGU{Yy=k9}Z%ZvD%f1l3M-5DRo$EO!9^J<{fZ zhSoMUV{$9$JEI&>V-Vb82}+#pZ1CqxtG0XIALYIQSM2sg+Fw`fU92ECUK119`WJkX z1VRcZp!sa|sV4{(AS595C6 zF7ny>`=f&Hp9Pneb2>rbVu+ewNGc~!G*c+jZ_s}TKC9vdkU@$r61Or5ut&F91JR>nYZ^CPFCf+k4LCVCCenRij>Y?asskuw#*F2#y;4Ut#d2N6L4)_G)Q#Q`ee z!I6=N&GQRvsbePc#tfWKk2)C9Z&-lhxbT@k9vob@KU?=u$FZZ_Y*`FrsO%tt#R9?LX6oJ1$Yd zW>cz7Ow?YpAL&7K&FZ3qX5roe+V^c2Py=#L1%JuY{lqgtboXcZTX{y+hUYJ)Q?$(j zSW?ld`+$iz=9&r4HyM)mioF*elJ5XfAQRjtIniV3@y+l1I#~Uf zH(AG8^(jr9X|RYQZ@wjE#%qjoTI6>V&7P=NX-@?HY$=m|`}GsB1Vuq$(H5KccSN2u z>5Is#^z^=}FtV|U^Kk=eIlFJ@j-C6-kGVxGpL23h;TKnW?R8S7JbaU+B1y!E{;7oxM)cy!~%ez5|@S zNhJ`PoZ`b=S4><_aHC+GQY{@K8<~o%3jR<=`e^Y7l*5J5-rB$R(w==Kh+juqpxFZc za)xRwAW8v-jpVWaO}u{<{(pJk{5@4-*yMhz@GP?&w$ZWNCo}m2TA-31ytE`s6;F2W zT`xKE^Fo;F!(?fp%10lN;%Us;maTDD#GoXF&;#*^>r+@fm3`#=%V_aD(60=FZER>C z?~#U=F$Phs1++O4!GCr}H>rgCkGCv>2e0S+R&Lw9>dQ)`vCKGv?5!k^afv^&Wy>Co z);NM>!z*YKgzX0tk7=R2hV-~xuCCuY*ZJ@N{wb(J^VeV|xB;IHnNm>#N~LxslcNbU z;oFy(GUz5q{*|4Hkm12|ITjX6?#SC`+P>vQ00lExh$Y}kOo>5}k&(U6T0wi4O^5oL6P53lEU)&%ZKo7g_`#3$NMiT{ta4aM9*3+*Ub zGpo!Fr$3Hb5%*kG?aQ2P+s8DEjAmRn{61vb!Oqnq^J`tT-rpfxlJMAOcvYcOUM+~z zm%QP11bF;3Yo_Lo^#viHXoia|jTF&$#m*X!q_K4j#R4t7cPpKT^X?cXz6zKSy~l`z zbEd#k3bx1Ig8bd~dixY&)Vs<@(lFw(nFnr8+(r zUuz-loPk$1qCg6_r@BNXspbS5V)@2(5Ei!Ahy}Jhf@Z925*PHmC`ypBtrf0;{;#BP z=yKt95JT0zFE9*cKGeaRNeT$Qg z-CEK4P}ogJ=kbY_R_D|7&%GW4z_@VyN&cTEn6k=`iqCSJ$A)IX=j#k9asG#Q;TF;! zF~T*dtww4Z8Aw-!W?nbf^=v%ghv`k9I~rvC9{;?T0kQ=$I&AVO>yowM6w zGrj+DoGwQ=9cgmSEp?D2mqjj9E@Q1o>7?YAUtvf35hHf-BQp=8q$8YMhm<5t*2Y}Q z6w+LxO{-O+w3*x4q;x*t(eH82KflNM=l9qC_c}Z9-=UIejQd@zx z3Pe8~xewrn;4)BXsGoZSTm2Bme{BJzlR+x}jzT1X(;BqtR#(3O7`aAm^aWKA+cepx z=t*0hf)4!6{#fw{JL0^j$<5)kraqQ0Mjk5ro)=hQ&U8PnexeZ*423E!K)LXp{k`)H z1Adq^zTkeQn(IYwQ4l-qdERJ%XW1TgSGRl2(0OhVh;E@rXEWCt)nj(n!y$l`=ii`& zIz40VFrLhzBGOS`DSw41MlE%h%QiPuTtF6s%uRH?`Y;|*1LviK4LI=`#_oZPG%t$u zAYoJib?vsi0nm{uv)4PxXR9Mk&DN659gz0w}${xG?Nc;c(9OuThwIssn@r0g;rdw2 z5~#fQtN_sZjoH2|$Y&GO%&Hd|YDM;B<0L7UNcU}M)> zJ;afdQU+ko4s?>dk2%~$adB~t#XA~vtpOw?3Prw8uB=FjyHHolP68**x(KL$%NYp| zc|AWrzm=2(F!rf*$@OO<*+u#2GMRMMF0*X|m>iYpQp492JalT|NZ4a}$a z42H4RMGFBCSuCvGSPJw%Ql%ESU#I}<`#TMpHzG2xQVH}VGVNzUIEn`Mf>K z15klO-rg!$(O#>LR7U7wrH&^NBu)R_fB zy2B7&byHWl20je^kfjjTHan8>-&{mc53cpd464f6Z~tQo%?jZ9UBL}-ofA-Mp?FhH zSzQbxb-@!+a<)Rg|L`pq-}R^P-ZLP0P&&!F7rT@=Rqh}hP^IQ6<=a5PkdvWclUlDM zWubb`7b!swe0)s>gx)O{&C;2^+^%9yO1Y$oPD(PGt70xl8u;Gc z{)2Sf|4PpN|F30}NTiF6UX#g6h`oM(+Mlp?fK9yl@<7CWUc3o3?c$z~Xe2l*%o;Ld zZ1``XbulsU?c29Q$h&(sMc@?Sv<~von36|kwc#%;D3z7JNrbQ2$|FBbs*kd}JcV7SG|GC2Dz=JMsrq$vHlWH&U^sdggWWTlpli`z7%g zu-d2=oGT~hcQwDv#Olk*=yTXhe;iq7kN~~|`yd9#@3@>fRlq|o(uKoQjc)~&{a1FP zLSN-418_fh`pqNR?p>TbMOu&kRDC1weomP-nvDH2`$la#;c@1x8Uj)yzyay;^p|Za zDW$v+uzrLiTz`0fTGuQi;MZv=!~h%!uY?uWPZGtaeG5#)-EX_Qy+~L${-{wuP@ef^EyM$Dq4)2|`@y-|Oyj5>u z`s$P|sI|7lFEEnWI^y~KkjRgZsQaE$)tHeL-CuWqTOUuH#rS~VfK0{#Bta%o|LbDw zJ?vaMVXnwx-`!lUE8QL~>aiR{?AxN}=ID_}AnDS(;iNGWQ4Z$3Pw<}4kqXTc!N-ty zEn`fvZ`K8~)czA|$qy-n;kIad@>NFj(xGA#ILcp5ipJ@(o-IY# zW%3u@t-H!IsAj7Q0*Kmonf}Q|c@1_c=fZ;>QhM()clP0|tlu6oCi+6Lk2;N5$z`-z z>$z^CBF<4eEqjg4jS2XN5;+++yUclvxPRByEi|!p+)y+2p%=5yu0vVqLvq>rR^Wv; zjk~XXuHfC}cQ07=XI$rogLKt;0LggpW@M^T6WRub4b1mEB z^W}~4mrql<242}Yc#mcyaRB2Z_#YP6Qf9F;pwPjPnwpvlG)(`gd&M-`aqfXKA4Z&c z^)I&@!9#4KLhYKA-7Wq13Uu?Gz0pK~1c`iGK^W9H@+(>m*IZQNCv{$cQVJgT=YrY)eFWVi;yO;f-aJ0+I;kBEl zoaC#@G4g|#29os`l`pSe)x@L2O_5*PsLt|wZm(Wp_5QcL20Q(+c=gI)R6+WSwzu&K zV$Y3SA+suSa~h#LcU*b0v}9D*u+osdD7HCvynK@3X0?0m4y*rrh+F=t_D{nr1)GWn z8vHLmCg|j`kx;QoL+A7@C8Lovt!eQ7xvkK?#zoRNuBFr`MXdJLzcw&lxp3d8qP0s=-h z{gkdFfrh7S%V&aBn_yEvKeFh4NmW&>ot>R;0DhIOhebV|P|v0?&mEX{-p<2=?>#VZ z%bO!&x|HH$MO|HTd%I}4PT8G#*P}$}X22${=89=0jH~-;D5lb>w$9`?JMiw&0d>-G z>c0I~QkuevMDX&SYYX2O5MySGmVRvap_S4$jJ%27-LbIdBQSV?Yk8BnNf0P1`_!>Q zldDwS(4a8CB^s6EX-eR;eZg(Co0wZwDRmK23xfhvQ#A|>%53N8pBOZ+n;cj2;oWc9 z07^P&lD#c0E!Q_U9UK2T3PTp8=CF$Lt8Z^TcXooL4IQbyz}z2d5kj|Myh5E1gIxvI zx$cjWaQ?C|m6Gg|7V6Mc5BtZRwTYc-GpM`yKKqud^gDD0=1ZeQ#j36b#k@5SI2sh^ zYq;l3@k-R9L-z*mqn9W_#Hd4`O(UB1$4{>=f=nSheRefbS7+_)Y8N3iEH{&r8l%ig z7iSsxh$1B*6CR~WCOqgPz`;QoV$Tnfq3O;b@tfiV7DqF&Wj~W`DGWr8N&r z&;3Ot>hA6yA4hI(Zq*LdZQU3gN3m59FWv(N%UzxG@bFkzS=ohAiRNuae4berF{odT zXK3+NgcBHw%46pJF;{s6VTiSasC|FOL!<3gg|X$gBRH8_;9wy4t8;U-_a_04FZ2)O zV#hp}7vna7(nsx_I#62swvJy*ZodH~(@LO^3*;^;zinm1H9#To0hc3rmh=6b>6^}n z^uE8S4ulj60+mTn&VxvX)mG2`M`?Wx4V;1=Hy9>{^Os-nnNqDRF~77eaJ<{)w!%pp z90L=*mpT8(BkYXlNte21G)4E_zDQXWm1WN+*rX!<(axGL=0OXiMUWMyx6_nWa^cl*&sNmvAs z0dk&n!gH`#v_RialjU@_uMW7QLT%p2gF zE4}$CQ}gEF&u19!szj-_t5rz3jeLxfoEy0 z&)o4-{XufgiLs!<4UOA!4X~KA>=lyhbNT59;Z8V@5;cz3^$gO`#OyfGlh{HlXiu zrO-vkIO1n*wSDMr19SH?ioA*toH&=tNyb}9Mf|G{Pp~Z-za-) z`$fR`9YSE*AKwSC_6nMhJ$9Uvtk zZFGZYCla0k%Kp;j(AUur3vE4`nlq>NY)zGct&+Rc=X|PgqrBx%iy0N8dKt^}gIb z2Zvar0|>Hr&Rq`NKqR4T6b)L+8FgX%c&cKlsN~2yKVKNvb{MOxe>O78zJ&*z`BZMo z_3Ze3^E{uNPt(p%D1dj)?*>}7jqf+ec{(kCI99!fu^Vq%j|Jj=vQ3Il2mRUAc+8cq zcH1~kj=d=_m;%YSY57Z>XE;nc{psPu!-}IZ%+kMUSKT@rl=HlfoY&Zi-*R#v8pHy$PTd-Bb)Q$7h;q3`MK+;tnl02&G!^aC@ zR*gJ$|8%Bq7JNArdhf7Jlsu#=b#m&Ib%H3`P~3J&mr(U!U8Xo}#dEa&ablpo76|mY zQ;c~nW=6cfQCWhu3%C7q=bzuJ=(_Obgr2mK3WI|wGTSV}66>)-07SR0LAT%o`Bg7k zAgtp!>tljMlB5eZvsFm{7$V2gER$&H zNUy-bCfl#Wp03l7rpqsbAQ{0niQetUAe~Y*-{}i9CNpXXGLyPszI6f&0mgs2S^yhO zL>#@O=T2cAJv%ObJ9fr1j%J>U5#5afxDIPAAJSOjaUORXOJclGzxdjq*O`a7w>uOS zwS?v5*^cDnE-bOA8Mc+b+U7gy@ZEO>)Ya8}v6rD%39*-ShCYs+FweUM(45S2<=7`3 z_|3a>syO> zv>U$_P6maOu#jYzS$I`U)$RZG{vI;;;G;8oE0FaI2p_lLC;et*06x4%v6w!v(1O3O zFL(Nm|6Y)S5W~K#9JEar?pmznQjPs*+SD?@g93y-)(9DUSR&P-t_3i&CptlQCZIgBZAk-b4&zAB z!jigjk?ph>&{WM3Xdcn55x?Kx!Y{Um3q)nV0AF2E*V(C)Ga7XyWA|~=lLvNDSCrurd;uRld;NmCEwcm@Jvl79fmaZ1)>_(*vU*27c z{dl^)LeIbcZa%<#-^jCeC#ytID9PDG!ZAX6lDAbhNCuR3B-z}gB!?_ z@avUR&($H9unZIK=IF`4Z?(8P$1QwxGVk$7$yF7N>3vebH2Vy%f9CbA1cmFNxctpaKzU}uHGhNpYsAuUrhdxZ1iis8a(G$l;aYF=Kfj{qd;pKKO(TZrdlq2 z@9|tD@aKCioo`g@x1}7v{_b1mm1J>92+1m_QC)za`%(*iE5uea10&arqs(m_F?$bd z|BQxuiJrfeFI*R>cq^SN9OfmD;WCMeulwP9{S7vJmVgoyonyr)CU?!~U5UQT&e!2g z`^Yxx8s0HXEqrXn*+IHguBHpMzunYh0Fwji4J#%-8ttA2c}HpDuN}hg%CwqY+1`iz z5~Y+zN}T4DpWmXZ^j+)S=hVRW|4Bx|Av2iAc#$cgRsysg!|G?;czxAsyy+M*+5WqL z2F~z>gss`3pwAGiI8bAJD;*_DI{lq*Z9uQjxhn>w_?PC)om)S<>uJc74{c-xWF(_0 zonQmrZW+z@$j|1>LKL>36v6+k7Ih_R_@gRHr`q-!h^N~gaL#OsgF{D0?|Q11i;eLA z(VdAuW~gBKWz5{dXsTGIlOj>E_E|<;oIvq;fOiNuo^WsWMS4eeLr)F6M5Wv&|KmPQ z_q|67IBRbn?39p0+is&tZW}LBJ<`5Cr%P07hz-~~C?S7bacU+^Rn@T~om6Py{m_)` zdP811n{$x8EGwpVR^t3PK<{WRP5=m690**YT$DmvH9aJFgT8z4kDNs*nm9(a{3x4l zCCYr_-EOI!ng1!Ld5vZ^Q2HW7atcj`z8?+tw)YQZ_PVVmFU-W70d4^DR~;E8$AQR0 z^Y~c8;3e7a;t={B)dK!2Whq$p0k9cr16Odn$#K`d;xiv=qNBp@KPx=rEIv?6wUO%W>IBZMY@_HnHbh@qAC$@ zVHJYDgJ8!k2w{>wFntMo)%aJSLg<;85|%hy&PS#C4Pcn+?uesRe3Khi>1Odz1CYTSQ=0hVUP>Y`JJT7AEQOL zDi{%@gE`HUi~);0AlP01kOpA*Wssi0K>jSpQs~{FD=Q|e@#PRuW z2EBO;M-Ml-Rg*V@vg@N{-wE4G*U==!EqY9{D$ua*GSf8Y3$yw*69dhm`8HOo&Yfcj=lnGnQqU8y2 zKAV1d4I3uYtFD~7{+*JB6JleaZ~K?E(DA#mf79emeb(@)9!~g@O%>}$F7}*45t-}J zOlK#deLo#(vrD4~V3wc^va^7Yf;3SnKAM`XOVy;hF^wo9lL?{_1jw`LWk^cMVHq}= z@U8gp^a>Yfj}CPEnF(i2O4gEBP+UI^(YKl>wDr*hm*2fEhe*wc+gg706t0ZBTHl~` zC#wx{2^lNl%a_!*YOYDU-Jmzz6oo=Xz7d82N8_l4H?Y)1O6nEf?mc9dI+2=dZIjrS z7MNS4W&8y4Hmd&GGK}G!o8Yget6XCZjJCdxrV6n33(E{<#ofXi{6x@I1|3lkz|1Zv zqgfLMd>tpjIpNy!<~h*wsgJZ$_p-!MaB;yGJ?papW|un2@kbRz8hVx+>@ybz8XqMJ zIG#>4eBT?!qX>ZHoIF`F=G~yzkD@p=k&P*}*GFQ&KsFqwqXE0N9Oy8ge~pH1wj#B^ zJB+!R2M}v(qGQvx=dE%op5o-9P4;KR6@~{=Fla8wew>GBIS}gVH)WyLV%VvW>8aKx zz?6%2Y&0E*lF={|aIzy%HO`ZpRU$nO>wxDt}VvM)b4nG5fOT&y%N%6mbp ziOZ9CjE?P~nOyXU4xvJIav9!9r)Yk7cxvs>uO}DKjvRthk188zt+nx~!~(|kr+vL` z>o4(6w#yjwcfV(mOrp4dfz(Qf7(7u`GU6L5#!g zkY&fS;nlc(zr&vO&!M={yj4LhrcQL@T8(G7FCNCs_-mMaqB?nklK&qK&TvS`W2D%y{(*^L)bAw*|$`kT#UC{SGwe?UN^!8 zfV$!=!EtCjA=LSA?v^lINUXI~L;uiDW+@4!rRm$O8dsYQzVSl$R7{*+09B-Q*QV=T zfNUD8P4eecECe@(7!7}+d=_$=TxO}ybfI z<;2Pu?GEcQ0UyeGwjkUe+FplAUI4c)&bcypmB%KMsKyT5O!0RI+%7cnjmXG&E#JF5 zIG`XSGv6)8SIt51<-F0`L)&MnKw7l0&2d_!E>&B3?r|_fXmf%5 z;K@0P(7*WnMjXfjR=8RttX;!Z=?H2ruigm6RWjXD>f zV2mo9@m92ZUo!-@vZuwjqK~%u7!lka5yLK`M^Vbi$Ct6x-LFgCk)`hgCuJU)LSz zP#WjmReBE*sADb7rInl+A?=B(p;1>jIC~k~j-9C-)UVx^{6hPJ5@@iUoLhAqyBcWf zYE<;L5pkIK?-sh~1$lXxAXI2u9PG#Xs1>^TRu zK3cyTPU(RjP9Bwi>rO*+^T^Tgx79-=95*doKgFscRRl%sAhd-*JlYTOq}y=F*^bg zus}E`+br``}*!fX+{n^yZv-78|ZGoOqvy=O?id6 zr+Mb={p(>!bdq8Y@af|s^rq(e$xZju(e+= zi}KKRh!R-tv#FL}VP?UhU0Y`Qn9({!14E8EVgbXyjXPVq#yw>tMF>^F;RMrg(bh@h zzxIb2Ttd}so!iSnaFH=OZElY7uDUOBMwqdHomadDqSdYblhrTjZejOa-qs zsW^mF5QOUO)xlZi?()zUk{shCC*5Avnze{Ioj(Eiw_naNUeLN1<6BhNQ`DM?e1Uy+ z(pJz_6TA&Z@h6VXCxyfy(=$-RmR&hRPpk$GO4m)NXSF~#v9QyY`2wU==EqJ403f!$o_A?! zsowfjEv&B15614)JkaVxMBk*M>qYo-Us6 z1l|<%G;fmiC9je^K+ZPtPJ0=j?{Ic}DT)Nfu+wREHTm!^+u{dLl1=ZiSDyn@fmULX zUC;E-w~5%JrQhwYT~+o7RS^(^haLLw^fziBfa#)6X$oEO$Ru#4;Gw5)A7YaIG&MVbQ$lgiOueA)i`$ z>N+ANo=!Lf@rdKP&e${$p-ip5X`zg+lQtBnIC;y=rB?Dy%JYcNPeIQYA427j-Gu0V z4o8WO1>(~wuGD#Oj~nK>e)I9NJR8^1(HR^brlF?(`2yq38+?~mgIwgDvJZ;$yn;!L zI6AGx882iC&dYO?&;8M+6~*R*d~!S^Cag?e=RQM)KD8C6>wW7Mz5Ed$=oL#;P=OS1 zNkE^Bx7_8OxlpU1afp`K^_bLSOl_x2w3;9Du)}2*$Z?5zaupc&`i;a>OuaR23@>oY zODyyCGGCIY@+lcs4Nr?k1`9zzJ~QLO|6N{v2ECcl`@Df~T8M zF`jleGxhdAT;kKe4%*!Mhe66xD|9R8MEbsDg00e2&t_{5Xp*co4Zf?#inh_<-TZ!eX*CxsmpesadeMxugZ|VwH_08%m4O6z93bWGI-(R?4qTz zotuuslwRW4EDupLc8dH_A8RbEr{xe&(!qhB*6OX1rOb?{l1qk9PO92oR`bWEoTzI} z23D5#!B4;aOoMb+u)-5dkc>Nf8p!fyVxu#|EI&<0lkincrRaW!)(hT`y@{Vf%cyqb zF(%*`%l!vwB%I0*Xm$~cYn|MATeZusm3);YG-+Wo#PQownST!GlN4x@(l+X@@`=d5 z*A~=dPFdLpgVVEnzt>#mQ&!XC@(3Zz=blO=*aKTmE_=`9floZugE-P4G=*KnfpRQh zUSpN@Qa8g7gaoNyq!k+S4Rh5=87;N*l@~u)uN~GbkFOqYS$Lalz>8ike4w&gSF=Jf zoI-)=SF*mhOXyh_z|3Uu&ch;i6)4+8ZTtV zk2(vTAB#&ooN1f69aHeTvz=p7qlK;NR>Cmf}z z8lk;BOZUzp>`9I8lpg3)*Ia|x7*Ncr^nUY#Ec1Z!Krz29e05fSQRrt*PssSHlOm1` zj@xTS>R%qzelji7mpEfMnVkYFf?|4@roh~Mge25kMjpUN8fI;5>FS=c$XCYT8+|@Q z#sIQ;C)h8((GttEz1lAC{J=vv>bgV1tp9cghZ?`hn*p4#i%q%48V6ln?RGS@uy77V z4#A5k4M}G)t*?r*vjmSF&mc)j%}SfZ=0U5G|Kt0xp3a6HD;o}_?WI(8&CQTzZ(z><4c+=FyDxK%y2GD_#2+A7x%-HA)lbA@sTb{1Qmi~HQM81naZoD0xI^=YQRWc(>ZLJ4e)dy%n{NNZ9KyDYtPNxon{sOH35Jql&cmi?%Z2p?L;ki7aS@vp);^BPOG$kX%8S_JO z`pw&y+n(7W)YN#uu0+G%64g0pR%sGhJU4_OKw+lKNVbXDbN~q6g;BRbn zIW$m+LZhY=`nmASU3kx7q%jqb7+3b|TJRioQwJcwg~`O*NsJu;?U%lyRiyIVTv-0i zFOjkS1G|SO>VV`ZkGmvB0m4uhi+qP1n&i!xdcCcg*w}xY7k~L8Oq*oLL<<(<9K8IF zJwCK1XW~x=a-8~f7;cZ1IaOJwE_#()e;j)2@oF~Kakuub}c1y*Rw|!drq`rgRWHL3|JEtA$$XI;Cr!yb(*mjV@ zEHdS(ey+*k+8w#R{yixv>CfN4G$JAf2A0Kt{``9XfztW$7VzoQr+TJ!Bh#n6V1Yve z)GptjZj4>K(H#CZ_zi=fW!}|9jwkZ;Q`KKsFG`1|T~*#NjJ6U<8!MPP2L!h@U1qH4 z(Y*{!fA~Ggm=Qdw?hCSV<(W8_dhsy|^h-BH-6*QFwKcQ3S#THv;ZsynY8Z#Z{n-M> zp;Y+yfvdRf5Kj4#`>Mh6c|Tlo{t;DLm4hemNAb^!uF#sX?FC|VvZW||j8xgJ0HKgX zAqqG9wrKX3>gWMDV}`>%HPoXvxGV^x(9JTT*&H>-Ne%aLl9@7Q7z^dC7klou>$RJ) zms~r4lWj)i1y`Eqc$1}pxnx$*(6(Nmncwu#T(Rid$SEJdTmWPugUMfpS*BI-lz&`e zjur7X`CycO660M6!*NJ3P;F|vb8kG?J`etqG-&ZU%h{bQDggp~T9#*KW>Ol_M!zJU z?Xh)p)UWRJ?2xMzx*4f8CmZw({!MKxSARpBQ(Ggf+>^?aB@~dURZ?~7;1G7usx6Kv zQM*9B+?^gD*Icb*#TOq%EI`N~FIgmOO=>&O-{g)<9T0*U1G0PiZt){?a|`^{4Y&!- zohWrV9=LrzU$B^Spb+e}c|Wmw9Wii@dE)wfkw)HxGmfw3yDa^a_TV-_?ApWMR$eYf zl4;KhB;u3&cRd2k6dKF&80CbMz7z2 zyLPYK6fllJsPN6->*H777pAel*|KamdaP|Q=tfzle;FQvu{TO zXCF+GRz*uYa198&WA8QzNnLwUqrBxI1$1LsDwZo&CA5qw3!9&WS@7;QlTp&L`HlPh z@RF|*Y75ob`$h8$gu|kB@%D57ZtcoFx&KNK->vf=kd~)HJaP=38RBigkVtgd6V$nf zK|h-k{Z{wRA(h}y_To*^b0VEtElJUOSzXNixrM{ROAs@`0AY%PX-UfkWJ&Zk7qW1a zHjnJOy-E$fiidY82ws)SGxm;p!=r$4tfTtmxuvDQU$LhkN>%K=cXIKRwm3hoAExS* zM3^oSEZy{RQ@br0gm`tG<6FRRTt zwr#c}GeX>EEY{IX1|>++5R$!D;Gnlta)6sA}2$o0S(KK^`QK4yf-j1DAG zP8XXUh%V%~zw?~SL%ex0>)I!r-4Hp)yu3VpwGRB&<7}Yh_|VW~U*O3SA6q%_w!-f~ ztnunErc3&e&>_|m)D~K`*8R)?^|Y_s^>zr!6mhv*ZXZ%{qQZ4MS~O7BM|ZjvBU}yw zWP3@q_*{t}RZV$1>$`8uJZep{%zNqjtCXzrLkY-kJNjnv zo{?%66SL#-L^P7Y`~LJntyj99lW*X2kOU&&t@21YCCbGIfnMz{o3@emO2bB03|<|L zd9lYm7IMqsK(ofKX-Bm`n%7g?Zl1Q>?Z6zP+u^gKA^GoJ%^!O0O-&&tzL_`!MvmDUoU?= zA(!3Woz8X|ogEYAH&m%=_rL`jdpk ziQ#L|&7dibYw{IUea?WticJss3kRA-WGf z*8eR^og5b?bjZDh60Q7k`T5-!ubi6GCA^M3v93GH+Rl=b3bZi>#9Me_9 zk8AUwe}63U-#Fp20ZUyEhAvQxWW5m$dHPeEn~~^B_^CY@pv}9GfLAYgw517|CnIv33f~}7N;A|CfZ?`EElL*K>YY^@Kd_KFG*`9Rw^z~ zEB^S!gi&xrjP|v57siUL%`xOJbhp+F96BzIrCImI&1Ez`w$M)*2L&N&yA05K!|wWB z*g-a+Do{Px$1!n@ZR#WqQ{(ggMrbziI0ak86J6$s#z@=}iATittPHpbd3 z=6~|El&Q3D=(!vwwYPtVskZcqrObW_)}fZs_%{G}GeepJJ~Y3U4@v1vS*v(->h;)7ID8!Q#Q%Q#=`3%!Cp>zrcLZ- ze5z67mHg0(EuA)!YzZZTl%YERl(Dz)X)g*240E-M((8 z{7{K+ZkO_VQ-$q4EZSef2{jO6TX^coYmy@vV)E^^ZwVPB(-cOqcs4!9blKJuxm6+* zu*u!ukyGTh$~`R@MuLn{0kMAf9)N!J_ciGFc8$--O?*`&6zlY+%V zFpqZPHli#k&27bC6f!e@!E0Ib7p4WH|U1c*>$rgo@u0w_E-ZhCP zWqfDnnTk(G)!EWXy*=%F6J4p{ln#F8_K-%P=bv?%ba4UlVE@p{4ic89XNX#)ZG*wp zX?U2T2D72LBNO(te#C;8x6O0fx{8)si#wd$R>p%BJ$sC=@Jj7@<@uk4v2vk#KB(NS zL~t)#61q)?2}^GFym*(#C9}^8#n3I;blV~?L7djwHEHDH4zQhlbj*3m0$m7g7t3DW zFb1v5T9K(*MIzw19+&-Mezbc*9+8x7(HTqbV*xJVJdE*)ci+h$Z^^v^@b*`~<#TWI z&$lB4J&x6uhwZ;_TT(vXK0i_2KUPWEnmC)1{iR{JWj9($4F!wi%tFvez(GULPhC&J zWX#NW@T;|pE=0+iTS-fNW4XL)kZ^20sSpk7S4FJS3JTEwxDiVLyPaZ}1M`@h7vRA;U`<{y|*>hyx z(zDdLK+FunI@|=#2I28)xJ7)aQNjZBIKOm(}=daW*%dCe@i{_Fr9LDdJ z8f7wUsCDgFb#tB;>Bl`ElAaVZYj%lBec8+0Rk1%^9uQ~?N$po{AJ5uaOiC^yD(&gD zwlVHtnE4Y_ zXXMuu@Ynr5YdizE#p^=!V0>3bh}Q7F#N_<3s>=204s570#JGE_P+2lpYbrG1u3qdg z{{rjkcXDR+kV;akBx*!VreJ|i2 zlgy7Y+Oi+yq&=*THfFluFX4iYk%l9p`&kM2n>9GYh19m|mSN=!OH%3#Gx?S+idQtx zfvwB?sgp$E0nTbEg@?Z;s9@8bE)S*Ni?SiaP}MvPNa|y}qOlDpdwwGL;GwmhEh48{ zpfMAJ$a?5^imbG;_V-Qm0O0`_KO2}bfjQRl{Zx=|hOq?h%~bUCG8&1uSIrWqu$+3M z<^3k&hi=Ei{ym@RzL6M|M?3>FPPI0POKkp9uDZa|)&V>mD4qVXJVfKppFU=n&hq6T zQEYk5Z`o4W(QO!VuDB*-d4wxd;k?OhAD*q&Hok4d^E<&O3J4GiI!5D2+m%K+mD5SR zjZ+Z49F`z1`G3|_(zA_@q7Pvmtq%vw(n>bx0^ZYUa%vwg`$uAJ!o>BlaC z81lTM#vGenNVbNV>ceEGV-f$RU{Et97*YSqpGVX*iFe*vXjQ}Z+MRzBQDU;Gy99RV z73+j=$yxAXor`iOxNCI&u^j9;mYNNXvrZJa+q}9m?7S8Nh75B$oK!ER3=o<4kB0py z^M3%<qm`cyuIs~bYUNa>*wspjk`u2U3`DxME_W>YYD zljNv*Wt*)pKFgsmMs)t{&U{YokF;Q;NdQyew{0gR1fp;Cu=3E1)|d5DG-dSDZ}jh6 zs*I_qz0a`1O1tJ68DTR`^r1K!8RIa=$x)5VaT9^9rml7GIB(7ohTs#?b~>FnoGJND zW1#z9qrZ!vU&}Q~+0*lO7aoa8+`GGYzD^~>w z?69W$gWEDFFF=a}5mlMcvhkfzwxBcdKGAt{jOwZJ<@c|=7;ou>zS?Wbu}Cg`3Ut~y zkvRRhWR^5NE!-&W_+k~dxasSeblt?6?#|GwZf-Ci9^R8sWUIdz4{pt9SLssn;K20n zT`VtMZ~>p}UOE)imVkTv9s+*rI3l050f+-s~?GPSAcPPU1-M#9mz$p|= zcGi62P)Bd#lI#Sg7-{pWL_ic;8XM|$m2_VSIxHCIk2!8Tt)X=pTEoW(HK?LO9pq|o ztc%~eHEO-_CK1Kv0F=BrQu z8c%#K0v_un3&H%FI2%&#bLiL6N~{%t@2i~EjgD4Yop?-cXPtYe%ljWMvm|2Os@_M7 z2rNa4Qb`}xHm3&=#J%Y5z{RHYLPOl)Fc&y-q-kPk1hU4`pWc=RML_r4ce-i&Ol@*j z4-_<#Z{&R@;=X!YHxb5x{wSg_upg0_PfxRwjqbOaobS63HL9?zs71@iC={P-h@ck- zYC1_z)$E!D{K;&dkjLfDzvv)Wfo1xeetUW8-X`8iPA820BH_SRoZ0|Kt(m{(fJ1th zZd2azxa0uZUrQOxXb6st+}Rk@+BE2%l#OwH&Z zP89(_=mv2_k?8}l~3yfXHl+TxvQ1+BZQEMfd^uTx7lH|;`cSZW~68OL9JfC{w zXh5YQu(_Xk$Z)F3S#iWYDXHQAs$trE`!P(b1hoA{nt11=%xsM{L*Pv*Woop-aJ%d$ zbEy`XzrG_1)bxCa0yjB#;%iX+m35EJK9!ZbT~2BdR+3ZyH2)kmO(bveb{{7h2Wkho zWpsN+mJ!&nD!3o$8Au@<6V!rxsjzAP1P&h-T<*{EVb%S zORxdH5F2Q;8z7UI>eO7j>xZeORN;`hSjhXAIM{H`yZ>B&zmhl+Fn*+*kXkw!HF+`0 zh^%GnpIw@*Rk^1K;l+ws-Px{RHbZCcn&Qg! zxGm-N{xb*w?Q6gQnNV-|e`bS}K}J-N_cCoACK(u+b?6%y_%uuGo1q(cIvzx9rz-Lv zHraI+42z%&kUIA>OS7GNS->;cA5ElZU@+ZSbPl`hM5p)u@2nTZ&FybuV!|jUmK_y^ zhI6ZykeVv|9R&6F;=deM6g3to$je)h1h z(5|s?yUU9%m+UfuUZ3S6a3`1Y4nlu>`YEK*mV8TL6>oQ~B{ccxLvTb)dDYb8pPvc^ ziIL&BZ>oIVuF?cFV0K&-t_J8GH`?N-k`+V+$Kb{nY@F+ANy>aVF`-^tSJxoJ`nRmC zfxtmiOA8eZEx)4!0YzV}(T1e1k#8B#vo=bU!7!s7-(NO+nR`7#7x<;(^IP2c;)uio z+`DXk*Ax(~IHT?bZiN`)?f%wtPc^)A0v-hIdBozukMKmZD0favv2(_iH8#95c?3?U zR;E?qsa1DNRlyw^Z~f-7iTYZ2GQ8(9UXnK8x{%;CvWTgJ(JtxBvYe~I@aZ(hF^ z?@Qt5T-E70dj>y8c0GPd8MvUfJH_W@s_z1j7c%X<_64JToix3>FPRW+8MJ`D4R*Ru zcHR?Gn$wW%6b0AmY3w(~(zi^6*+=hm;#N{be;y~~iG?+LK9K7=K4?g6^Z{(PZ=K@t z$caY7q6NH{t*zwlw!`z~{ow(T949>u6^yr65bXm7rdCp=j=*q6Q-Mg9eb{n-UM$EY ztI%|@uEiSBKZ%4S_T8#B{P#L5gNmn>0i^%tU6Lz|qyke0U+ipdmB%$kH~HJQc*ASO zjM}aNiwF?L_|Dc#32TuZCK~6~f!l5vuIPLE=MU3pvP8gxHr~U(f3NsKDh8m_C&XhW zI6g9xT0j8yRc_Vf;~m?s@TJhrQFOZ@Pn~-1fYZ}e(AHF{8~k#bN0id5W=f;4C_48% zJ_tt7$0J1-cq=sdA)HI(;o`9gPkTICGwpB7<&lAV#fMEqbP)a1BJaQjn)vhMbI&}G z{i;N*i>|ovS3-1T47wA5qHN<*?z!jShuHx5@5Q5V+-oBqfhXrE@Y`vOuMDm`I=mrn z6j_L_ASGX;T$K4OK`*vo*g1&wx;Q(a3hO6a!p_0L` z+TXuF8~l8c`a03C3QL~RkvaU~$^@&q$h0?{#@pZ^L={n~on29bQ%e7n9h^f>>oz~M zj4F4d{}4l{gk6J+>G?=)^Esl*HVyqtrMBB?e4sm7c;v%28i!2Q@&3@R$@7)<_B{Xe zl-S_rA+}4dzl?94A%}kG(Mij)V6x=dA3^YZjr+CO?#!9wgz=Hd8BNbKX>2dx&fi}u zOJSkM&A)1;u(ToW2bqhgcMiB;t!ko`OJCn*$IW|sYSmUux4<26RuTj~Oxmc}R zv2={_BK3%f2rsQwtOLuE?E17cLbenj_+`=GOYvn=_{bt6XVwmJ3-I3BHm^hs)~4X zWWZK4&gc^5$ra{}Gn0^; ztEQ+JWn=Z;CNa1^gl28rAb4$ww~n{1=TtayVZPY5v|>ZXK*BzbEoNG5rs4(IH%@)qvFdBA+CCQpg(oE6&s=H~t0F(HoL@c&w8Z1DJu~ZfNtWX+3n6rmm$QG{Lb|w?hO+W|2O!VpHN{$Ghaw3dZmnSGNZrD<5DrtU8X{ z%|h1i-=W+?reX%K9k)Me&SQnU32q zpDc8#)ww>Vr zs_T7I_bU=zt=sO|?sqk>@SyXl;w>)oJ)N&85&=qUr^fwzLB#x_VJVTF^Yp>_7?Em= zaaRNo@|#pUKYekN#hB?1Qr+uii?&cg4sq>(8K4HxF zslJ(|=(xrZPZST{`;JzS$KY5Z=kub4d}2K%9?3;E$xbA{&u*vfy-6$ zk0!@P3FUMq=jMvxX@c%9E(xj9S8;+JmchJ)^q2GmHByFn?F4B&RLq@h%(~SdDgD*nC$^n{0E$sRq0~&xgA?*0`w{S=4cooIRLUf_xU~ z6LicSOH@?ECBT!-a7zJp@qD_vL^6b(^&ah{eyku=;6s|^PC-I&*g0pQp|38!CEcmI zfgu)i_NHbcZ)Zf6y9%>maTHZb|7h%&I+S+&jiKr(Qpjw9rXuKlo_+4!$f%-oL*! zHH*#3$2;~*5_HvY<@wGVuG#=>-;ag@Di$xv?p9^~93Iw;KpS5J>Z743T&664b$x%* z*LvK{o3;pJo6TtCl5y`Sh0i9~CYmHMD?U$CpfI)EH#q%qk??%eQg1Gj`z7+Wq z`2*4C2+iJelb@Eg&$=&i6w^H$;^~S4{Ys{VPlQNU=35Ej^wi9r`kq=;Io)<**6KhXC>Z$cFGn-#=JlRQMH)BRkc#*yn;_VRy3)LDm%d$V~7k~ zTC5#&c6wvDq!!+b2_=)zYyE}W_i?npHVS@Byy4t9Mh!-O{TXnjx2m_nF1~XdtI{tPqx^GwkHy}`MsF`@`OQ(K zik$#!Ha{1jbbjsozM~No8m$7V+gY|*bnmaN`?VYMc+Ca>7#jJdz=SG2k{&wqWs65g z8ZecDEA_iT;(VltbcaTPC=2WhgnfSwxY^er{2(Aa{&_R`_z=dRySjI>(XHF|yOVPr zdS#>Gd)qU1vizCvTw`0c^o^vT*3TPzK7A5>kRU;iifA!!)UYeV90)h6CMFwtd-^=E zi1k90oCwn!^}_SK;vXq2CkC|ht#CdQZoU^OUNNydKDtwD+AGGF=xgRRSQr+h9gW2A z*4P`NS%rn30pK_0alj57s+&E>sjrA-fulCki$v7ZQmZd6oTZjZ()XXIddaXtZ^h`9 z+7UMAyTUn`-^__>9}|34t)!l04!C=x1W?|K2aWO1>9cG0%v#e z^?fJV&9lnUHxDS@W*?dx(iN(>WAUeP2ld0TtY&PBrX(7)PU~F5yaVUe`iPM|GP5$B zQ0^XtnylwL;?L9k8s#5fl=&=GviPP}A@B{>S_oh{%w=6BVZX8%=h$0bklZMydE$@- zih72Y5PC8rWTv=&Ogb?PnseLqHzHx6n$r@%I-%ZDCF-7QZlkGDDtj{flo0tHA=Tdb#- z?yl2a@Tbd#?DTESG&R;)cJnSmoj!jq_I0$Rs@|gW-oGF<7q}-_$vgqBOX9P#%3Ok< zy7&;Te(Y-sYD?dK`?QV)Gisw%SmEn2k(cw)iQ(7P^0ru*N=v3x*_N)@YA;lN}FnYwk zc-EJRvduNQRY!MWOQO^EG zlAZIdeUsw~@^DeGh6+sc*G~CPwcyi9P^{@kf9y%pXX1Vnv*s_O8h!X(5Xw+jIcA6H zFwL{9HQ!BOJkNv${K>C+LaOplK>f&wI5n-V!|n(AEa-aa6_tRO9O+%@Wz_Et^xv!9 z_CSfjyqxeV`g-m+8l!Aol$@)2FE0dzR}8dsJnw%`zm4jp5w#OJ;S(#wN7*$%E0AGg z@7G&o1rXZhkoo$H{Hf=1Kc_Q81CYQ*yBAODGh2+MekD{FoDHa!{M)0(TVq#TlYsdF z+eYz%hY7fZ!13FE_AoIt!Y>e7e!Bng?3gyBR6E!7d*Rk5E_aId-= zUUHf|n4Y)-L$}u}=#m~j`d2VhcUd_^>DlZ}9lg1U9m3wu&L-!2AU`m6OIX~K8s&O< z3sX5lNA?6xiYh8ukbfb=Dq=cRz1=v*D!2j}nUrl#|2c2<=GgLgnbl`Ts*FwvkF}kY z9FLoRkX2j~1x90vIHS_AEd}M>f$h;bTu|a1H`2}EAp%%|){4ETFrHdc)xr*cDzXlZ z`}3W$v5BfGo{Gx8cxuKYf984j7N^c3tZBVediXF_EH^ZB@Ui1cVecfX{Bgza=2THB z69pSNDx(PxmMVer>q`RdQ|im(kM{u_8FU2$rJFoJ_xks7wK67whdD4M>P<)8>2j5o zt2oI|BsumKd2UVh>T0^;slm|i{Qemf^m1OZ1Q1@Vx7(sK&*)a)59+Y7p}-@zC=`=q z@&icFj6|ZjMbBtT;=~m`HBK|qTCcPf8OjVA`%|l_pGK~xs!QMS4CvMtWUu(W` za*mjuo^CGG>3l-7o6U{_8<@bNouujqhL(MG^3LWZMYQ&spm2;&D2!Vnpq7wT4Lmkz zDZXzs^OGTtoEiOTnR~iMD{O-r+XMk$mSOD2E|*)^CF~Eex66L+MR9Zuhi113oXQ3W z8LP2uZc(0};QC7FUZNrK&}uUsvT`PF^g9@xL&=U=DFUF~!*IoWSxF}U(2-C0EB-g6Zw2

YR^Ox=}f#`1Y7rmoF#Tg-19zu#o*l#8KE}vjIkGqrU_w-gMUs# zuNArp661=C1ODmHw&-%U4vK)=t)t;V3+N`Xiwl&8r^kYcwP}w)^_0k>R?wiR}pD1=E?jRtG%}TiejBFWcSuqyZHX*K5mD33Z zWrNLs)-^z8%&}CkkzI*c@ks|Jj%i$?LpXKNnA7O9u4`O@NEi`-hG$PieK3E+0G+C~ zth~5MV=&@Q5rIV*}w@pD{Ydr@E&Nt!b9M! zvMLElnMZKal6e;n+biMz^K3O5-ZJxopp9Iy=fi7$V~b#)ZrdP{+WmPN;EI*HixDAQ zN~h(*;U7dLLz2wcQvq1gXbk0OSU`3TMztCfb-L}q?@jA*V*i^BT%XmkU>nifI$Dre z4sl4T%4hvOLNEp&mtCJ$Y^#99wNffhMNlpJ5qYnsH|#x(2Zu$U+!`7!9Io+mHF}bwlX_2X~5_W+Vwsyn;%5s&z@U-42`l* zcqQ8;JJUaql!*L9eqE5S=a;Bax9b|$UihQ7Ds?G#rA6MG`dH*w2gi9i*`Q%q+>B`C zbng8tCGy9V>e_w(%oLbDlBaYkELc^uWqgd%kXL$f@UTx;hKIUsZDuo}r%u%{fbDid zxL5Hq{k?JQhqsRYpCcBntvJ|tV*dPa+EoVCnkx?`69ZmGovB@>ZZo)Zay`$&TO<9M z{zN|XCQc)P3;eu5uRgcJd$u+GrTm>PyXi;MBr`3)EKZ6)v`R0XqM9`8M_r;ziR$0Wm9Pk^*ehIYt(Nm<$YMx8OZ;+5`8cbhlKnHH7otkQDhac`V%%2TQ|d(Tt)s3P|P zoU~xTwMjy89?9s@+PzO>)WcSfTZt;=X;@?qcXUOFFWaKmGoAO{>O=(DqLdis)J_Nyi{TmJ!C0Z*h}$-{N6U;Gp-W4`#2;#xj4~R!e)2FN!B*bC$$5#J%90 zq!bux6drtn`ZRb`BV5JMxKo>?9JJqt@vtpx&Dl-^14eu%2WR^S)*P%d^|E=ogG~zq zZvBD1$8|}`i03^q&ek2th*0DU>jI7J*2BEsj%0Y~qdtm0p8+#UD_-8pv$Hc%sc#{k znp#YS)j0f^tkf>^jEvZ|@Rp=YkdIFo)Lb0BGW4DYE2LR{DJ1^;k3_QJ`ae}Q!)Oz^ ziYYx4lh0BqnefxUmxjp#;6Vxs_aFyL6Ve6-20JL!bbAM5-`D#!n)CQh2n2$Ii%V2O zqBXxnKGZkk$DXO z02;<{=2a0<(Vg90ecWrU;~m+vGN>(X`9F!Hvg;r9>1{<)hN`IG$H&JnpnE)UN|W7C zuL$G63I)gKT{TIaW&S2LgbyqQ&(%gfV(9~)_f7e z;Jos)sX)W8!^v*&dw%VT+v#)ZN)^xJs@!t{+*9LTP6lAktFeLoQi{9o{UvK2n&SMt zK4Y$Xd`axjjwSXS?`c#5E_6a_YQ_K@pQDlys9B&!(23N~Dp>x?D*-$F^ey>qp?$zU z2%5Fg=TKFRKtXc;k`EqWR=EF%hPUJkvt8mQa<*0s!=fnat@1RZvm{Pf{%Su&b;RA? zJ3MOhrwuEHACSkA_C7=1#e(j5a7@U3-uX?2hFJjf7?vx?WpoWAa~D%duVd2GYQp5f z0g_6bOn_D$`pyX$kf|+yCyq<#lAhI0af;;va88Lk+TpV^O(*|a&?J_sEW<6?F&Ob7 z>(RI{&%FQquz~HO`24lP>fy@7uTu*Eo&F!*?Se%PLBhEDq+cKp z)s$B61I`JF7g`P{{GpVENE*)M-?WtbA-P1e2K3_)3dkxA8GdwZp(n@ox^622u2+)`kc&E-Ak5MvX~;-7Qk|~qYVOhv2KnXR@~%hWr<>@is=pPyls`hG`l^n(zG_*+^PV6xHDzJ%xdR@f z6V#3*baeF}MdwcZpGgi`y)PKexSJAxs>CrjrnDqumj&@q&jZ#kIci6YGrBwb{~; z+@rmE2h5X<&9OA0F=rb?lJi?en%K8{nwQmod3QCV?e}mjh>v&5ui>nY6|J!-GI>`Q zU*3eLX(qWcSQn0$=t|Xd4Q2rH_~2_?1*MnE2S+0pYjuuf-eg7?3odh0ul&O2#1s}X zQ&Lhk9L6H5&aSa3#LkZ0+}wPn zr8*n;@ylZQU}yDxTp&|g-lG41A!Gx!gLf<^-rNotes

  • The watershed layer included in RSA database was not used due to issue found. For example, layer “Watersheds_Notes_RSA2021” has missing area in Rock Creek. Instead used: DC -Subwatersheds  Open Data Layer.

  • +Subwatersheds Open Data Layer.


    @@ -96,9 +96,10 @@

    Steps for processi
  • Intersect DC -Subwatersheds  Open Data Layer with StreamPoints layer

  • +Subwatersheds Open Data Layer with StreamPoints layer

  • -

    Run the ‘Add Attribute Index’ tool:

    +

    Run the ‘Add Attribute Index’ tool on the Intersected layer so +that ID field is indexed:

    • Field to Index: GlobalID_copy

    • Index Name: GlobalID_index

    • @@ -107,19 +108,24 @@

      Steps for processi
    • Join StreamPoints attributes table to StreamPoints feature layer. See settings in screenshot below:

      +
        +
      • If ’One to Many” Join Operation not an option, make sure attribute +table and feature layer are in the same geodatabase
      • +

    • -

      Export joined StreamPoints table to Excel

      +

      Export joined StreamPoints by using ‘Table to Excel’ +geoprocessing tool

        -
      • Save file to R package data folder
      • -
      • Update file name input variable in workflow.R file
      • +
      • Save file to R package data folder

      • +
      • Update file name input variable in workflow.R file

    • In the input variables section of workflow.R, update the field name prefix variables (if needed)

    • -

      Repeat steps 1, 3,4, 5 and 6 using StreamReaches feature layer +

      Repeat steps 1, 3, 4, 5 and 6 using StreamReaches feature layer and attribute table

      • Skip Step 2 since subshed information is already included as a field diff --git a/pkgdown.yml b/pkgdown.yml index 61091ed..81b0e5b 100644 --- a/pkgdown.yml +++ b/pkgdown.yml @@ -4,7 +4,7 @@ pkgdown_sha: ~ articles: faqs: faqs.html pre-processing: pre-processing.html -last_built: 2025-02-05T20:30Z +last_built: 2025-02-06T23:18Z urls: reference: https://limnotech.github.io/SCI/reference article: https://limnotech.github.io/SCI/articles diff --git a/reference/assess_dumpsites.html b/reference/assess_dumpsites.html index 488af1c..8e74677 100644 --- a/reference/assess_dumpsites.html +++ b/reference/assess_dumpsites.html @@ -46,8 +46,6 @@

        Usage
        assess_dumpsites(
           df_point,
           df_reach,
        -  start_date,
        -  end_date,
           reach_prefix_from_table,
           reach_prefix_from_layer,
           point_prefix_from_table,
        @@ -67,16 +65,6 @@ 

        Argumentsstart_date -

        Beginning of timeframe being evaluated (ex. first day of -five-year assessment period).

        - - -
        end_date
        -

        End of timeframe being evaluated (ex. last day of -five-year assessment period).

        - -
        reach_prefix_from_table

        Field name prefix sourced from RSA reach attribute table. Should be consistent with the name of the stream reach @@ -110,24 +98,24 @@

        Value

        Examples

        assess_dumpsites(df_point_example,
                                    df_reach_example,
        -                           "2020-01-01",
        -                           "2024-12-31",
                                    "StreamReachAttributes",
                                    "StreamReaches",
                                    "StreamPointAttributes",
                                    "StreamPoints_Intersect" )
         #> $summary
        -#> # A tibble: 7 × 8
        +#> # A tibble: 9 × 8
         #> # Rowwise: 
         #>   sci_subshed            dumpsite_count applied_influence subshed_length_meters
         #>   <chr>                           <int>             <dbl>                 <dbl>
        -#> 1 Broad Branch                        1                -1                 9090.
        -#> 2 Fort Davis Tributary                1                -1                 1198.
        -#> 3 Fort Dupont Tributary               1                -2                 8637.
        -#> 4 Fort Stanton Tributary              1                -2                 2214.
        -#> 5 Oxon Run                            1                -1                 9965.
        -#> 6 Rock Creek                          2                -2                32011.
        -#> 7 Watts Branch                        1                -2                 6013.
        +#> 1 Broad Branch                        2                -1                 9840.
        +#> 2 Fort Davis Tributary                3                -2                 1469.
        +#> 3 Fort Dupont Tributary               6                -2                 9175.
        +#> 4 Fort Stanton Tributary              1                -2                 2507.
        +#> 5 Normanstone Creek                   2                -1                 3527.
        +#> 6 Oxon Run                            1                -1                10144.
        +#> 7 Pinehurst Branch                    1                -1                 6001.
        +#> 8 Rock Creek                          4                -2                38328.
        +#> 9 Watts Branch                        2                -2                 6035.
         #> # ℹ 4 more variables: subshed_length_miles <dbl>, sites_per_mile <dbl>,
         #> #   score <dbl>, score_weighted <dbl>
         #> 
        @@ -136,16 +124,16 @@ 

        Examples#> # Rowwise: #> sci_subshed Dumpsites #> <chr> <dbl> -#> 1 Broad Branch 8 -#> 2 Fort Davis Tributary 2 -#> 3 Fort Dupont Tributary 7 +#> 1 Broad Branch 7 +#> 2 Fort Davis Tributary 1 +#> 3 Fort Dupont Tributary 2 #> 4 Fort Stanton Tributary 4 -#> 5 Oxon Run 8 -#> 6 Rock Creek 7 -#> 7 Watts Branch 6 -#> 8 Battery Kemble Creek 10 -#> 9 Dalecarlia Tributary 10 -#> 10 Dumbarton Oaks 10 +#> 5 Normanstone Creek 4 +#> 6 Oxon Run 8 +#> 7 Pinehurst Branch 7 +#> 8 Rock Creek 7 +#> 9 Watts Branch 5 +#> 10 Battery Kemble Creek 10 #> # ℹ 15 more rows #>

        diff --git a/reference/assess_trash.html b/reference/assess_trash.html index 7c4da6b..3f35b3c 100644 --- a/reference/assess_trash.html +++ b/reference/assess_trash.html @@ -43,13 +43,7 @@

        Calculate trash scores.

        Usage

        -
        assess_trash(
        -  df_reach,
        -  start_date,
        -  end_date,
        -  reach_prefix_from_table,
        -  reach_prefix_from_layer
        -)
        +
        assess_trash(df_reach, reach_prefix_from_table, reach_prefix_from_layer)
        @@ -60,16 +54,6 @@

        Argumentsstart_date -

        Beginning of timeframe being evaluated (ex. first day of -five-year assessment period).

        - - -
        end_date
        -

        End of timeframe being evaluated (ex. last day of -five-year assessment period).

        - -
        reach_prefix_from_table

        Field name prefix sourced from RSA reach attribute table. Should be consistent with the name of the stream reach @@ -89,41 +73,12 @@

        Value

        Examples

        -
        assess_trash(df_reach_example, "2020-01-01", "2024-12-31",
        -    "StreamReachAttributes", "StreamReaches")
        -#> $summary
        -#> # A tibble: 733 × 7
        -#>    location_name        id    date                trash_extent score sci_subshed
        -#>    <chr>                <chr> <dttm>              <chr>        <dbl> <chr>      
        -#>  1 Hickey Run           {000… 2024-10-04 15:09:00 extensive        1 Hickey Run 
        -#>  2 Piney Branch         {001… 2024-06-11 16:37:00 absent          10 Piney Bran…
        -#>  3 NA                   {003… 2022-06-08 17:06:00 absent          10 NA         
        -#>  4 Pinehurst Branch     {004… 2024-05-31 14:58:00 absent          10 Pinehurst …
        -#>  5 Rock Creek           {008… 2021-07-08 13:42:00 absent          10 Rock Creek 
        -#>  6 NA                   {015… 2021-07-29 16:44:00 minor            7 NA         
        -#>  7 Broad Branch         {024… 2023-08-25 19:20:00 minor            7 Broad Bran…
        -#>  8 Broad Branch         {027… 2023-08-25 19:43:00 minor            7 Broad Bran…
        -#>  9 Fort Davis Tributary {02A… 2024-09-17 16:35:00 absent          10 Fort Davis…
        -#> 10 Klingle Valley Run   {031… 2024-06-18 14:50:00 minor            7 Klingle Va…
        -#> # ℹ 723 more rows
        -#> # ℹ 1 more variable: ecoregion <chr>
        -#> 
        -#> $score
        -#> # A tibble: 22 × 2
        -#>    sci_subshed            Trash
        -#>    <chr>                  <dbl>
        -#>  1 Battery Kemble Creek    8.43
        -#>  2 Broad Branch            8.88
        -#>  3 Dumbarton Oaks         10   
        -#>  4 Fenwick Branch          7.43
        -#>  5 Fort Davis Tributary    8.33
        -#>  6 Fort Dupont Tributary   7.64
        -#>  7 Fort Stanton Tributary  6.47
        -#>  8 Foundry Branch          8.83
        -#>  9 Hickey Run              8.04
        -#> 10 Klingle Valley Run      6.25
        -#> # ℹ 12 more rows
        -#> 
        +    
        assess_trash(df_reach_example,
        +                       "StreamReachAttributes",
        +                       "StreamReaches")
        +#> Error in dplyr::filter(., date >= start_date & date <= end_date):  In argument: `date >= start_date & date <= end_date`.
        +#> Caused by error:
        +#> ! object 'start_date' not found
         

        diff --git a/search.json b/search.json index 1714e9d..ca2bf64 100644 --- a/search.json +++ b/search.json @@ -1 +1 @@ -[{"path":"https://limnotech.github.io/SCI/LICENSE.html","id":null,"dir":"","previous_headings":"","what":"MIT License","title":"MIT License","text":"Copyright (c) 2024 SCI authors Permission hereby granted, free charge, person obtaining copy software associated documentation files (“Software”), deal Software without restriction, including without limitation rights use, copy, modify, merge, publish, distribute, sublicense, /sell copies Software, permit persons Software furnished , subject following conditions: copyright notice permission notice shall included copies substantial portions Software. SOFTWARE PROVIDED “”, WITHOUT WARRANTY KIND, EXPRESS IMPLIED, INCLUDING LIMITED WARRANTIES MERCHANTABILITY, FITNESS PARTICULAR PURPOSE NONINFRINGEMENT. EVENT SHALL AUTHORS COPYRIGHT HOLDERS LIABLE CLAIM, DAMAGES LIABILITY, WHETHER ACTION CONTRACT, TORT OTHERWISE, ARISING , CONNECTION SOFTWARE USE DEALINGS SOFTWARE.","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"effective-impervious-area-eia","dir":"Articles","previous_headings":"","what":"Effective Impervious Area (EIA)","title":"FAQs","text":"Rock Creek mainstem includes MS4 CSS. Rock Creek watershed also includes additional area beyond Rock Creek SCI subwatersheds. area Rock Creek used SCI analysis? entire Rock Creek watershed used SCI analysis. See “eia_rock_creek” tab lookup spreadsheet full list subsheds included.","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"water-quality-index","dir":"Articles","previous_headings":"","what":"Water Quality Index","title":"FAQs","text":"non-detects handled dataset? reporting limit used non-detect values. example, results reported <0.010, got incorporated analysis 0.010. water quality criteria used various parameters? See Section 1 (Indicators Metrics) Stream Condition Index Methodology (2022) Among non-nutrient WQ indicators, maximum criteria, minimum criteria, maximum minimum criteria. assess_wq() function handles three cases. nutrient WQ indicatiors (TN TP), assess_wq_nutrients() requires criteria input Piedmont Coastal Plain Level III Ecoregions","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"aquatic-life-index","dir":"Articles","previous_headings":"","what":"Aquatic Life Index","title":"FAQs","text":"aquatic life metrics calculated?  four metrics get calculated outside R package get incorporated lookup tables data/lookup_tables.xlsx. Connection Mainstem unlikely change year year. Subsheds assigned score either 10 (connected mainstem) 1 (connected mainstem). option calculating fish IBI, habitat PHI, macroinvertebrate HBI MBSStools R package.  IBI, PHI, HBI scores get translated SCI score added relevant lookup table. IBI, PHI, HBI scores translate SCI scores? See conversion tables :","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"stations","dir":"Articles","previous_headings":"","what":"Stations","title":"FAQs","text":"stations included tributary? station locations within given tributary incorporated analysis included location_id lookup table. Rock Creek, include stations specifically Upper Rock Creek Lower Rock Creek. separate scores calculated upper lower segments Rock Creek Watts Branch? , single score calculated indicator Rock Creek Watts Branch","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"timeframe","dir":"Articles","previous_headings":"","what":"Timeframe","title":"FAQs","text":"timeframe used analysis? goal use five recent calendar years data (e.g. 01/01/2020 – 12/31/2024). exceptions may made extend timeframe beyond five years found inadequate data.","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"troubleshooting","dir":"Articles","previous_headings":"","what":"Troubleshooting","title":"FAQs","text":"try run R code, get warning. ? warnings expected issue running R code. able proceed steps normal. still questions concerns, contact package maintainer (see ‘Contacts’ subsection FAQs) try run R code, get error. ? Possible things try include: Read error message see can pinpoint error occurred. ·sure input variables specified correct. Confirm file names pathways correct. Check whether datasets format changes. Changes field names, changes waterbody names, new sampling IDs, etc. can disrupt code. ’re unable resolve error, still questions concerns, contact package maintainer (see ‘Contacts’ subsection FAQs)","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"contacts","dir":"Articles","previous_headings":"","what":"Contacts","title":"FAQs","text":"contact questions R package? R package currently maintained Shanna Rucker Blount (LimnoTech). can contact srucker@limno.com","code":""},{"path":[]},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"notes","dir":"Articles","previous_headings":"Trash and Dumpsites","what":"Notes","title":"Pre-Processing","text":"R functions assess_trash() assess_dumpsites() require dataframe inputs processed RSA reach data RSA point data (dumpsites ) watershed layer included RSA database used due issue found. example, layer “Watersheds_Notes_RSA2021” missing area Rock Creek. Instead used: DC Subwatersheds  Open Data Layer.","code":""},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"steps-for-processing-rsa-stream-reach-and-point-data","dir":"Articles","previous_headings":"Trash and Dumpsites","what":"Steps for processing RSA stream reach and point data","title":"Pre-Processing","text":"StreamPoints feature layer, copy GlobalID new field Set data type new GUID New field name suggestion: GlobalID_copy step needed GlobalID information included intersected table Intersect DC Subwatersheds  Open Data Layer StreamPoints layer Run ‘Add Attribute Index’ tool: Field Index: GlobalID_copy Index Name: GlobalID_index Join StreamPoints attributes table StreamPoints feature layer. See settings screenshot : Export joined StreamPoints table Excel Save file R package data folder Update file name input variable workflow.R file input variables section workflow.R, update field name prefix variables (needed) Repeat steps 1, 3,4, 5 6 using StreamReaches feature layer attribute table Skip Step 2 since subshed information already included field StreamReaches feature layer","code":""},{"path":[]},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"notes-1","dir":"Articles","previous_headings":"Effective Impervious Area (EIA)","what":"Notes","title":"Pre-Processing","text":"Subshed area calculations done external R package. Lookup tables must updated manually. final EIA percentage score calculations done R. BMP Data timeframe always begin January 1, 2000. two lookup tables need manually updated annual SCI run: eia_subsheds eia_rock_creek addition sci_subshed field, main columns lookup tables get used R functions prefix “analysis_” Analysis_subshed_m2 Analysis_impervious_m2 Analysis_impervious_treated_m2 columns particular years reference GIS analysis repeated numbers needed later Special subshed notes: Piney Branch sci_subshed MS4 CSS. eia_subsheds lookup table includes MS4 CSS subsheds two separate segments. Reservation 630 used sci_subshed, Melvin Hazen Valley Branch used lookup subwatershed data sources","code":""},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"steps-for-impervious-area-calculation","dir":"Articles","previous_headings":"Effective Impervious Area (EIA)","what":"Steps for impervious area calculation","title":"Pre-Processing","text":"perform following steps new impervious surface available Open Data DC Start latest impervious surface layer Open Data DC (e.g., 2023 Impervious Surface) Intersect impervious surface layer Open Data DC Subwatersheds layer Dissolve intersected layer summarize reduce number features (See screenshot Dissolve settings) Export Dissolve table excel. Update lookup tables R package spreadsheet. data/lookup_tables.xlsx, modify eia_subsheds eia_rock_creek spreadsheet tabs: Add new column called “20XX_impervious_m2” bring data Dissolve table Update “analysis_impervious_m2” column referencing new column created Save spreadsheet workflow.R script, uncomment code eia_subsheds eia_rock_creek. Run code update two .rda files. Re-comment code.","code":""},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"steps-for-impervious-treated-calculation","dir":"Articles","previous_headings":"Effective Impervious Area (EIA)","what":"Steps for impervious treated calculation","title":"Pre-Processing","text":"Export BMP data SGS using IPMT query - Start date: January 1, 2000 - End date: December 31, 20XX (use last year target 5-year period) Excel, remove BMP records meet Chesapeake Bay inspection criteria. BMPs installed 10 years ago, BMPs must passing inspection within last 10 years. Use following three fields identify remove BMP records older 10 years passing inspection: (1) Installation Date; (2) recent Inspection (3) Passed failed recent maintenance inspection? example, SCI timeframe 2020-2024, BMPs installed prior January 1, 2015 must passing inspection. ArcGIS Pro, bring table SGS export display XY data Intersect BMP point data Open Data DC Subwatersheds layer Dissolve using settings screenshot Export Dissolve table Excel Convert impervious area sum ft2 m2 Update lookup tables R package spreadsheet. data/lookup_tables.xlsx, modify eia_subsheds eia_rock_creek spreadsheet tabs: Add new column called “20XX_impervious_treated_m2” bring data Dissolve table Update “analysis_impervious_treated_m2” column referencing new column created Save spreadsheet workflow.R script, uncomment code eia_subsheds eia_rock_creek. Run code update two .rda files. Re-comment code.","code":""},{"path":"https://limnotech.github.io/SCI/authors.html","id":null,"dir":"","previous_headings":"","what":"Authors","title":"Authors and Citation","text":"Shanna Rucker. Author, maintainer.","code":""},{"path":"https://limnotech.github.io/SCI/authors.html","id":"citation","dir":"","previous_headings":"","what":"Citation","title":"Authors and Citation","text":"Rucker S (2025). SCI: Stream Condition Index. R package version 0.0.0.9000, https://limnotech.github.io/SCI/, https://github.com/LimnoTech/SCI.","code":"@Manual{, title = {SCI: Stream Condition Index}, author = {Shanna Rucker}, year = {2025}, note = {R package version 0.0.0.9000, https://limnotech.github.io/SCI/}, url = {https://github.com/LimnoTech/SCI}, }"},{"path":"https://limnotech.github.io/SCI/index.html","id":"sci","dir":"","previous_headings":"","what":"Stream Condition Index","title":"Stream Condition Index","text":"goal SCI calculate scores hydrology, water quality, human health, aquatic biology metrics assess streams within District Columbia.","code":""},{"path":"https://limnotech.github.io/SCI/index.html","id":"installation","dir":"","previous_headings":"","what":"Installation","title":"Stream Condition Index","text":"can install development version SCI GitHub :","code":"# install.packages(\"devtools\") devtools::install_github(\"LimnoTech/SCI\")"},{"path":"https://limnotech.github.io/SCI/index.html","id":"using-the-package","dir":"","previous_headings":"","what":"Using the Package","title":"Stream Condition Index","text":"Produce SCI scores just simple steps: Complete pre-processing steps, needed. RStudio, open workflow.R script. Update “Input Variables” section script confirm latest data files saved data folder. Update lookup_tables.xlsx run code “Lookup Tables” section script, needed. Note: step typically necessary. Click Source button run workflow.R script. View results all_scores.csv.","code":""},{"path":"https://limnotech.github.io/SCI/reference/SCI-package.html","id":null,"dir":"Reference","previous_headings":"","what":"SCI: Stream Condition Index — SCI-package","title":"SCI: Stream Condition Index — SCI-package","text":"Calculates scores hydrology, water quality, human health, aquatic biology metrics assess streams within District Columbia.","code":""},{"path":[]},{"path":"https://limnotech.github.io/SCI/reference/SCI-package.html","id":"author","dir":"Reference","previous_headings":"","what":"Author","title":"SCI: Stream Condition Index — SCI-package","text":"Maintainer: Shanna Rucker srucker@limno.com","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate dumpsite scores — assess_dumpsites","title":"Calculate dumpsite scores — assess_dumpsites","text":"Calculate dumpsite scores","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate dumpsite scores — assess_dumpsites","text":"","code":"assess_dumpsites( df_point, df_reach, start_date, end_date, reach_prefix_from_table, reach_prefix_from_layer, point_prefix_from_table, point_prefix_from_layer )"},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate dumpsite scores — assess_dumpsites","text":"df_point Point data Rapid Stream Assessment (RSA) dataset. df_reach Reach data Rapid Stream Assessment (RSA) dataset. start_date Beginning timeframe evaluated (ex. first day five-year assessment period). end_date End timeframe evaluated (ex. last day five-year assessment period). reach_prefix_from_table Field name prefix sourced RSA reach attribute table. consistent name stream reach attribute table RSA geodatabase. reach_prefix_from_layer Field name prefix sourced RSA stream reach layer. consistent name stream reach polyline layer RSA geodatabase. point_prefix_from_table Field name prefix sourced RSA point attribute table. consistent name point assessment attribute table RSA geodatabase. point_prefix_from_layer Field name prefix sourced RSA point layer. consistent name intersected stream point layer RSA geodatabase.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate dumpsite scores — assess_dumpsites","text":"List 2 dataframes: summary score.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate dumpsite scores — assess_dumpsites","text":"","code":"assess_dumpsites(df_point_example, df_reach_example, \"2020-01-01\", \"2024-12-31\", \"StreamReachAttributes\", \"StreamReaches\", \"StreamPointAttributes\", \"StreamPoints_Intersect\" ) #> $summary #> # A tibble: 7 × 8 #> # Rowwise: #> sci_subshed dumpsite_count applied_influence subshed_length_meters #> #> 1 Broad Branch 1 -1 9090. #> 2 Fort Davis Tributary 1 -1 1198. #> 3 Fort Dupont Tributary 1 -2 8637. #> 4 Fort Stanton Tributary 1 -2 2214. #> 5 Oxon Run 1 -1 9965. #> 6 Rock Creek 2 -2 32011. #> 7 Watts Branch 1 -2 6013. #> # ℹ 4 more variables: subshed_length_miles , sites_per_mile , #> # score , score_weighted #> #> $score #> # A tibble: 25 × 2 #> # Rowwise: #> sci_subshed Dumpsites #> #> 1 Broad Branch 8 #> 2 Fort Davis Tributary 2 #> 3 Fort Dupont Tributary 7 #> 4 Fort Stanton Tributary 4 #> 5 Oxon Run 8 #> 6 Rock Creek 7 #> 7 Watts Branch 6 #> 8 Battery Kemble Creek 10 #> 9 Dalecarlia Tributary 10 #> 10 Dumbarton Oaks 10 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate scores for effective impervious area (EIA) — assess_eia","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"Calculate scores effective impervious area (EIA)","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"","code":"assess_eia()"},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"List 2 dataframes: summary score.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"","code":"assess_eia() #> $summary #> # A tibble: 25 × 7 #> # Rowwise: #> sci_subshed subshed impervious impervious_treated effective_impervious #> #> 1 Battery Kemble Cr… 9.37e5 201840. 26487. 175353. #> 2 Broad Branch 4.65e6 1627836. 76874. 1550962. #> 3 Dalecarlia Tribut… 1.09e6 295100. 23139. 271960. #> 4 Dumbarton Oaks 5.51e5 153336. 11829. 141508. #> 5 Fenwick Branch 8.94e5 284547. 6249. 278298. #> 6 Fort Chaplin Trib… 1.18e6 423427. 30927. 392500. #> 7 Fort Davis Tribut… 9.47e5 328257. 12328. 315929. #> 8 Fort Dupont Tribu… 1.66e6 154824. 15534. 139289. #> 9 Fort Stanton Trib… 1.12e6 488751. 64799. 423952. #> 10 Foundry Branch 4.59e6 1725013. 198108. 1526905. #> # ℹ 15 more rows #> # ℹ 2 more variables: percent_effective_impervious , score #> #> $score #> # A tibble: 25 × 2 #> # Rowwise: #> sci_subshed `Effective Impervious Area` #> #> 1 Battery Kemble Creek 6 #> 2 Broad Branch 2 #> 3 Dalecarlia Tributary 4 #> 4 Dumbarton Oaks 4 #> 5 Fenwick Branch 3 #> 6 Fort Chaplin Tributary 2 #> 7 Fort Davis Tributary 2 #> 8 Fort Dupont Tributary 8 #> 9 Fort Stanton Tributary 1 #> 10 Foundry Branch 2 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate trash scores. — assess_trash","title":"Calculate trash scores. — assess_trash","text":"Calculate trash scores.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate trash scores. — assess_trash","text":"","code":"assess_trash( df_reach, start_date, end_date, reach_prefix_from_table, reach_prefix_from_layer )"},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate trash scores. — assess_trash","text":"df_reach Reach data Rapid Stream Assessment (RSA) dataset. start_date Beginning timeframe evaluated (ex. first day five-year assessment period). end_date End timeframe evaluated (ex. last day five-year assessment period). reach_prefix_from_table Field name prefix sourced RSA reach attribute table. consistent name stream reach attribute table RSA geodatabase. reach_prefix_from_layer Field name prefix sourced RSA stream reach layer. consistent name stream reach polyline layer RSA geodatabase.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate trash scores. — assess_trash","text":"List 2 dataframes: summary score.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate trash scores. — assess_trash","text":"","code":"assess_trash(df_reach_example, \"2020-01-01\", \"2024-12-31\", \"StreamReachAttributes\", \"StreamReaches\") #> $summary #> # A tibble: 733 × 7 #> location_name id date trash_extent score sci_subshed #> #> 1 Hickey Run {000… 2024-10-04 15:09:00 extensive 1 Hickey Run #> 2 Piney Branch {001… 2024-06-11 16:37:00 absent 10 Piney Bran… #> 3 NA {003… 2022-06-08 17:06:00 absent 10 NA #> 4 Pinehurst Branch {004… 2024-05-31 14:58:00 absent 10 Pinehurst … #> 5 Rock Creek {008… 2021-07-08 13:42:00 absent 10 Rock Creek #> 6 NA {015… 2021-07-29 16:44:00 minor 7 NA #> 7 Broad Branch {024… 2023-08-25 19:20:00 minor 7 Broad Bran… #> 8 Broad Branch {027… 2023-08-25 19:43:00 minor 7 Broad Bran… #> 9 Fort Davis Tributary {02A… 2024-09-17 16:35:00 absent 10 Fort Davis… #> 10 Klingle Valley Run {031… 2024-06-18 14:50:00 minor 7 Klingle Va… #> # ℹ 723 more rows #> # ℹ 1 more variable: ecoregion #> #> $score #> # A tibble: 22 × 2 #> sci_subshed Trash #> #> 1 Battery Kemble Creek 8.43 #> 2 Broad Branch 8.88 #> 3 Dumbarton Oaks 10 #> 4 Fenwick Branch 7.43 #> 5 Fort Davis Tributary 8.33 #> 6 Fort Dupont Tributary 7.64 #> 7 Fort Stanton Tributary 6.47 #> 8 Foundry Branch 8.83 #> 9 Hickey Run 8.04 #> 10 Klingle Valley Run 6.25 #> # ℹ 12 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate water quality scores — assess_wq","title":"Calculate water quality scores — assess_wq","text":"Calculate water quality scores","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate water quality scores — assess_wq","text":"","code":"assess_wq(df, parameter_name, unit, max_criteria = NULL, min_criteria = NULL)"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate water quality scores — assess_wq","text":"df DOEE Ambient Water Quality Data processed formatted parameter_name String indicate target parameter: \"Conductivity\", \"Dissolved oxygen ()\", \"Escherichia coli\", \"pH\", \"Temperature, water\", \"Turbidity\" unit String indicate unit water quality parameter max_criteria DC water quality criteria maximum value min_criteria DC water quality criteria minimum value","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate water quality scores — assess_wq","text":"List three dataframes: results, summary, score","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate water quality scores — assess_wq","text":"","code":"assess_wq(df_wq_processed_example, parameter_name = \"Temperature, water\", unit = \"deg C\", max_criteria = 24) #> $results #> # A tibble: 586 × 9 #> sample_id location_id date parameter unit result qualifier excursion #> #> 1 RCR01_1507… RCR01 2015-07-14 Temperat… deg C 23.2 \"\" NA #> 2 RCR01_1508… RCR01 2015-08-11 Temperat… deg C 22.8 \"\" NA #> 3 RCR01_1509… RCR01 2015-09-21 Temperat… deg C 19.4 \"\" NA #> 4 RCR01_1510… RCR01 2015-10-20 Temperat… deg C 9.10 \"\" NA #> 5 RCR01_1511… RCR01 2015-11-03 Temperat… deg C 12.5 \"\" NA #> 6 RCR01_1512… RCR01 2015-12-15 Temperat… deg C 12.4 \"\" NA #> 7 RCR01_1601… RCR01 2016-01-05 Temperat… deg C 1.85 \"\" NA #> 8 RCR01_1602… RCR01 2016-02-02 Temperat… deg C 3.21 \"\" NA #> 9 RCR01_1603… RCR01 2016-03-15 Temperat… deg C 10.8 \"\" NA #> 10 RCR01_1604… RCR01 2016-04-05 Temperat… deg C 10.3 \"\" NA #> # ℹ 576 more rows #> # ℹ 1 more variable: outside_limit #> #> $summary #> # A tibble: 25 × 9 #> sci_subshed n_tests n_tests_failed excursion nse f2 f3 wqi score #> #> 1 Battery Kem… 14 0 0 0 0 0 100 10 #> 2 Broad Branch 14 0 0 0 0 0 100 10 #> 3 Dalecarlia … 14 1 0.0904 0.00646 7.14 0.642 94.9 9.49 #> 4 Dumbarton O… 14 0 0 0 0 0 100 10 #> 5 Fenwick Bra… 14 1 0.0179 0.00128 7.14 0.128 94.9 9.49 #> 6 Fort Chapli… 13 0 0 0 0 0 100 10 #> 7 Fort Davis … 14 0 0 0 0 0 100 10 #> 8 Fort Dupont… 10 0 0 0 0 0 100 10 #> 9 Fort Stanto… 14 1 0.0321 0.00229 7.14 0.229 94.9 9.49 #> 10 Foundry Bra… 14 0 0 0 0 0 100 10 #> # ℹ 15 more rows #> #> $score #> # A tibble: 25 × 2 #> sci_subshed `Temperature, water` #> #> 1 Battery Kemble Creek 10 #> 2 Broad Branch 10 #> 3 Dalecarlia Tributary 9.49 #> 4 Dumbarton Oaks 10 #> 5 Fenwick Branch 9.49 #> 6 Fort Chaplin Tributary 10 #> 7 Fort Davis Tributary 10 #> 8 Fort Dupont Tributary 10 #> 9 Fort Stanton Tributary 9.49 #> 10 Foundry Branch 10 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate water quality scores for nutrients — assess_wq_nutrients","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"Calculate water quality scores nutrients","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"","code":"assess_wq_nutrients( df, parameter_name, piedmont_criteria, coastal_plain_criteria )"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"df DOEE Ambient Water Quality Data processed formatted parameter_name String indicate target parameter: \"Nitrogen\", \"Phosphorus, Total (P)\" piedmont_criteria Nutrient criteria recommendation Piedmont Ecoregion 64 mg/L coastal_plain_criteria Nutrient criteria recommendation Coastal Plain Ecoregion 65 mg/L","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"List three dataframes: results, summary, score","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"","code":"assess_wq_nutrients(df_wq_processed_example, parameter_name = \"Phosphorus, Total (as P)\", piedmont_criteria = 0.04, coastal_plain_criteria = 0.0225) #> $results #> # A tibble: 28 × 12 #> sample_id location_id date parameter unit result qualifier #> #> 1 20C0524-02 TNS01 2020-03-03 Phosphorus, Total (… mg/l 0.016 \"\" #> 2 20C0524-03 TDA01 2020-03-03 Phosphorus, Total (… mg/l 0.12 \"\" #> 3 20C0524-04 TBK01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 4 20C0524-06 TFB02 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 5 20C0524-07 TDO01 2020-03-03 Phosphorus, Total (… mg/l 0.046 \"\" #> 6 20C0524-10 RCR09 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 7 20C0524-11 TKV01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 8 20C0524-12 TPY01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 9 20C0524-13 TMH01 2020-03-03 Phosphorus, Total (… mg/l 0.021 \"\" #> 10 20C0524-14 TBR01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> # ℹ 18 more rows #> # ℹ 5 more variables: location_name , sci_subshed , ecoregion , #> # excursion , outside_limit #> #> $summary #> # A tibble: 25 × 9 #> sci_subshed n_tests n_tests_failed excursion nse f2 f3 wqi score #> #> 1 Battery Kembl… 1 0 0 0 0 0 100 10 #> 2 Broad Branch 1 0 0 0 0 0 100 10 #> 3 Dalecarlia Tr… 1 1 2 2 100 66.7 15.0 1.50 #> 4 Dumbarton Oaks 1 1 0.150 0.150 100 13.0 28.7 2.87 #> 5 Fenwick Branch 1 0 0 0 0 0 100 10 #> 6 Fort Chaplin … 1 0 0 0 0 0 100 10 #> 7 Fort Davis Tr… 1 1 1.44 1.44 100 59.1 17.9 1.79 #> 8 Fort Dupont T… 1 0 0 0 0 0 100 10 #> 9 Fort Stanton … 1 0 0 0 0 0 100 10 #> 10 Foundry Branch 1 0 0 0 0 0 100 10 #> # ℹ 15 more rows #> #> $score #> # A tibble: 25 × 2 #> sci_subshed `Phosphorus, Total (as P)` #> #> 1 Battery Kemble Creek 10 #> 2 Broad Branch 10 #> 3 Dalecarlia Tributary 1.50 #> 4 Dumbarton Oaks 2.87 #> 5 Fenwick Branch 10 #> 6 Fort Chaplin Tributary 10 #> 7 Fort Davis Tributary 1.79 #> 8 Fort Dupont Tributary 10 #> 9 Fort Stanton Tributary 10 #> 10 Foundry Branch 10 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":null,"dir":"Reference","previous_headings":"","what":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"Assign score based lookup maximum number. Used dumpsites EIA.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"","code":"assign_score_from_maximum(value, lookup_table, lookup_field)"},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"value Number dumpsites per mile given subwatershed. lookup_table Lookup table number dumpsites per mile allowed per score (1-10). lookup_field Column name containing maximum number sites per mile.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"Numerical score (1-10)","code":""},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"","code":"assign_score_from_maximum(0.5, dumpsite_score, \"max_sites_per_mile\") #> [1] 7"},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":null,"dir":"Reference","previous_headings":"","what":"Format downloaded water quality data — format_wq","title":"Format downloaded water quality data — format_wq","text":"Format downloaded water quality data","code":""},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Format downloaded water quality data — format_wq","text":"","code":"format_wq(df)"},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Format downloaded water quality data — format_wq","text":"df DOEE Ambient Water Quality data downloaded https://dcdoeepub.equisonline.com/PUBLIC.html","code":""},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Format downloaded water quality data — format_wq","text":"Formatted dataframe","code":""},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Format downloaded water quality data — format_wq","text":"","code":"format_wq(df_wq_example) #> Warning: There was 1 warning in `dplyr::mutate()`. #> ℹ In argument: `result_format = as.numeric(result_format)`. #> Caused by warning: #> ! NAs introduced by coercion #> # A tibble: 98,073 × 14 #> `Sample ID` Location Date Analyte Unit SAMPLETIME Result `Lab Qualifiers` #> #> 1 20C0524-01 PTB01 03 Mar… Nitrat… mg/l 10:03 1.23 M4 #> 2 20C0524-01 PTB01 03 Mar… Nitrog… mg/l 10:03 2.91 NA #> 3 20C0524-01 PTB01 03 Mar… Phosph… NA 10:03 < 0.0… NA #> 4 20C0524-01 PTB01 03 Mar… Total … mg/l 10:03 1.68 NA #> 5 20C0524-02 TNS01 03 Mar… Nitrat… mg/l 10:29 0.489 NA #> 6 20C0524-02 TNS01 03 Mar… Nitrog… mg/l 10:29 1.89 NA #> 7 20C0524-02 TNS01 03 Mar… Phosph… mg/l 10:29 0.016 NA #> 8 20C0524-02 TNS01 03 Mar… Total … mg/l 10:29 1.40 NA #> 9 20C0524-03 TDA01 03 Mar… Nitrat… mg/l 10:53 2.10 NA #> 10 20C0524-03 TDA01 03 Mar… Nitrog… mg/l 10:53 3.50 NA #> # ℹ 98,063 more rows #> # ℹ 6 more variables: `Validator Qualifiers` , #> # `Interpreted Qualifiers` , result_format , #> # qualifier_format , location_format , date_format "},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":null,"dir":"Reference","previous_headings":"","what":"Pipe operator — %>%","title":"Pipe operator — %>%","text":"See magrittr::%>% details.","code":""},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Pipe operator — %>%","text":"","code":"lhs %>% rhs"},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Pipe operator — %>%","text":"lhs value magrittr placeholder. rhs function call using magrittr semantics.","code":""},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Pipe operator — %>%","text":"result calling rhs(lhs).","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":null,"dir":"Reference","previous_headings":"","what":"Process formatted water quality data — process_wq","title":"Process formatted water quality data — process_wq","text":"Process formatted water quality data","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Process formatted water quality data — process_wq","text":"","code":"process_wq(df, start_date, end_date)"},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Process formatted water quality data — process_wq","text":"df Formatted dataframe returned format_wq start_date Beginning timeframe evaluated (ex. first day five-year assessment period). end_date End timeframe evaluated (ex. last day five-year assessment period).","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Process formatted water quality data — process_wq","text":"Processed dataframe","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Process formatted water quality data — process_wq","text":"","code":"process_wq(df_wq_formatted_example, \"2015-07-01\", \"2020-06-30\") #> # A tibble: 3,734 × 7 #> sample_id location_id date parameter unit result qualifier #> #> 1 20C0524-02 TNS01 2020-03-03 Nitrogen mg/l 1.89 \"\" #> 2 20C0524-02 TNS01 2020-03-03 Phosphorus, Total (… mg/l 0.016 \"\" #> 3 20C0524-03 TDA01 2020-03-03 Nitrogen mg/l 3.5 \"\" #> 4 20C0524-03 TDA01 2020-03-03 Phosphorus, Total (… mg/l 0.12 \"\" #> 5 20C0524-04 TBK01 2020-03-03 Nitrogen mg/l 5.12 \"\" #> 6 20C0524-04 TBK01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 7 20C0524-06 TFB02 2020-03-03 Nitrogen mg/l 2.65 \"\" #> 8 20C0524-06 TFB02 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 9 20C0524-07 TDO01 2020-03-03 Nitrogen mg/l 4.46 \"\" #> 10 20C0524-07 TDO01 2020-03-03 Phosphorus, Total (… mg/l 0.046 \"\" #> # ℹ 3,724 more rows"}] +[{"path":"https://limnotech.github.io/SCI/LICENSE.html","id":null,"dir":"","previous_headings":"","what":"MIT License","title":"MIT License","text":"Copyright (c) 2024 SCI authors Permission hereby granted, free charge, person obtaining copy software associated documentation files (“Software”), deal Software without restriction, including without limitation rights use, copy, modify, merge, publish, distribute, sublicense, /sell copies Software, permit persons Software furnished , subject following conditions: copyright notice permission notice shall included copies substantial portions Software. SOFTWARE PROVIDED “”, WITHOUT WARRANTY KIND, EXPRESS IMPLIED, INCLUDING LIMITED WARRANTIES MERCHANTABILITY, FITNESS PARTICULAR PURPOSE NONINFRINGEMENT. EVENT SHALL AUTHORS COPYRIGHT HOLDERS LIABLE CLAIM, DAMAGES LIABILITY, WHETHER ACTION CONTRACT, TORT OTHERWISE, ARISING , CONNECTION SOFTWARE USE DEALINGS SOFTWARE.","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"effective-impervious-area-eia","dir":"Articles","previous_headings":"","what":"Effective Impervious Area (EIA)","title":"FAQs","text":"Rock Creek mainstem includes MS4 CSS. Rock Creek watershed also includes additional area beyond Rock Creek SCI subwatersheds. area Rock Creek used SCI analysis? entire Rock Creek watershed used SCI analysis. See “eia_rock_creek” tab lookup spreadsheet full list subsheds included.","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"water-quality-index","dir":"Articles","previous_headings":"","what":"Water Quality Index","title":"FAQs","text":"non-detects handled dataset? reporting limit used non-detect values. example, results reported <0.010, got incorporated analysis 0.010. water quality criteria used various parameters? See Section 1 (Indicators Metrics) Stream Condition Index Methodology (2022) Among non-nutrient WQ indicators, maximum criteria, minimum criteria, maximum minimum criteria. assess_wq() function handles three cases. nutrient WQ indicatiors (TN TP), assess_wq_nutrients() requires criteria input Piedmont Coastal Plain Level III Ecoregions","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"aquatic-life-index","dir":"Articles","previous_headings":"","what":"Aquatic Life Index","title":"FAQs","text":"aquatic life metrics calculated?  four metrics get calculated outside R package get incorporated lookup tables data/lookup_tables.xlsx. Connection Mainstem unlikely change year year. Subsheds assigned score either 10 (connected mainstem) 1 (connected mainstem). option calculating fish IBI, habitat PHI, macroinvertebrate HBI MBSStools R package.  IBI, PHI, HBI scores get translated SCI score added relevant lookup table. IBI, PHI, HBI scores translate SCI scores? See conversion tables :","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"stations","dir":"Articles","previous_headings":"","what":"Stations","title":"FAQs","text":"stations included tributary? station locations within given tributary incorporated analysis included location_id lookup table. Rock Creek, include stations specifically Upper Rock Creek Lower Rock Creek. separate scores calculated upper lower segments Rock Creek Watts Branch? , single score calculated indicator Rock Creek Watts Branch","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"timeframe","dir":"Articles","previous_headings":"","what":"Timeframe","title":"FAQs","text":"timeframe used analysis? goal use five recent calendar years data (e.g. 01/01/2020 – 12/31/2024). exceptions may made extend timeframe beyond five years found inadequate data.","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"troubleshooting","dir":"Articles","previous_headings":"","what":"Troubleshooting","title":"FAQs","text":"try run R code, get warning. ? warnings expected issue running R code. able proceed steps normal. still questions concerns, contact package maintainer (see ‘Contacts’ subsection FAQs) try run R code, get error. ? Possible things try include: Read error message see can pinpoint error occurred. ·sure input variables specified correct. Confirm file names pathways correct. Check whether datasets format changes. Changes field names, changes waterbody names, new sampling IDs, etc. can disrupt code. ’re unable resolve error, still questions concerns, contact package maintainer (see ‘Contacts’ subsection FAQs)","code":""},{"path":"https://limnotech.github.io/SCI/articles/faqs.html","id":"contacts","dir":"Articles","previous_headings":"","what":"Contacts","title":"FAQs","text":"contact questions R package? R package currently maintained Shanna Rucker Blount (LimnoTech). can contact srucker@limno.com","code":""},{"path":[]},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"notes","dir":"Articles","previous_headings":"Trash and Dumpsites","what":"Notes","title":"Pre-Processing","text":"R functions assess_trash() assess_dumpsites() require dataframe inputs processed RSA reach data RSA point data (dumpsites ) watershed layer included RSA database used due issue found. example, layer “Watersheds_Notes_RSA2021” missing area Rock Creek. Instead used: DC Subwatersheds Open Data Layer.","code":""},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"steps-for-processing-rsa-stream-reach-and-point-data","dir":"Articles","previous_headings":"Trash and Dumpsites","what":"Steps for processing RSA stream reach and point data","title":"Pre-Processing","text":"StreamPoints feature layer, copy GlobalID new field Set data type new GUID New field name suggestion: GlobalID_copy step needed GlobalID information included intersected table Intersect DC Subwatersheds Open Data Layer StreamPoints layer Run ‘Add Attribute Index’ tool Intersected layer ID field indexed: Field Index: GlobalID_copy Index Name: GlobalID_index Join StreamPoints attributes table StreamPoints feature layer. See settings screenshot : ’One Many” Join Operation option, make sure attribute table feature layer geodatabase Export joined StreamPoints using ‘Table Excel’ geoprocessing tool Save file R package data folder Update file name input variable workflow.R file input variables section workflow.R, update field name prefix variables (needed) Repeat steps 1, 3, 4, 5 6 using StreamReaches feature layer attribute table Skip Step 2 since subshed information already included field StreamReaches feature layer","code":""},{"path":[]},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"notes-1","dir":"Articles","previous_headings":"Effective Impervious Area (EIA)","what":"Notes","title":"Pre-Processing","text":"Subshed area calculations done external R package. Lookup tables must updated manually. final EIA percentage score calculations done R. BMP Data timeframe always begin January 1, 2000. two lookup tables need manually updated annual SCI run: eia_subsheds eia_rock_creek addition sci_subshed field, main columns lookup tables get used R functions prefix “analysis_” Analysis_subshed_m2 Analysis_impervious_m2 Analysis_impervious_treated_m2 columns particular years reference GIS analysis repeated numbers needed later Special subshed notes: Piney Branch sci_subshed MS4 CSS. eia_subsheds lookup table includes MS4 CSS subsheds two separate segments. Reservation 630 used sci_subshed, Melvin Hazen Valley Branch used lookup subwatershed data sources","code":""},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"steps-for-impervious-area-calculation","dir":"Articles","previous_headings":"Effective Impervious Area (EIA)","what":"Steps for impervious area calculation","title":"Pre-Processing","text":"perform following steps new impervious surface available Open Data DC Start latest impervious surface layer Open Data DC (e.g., 2023 Impervious Surface) Intersect impervious surface layer Open Data DC Subwatersheds layer Dissolve intersected layer summarize reduce number features (See screenshot Dissolve settings) Export Dissolve table excel. Update lookup tables R package spreadsheet. data/lookup_tables.xlsx, modify eia_subsheds eia_rock_creek spreadsheet tabs: Add new column called “20XX_impervious_m2” bring data Dissolve table Update “analysis_impervious_m2” column referencing new column created Save spreadsheet workflow.R script, uncomment code eia_subsheds eia_rock_creek. Run code update two .rda files. Re-comment code.","code":""},{"path":"https://limnotech.github.io/SCI/articles/pre-processing.html","id":"steps-for-impervious-treated-calculation","dir":"Articles","previous_headings":"Effective Impervious Area (EIA)","what":"Steps for impervious treated calculation","title":"Pre-Processing","text":"Export BMP data SGS using IPMT query - Start date: January 1, 2000 - End date: December 31, 20XX (use last year target 5-year period) Excel, remove BMP records meet Chesapeake Bay inspection criteria. BMPs installed 10 years ago, BMPs must passing inspection within last 10 years. Use following three fields identify remove BMP records older 10 years passing inspection: (1) Installation Date; (2) recent Inspection (3) Passed failed recent maintenance inspection? example, SCI timeframe 2020-2024, BMPs installed prior January 1, 2015 must passing inspection. ArcGIS Pro, bring table SGS export display XY data Intersect BMP point data Open Data DC Subwatersheds layer Dissolve using settings screenshot Export Dissolve table Excel Convert impervious area sum ft2 m2 Update lookup tables R package spreadsheet. data/lookup_tables.xlsx, modify eia_subsheds eia_rock_creek spreadsheet tabs: Add new column called “20XX_impervious_treated_m2” bring data Dissolve table Update “analysis_impervious_treated_m2” column referencing new column created Save spreadsheet workflow.R script, uncomment code eia_subsheds eia_rock_creek. Run code update two .rda files. Re-comment code.","code":""},{"path":"https://limnotech.github.io/SCI/authors.html","id":null,"dir":"","previous_headings":"","what":"Authors","title":"Authors and Citation","text":"Shanna Rucker. Author, maintainer.","code":""},{"path":"https://limnotech.github.io/SCI/authors.html","id":"citation","dir":"","previous_headings":"","what":"Citation","title":"Authors and Citation","text":"Rucker S (2025). SCI: Stream Condition Index. R package version 0.0.0.9000, https://limnotech.github.io/SCI/, https://github.com/LimnoTech/SCI.","code":"@Manual{, title = {SCI: Stream Condition Index}, author = {Shanna Rucker}, year = {2025}, note = {R package version 0.0.0.9000, https://limnotech.github.io/SCI/}, url = {https://github.com/LimnoTech/SCI}, }"},{"path":"https://limnotech.github.io/SCI/index.html","id":"sci","dir":"","previous_headings":"","what":"Stream Condition Index","title":"Stream Condition Index","text":"goal SCI calculate scores hydrology, water quality, human health, aquatic biology metrics assess streams within District Columbia.","code":""},{"path":"https://limnotech.github.io/SCI/index.html","id":"installation","dir":"","previous_headings":"","what":"Installation","title":"Stream Condition Index","text":"can install development version SCI GitHub :","code":"# install.packages(\"devtools\") devtools::install_github(\"LimnoTech/SCI\")"},{"path":"https://limnotech.github.io/SCI/index.html","id":"using-the-package","dir":"","previous_headings":"","what":"Using the Package","title":"Stream Condition Index","text":"Produce SCI scores just simple steps: Complete pre-processing steps, needed. RStudio, open workflow.R script. Update “Input Variables” section script confirm latest data files saved data folder. Update lookup_tables.xlsx run code “Lookup Tables” section script, needed. Note: step typically necessary. Click Source button run workflow.R script. View results all_scores.csv.","code":""},{"path":"https://limnotech.github.io/SCI/reference/SCI-package.html","id":null,"dir":"Reference","previous_headings":"","what":"SCI: Stream Condition Index — SCI-package","title":"SCI: Stream Condition Index — SCI-package","text":"Calculates scores hydrology, water quality, human health, aquatic biology metrics assess streams within District Columbia.","code":""},{"path":[]},{"path":"https://limnotech.github.io/SCI/reference/SCI-package.html","id":"author","dir":"Reference","previous_headings":"","what":"Author","title":"SCI: Stream Condition Index — SCI-package","text":"Maintainer: Shanna Rucker srucker@limno.com","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate dumpsite scores — assess_dumpsites","title":"Calculate dumpsite scores — assess_dumpsites","text":"Calculate dumpsite scores","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate dumpsite scores — assess_dumpsites","text":"","code":"assess_dumpsites( df_point, df_reach, reach_prefix_from_table, reach_prefix_from_layer, point_prefix_from_table, point_prefix_from_layer )"},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate dumpsite scores — assess_dumpsites","text":"df_point Point data Rapid Stream Assessment (RSA) dataset. df_reach Reach data Rapid Stream Assessment (RSA) dataset. reach_prefix_from_table Field name prefix sourced RSA reach attribute table. consistent name stream reach attribute table RSA geodatabase. reach_prefix_from_layer Field name prefix sourced RSA stream reach layer. consistent name stream reach polyline layer RSA geodatabase. point_prefix_from_table Field name prefix sourced RSA point attribute table. consistent name point assessment attribute table RSA geodatabase. point_prefix_from_layer Field name prefix sourced RSA point layer. consistent name intersected stream point layer RSA geodatabase.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate dumpsite scores — assess_dumpsites","text":"List 2 dataframes: summary score.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_dumpsites.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate dumpsite scores — assess_dumpsites","text":"","code":"assess_dumpsites(df_point_example, df_reach_example, \"StreamReachAttributes\", \"StreamReaches\", \"StreamPointAttributes\", \"StreamPoints_Intersect\" ) #> $summary #> # A tibble: 9 × 8 #> # Rowwise: #> sci_subshed dumpsite_count applied_influence subshed_length_meters #> #> 1 Broad Branch 2 -1 9840. #> 2 Fort Davis Tributary 3 -2 1469. #> 3 Fort Dupont Tributary 6 -2 9175. #> 4 Fort Stanton Tributary 1 -2 2507. #> 5 Normanstone Creek 2 -1 3527. #> 6 Oxon Run 1 -1 10144. #> 7 Pinehurst Branch 1 -1 6001. #> 8 Rock Creek 4 -2 38328. #> 9 Watts Branch 2 -2 6035. #> # ℹ 4 more variables: subshed_length_miles , sites_per_mile , #> # score , score_weighted #> #> $score #> # A tibble: 25 × 2 #> # Rowwise: #> sci_subshed Dumpsites #> #> 1 Broad Branch 7 #> 2 Fort Davis Tributary 1 #> 3 Fort Dupont Tributary 2 #> 4 Fort Stanton Tributary 4 #> 5 Normanstone Creek 4 #> 6 Oxon Run 8 #> 7 Pinehurst Branch 7 #> 8 Rock Creek 7 #> 9 Watts Branch 5 #> 10 Battery Kemble Creek 10 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate scores for effective impervious area (EIA) — assess_eia","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"Calculate scores effective impervious area (EIA)","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"","code":"assess_eia()"},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"List 2 dataframes: summary score.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_eia.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate scores for effective impervious area (EIA) — assess_eia","text":"","code":"assess_eia() #> $summary #> # A tibble: 25 × 7 #> # Rowwise: #> sci_subshed subshed impervious impervious_treated effective_impervious #> #> 1 Battery Kemble Cr… 9.37e5 201840. 26487. 175353. #> 2 Broad Branch 4.65e6 1627836. 76874. 1550962. #> 3 Dalecarlia Tribut… 1.09e6 295100. 23139. 271960. #> 4 Dumbarton Oaks 5.51e5 153336. 11829. 141508. #> 5 Fenwick Branch 8.94e5 284547. 6249. 278298. #> 6 Fort Chaplin Trib… 1.18e6 423427. 30927. 392500. #> 7 Fort Davis Tribut… 9.47e5 328257. 12328. 315929. #> 8 Fort Dupont Tribu… 1.66e6 154824. 15534. 139289. #> 9 Fort Stanton Trib… 1.12e6 488751. 64799. 423952. #> 10 Foundry Branch 4.59e6 1725013. 198108. 1526905. #> # ℹ 15 more rows #> # ℹ 2 more variables: percent_effective_impervious , score #> #> $score #> # A tibble: 25 × 2 #> # Rowwise: #> sci_subshed `Effective Impervious Area` #> #> 1 Battery Kemble Creek 6 #> 2 Broad Branch 2 #> 3 Dalecarlia Tributary 4 #> 4 Dumbarton Oaks 4 #> 5 Fenwick Branch 3 #> 6 Fort Chaplin Tributary 2 #> 7 Fort Davis Tributary 2 #> 8 Fort Dupont Tributary 8 #> 9 Fort Stanton Tributary 1 #> 10 Foundry Branch 2 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate trash scores. — assess_trash","title":"Calculate trash scores. — assess_trash","text":"Calculate trash scores.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate trash scores. — assess_trash","text":"","code":"assess_trash(df_reach, reach_prefix_from_table, reach_prefix_from_layer)"},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate trash scores. — assess_trash","text":"df_reach Reach data Rapid Stream Assessment (RSA) dataset. reach_prefix_from_table Field name prefix sourced RSA reach attribute table. consistent name stream reach attribute table RSA geodatabase. reach_prefix_from_layer Field name prefix sourced RSA stream reach layer. consistent name stream reach polyline layer RSA geodatabase.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate trash scores. — assess_trash","text":"List 2 dataframes: summary score.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_trash.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate trash scores. — assess_trash","text":"","code":"assess_trash(df_reach_example, \"StreamReachAttributes\", \"StreamReaches\") #> Error in dplyr::filter(., date >= start_date & date <= end_date): ℹ In argument: `date >= start_date & date <= end_date`. #> Caused by error: #> ! object 'start_date' not found"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate water quality scores — assess_wq","title":"Calculate water quality scores — assess_wq","text":"Calculate water quality scores","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate water quality scores — assess_wq","text":"","code":"assess_wq(df, parameter_name, unit, max_criteria = NULL, min_criteria = NULL)"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate water quality scores — assess_wq","text":"df DOEE Ambient Water Quality Data processed formatted parameter_name String indicate target parameter: \"Conductivity\", \"Dissolved oxygen ()\", \"Escherichia coli\", \"pH\", \"Temperature, water\", \"Turbidity\" unit String indicate unit water quality parameter max_criteria DC water quality criteria maximum value min_criteria DC water quality criteria minimum value","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate water quality scores — assess_wq","text":"List three dataframes: results, summary, score","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate water quality scores — assess_wq","text":"","code":"assess_wq(df_wq_processed_example, parameter_name = \"Temperature, water\", unit = \"deg C\", max_criteria = 24) #> $results #> # A tibble: 586 × 9 #> sample_id location_id date parameter unit result qualifier excursion #> #> 1 RCR01_1507… RCR01 2015-07-14 Temperat… deg C 23.2 \"\" NA #> 2 RCR01_1508… RCR01 2015-08-11 Temperat… deg C 22.8 \"\" NA #> 3 RCR01_1509… RCR01 2015-09-21 Temperat… deg C 19.4 \"\" NA #> 4 RCR01_1510… RCR01 2015-10-20 Temperat… deg C 9.10 \"\" NA #> 5 RCR01_1511… RCR01 2015-11-03 Temperat… deg C 12.5 \"\" NA #> 6 RCR01_1512… RCR01 2015-12-15 Temperat… deg C 12.4 \"\" NA #> 7 RCR01_1601… RCR01 2016-01-05 Temperat… deg C 1.85 \"\" NA #> 8 RCR01_1602… RCR01 2016-02-02 Temperat… deg C 3.21 \"\" NA #> 9 RCR01_1603… RCR01 2016-03-15 Temperat… deg C 10.8 \"\" NA #> 10 RCR01_1604… RCR01 2016-04-05 Temperat… deg C 10.3 \"\" NA #> # ℹ 576 more rows #> # ℹ 1 more variable: outside_limit #> #> $summary #> # A tibble: 25 × 9 #> sci_subshed n_tests n_tests_failed excursion nse f2 f3 wqi score #> #> 1 Battery Kem… 14 0 0 0 0 0 100 10 #> 2 Broad Branch 14 0 0 0 0 0 100 10 #> 3 Dalecarlia … 14 1 0.0904 0.00646 7.14 0.642 94.9 9.49 #> 4 Dumbarton O… 14 0 0 0 0 0 100 10 #> 5 Fenwick Bra… 14 1 0.0179 0.00128 7.14 0.128 94.9 9.49 #> 6 Fort Chapli… 13 0 0 0 0 0 100 10 #> 7 Fort Davis … 14 0 0 0 0 0 100 10 #> 8 Fort Dupont… 10 0 0 0 0 0 100 10 #> 9 Fort Stanto… 14 1 0.0321 0.00229 7.14 0.229 94.9 9.49 #> 10 Foundry Bra… 14 0 0 0 0 0 100 10 #> # ℹ 15 more rows #> #> $score #> # A tibble: 25 × 2 #> sci_subshed `Temperature, water` #> #> 1 Battery Kemble Creek 10 #> 2 Broad Branch 10 #> 3 Dalecarlia Tributary 9.49 #> 4 Dumbarton Oaks 10 #> 5 Fenwick Branch 9.49 #> 6 Fort Chaplin Tributary 10 #> 7 Fort Davis Tributary 10 #> 8 Fort Dupont Tributary 10 #> 9 Fort Stanton Tributary 9.49 #> 10 Foundry Branch 10 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":null,"dir":"Reference","previous_headings":"","what":"Calculate water quality scores for nutrients — assess_wq_nutrients","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"Calculate water quality scores nutrients","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"","code":"assess_wq_nutrients( df, parameter_name, piedmont_criteria, coastal_plain_criteria )"},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"df DOEE Ambient Water Quality Data processed formatted parameter_name String indicate target parameter: \"Nitrogen\", \"Phosphorus, Total (P)\" piedmont_criteria Nutrient criteria recommendation Piedmont Ecoregion 64 mg/L coastal_plain_criteria Nutrient criteria recommendation Coastal Plain Ecoregion 65 mg/L","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"List three dataframes: results, summary, score","code":""},{"path":"https://limnotech.github.io/SCI/reference/assess_wq_nutrients.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Calculate water quality scores for nutrients — assess_wq_nutrients","text":"","code":"assess_wq_nutrients(df_wq_processed_example, parameter_name = \"Phosphorus, Total (as P)\", piedmont_criteria = 0.04, coastal_plain_criteria = 0.0225) #> $results #> # A tibble: 28 × 12 #> sample_id location_id date parameter unit result qualifier #> #> 1 20C0524-02 TNS01 2020-03-03 Phosphorus, Total (… mg/l 0.016 \"\" #> 2 20C0524-03 TDA01 2020-03-03 Phosphorus, Total (… mg/l 0.12 \"\" #> 3 20C0524-04 TBK01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 4 20C0524-06 TFB02 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 5 20C0524-07 TDO01 2020-03-03 Phosphorus, Total (… mg/l 0.046 \"\" #> 6 20C0524-10 RCR09 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 7 20C0524-11 TKV01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 8 20C0524-12 TPY01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 9 20C0524-13 TMH01 2020-03-03 Phosphorus, Total (… mg/l 0.021 \"\" #> 10 20C0524-14 TBR01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> # ℹ 18 more rows #> # ℹ 5 more variables: location_name , sci_subshed , ecoregion , #> # excursion , outside_limit #> #> $summary #> # A tibble: 25 × 9 #> sci_subshed n_tests n_tests_failed excursion nse f2 f3 wqi score #> #> 1 Battery Kembl… 1 0 0 0 0 0 100 10 #> 2 Broad Branch 1 0 0 0 0 0 100 10 #> 3 Dalecarlia Tr… 1 1 2 2 100 66.7 15.0 1.50 #> 4 Dumbarton Oaks 1 1 0.150 0.150 100 13.0 28.7 2.87 #> 5 Fenwick Branch 1 0 0 0 0 0 100 10 #> 6 Fort Chaplin … 1 0 0 0 0 0 100 10 #> 7 Fort Davis Tr… 1 1 1.44 1.44 100 59.1 17.9 1.79 #> 8 Fort Dupont T… 1 0 0 0 0 0 100 10 #> 9 Fort Stanton … 1 0 0 0 0 0 100 10 #> 10 Foundry Branch 1 0 0 0 0 0 100 10 #> # ℹ 15 more rows #> #> $score #> # A tibble: 25 × 2 #> sci_subshed `Phosphorus, Total (as P)` #> #> 1 Battery Kemble Creek 10 #> 2 Broad Branch 10 #> 3 Dalecarlia Tributary 1.50 #> 4 Dumbarton Oaks 2.87 #> 5 Fenwick Branch 10 #> 6 Fort Chaplin Tributary 10 #> 7 Fort Davis Tributary 1.79 #> 8 Fort Dupont Tributary 10 #> 9 Fort Stanton Tributary 10 #> 10 Foundry Branch 10 #> # ℹ 15 more rows #>"},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":null,"dir":"Reference","previous_headings":"","what":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"Assign score based lookup maximum number. Used dumpsites EIA.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"","code":"assign_score_from_maximum(value, lookup_table, lookup_field)"},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"value Number dumpsites per mile given subwatershed. lookup_table Lookup table number dumpsites per mile allowed per score (1-10). lookup_field Column name containing maximum number sites per mile.","code":""},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"Numerical score (1-10)","code":""},{"path":"https://limnotech.github.io/SCI/reference/assign_score_from_maximum.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Assign score based on lookup of maximum number. Used for dumpsites and EIA. — assign_score_from_maximum","text":"","code":"assign_score_from_maximum(0.5, dumpsite_score, \"max_sites_per_mile\") #> [1] 7"},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":null,"dir":"Reference","previous_headings":"","what":"Format downloaded water quality data — format_wq","title":"Format downloaded water quality data — format_wq","text":"Format downloaded water quality data","code":""},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Format downloaded water quality data — format_wq","text":"","code":"format_wq(df)"},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Format downloaded water quality data — format_wq","text":"df DOEE Ambient Water Quality data downloaded https://dcdoeepub.equisonline.com/PUBLIC.html","code":""},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Format downloaded water quality data — format_wq","text":"Formatted dataframe","code":""},{"path":"https://limnotech.github.io/SCI/reference/format_wq.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Format downloaded water quality data — format_wq","text":"","code":"format_wq(df_wq_example) #> Warning: There was 1 warning in `dplyr::mutate()`. #> ℹ In argument: `result_format = as.numeric(result_format)`. #> Caused by warning: #> ! NAs introduced by coercion #> # A tibble: 98,073 × 14 #> `Sample ID` Location Date Analyte Unit SAMPLETIME Result `Lab Qualifiers` #> #> 1 20C0524-01 PTB01 03 Mar… Nitrat… mg/l 10:03 1.23 M4 #> 2 20C0524-01 PTB01 03 Mar… Nitrog… mg/l 10:03 2.91 NA #> 3 20C0524-01 PTB01 03 Mar… Phosph… NA 10:03 < 0.0… NA #> 4 20C0524-01 PTB01 03 Mar… Total … mg/l 10:03 1.68 NA #> 5 20C0524-02 TNS01 03 Mar… Nitrat… mg/l 10:29 0.489 NA #> 6 20C0524-02 TNS01 03 Mar… Nitrog… mg/l 10:29 1.89 NA #> 7 20C0524-02 TNS01 03 Mar… Phosph… mg/l 10:29 0.016 NA #> 8 20C0524-02 TNS01 03 Mar… Total … mg/l 10:29 1.40 NA #> 9 20C0524-03 TDA01 03 Mar… Nitrat… mg/l 10:53 2.10 NA #> 10 20C0524-03 TDA01 03 Mar… Nitrog… mg/l 10:53 3.50 NA #> # ℹ 98,063 more rows #> # ℹ 6 more variables: `Validator Qualifiers` , #> # `Interpreted Qualifiers` , result_format , #> # qualifier_format , location_format , date_format "},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":null,"dir":"Reference","previous_headings":"","what":"Pipe operator — %>%","title":"Pipe operator — %>%","text":"See magrittr::%>% details.","code":""},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Pipe operator — %>%","text":"","code":"lhs %>% rhs"},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Pipe operator — %>%","text":"lhs value magrittr placeholder. rhs function call using magrittr semantics.","code":""},{"path":"https://limnotech.github.io/SCI/reference/pipe.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Pipe operator — %>%","text":"result calling rhs(lhs).","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":null,"dir":"Reference","previous_headings":"","what":"Process formatted water quality data — process_wq","title":"Process formatted water quality data — process_wq","text":"Process formatted water quality data","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Process formatted water quality data — process_wq","text":"","code":"process_wq(df, start_date, end_date)"},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Process formatted water quality data — process_wq","text":"df Formatted dataframe returned format_wq start_date Beginning timeframe evaluated (ex. first day five-year assessment period). end_date End timeframe evaluated (ex. last day five-year assessment period).","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Process formatted water quality data — process_wq","text":"Processed dataframe","code":""},{"path":"https://limnotech.github.io/SCI/reference/process_wq.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Process formatted water quality data — process_wq","text":"","code":"process_wq(df_wq_formatted_example, \"2015-07-01\", \"2020-06-30\") #> # A tibble: 3,734 × 7 #> sample_id location_id date parameter unit result qualifier #> #> 1 20C0524-02 TNS01 2020-03-03 Nitrogen mg/l 1.89 \"\" #> 2 20C0524-02 TNS01 2020-03-03 Phosphorus, Total (… mg/l 0.016 \"\" #> 3 20C0524-03 TDA01 2020-03-03 Nitrogen mg/l 3.5 \"\" #> 4 20C0524-03 TDA01 2020-03-03 Phosphorus, Total (… mg/l 0.12 \"\" #> 5 20C0524-04 TBK01 2020-03-03 Nitrogen mg/l 5.12 \"\" #> 6 20C0524-04 TBK01 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 7 20C0524-06 TFB02 2020-03-03 Nitrogen mg/l 2.65 \"\" #> 8 20C0524-06 TFB02 2020-03-03 Phosphorus, Total (… NA 0.01 \"<\" #> 9 20C0524-07 TDO01 2020-03-03 Nitrogen mg/l 4.46 \"\" #> 10 20C0524-07 TDO01 2020-03-03 Phosphorus, Total (… mg/l 0.046 \"\" #> # ℹ 3,724 more rows"}]