The dataset provides growth and toxicokinetic parameters for the simulations presented in the manuscript by Pollesch et al. Developing Integral Projection Models for Aquatic Ecotoxicology. This dataset is associated with the following publication: Pollesch, N., K. Flynn, S. Kadlec, J. Swintek, S. Raimondo, and M. Etterson. Developing integral projection models for ecotoxicology. ECOLOGICAL MODELLING. Elsevier Science BV, Amsterdam, NETHERLANDS, 464: 109813, (2022).

These are the data associated with each of the figures in the publication. This dataset is associated with the following publication: Thursby, G., K. Sappington, and M. Etterson. Coupling Toxicokinetic-Toxicodynamic and Population Models for Assessing Aquatic Ecological Risks to Time-Varying Pesticide Exposures. ENVIRONMENTAL TOXICOLOGY AND CHEMISTRY. Society of Environmental Toxicology and Chemistry, Pensacola, FL, USA, 37(10): 2633-2644, (2018).

These are the data associated with each of the figures in the publication. This dataset is associated with the following publication: Thursby, G., K. Sappington, and M. Etterson. Coupling Toxicokinetic-Toxicodynamic and Population Models for Assessing Aquatic Ecological Risks to Time-Varying Pesticide Exposures. ENVIRONMENTAL TOXICOLOGY AND CHEMISTRY. Society of Environmental Toxicology and Chemistry, Pensacola, FL, USA, 37(10): 2633-2644, (2018).

Data sources for Pesticide in Water Calculator scenario simulations. This dataset is associated with the following publication: Awkerman, J., T. Purucker, S. Raimondo, and L. Oliver. Long-term, landscape-level assessment of aquatic pesticide exposure to identify amphibian ontological traits affecting vulnerability. Integrated Environmental Assessment and Management. Allen Press, Inc., Lawrence, KS, USA, 20(5): 1667-1676, (2024).

