Methods

This is a compilation of three data sets

North Temperate Lakes LTER: Chemical Limnology of Primary Study Lakes: Nutrients, pH and Carbon 1981 - current
https://lter.limnology.wisc.edu/dataset/north-temperate-lakes-lter-chemical-limnology-primary-study-lakes-nutrients-ph-and-carbon-19

North Temperate Lakes LTER: Physical Limnology of Primary Study Lakes 1981 - current
https://lter.limnology.wisc.edu/dataset/north-temperate-lakes-lter-physical-limnology-primary-study-lakes-1981-current

North Temperate Lakes LTER: Ice Duration - Trout Lake Area 1981 - current
https://lter.limnology.wisc.edu/dataset/north-temperate-lakes-lter-ice-duration-trout-lake-area-1981-current

Methods

This dataset is a compilations of:
USGS hydrologic model estimates (Hunt et al. in preparation)
and
North Temperate Lakes LTER: Ice Duration - Trout Lake Area 1981 - current
https://lter.limnology.wisc.edu/dataset/north-temperate-lakes-lter-ice-duration-trout-lake-area-1981-current

Methods

Each ecosystem service was quantified and mapped by using empirical estimates and spatially explicit model for the terrestrial landscape of the Yahara Watershed for 2006. Crop production (expected annual crop yield, bu per yr) Crop yield was estimated for the four major crop types (corn, soybean, winter wheat and oats) that account for 98.5 percent of the cultivated land in the watershed by overlaying maps of crop types and soil-specific crop yield estimates. The spatial distribution of each crop was obtained from the 2006 Cropland Data Layer (CDL) from the National Agricultural Statistics Service (NASS) and soil productivity data were extracted from Soil Survey Geographic (SSURGO) database. Crop and soil data were converted to 30 m resolution and the two maps were overlain to estimate crop yield in each cell. For each crop-soil combination, crop area was multiplied by the estimated yield per unit area. Estimates for each crop type were summed to map estimated crop yield for 2006. Pasture production (expected annual forage yield, animal-unit-month per year ) As for crop production, forage yield was estimated by overlaying the distribution of all forage crops (alfalfa, hay and pasture/grass) and soil specific yield estimates. The spatial distribution of each forage crop was also derived from 2006 CDL, and rescaled to 30 m grid prior to calculation. The SSURGO soil productivity layer provided estimates of potential annual yield per unit area for each forage crop. Overlay analyses were performed for each forage-soil combination, as done for crops, and summed to obtain the total expected forage yield in the watershed for 2006. Freshwater supply (annual groundwater recharge, cm per year) . Groundwater recharge was quantified and mapped using the modified Thornthwaite-Mather Soil-Water-Balance (SWB) model. SWB is a deterministic, physically based and quasi three-dimensional model that accounts for precipitation, evaporation, interception, surface runoff, soil moisture storage and snowmelt. Groundwater recharge was calculated on a grid cell basis at a daily step with the following mass balance equation<p align="center">Recharge= (precipitation + snowmelt + inflow) &ndash;<p align="center">(interception + outflow + evapotranspiration) &ndash; delta soil moisture<p align="center"> We ran the model for three years (2004 to 2006) at 30m resolution, with the first two years as spin up of antecedent conditions (e.g. soil moisture and snow cover) that influence groundwater recharge for the focal year of 2006. Carbon storage (metric tonsper ha) We estimated the amount of carbon stored in each 30 m cell in the Yahara Watershed by summing four major carbon pools: aboveground biomass, belowground biomass, soil carbon and deadwood/litter. Our quantification for each pool was based mainly on carbon estimates from the IPCC tier-I approach and other published field studies of carbon density and was estimated by land-use/cover type.Groundwater quality (probability of groundwater nitrate concentration greater than 3.0 mg per liter, unitless 0 to1) Groundwater nitrate data were obtained from Groundwater Retrieve Network (GRN), Wisconsin Department of Natural Resources (DNR). A total of 528 shallow groundwater well (well depth less than the depth from surface to Eau Claire shale) nitrate samples collected in 2006 were used for our study. We performed kriging analysis to interpolate the spatial distribution of the probability of groundwater nitrate concentration greater than 3 mgper liter. We mapped the interpolation results at a 30m spatial resolution using Geostatistical Analyst extension in ArcGIS (ESRI). In this map, areas with lower probability values provided more groundwater quality service, and vice versa. Surface water quality (annual phosphorus loading, kg per hectare). We adapted a spatially explicit, scenario-driven modeling tool, Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) to simulate discharge of nonpoint-source phosphorus. A grid cells phosphorus contribution was quantified as a function of water yield index, land use/cover, export coefficient, and downslope retention ability with the following equation:Expx = ALVx * sum of the products from y=x+1 to X for (1-Ey)where ALVx is the adjusted phosphorus export from pixel x , Ey is the filtration efficiency of each downstream pixel y , and X represents phosphorus transport route from where it originated to the downstream water bodies. Filtration efficiency was assigned by cover type: natural vegetation was assigned a high value, semi-natural vegetation an intermediate value, and developed or impervious covers were assigned low values. We ran the model for 2006 and mapped estimated phosphorus loading across the watershed. The ecosystem service of providing high quality surface water was the inverse of phosphorus loading. Therefore, areas with lower phosphorus loading values delivered more surface water quality, and areas with higher phosphorus loading values supplied less surface water quality.Soil retention (annual sediment yield, metric tons per hectare). We quantified annual sediment yield as the (inverse) indicator for soil retention by using the Modified Universal Soil Loss Equation (MUSLE). MUSLE is a storm event based model that estimates sediment yield as a function of runoff factor, soil erodibility, geomorphology, land use/cover and land management. Specifically, a grid cells contribution of sediment for a given storm event is calculated as:Sed= 11.8*(Q*q_p)^0.56 * K * LS * C * Pwhere Sed represents the amount of sediment that is transported downstream network (metric tons), Q is the surface runoff volume (m³), q_pis the peak flow rate (cubic meters per s), K is soil erodibility which is based on organic matter content, soil texture, permeability and profiles, LS is combined slope and steepness factor, and C* P is the product of plant cover and its associated management practice factor. We used the ArcSWAT interface of the Soil and Water Assessment Tool (SWAT) to perform all the simulations. We ran this model at a daily time step from 2004 to 2006, with the first two years as spin up , then mapped total sediment yield for 2006 across the watershed. Similar to surface water quality, the ecosystem service of soil retention was the inverse of sediment yield. In this map, areas with lower sediment yield provided more of this service, and areas with higher sediment yield delivered less. Flood regulation (flooding regulation capacity, unitless, 0 to 100) We used the capacity assessment approach to quantify the flood regulation service based on four hydrological parameters: interception, infiltration, surface runoff and peak flow. We first applied the Kinematic Runoff and Erosion (KINEROS) model to derive estimates of three parameters (infiltration, surface runoff and peak flow) for six sampled sub basins in this watershed. KINEROS is an event-oriented, physically based, distribution model that simulates interception, infiltration, surface runoff and erosion at sub-basin scales. In each simulation, a sub basin was first divided into smaller hydrological units. For the given pre-defined storm event, the model then calculated the amount of infiltration, surface runoff and peak flow for each unit. Second, we classified these estimates into 10 discrete capacity classes with range from 0 to 10 (0 indicates no capacity and 10 indicates the highest capacity) and united units with the same capacity values and overlaid with land cover map. Third, we calculated the distribution of all land use/cover classes within every spatial unit (with a particular capacity). We then assigned each land use/cover a capacity parameter based on its dominance (in percentage) within all capacity classes. As a result, every land use/cover was assigned a 0 to 10 capacity value for infiltration, surface runoff and peak flow. This procedure was repeated for six sub basins, and derived capacity values were averaged by cover type. We applied the same procedure to soil data and derived averaged capacity values for each soil type with the same set of three parameters. In addition, we obtained interceptions from published studies for each land use/cover and standardized to the same 0 to 10 range. Finally, the flood regulation capacity (FRC) for each 30m cell was calculated with the equation below:FRC= for each land use and land cover class the sum of (interception + infilitration + runoff + peakflow) + for each soil class the sum of (infiltration + runoff + peakflow).To simplify interpretation, we rescaled original flood regulation capacity values to a range of 0 to100, with 0 representing the lowest regulation capacity and 100 the highest. Forest recreation (recreation score, unitless, 0 to 100). We quantified the forest recreation service as a function of the amount of forest habitat, recreational opportunities provided, proximity to population center, and accessibility of the area for each 30m grid cell with the equation below:FRSi= Ai * sum of (Oppti + Popi + Roadi)where FRS is forest recreation score, A is the area of forest habitat, Oppt represents the recreation opportunities, Pop is the proximity to population centers, and Road stands for the distance to major roads. To simplify interpretation, we rescaled the original forest recreation score (ranging from 0 to 5200) to a range of 0 to 100, with 0 representing no forest recreation service and 100 representing highest service. Several assumptions were made for this assessment approach. Larger areas and places with more recreational opportunities would provide more recreational service, areas near large population centers would be visited and used more than remote areas, and proximity to major roads would increase access and thus recreational use of an area. Hunting recreation (recreation score, unitless 0 to100) We applied the same procedure used for forest recreation to quantify hunting service. Due to limited access to information regarding private land used for hunting, we only included public lands, mainly state parks, for this assessment. The hunting recreation service was estimated as a function of the extent of wildlife areas open for hunting, the number of game species, proximity to population center, and accessibility for each 30m grid cell with the following equation:<br />HRSi= Ai * sum of (Spei + Popi + Roadi)where HRS is hunting recreation score, A is the area of public wild areas open for hunting and fishing, Spe represents the number of game species, Pop stands for the proximity to population centers, and Road is the distance to major roads. To simplify interpretation, we rescaled the original hunting recreation score (ranging from 0 to 28000) to a range of 0 to100, with 0 representing no hunting recreation service and 100 representing highest service. Similar assumptions were made for this assessment. Larger areas and places with more game species would support more hunting, and areas closer to large population centers would be used more than remote areas. Finally, proximity to major roads would increase access and use of an area.

Methods

Data were collected from two lakes, Sparkling Lake (46.008, -89.701) and Peter Lake (46.253, -89.504), both located in the northern highlands Lake District of Wisconsin and the Upper Peninsula of Michigan over a 10 day period on each lake in July and August of 2007. Refer to Van de Bogert et al. 2011 for limnological characteristics of the study lakes.Measurements of dissolved oxygen and temperature were made every 10 minutes using multiple sondes dispersed horizontally throughout the mixed-layer in the two lakes (n=35 sondes for Sparkling Lake and n=27 sondes for Peter Lake). Dissolved oxygen measurements were corrected for drift.Conductivity, temperature compensated specific conductivity, pH, and oxidation reduction potential were also measured by a subset of sensors in each lake. Of the 35 sondes in Sparkling Lake, 31 were from YSI Incorporated: 15 of model 600XLM, 14 of model 6920, and 2 of model 6600). The remaining sondes placed in Sparkling Lake were 4 D-Opto sensors, Zebra-Tech, LTD. In Peter Lake, 14 YSI model 6920 and 13 YSI model 600XLM sondes were used.Sampling locations were stratified randomly so that a variety of water depths were represented, however, a higher density of sensors were placed in the littoral rather than pelagic zone. See Van de Bogert et al. 2012 for the thermal (stratification) profile of Sparkling Lake and Peter Lake during the period of observation, and for details on how locations were classified as littoral or pelagic. In Sparkling Lake, 11 sensors were placed within the shallowest zone, 12 in the off-shore littoral, and 6 in each of the remaining two zones, for a total of 23 littoral and 12 pelagic sensors. Similarly, 15 sensors were placed in the two littoral zones, and 12 sensors in the pelagic zone.Sensors were randomly assigned locations within each of the zones using rasterized bathymetric maps of the lakes and a random number generator in Matlab. Within each lake, one pelagic sensor was placed at the deep hole which is used for routine-long term sampling.Note that in Sparkling Lake this corresponds to the location of the long-term monitoring buoy. After locations were determined, sensors were randomly assigned to each location with the exception of the four D-Opto sensor is Sparkling Lake, which are a part of larger monitoring buoys used in the NTL-LTER program. One of these was located near the deep hole of the lake while the other three were assigned to random locations along the north shore, south shore and pelagic regions of the lake. Documentation: Van de Bogert, M.C., Bade, D.L., Carpenter, S.R., Cole, J.J., Pace, M.L., Hanson, P.C., Langman, O.C., 2012. Spatial heterogeneity strongly affects estimates of ecosystem metabolism in two north temperate lakes. Limnology and Oceanography 57, 1689-1700.