This dataset shows the effects of of the insecticide imidacloprid and the fungicide propiconazole on larval development of the acorn barnacle Amphibalanus amphitrite experiments conducted in 2018 and 2019. The aim of this project was to apply standard ecotoxicology protocols to determine the effects of the insecticide imidacloprid and the fungicide propiconazole on larval development rate of the acorn barnacle Amphibalanus Amphitrite. Larval development bioassays(4-d exposures) were conducted using a fungicide and insecticide that have been detected in the Great Barrier Reef catchment area (O'Brien et al., 2016). These toxicity data will enable improved assessment of the risks posed by pesticides to marine crustaceans for both regulatory purposes and for comparison with other taxa. Methods: Pesticide stock solutions were prepared using PESTANAL (Merck) analytical grade products (purity greater than or equal to 98%): imidacloprid (CAS 138261-41-3) and propiconazole (CAS 60207-90-1). This selection was based on application rates and detection in coastal waters of the GBR (O’Brien et al., 2016; Grant 2017). Pesticide stock solutions (100 – 1,000 mg L-1) were prepared by dissolving aliquots of the pure compounds in ultrapure water using clean, acid-washed (5% nitric acid) glass screw-top containers. Acetone was used to dissolve the imidacloprid and propiconazole (less than or equal to 0.01 % (v/v) in exposure solutions). Stock solutions were stored refrigerated and in the dark. Broodstock barnacles had been grown for several generations in the AIMS-NT aquaria facility (originally sourced from Darwin Harbour – 12°26'57.48"S, 130°51'7.51"E). Broodstock were fed freshly hatched brine shrimp (Artemia salina) and live rotifers daily. Broodstock were spawned as previously described (van Dam et al., 2016) and nauplii collected. Tests were conducted as previously described (van Dam et al., 2016). Nauplii were exposed in a custom-designed experimental test system that allowed for constant movement of the exposure media. The system consisted of a series of silanized glass funnels in which nauplii were exposed to increasing concentrations of imidacloprid or propiconazole and tested against control nauplii. Generally, a total of 24 funnels were used for 7 treatment concentrations and a control group, thus allowing for 3 replicate funnels per treatment. Each treatment vessel contained 100 mL exposure media, 50 newly released stage II nauplii and 1 x 107 cells of rinsed Chaetoceros muelleri. Every 24 h, 1 x 107 cells of rinsed C. muelleri were added to each funnel. After 96 h exposure, funnel contents were drained over a 150 µm nitrile mesh. The mesh was examined under a stereomicroscope and the number of cyprids and settled larvae scored. Quality control criteria (> 70% survival in control group) for test acceptability were met for each test used to derive toxicity estimates. Treatment effects were quantified by the percentage successful transition to cyprid in treatment groups relative to controls. Following prescribed statistical procedures (OECD 2006) the R package DRC (R-project 2015, Ritz & Streibig 2005), was used to model the test data and calculate toxicity estimates. Regression models evaluated included log-logistic and Weibull models of different levels of parametrisation. Model comparisons were conducted using the Akaike Information Criterion (AIC) and models that best described the data were applied to approximate pesticide concentrations eliciting 10 and 50% inhibition of successful transition relative to control animals (EC10 and EC50, respectively). The associated 95% confidence limits were estimated using the delta method. Format: The dataset is summarised in one file named ‘Amphibalanus amphitrite pesticide toxicity data_eAtlas.xlsx’ Data Dictionary: The excel spreadsheet has one tab for each pesticide. The last tab of the dataset shows the measured (start and end of test) water quality (WQ) parameters (pH, salinity, dissolved

This dataset shows the effects of herbicides and one fungicide (detected in Great Barrier Reef catchments) on the specific growth rates (from cell density data) of the microalgae Tisochrysis lutea during laboratory experiments conducted from 2018-2019. The aim of this project was to apply standard ecotoxicology protocols to determine the effects of Photosystem II (PSII), alternative herbicides and one fungicide on the growth of the marine microalgae Tisochrysis lutea. Growth bioassays were performed over 3-day exposures using pesticides that have been detected in the Great Barrier Reef catchment area (O'Brien et al., 2016). These toxicity data will enable improved assessment of the risks posed by PSII and alternative herbicides as well as the fungicide propiconazole to microalgae for both regulatory purposes and for comparison with other taxa. Methods: The haptophyte Tisochrysis lutea (formerly known as Isochrysis galbana)(Grant etal. 2017) (strain CS-177) was purchased from the Australian National Algae Supply Service, Hobart (CSIRO). Cultures of T. lutea were established in EDTA-free Guillard’s f/2 marine medium (Trenfield et al. 2015) (1 ml L-1 of f/2 medium in autoclaved natural seawater). Cultures were maintained in sterile 500 ml Erlenmeyer flasks as batch cultures in exponential growth phase with weekly aseptically transfers of 10 ml T. lutea suspension to 300 ml f/2 medium. Culture were maintained at 28 ± 1°C, 33 ± 1.5 psu and under a 12:12 h light:dark cycle (80 – 100 µmol photons m–2 s–1). Pesticide stock solutions were prepared using PESTANAL (Merck) analytical grade products (purity greater than or equal to 98%): diuron (CAS 330-54-1), metribuzin (CAS 21087-64-9), tebuthiuron (CAS 34014-18-1), bromacil (CAS 314-40-9), propazine (CAS 139-40-2), simazine (122-34-9), imazapic (CAS 104098-48-8), haloxyfop-p-methyl (CAS 72619-32-0), 2,4-D (CAS 94-75-7), MCPA (CAS 94-74-6), fluroxypyr (CAS 69377-81-7) and propiconazole (CAS 60207-90-1). The selection of pesticides was based on application rates and detection in coastal waters of the GBR (Grant et al. 2017, O’Brien et al. 2016). Pesticide stock solutions (100 – 1,000 mg L-1) were prepared by dissolving aliquots of the pure compounds in ultrapure water using clean, acid-washed (5% nitric acid) glass screw-top containers. Simazine, tebuthiuron and haloxyfop-p-methyl were dissolved using the carrier dimethyl sulfoxide (DMSO) (less than or equal to 0.02 % (v/v) in exposure solutions). Diuron, imazapic, metribuzin, bromacil, 2,4-D, propazine, MCPA, fluroxypyr and propiconazole were dissolved in acetone (less than or equal to 0.01 % (v/v) in exposure). Stock solutions were stored refrigerated and in the dark. Tests were conducted as previously described (Trenfield et al. 2015). Cultures of T. lutea were exposed to increasing concentrations of individual pesticides over a period of 72 h. Inoculum was taken from cultures in exponential growth phase (4-d old culture) and starting cell density assessed using a haemocytometer. For each treatment, a total volume of 250 mL test media was prepared in a clean, acid-washed 500 mL Schott bottle. Test media consisted of filtered (0.45 µm) seawater spiked with the respective pesticide stock, quarter strength EDTA-free f/2 media as nutrient source and T. lutea at a starting density of 3x103 or 1x104 cells mL-1. In each toxicity test, the response (specific growth rate of the culture) of the treatments exposed to pesticide were assessed against a seawater control group (no herbicide). For each test, 2 – 3 replicate 125 mL Erlenmeyer flasks (50 mL test volume) were assessed. Flasks were incubated at 27 – 29.0°C under a 12:12 h light:dark cycle (80 – 100 µmol photons m–2 s–1). After 72h, sub-samples (7 ml) were taken from each flask and cell densities measured using a flow cytometer (BD Accuri C6, BD Biosciences, CA, USA). Specific growth rates (SGR) were expressed as the logarithmic increase in cell density from day i (ti) to day j (tj) as per

This dataset shows the effects of the insecticide imidacloprid on larval development of the hermit crab Coenobita variabilis. Experiments were conducted in 2017. The aim of this project was to apply a standard ecotoxicology protocol to determine the effects of the insecticide imidacloprid (that has been detected in the Great Barrier Reef catchment area (O'Brien et al., 2016)), on larval development (6-d exposures) of the hermit crab Coenobita variabilis. These toxicity data will enable improved assessment of the risks posed to marine crustaceans for both regulatory purposes and for comparison with other taxa. Methods: Imidacloprid (CAS 138261-41-3) stock solutions were prepared using PESTANAL (Merck) analytical grade product (purity greater than or equal to 99.9%). Stock solutions (100 – 1,000 mg L-1) were prepared by dissolving aliquots of the pure compound in ultrapure water using clean, acid-washed (5% nitric acid) glass screw-top containers. Acetone was used to dissolve the imidacloprid (less than or equal to 0.01 % (v/v) in exposure solutions). Stock solutions were stored refrigerated and in the dark. Tests were conducted as previously described (in van Dam et al, 2018). Broodstock crabs were collected from the Nightcliff seashore (Darwin, Australia – 12°23'8.70"S, 130°50'34.59"E) and maintained in custom-built, flat-bottomed enclosures. Spawning was left to occur naturally and toxicity tests initiated immediately following collection of larvae. Transparent polystyrene cell culture plates (Nunc; Thermo Scientific) were employed as test chambers. Each replicate plate contained six wells with a volume of 13 mL each. Exposures were conducted in a high-precision environmental chamber maintained at 29 ± 1°C, under 80 – 100 µmol quanta m-2s-1 PAR irradiance and a 12h:12h diurnal light:dark cycle. Zoeae were exposed to increasing concentrations of imidacloprid and tested against control (no toxicant) larvae. Zoeae were allocated individually to a well as larvae became cannibalistic once transitioned to megalopae. Five wells within a discrete plate contained analogous treatment solutions. Per test, a total of 18 plates were employed for 5 treatment concentrations and a control group, allowing for 3 replicate plates per treatment. Ten mL of exposure media was added to individual wells before the tests were started by randomly placing a larva from the common pool into each well. Larvae were transferred every 48 h to fresh exposure solutions in clean plates. After 6 d exposure, tests were terminated and individuals scored under a stereo microscope. Quality control criteria (> 70% survival in control group) for test acceptability were met for each test. Treatment effects were quantified by the percentage successful transition to megalopae in treatment groups relative to controls. Following prescribed statistical procedures (OECD, 2006), the R package DRC (R project., 2015, Ritz and Stribig., 2005) was used to model the test data and calculate toxicity estimates. Regression models evaluated included log-logistic and Weibull models of different levels of parametrisation. Model comparisons were conducted using the Akaike Information Criterion (AIC) and models that best described the data were applied to approximate pesticide concentrations eliciting 10 and 50% inhibition of successful transition relative to control animals (EC10 and EC50, respectively). The associated 95% confidence limits were estimated using the delta method. Format: The dataset is summarised in one file named ‘Coenobita variabilis pesticide toxicity data_eAtlas.xlsx’ Data Dictionary: The excel spreadsheet has one tab for each pesticide. The last tab of the dataset shows the measured (start and end of test) water quality (WQ) parameters (pH, salinity, dissolved oxygen (DO), and temperature) for each test. For the ‘Imidacloprid_Development tab: Nominal (µg/L) = nominal herbicide concentrations used in the bioassays Measured (µg/L) = measured concentrations analysed by

This dataset shows the effects of herbicides (detected in the Great Barrier Reef catchments) on the growth rates (from stem length and biomass) on the stonewort Ceratophyllum demersum during laboratory experiments conducted in 2019. The aims of this project were to develop and apply standard ecotoxicology protocols to determine the effects of non-PSII herbicides on the growth of the stonewort Ceratophyllum demersum. Growth bioassays were performed over 7-day exposures using herbicides that have been detected in the Great Barrier Reef catchment area (O’Brien et al. 2016). This toxicity data will enable improved assessment of the risks posed by non-PSII herbicides to aquatic macrophytes for both regulatory purposes and for comparison with other taxa. Methods: The stonewort Ceratophyllum demersum was supplied by Watergarden Paradise Aquatic Nursery, Bass Hill, NSW. Cultures were maintained in 500 L outdoor plastic tanks in recirculating dechlorinated tap water, aerated and maintained at ambient outdoor temperature and lighting. Test replicates selected 48 h in advance and acclimated in dechlorinated tap water, 26 ± 2 °C, under a 12:12 hr light:dark cycle (102 ± 9 µmol photons m–2 s–1). Herbicide stock solutions were prepared using PESTANAL (Sigma-Aldrich) analytical grade products (HPLC less than or equal to 98%): haloxyfop-p-methyl (CAS 72619-32-0), imazapic (CAS 104098-48-8) and triclopyr (CAS 5535-06-3). The selection of herbicides was based on application rates and detection in coastal waters of the GBR (Grant et al. 2017, O’Brien et al. 2016). Stock solutions were prepared in 500 mL glass volumetric flasks using milli-Q water. Haloxyfop-p-methyl was dissolved using analytical grade acetone (< 0.01% (v/v) in exposures). Imazapic was dissolved in methanol (less than or equal to 0.01% (v/v) in exposure). No solvent carrier was used for the preparation of triclopyr. Cultures of Ceratophyllum demersum were exposed to a range of herbicide concentrations over a period of 7 days. Plants were sourced from actively growing cultures free of overt disease or deformity. Individual plants approximately 35 mm long with 5 whorls and an apical tip were added to 150 mL of each herbicide solution concentration and control treatment. In each toxicity test, a control (no herbicide) and solvent control (if used) treatments were added to support the validity of the test protocols and to monitor continued performance of the assays. Experiments were conducted in autoclaved, recirculating dechlorinated tap water. Five replicates of each treatment solution and control were prepared and incubated at 26.6 ± 0.5 °C under a 12:12 h light:dark cycle (90 ± 6 µmol photons m–2 s–1). Each replicate treatment was photographed at a standard height to measure stem length at Day 0 and Day 7. Biomass of a representative numbers of fronds were weighed at Day 0 to 3 significant figures using an analytical balance after blotting for 15 seconds to remove excess moisture. Fronds from each treatment replicate were weighed at Day 7 using the same technique. Specific growth rates (SGR) were expressed as the logarithmic increase in stem length (mm) or biomass (g) from day i (ti) to day j (tj) as per equation (1), where SGRi-j is the specific growth rate from time i to j; Xj is the stem length or biomass at day j and Xi is the stem length or biomass at day i (OECD 2006). SGRi-j =[(ln Xj - ln X i) / (tj - ti)] (day-1) (1) SGR relative to the control / solvent control treatment was used to derive chronic effect values for growth inhibition. A test was considered valid, if the SGR for frond number or surface area of control replicates was greater than or equal to 0.0495 day-1 determined from (OECD 2006 and Riethmuller et al 2003). Physical and chemical characteristics of each treatment were measured at 0 and 7 days including pH, electrical conductivity and temperature. Temperature was also logged in 15-min intervals over the total test duration. Analytical samples were taken at

This dataset shows the effects of herbicides (detected in the Great Barrier Reef catchments) on growth rates (from surface area and biomass) and photosynthesis (effective quantum yield) on the aquatic fern Azolla pinnata during laboratory experiments conducted in 2019. The aims of this project were to develop and apply standard ecotoxicology protocols to determine the effects of Photosystem II (PSII) and alternative herbicides on the growth and photosynthetic efficiency of the aquatic fern Azolla pinnata. Growth bioassays were performed over 14-day exposures using herbicides that have been detected in the Great Barrier Reef catchment area (O’Brien et al. 2016). Chronic effects of herbicides on the photophysiology of A. pinnata, measured by chlorophyll fluorescence as the effective quantum yield (Delta F/Fm’) were investigated using PAM fluorometry after 14-day herbicide exposure. These toxicity data will enable improved assessment of the risks posed by PSII and alternative herbicides to aquatic macrophytes for both regulatory purposes and for comparison with other taxa. Methods: The aquatic fern Azolla pinnata was sourced from Watergarden Paradise Nursery, NSW. Cultures were established in IRRI2 medium (Pereira & Carrapiço 2009). Cultures were maintained in 10 L tubs containing 3–5 L IRRI2 as batch cultures with weekly transfers to fresh medium. Clean culture solutions were maintained at 26 ± 1 °C, under a 12:12 hr light:dark cycle (65-77µmol photons m–2 s–1). Herbicide stock solutions were prepared using PESTANAL (Sigma-Aldrich) analytical grade products (HPLC greater than or equal to 98%): diuron (CAS 330-54-1), fluometuron (CAS 2164-17-2), fluroxypyr (CAS 69377-81-7), haloxyfop-p-methyl (CAS 72619-32-0), imazapic (CAS 104098-48-8), isoxaflutole (CAS 141112-29-0) and triclopyr (CAS 5535-06-3). The selection of herbicides was based on application rates and detection in coastal waters of the GBR (Grant et al. 2017, O’Brien et al. 2016). Stock solutions were prepared in 100 mL glass volumetric flasks using milli-Q water. Diuron, haloxyfop-p-methyl and isoxaflutole were dissolved using analytical grade acetone (< 0.01% (v/v) in exposures). Imazapic was dissolved in methanol (less than 0.01% (v/v) in exposure). No solvent carrier was used for the preparation of the remaining herbicide stock solutions. Cultures of A. pinnata were exposed to a range of herbicide concentrations over a period of 14 days. Fronds were selected from actively growing cultures free of overt disease or deformity. Four triplicate fronds each comprising eight ramets were added to 100 mL of each herbicide solution concentration and control treatment. In each toxicity test, control (no herbicide) and solvent control (if used) treatments were added to support the validity of the test protocols and to monitor continued performance of the assays. Experiments were conducted in IRRI2 medium (Pereira & Carrapiço 2009) with solutions replaced at Day 7. Three replicates of each treatment solution and control were prepared and incubated at 26.6 ± 0.5 °C under a 12:12 h light:dark cycle (90 ± 6 µmol photons m–2 s–1). Each replicate treatment was photographed at a standard height to estimate surface area at Day 0 and Day 14. Biomass of a representative numbers of fronds were weighed to 4 significant figures using an analytical balance after blotting for 15 seconds to remove excess moisture. Fronds from each treatment replicate were weighed at Day 14 using the same technique. Specific growth rates (SGR) were expressed as the logarithmic increase in surface area or biomass from day i (ti) to day j (tj) as per equation (1), where SGRi-j is the specific growth rate from time i to j; Xj is the surface area or biomass at day j and Xi is the surface area or biomass at day i (OECD 2006). SGR i-j = [(ln Xj - ln Xi )/(tj - ti )] (day-1) SGR relative to the control / solvent control treatment was used to derive chronic effect values for growth inhibition. A test was considered