Methods

I. Field MethodsThe site where this manipulative field experiment was conducted on the Kalfastrond peninsula at Lake Myvatn is approximately 150 meters long and 75 meters wide. The vegetation consists of grasses Deschampsia spp., Poa spp., and Agrostis spp.), sedges (Carex spp.), and forbs (Ranunculus acris, Geum rivale,and Potentilla palustris). The experimental midge exclusions occurred from the middle or end of May to the beginning of August in each year, the period corresponding to the active growing season of plants and the flight activity of midges at this site. 2 by 2 meter plots were established in 6 blocks spread across the site. Each block included 3 treatment levels, an open control plot, a full exclusion cage and a partial exclusion cage, for a total of 18 experimental plots. Control plots were open to allow complete access to all arthropods. Experimental midge exclusion cages were 1 meter high and constructed from white PVC tubing affixed to rebar posts on each corner of the plot, Plate 1. Full exclusion cages were entirely covered with white polyester netting, 200 holes per square inch, Barre Army Navy Store, Barre VT, USA, to prevent midges from entering the plot. The mesh netting completely enclosed the 2 by 2 by 1 meter frame to prevent flying insects from entering, however the mesh was not secured to the ground in order to allow non flying,ground crawling, arthropods to freely enter and exit the cages. Partial exclusion cages had one 0.5 meter strip of mesh stretched around the outside of the frame and another 0.75 meter strip draped over the top. Partial cages were designed to examine any effects of cages themselves on terrestrial responses while minimally affecting midge inputs into the plots and arthropod movement.The partial exclusion treatment was discontinued in 2011. Each plot contains a pitfall and an infall trap that are continuously sampled during the summer, while the cages are up. Vacuum samples were taken from the plots about once per month in 2008 through 2010 and only once per summer for subsequent summers.Midge activity was measured in all plots to document changes in midge abundance over the course of a season and between years and to assess the degree to which cages excluded midges. Midge abundance in the plots was continuously measured using passive aerial infall traps consisting of a 1000 milliliter clear plastic cup, 95 square centimeter opening, attached to a post 0.5 meters high and filled with 250 milliliters of a 1 to 1 ethylene glycol to water solution and a small amount of unscented detergent to capture and kill insects that alighted upon the surface. Infall traps were emptied about every 10 days.II. AnalysisMidges were counted and identified to morphospecies, small and large. The midge (Diptera,Chrionomidae) assemblage at Myvatn is dominated by two species,Chironomus islandicus (Kieffer)(large, 1.1 mg dw) and Tanytarsus gracilentus(Holmgren)(small, 0.1 mg dw), together comprising 90 percent of total midge abundance (Lindegaard and Jonasson 1979). First, the midges collected in the infall traps were spread out in trays, and counted if there were only a few. Some midges were only identified to the family level of Simuliidae,and other arthropods were counted and categorized as the group, others. Arthropods only identified to the family level Simuliidae or classified as others were not dually counted as Chironomus islandicus or Tanytarsus gracilentus. If there were many midges, generally if there were hundreds to thousands, in an infall trap,subsamples were taken. Subsampling was done using plastic rings that were dropped into the tray. The rings were relatively small compared to the tray, about 2 percent of the area of a tray was represented in a ring. The area inside a ring and the total area of the trays were also measured. Note that different sized rings and trays were used in subsample analysis. These are as follows, Trays, small (area of 731 square centimeters), large1 (area of 1862.40 square centimeters), and large2 (area of 1247 square centimeters). Rings, standard ring (diameter of 7.30 centimeters, subsample area is 41.85 square centimeters) and small ring (diameter of 6.5 centimeters, subsample area is 33.18 square centimeters). A small ring was only used to subsample trays classified as type large2.The fraction subsampled was then calculated depending on the size of the tray and ring used for the subsample analysis. If the entire tray was counted and no subsampling was done then the fraction subsampled was assigned a value of 1.0. If subsampling was done the fraction subsampled was calculated as the number of subsamples taken multiplied by the fraction of the tray that a subsample ring area covers (number of subsamples multiplied by (ring area divided by tray area)). Note that this is dependent on the tray and ring used for subsample analysis. Finally, the number of midges in an infall trap accounting for subsampling was calculated as the raw count of midges divided by the fraction subsampled (raw count divided by fraction subsampled).Other metrics such as total insects in meters squared per day, and total insect biomass in grams per meter squared day can be calculated with these data. in addition to the estimated average individual midge masses in grams, For 2008 through 2010 average midge masses were calculated as, Tanytarsus equal to .0001104 grams, Chironomus equal to .0010837 grams. For 2011 average midge masses were, Tanytarsus equal to .000182 grams, Chironomus equal to .001268 grams.

Methods

DL is used to represent “detection limit” where known.NOTE (1): Each method listed below corresponds with a USGS Parameter Code, which is listed after the variable name. NOTE (2): If the NEMI method # is known, it is also specified at the end of each method description.NOTE (3): Some of the variables are calculated using algorithms within QWDATA. If this is the case see Appendix D of the NWIS User’s Manual for additional information. However, appendix D does not list the algorithm used by the USGS. If a variable is calculated with an algorithm the term: algor, will be listed after the variable name.anc: 99431, Alkalinity is determined in the field by using the gran function plot methods, see TWRI Book 9 Chapter A6.1. anc_1: 90410 and 00410, Alkalinity is determined by titrating the water sample with a standard solution of a strong acid. The end point of the titration is selected as pH 4.5. See USGS TWRI 5-A1/1989, p 57, NEMI method #: I-2030-89.2. c13_c12_ratio: 82081, Exact method unknown. The following method is suspected: Automated dual inlet isotope ratio analysis with sample preparation by precipitation with ammoniacal strontium chloride solution, filtration, purification, acidified of strontium carbonate; sample size is greater than 25 micromoles of carbon; one-sigma uncertainty is approximately ± 0.1 ‰. See USGS Determination of the delta13 C of Dissolved Inorganic Carbon in Water, RSIL Lab Code 1710. Chapter 18 of Section C, Stable Isotope-Ratio Methods Book 10, Methods of the Reston Stable Isotope Laboratory.3. ca, mg, mn, na, and sr all share the same method. The USGS parameter codes are listed first, then the method description with NEMI method #, and finally DL’s:ca- 00915, mg- 00925, mn- 01056, na- 00930, sr- 01080All metals are determined simultaneously on a single sample by a direct reading emission spectrometric method using an inductively coupled argon plasma as an excitation source. Samples are pumped into a crossflow pneumatic nebulizer, and introduced into the plasma through a spray chamber and torch assembly. Each analysis is determined on the basis of the average of three replicate integrations, each of which is background corrected by a spectrum shifting technique except for lithium (670.7 nm) and sodium (589.0 nm). A series of five mixed-element standards and a blank are used for calibration. Method requires an autosampler and emission spectrometry system. See USGS OF 93-125, p 101, NEMI Method #: I-1472-87.DL’s: ca- .02 mg/l, mg-.01 mg/l, mn-1.0 ug/l, na- .2 mg/l, sr- .5 ug/l4. cl, f, and so4 all share the same method. The USGS parameter codes are listed first, then the method description with NEMI method #, and finally DL’s:cl- 00940, f-00950, so4-00945All three anions (chloride, flouride, and sulfate) are separated chromatographically following a single sample injection on an ion exchange column. Ions are separated on the basis of their affinity for the exchange sites of the resin. The separated anions in their acid form are measured using an electrical conductivity cell. Anions are identified on the basis of their retention times compared with known standards. 19 The peak height or area is measured and compared with an analytical curve generated from known standards to quantify the results. See USGS OF 93-125, p 19, NEMI method #: I-2057.DL’s: cl-.2 mg/l, f-.1 mg/l, so4-.2 mg/lco2: 00405, algor, see NWIS User's Manual, QW System, Appendix D, Page 285.co3: 00445, algor.color: 00080, The color of the water is compared to that of the colored glass disks that have been calibrated to correspond to the platinum-cobalt scale of Hazen (1892), See USGS TWRI 5-A1 or1989, P.191, NEMI Method #: I-1250. DL: 1 Pt-Co colorconductance_field: 00094 and 00095, specific conductance is determined in the field using a standard YSI multimeter, See USGS TWRI 9, 6.3.3.A, P. 13, NEMI method #: NFM 6.3.3.A-SW.conductance_lab: 90095, specific conductance is determined by using a wheat and one bridge in which a variable resistance is adjusted so that it is equal to the resistance of the unknown solution between platinized electrodes of a standardized conductivity cell, sample at 25 degrees celcius, See USGS TWRI 5-A1/1989, p 461, NEMI method #: I-1780-85.dic: 00691, This test method can be used to make independent measurements of IC and TC and can also determine TOC as the difference of TC and IC. The basic steps of the procedure are as follows:(1) Removal of IC, if desired, by vacuum degassing;(2) Conversion of remaining inorganic carbon to CO₂ by action of acid in both channels and oxidation of total carbon to CO₂ by action of ultraviolet (UV) radiation in the TC channel. For further information, See ASTM Standards, NEMI method #: D6317. DL: n/adkn: 00623 and 99894, Organic nitrogen compounds are reduced to the ammonium ion by digestion with sulfuric acid in the presence of mercuric sulfate, which acts as a catalyst, and potassium sulfate. The ammonium ion produced by this digestion, as well as the ammonium ion originally present, is determined by reaction with sodium salicylate, sodium nitroprusside, and sodium hypochlorite in an alkaline medium. The resulting color is directly proportional to the concentration of ammonia present, see USGS TWRI 5-A1/1989, p 327, NEMI method #: 351.2. DL: .10 mg/Ldo: 0300, Dissolved oxygen is measured in the field with a standard YSI multimeter, NEMI Method #: NFM 6.2.1-Lum. DL: 1 mg/L.doc: 00681, The sample is acidified, purged to remove carbonates and bicarbonates, and the organic carbon is oxidized to carbon dioxide with persulfate, in the presence of an ultraviolet light. The carbon dioxide is measured by nondispersive infrared spectrometry, see USGS OF 92-480, NEMI Method #: O-1122-92. DL: .10 mg/L.don: 00607, algor, see NWIS User's Manual, QW System, Appendix D, page 291.dp: 00666 and 99893, All forms of phosphorus, including organic phosphorus, are converted to orthophosphate ions using reagents and reaction parameters identical to those used in the block digester procedure for determination of organic nitrogen plus ammonia, that is, sulfuric acid, potassium sulfate, and mercury (II) at a temperature of 370 deg, see USGS OF Report 92-146, or USGS TWRI 5-A1/1979, p 453, NEMI method #: I-2610-91. DL= .012 mg/L.fe: 01046, Iron is determined by atomic absorption spectrometry by direct aspiration of the sample solution into an air-acetylene flame, see USGS TWRI 5-A1/1985, NEMI method #: I-1381. DL= 10µg/L.h_ion: 00191, algor.h20_hardness: 00900, algor.h20_hardness_2: 00902, algor.hco3: 00440, algor.k: 00935, Potassium is determined by atomic absorption spectrometry by direct aspiration of the sample solution into an air-acetylene flame , see USGS TWRI 5-A1/1989, p 393, NEMI method #: I-1630-85. DL= .01 mg/L.n_mixed: 00600, algor.n_mixed_1: 00602, algor.n_mixed_2: 71887, algor.nh3_nh4: 00608, Ammonia reacts with salicylate and hypochlorite ions in the presence of ferricyanide ions to form the salicylic acid analog of indophenol blue (Reardon and others, 1966; Patton and Crouch, 1977; Harfmann and Crouch, 1989). The resulting color is directly proportional to the concentration of ammonia present, See USGS OF 93-125, p 125/1986 (mg/l as N), NEMI Method #: I-2525. DL= .01 mg/L.nh3_nh4_1: 71846, algor.nh3_nh4_2: 00610, same method as 00608, except see USGS TWRI 5-A1/1989, p 321. DL = .01 mg/L.nh3_nh4_3: 71845, algor.no2: 00613, Nitrite ion reacts with sulfanilamide under acidic conditions to form a diazo compound which then couples with N-1-naphthylethylenediamine dihydrochloride to form a red compound, the absorbance of which is measured colorimetrically, see USGS TWRI 5-A1/1989, p 343, NEMI method #: I-2540-90. DL= .01 mg/L.no2_2: 71856, algor.no3: 00618, Nitrate is determined sequentially with six other anions by ion-exchange chromatography, see USGS TWRI 5-A1/1989, P. 339, NEMI method #: I-2057. DL= .05 mg/L.no3_2: 71851, algor.no32: 00630, An acidified sodium chloride extraction procedure is used to extract nitrate and nitrite from samples of bottom material for this determination(Jackson, 1958). Nitrate is reduced to nitrite by cadmium metal. Imidazole is used to buffer the analytical stream. The sample stream then is treated with sulfanilamide to yield a diazo compound, which couples with N-lnaphthylethylenediamine dihydrochloride to form an azo dye, the absorbance of which is measured colorimetrically. Procedure is used to extract nitrate and nitrite from bottom material for this determination (Jackson, 1958), see USGS TWRI 5-A1/1989, p 351. DL= .1 mg/Lno32_2: 00631, same as description for no32, except see USGS OF 93-125, p 157. DL= .1 mg/L.o18_o16_ratio: 82085, Sample preparation by equilibration with carbon dioxide and automated analysis; sample size is 0.1 to 2.0 milliliters of water. For 2-mL samples, the 2-sigma uncertainties of oxygen isotopic measurement results are 0.2 ‰. This means that if the same sample were resubmitted for isotopic analysis, the newly measured value would lie within the uncertainty bounds 95 percent of the time. Water is extracted from soils and plants by distillation with toluene; recommended sample size is 1-5 ml water per analysis, see USGS Determination of the Determination of the delta (18 O or 16O) of Water, RSIL Lab Code 489.o2sat: Dissolved oxygen is measured in the field with a standard YSI multimeter, which also measures % oxygen saturation, NEMI Method #: NFM 6.2.1-Lum.ph_field: 00400, pH determined in situ, using a standard YSI multimeter, see USGS Techniques of Water-Resources Investigations, book 9, Chaps. A1-A9, Chap. A6.4 "pH," NEMI method # NFM 6.4.3.A-SW. DL= .01 pH.ph_lab: 00403, involves use of laboratory pH meter, see USGS TWRI 5-A1/1989, p 363, NEMI method #: I-1586.po4: 00660, algor, see NWIS User's Manual, QW System, Appendix D, Page 286.po4_2: 00671, see USGS TWRI 5-A1/1989, NEMI method #: I-2602. DL= .01 mg/L.s: 63719, cannot determine exact method used. USGS method code: 7704-34-9 is typically used to measure sulfur as a percentage, with an DL =.01 µg/L. It is known that the units for sulfur measurements in this data set are micrograms per liter.sar: 00931, algor, see NWIS User's Manual, QW System, Appendix D, Page 288.si: 00955, Silica reacts with molybdate reagent in acid media to form a yellow silicomolybdate complex. This complex is reduced by ascorbic acid to form the molybdate blue color. The silicomolybdate complex may form either as an alpha or beta polymorph or as a mixture of both. Because the two polymorphic forms have absorbance maxima at different wavelengths, the pH of the mixture is kept below 2.5, a condition that favors formation of the beta polymorph (Govett, 1961; Mullen and Riley, 1955; Strickland, 1952), see USGS TWRI 5-A1/1989, p 417, NEMI method #: I-2700-85. DL= .10 mg/L.spc: 00932, algor, see NWIS User's Manual, QW System, Appendix D, Page 289.tds: 70300 and 70301, A well-mixed sample is filtered through a standard glass fiber filter. The filtrate is evaporated and dried to constant weight at 180 deg C, see " Filterable Residue by Drying Oven," NEMI method #: 160.1, DL= 10 mg/l. Note: despite DL values occur in the data set that are less than 10 mg/l.tds_1: 70301, algor, see NWIS User's Manual, QW System, Appendix D, Page 289.tds_2: 70303, algor, see NWIS User's Manual, QW System, Appendix D, Page 290.tkn: 00625 and 99892, Block digester procedure for determination of organic nitrogen plus ammonia, that is, sulfuric acid, potassium sulfate, and Mercury (II) at a temperature of 370°C. See the USGS Open File Report 92-146 for further details. DL: .10 mg/L.toc: 00680, The sample is acidified, purged to remove carbonates and bicarbonates, and the organic carbon is oxidized to carbon dioxide with persulfate, in the presence of an ultraviolet light. The carbon dioxide is measured by nondispersive infrared spectrometry, see USGS TWRI 5-A3/1987, p 15, NEMI Method #: O-1122-92. DL=.10 mg/L.ton: 00605, algor, See NWIS User's Manual, QW System, Appendix D, page 286.tp: 00665 and 99891, This method may be used to analyze most water, wastewater, brines, and water-suspended sediment containing from 0.01 to 1.0 mg/L of phosphorus. Samples containing greater concentrations need to be diluted, see USGS TWRI 5-A1/1989, p 367, NEMI method #: I-4607. tp_2: 71886, algor.tpc: 00694, The basic steps of this test method are:1) Conversion of remaining IC to CO2 by action of acid, 2) Removal of IC, if desired, by vacuum degassing, 3) Split of flow into two streams to provide for separate IC and TC measurements, 4) Oxidation of TC to CO2 by action of acid-persulfate aided by ultraviolet (UV) radiation in the TC channel, 5) Detection of CO2 by passing each liquid stream over membranes that allow the specific passage of CO2 to high-purity water where change in conductivity is measured, and 6) Conversion of the conductivity detector signal to a display of carbon concentration in parts per million (ppm = mg/L) or parts per billion (ppb = ug/L). The IC channel reading is subtracted from the TC channel reading to give a TOC reading, see ASTM Standards, NEMI Method #: D5997. DL= .06 µg/L.tpn: 49570, A weighed amount of dried particulate (from water) or sediment is combusted at a high temperature using an elemental analyzer. The combustion products are passed over a copper reduction tube to covert nitrogen oxides to molecular nitrogen. Carbon dioxide, nitrogen, and water vapor are mixed at a known volume, temperature, and pressure. The concentrations of nitrogen and carbon are determined using a series of thermal conductivity detectors/traps, measuring in turn by difference hydrogen (as water vapor), carbon (as carbon dioxide), and nitrogen (as molecular nitrogen). Procedures also are provided to differentiate between organic and inorganic carbon, if desired, see USEPA Method 440, NEMI method #: 440. DL= .01 mg/L.