A benthic bioindicator reveals distinct land and ocean–Based influences in an urbanized coastal embayment

Biogeochemical maps of coastal regions can be used to identify important influences and inputs that define nearshore environments and biota. Biogeochemical tracers can also track animal movement and their diet, monitor human coastal development, and evaluate the condition of habitats and species. However, the beneficial applications of spatial biogeochemical analysis are hindered by a limited understanding of how tracer distribution is affected by different land and ocean–based influences. To help address these knowledge gaps, we determined the spatial trends of three stable isotopes (δ13C-carbon, δ15N-nitrogen, δ34S-sulfur) and 13 major and trace elements in an urbanized coastal embayment (Moreton Bay, Australia), as incorporated into the muscle tissue of a marine consumer, the eastern king prawn Melicertus plebejus. Results were used to identify unique biochemical regions within the bay and to discuss how spatial patterns in tracers could be used to indicate the relative importance of catchment, urban and offshore drivers in coastal bays. Discriminant analysis identified seven biogeochemical regions that were likely distinguished by variation in catchment, urban, and offshore input, and habitat type. δ13C and δ15N patterns suggested nearshore areas could be distinguished by increased sediment resuspension and higher wastewater inputs from catchments. High inshore lead (Pb) and copper (Cu) concentrations were likely the result of urban input. Arsenic (As) and cadmium (Cd) increased further from shore. This trend implied oceanic influences were a significant control over As and Cd bioavailability. Cobalt (Co) and rare earths were also used to differentiate some nearshore areas, but incongruent distribution patterns in Co suggested it may be less reliable. Overall, results indicated that δ15N, δ13C, Cd, Cu, Pb and rare earth elements were the most reliable tracers to differentiate nearshore and offshore environments, and catchment–based effects. We encourage future studies to consider using a similar multivariate approach in coastal spatial analysis, and to include unrelated tracers that reflect distinct coastal influences.


Introduction
Coastal and estuarine ecosystems are diverse, productive environments that contain a wide variety of habitats and species [1][2][3]. Healthy coastal ecosystems provide critical services for human communities, including waste-water purification, erosion control, and support for coastal fisheries [4]. However, coastal areas are also highly dynamic and often heavily developed environments that experience substantial and sometimes abrupt spatial changes across a range of abiotic and biotic factors, such as offshore water exchange [5], or catchment and urban inputs [6,7]. The condition of coastal habitats and their capacity to support biodiversity is largely determined by spatial variation in these and other environmental parameters [8,9]. Consequently, mapping physical, biological, and biogeochemical variations in coastal environments is essential for effective monitoring, management, and conservation.
Coastal zones and communities have traditionally been defined by structural factors such as exposure, tidal range, sediment supply, and climate. Biogeochemistry offers a different perspective on the configuration of coastal regions that incorporates both structural and processbased information. Spatial patterns in biogeochemical tracers, such as stable isotopes or trace metals, can be used to evaluate fine and large-scale abiotic and biotic variation in aquatic ecosystems [10][11][12]. For example, coastal rivers that drain catchments with unique geological and anthropogenic influences will release distinctive mixtures of sediments, nutrients, and pollutants into coastal basins [12][13][14]. Biogeochemical maps of these tracers can then be used to identify these influences, define distinct biogeochemical regions, and evaluate the condition of nearshore habitats and species. The information provided by tracers and biogeochemical mapping also enables landscape-level evaluations of changes in coastal communities. For example, changes in tracer concentrations and distribution patterns over time can be used to monitor the effects of coastal development or the success of management programs. Biogeochemical maps therefore allow us to better define input areas and processes within coastal systems.
The beneficial applications of biogeochemical mapping and monitoring programs are limited by our relatively poor understanding of how tracers disperse through coastal habitats and the spatial relationships that may exist between some tracers. As a result, most biogeochemical studies rely on only one or a few tracers (e.g. stable isotopes or trace metals), and there is little work that shows how different tracers are simultaneously transferred through coastal habitats. Using a small number of tracers can also limit the scope and precision of a particular study. For example, coastal ecosystems often have a larger number of sources than isotope tracers, which can lead to inaccurate model outputs [15,16]. The isotope values of distinct food sources or habitats may also overlap, making it difficult to distinguish different populations or environments using biochemical methods [17][18][19]. In dynamic and diverse coastal systems, additional and disparate tracers can improve spatial analysis by further distinguishing different habitats, populations, and environmental influences at smaller and more ecologically relevant scales [10,20]. Therefore, work that clarifies spatial biogeochemical trends in coastal areas, and includes a range of potential tracers, will advance a wide range of research fields and applications.
Bioindicators offer a valuable and practical way to measure biogeochemical variation in coastal environments [21,22]. A bioindicator is any organism that has the capacity to accumulate unique biogeochemical "fingerprints" that integrate and reflect local tracer concentrations from multiple potential sources (water, sediment, and diet) [20]. These fingerprints not only reflect tracer exposure or availability, but also the local environmental conditions which affect tracer bioavailability and uptake [23,24]. Secondary generalist consumers may be particularly valuable bioindicators because their diet, and resultant trace element values, represent an integration of disparate biogeochemical sources [25]. For example, lead (Pb) measured in the muscle tissue of the opportunistic consumer crab Callinectes arcuatus was used to determine Pb sources in the Gulf of California (Mexico) [26]. Increased metal accumulation in greentail prawns Metapenaeuss bennettae closely corresponded to sediment concentrations in contaminated estuaries in New South Wales, Australia [27]. Thus consumer tissue values can provide unique and useful insight into the distribution of bioavailable tracers.
The primary goal of this study was to create exploratory biogeochemical maps of an urbanized, coastal embayment (Moreton Bay, Australia) using stable isotope, major, trace, and rare earth element values as determined in the muscle tissue of the benthic invertebrate species, the eastern king prawn Melicertus plebejus. The results were used to investigate different tracers as identifiers of key costal influences. Melicertus plebejus was chosen as a potentially useful bioindicator because it has relatively fast isotope turnover (1-2 weeks) [28][29][30], and are one of the few benthic generalists that can be found in almost all areas of Moreton Bay [31][32][33][34]. Although M. plebejus is mobile, its residency time in Moreton Bay is relatively long in comparison to its tissue turnover rate [35]. As a result, M. plebejus tracer concentrations should reflect regional trends in bioavailable tracers. Where applicable, we discuss the utility of M. plebejus as a bioindicator, and how M. plebejus tracer accumulation can be used to better understand tracer dispersal. Our research will contribute to a better understanding of tracer dispersal patterns in urbanised bays, how tracers can be used to identify key coastal influences, and increase the number of potential tracers for future applications.

Ethics statement
Specimens were collected in cooperation with the Long-Term Monitoring Program (Fisheries Queensland), which conducts annual fishery-independent recruitment monitoring surveys of M. plebejus in southeast Queensland. Because samples were collected under this government program, no additional permits were required. No ethics permits were required from the Griffith University Office of Research because the study species was an invertebrate. Specimens were euthanized quickly and humanly by freezing individuals at -20˚C in on-board freezers.

Study site
Moreton Bay is a large (1845 km 2 ), shallow (mean depth = 8m), semi-enclosed estuary located on the southeast coast of subtropical Queensland, Australia (Fig 1). The bay is bounded by two barrier islands (Moreton Island and Stradbroke Island) that limit exchange with oceanic waters. Moreton Bay is adjacent to the city of Brisbane which has a steadily growing current population of approximately 2.25 million people. The surrounding areas are dominated by rural communities and are primarily used for agriculture and peri-urban activities. The bay contains a variety of habitats including rocky and sandy beaches, seagrass, mudflats, mangroves, river deltas, and coral communities. The bay also receives periodic terrestrial run-off from four primary rivers (Caboolture, Pine, Brisbane, Logan). Sewage treatment plant (STP) outlets are located along each of these primary river systems. The Brisbane River is the largest coastal river in Moreton Bay and drains both agricultural areas in the upper catchment as well as the urban centre of Brisbane near the coast. Average flow is low at 5 m 3 s -1 , and water residence time in the lower river near the coast is high at 50-100 days [36]. Summer dominated rainfall in combination with large storms can lead to increased river flow and periodic flood events that discharge water and associated sediment and contaminants into the bay [37]. The Brisbane River contains the largest and most active industrial port in the bay, however smaller industrial and recreational marinas are found all along the Moreton coastline. Freshwater inputs can eventually be transported out of the bay through the north passage by a tidal current that initially brings oceanic water southward along the eastern margin of the bay, but then turns northward along the western shore. Salinity in the bay ranges from 24.7 to 37.6 g/L (mean = 34.7 ± 1.69) and generally increases with distance from major coastal river mouths [38]. Salinity values are relatively consistent throughout the year but can decrease substantially (< 20 g/L) during periods of increased river flow or flooding. However, no major floods occurred during our sampling period (November-December 2013). Annual water temperature ranges from 13.8 to 30.5˚C (22.4± 3.5) [38]. Water temperatures between November and December typically range from 21.9 to 29.2˚C (24.8±1.4). Water quality monitoring studies over more than 10 years conclude that the nearshore areas (e.g. Bramble Bay, Deception Bay, Waterloo Bay) have lower water quality scores than areas further from the mainland (Central Bay and Eastern Bay) [39].

Study species
Melicertus plebejus is endemic to eastern Australia and is found from 20˚S (central Queensland) to 42˚S (north-eastern Tasmania) [40,41]. Juvenile M. plebejus use nearshore habitats as nurseries for two to four months in a one to four year life cycle [42]. As adults they eventually migrate offshore and do not return to estuarine waters [43,44]. Females spawn and release their eggs in oceanic waters. The eggs hatch, and a proportion of the free swimming post-larvae are trapped in tidal advective envelopes, dragging them into coastal estuaries to settle and grow throughout their juvenile phase [45][46][47]. Moreton Bay is considered an important nursery for M. plebejus at the northern end of their range [48]. In Moreton Bay prawns recruit from nearshore nursery habitats (which include bare substrate and seagrass habitats), to central areas of the bay in the austral spring, i.e. October and November [35, 48,49]. Melicertus plebejus ultimately migrate offshore, most likely via the northern passage, in January and February. Therefore all individuals collected in this study were juveniles. Melicertus plebejus is nocturnal and individuals burrow into the substrate during the day. In general, prawns have the ability to swim, although their primary mode of locomotion is walking. Melicertus plebejus is an opportunistic consumer. Its diet consists of small shellfish, other crustaceans, and worms, and it incidentally consumes plant material, detritus, sediment and micro-organisms [31][32][33].
Melicertus plebejus was selected as a bioindicator species because preliminary trawls and collections in Moreton Bay indicated that M. plebejus is one of the few benthic generalist consumers that can be found throughout Moreton Bay in a wide range of habitats [34,48]. Thus, by using M. plebejus as a study species, it was possible to examine tracer distribution across the majority of the bay without having to compare species with different tracer integration and accumulation patterns. Although M. plebejus is a mobile species, its migration patterns within Moreton Bay are well-understood and uni-directional (i.e. prawns do not exit and re-enter nearshore areas) [35]. Moreover, isotopic turnover of juvenile M. plebejus is relatively fast. Juveniles grow rapidly within nearshore habitats, increasing approximately 1000 times in biomass before they migrate offshore [45,50]. After 1-2 weeks of consuming a new food source, juvenile muscle tissue should reach full isotopic equilibrium [28][29][30]. Finally, previous research has shown prawns can be valuable bio-indicators of regional trends in anthropogenic pollution [27,51]. Therefore, although M. plebejus movement had the potential to affect our

Sample collection and analysis
Prawns were collected from eight areas in Moreton Bay in cooperation with the Long-Term Monitoring Program (Fisheries Queensland), which conducts annual fishery-independent recruitment monitoring surveys of M. plebejus in southeast Queensland. These areas were selected because they were thought to have distinct terrestrial and offshore inputs, and they contain different habitat types. We predicted these eight areas would exhibit distinct biogeochemical patterns which could be used to distinguish each area and identify different potential coastal influences. These eight areas included four inshore areas; (1) Deception Bay, which receives terrestrial input from the Caboolture River and northward flow from the Brisbane River; (2) Bramble Bay, which receives input from the Brisbane River and Pine River; (3) Waterloo Bay, which receives relatively limited input from the Brisbane River (as compared to Bramble Bay); and the (4) Southern Bay, which receives input from the Logan river, and oceanic water input via a small passage between the barrier islands (Jumpinpin). Bramble Bay primarily contains sandy bottom and mudflat habitats. Deception Bay and Waterloo Bay contain both seagrass beds and sandy/ mudflat benthic habitats. The Southern Bay primarily contains sandy habitats, coral, and seagrass. Samples were also collected from the: (5) Central Bay; (6) Eastern Bay; the transition zone between the Southern Bay and Central Bay, which was referred to as (7) Central South; and finally (8) Offshore sites, which were expected to be relatively pristine compared to inshore areas.
At each site, ten prawns representing the range of size classes present were collected using a beam trawl aboard the FRV Tom Marshal trawler. This occurred from November to December 2013. Prawns were kept frozen (-20º C) until they were processed for analysis. Prawns were rinsed with deionized H 2 O, blotted dry, weighed (g) and measured (total length, mm). Tail muscle tissue was dissected and the gut was removed. Muscle tissue was chosen for tracer analysis because muscle tracer values in crustaceans provide conservative estimates of trace metal accumulation [52]. Muscle tissue was expected to provide a more holistic representation of tracer distribution as it can be more strongly associated with diet (i.e. particulate sources) [53,54] as compared to the gills or hepatopancreas, where values may be more closely associated with tracer availability in the water column (dissolved sources) [54,55]. The isolated muscle tissue was rinsed three times to remove residual salts. Sub-samples were divided for isotope analysis and major and trace element analysis. Isotope values were determined for each individual, while 10 individuals were combined in composite samples to ensure there was enough sample for trace element analysis.
Stable isotope analysis. Isotope samples were oven dried at 60º C for 48 hours and ground into a fine powder using a stainless steel ball mill grinder. 9-10 mg of dried, ground tissue was weighed in tin capsules and an oxidizing agent (vanadium pentoxide) was added to aid combustion of S for δ 34 S analysis. Samples were analysed using a Sercon Europa EA-GSL elemental analyzer and a Hydra-20-20 continuous flow isotope ratio mass spectrometer. Stable isotope ratios were expressed in δ notation as deviations from standards in per mil (‰): δX = [(R sample /R standard )-1] x 1000, where X is 13 C, 15 N, or 34 S, R sample is the ratio ( 13 C / 12 C, 15 N/ 14 N, or 34 S/ 35 S) in the sample, and R standard is the ratio in the standard. Laboratory and National Institute of Standards and Technology (NIST) standards were analysed to estimate analytical precision and accuracy. Reference materials were Vienna Pee Dee Belemnite carbonate (VPDB) for δ 13 C, atmospheric N 2 for δ 15 N, and Vienna Cañon Diablo troilite (VCDT) for δ 34 S. The analytical precision (standard deviation) for NIST standard NBS 127 (n = 6) and NBS 1577 (bovine liver, n = 6) was 0.1‰ and 0.1‰ for δ 13 C, 0.1‰ and 0.3‰ for δ 15 N, and 0.6‰ and 0.7‰ for δ 34 S. Duplicate samples were run as an additional measure of analytical precision. The analytical precision for duplicates was calculated as the mean and standard deviation of the absolute difference between each set of duplicate samples. Differences were 0.2 ± 0.11‰, 0.3 ± 0.26 ‰, and 0.3 ± 0.34‰ for δ 13 C, δ 15 N, and 34 S, respectively.
Major and trace element analysis. Composite samples for trace element analyses were freeze-dried for 48 hours and powdered using a zirconium shatter box grinder. Samples were prepared using microwave-assisted digestion with laboratory grade distilled nitric acid (HNO 3 ) (70%), and analysed using an Agilent 7700 Series inductively coupled plasma mass spectrometer (ICP-MS). The primary certified reference materials (CRMs) included NIST 2976 (muscle tissue for trace elements), AGAL 3 (prawn tissue for trace elements), BCR 668 (muscle tissue for rare earth elements), and NCSZ 73034 (prawn tissue for trace metals and rare earths). A secondary reference material (SRM) was developed using locally sourced prawns, with means and standard deviations established with reference to CRM results of the same analytical run. The practical quantification limit (PQL) for individual elements were established by running 10 digest blanks with PQL calculated as: PQL = 3 × (mean + (3 standard deviations)). Spiked recovery tests were carried out to optimize the ICP-MS settings for analytes. Major (Al, Fe, Mn), trace (As, Cd, Co, Cu, Ni, Pb, V, Y), and rare earth (Ce, La) elements were simultaneously measured. Duplicate samples were run as an additional measure of analytical precision. The analytical precision for duplicates was evaluated by calculating the coefficient of variation, (100 x (standard deviation/mean)) for each set of duplicate samples (S1 Table).

Data analysis
Spatial interpolation. Spatial variation in prawn stable isotope, and major and trace elements values between sites were initially visualized using circle plots in ARC GIS. In addition, an interpolated raster layer was generated for each variable using the minimum curvature spline technique with barriers; whereby the interpolated surface was constrained by the geographic features of the bay (i.e. Moreton and Stradbroke Islands). The resultant isoscapes, and what will be referred to as "tracerscapes", were used to discuss broad trends in tracer dispersal across Moreton Bay.
Spatial isotope analysis. All remaining data analysis was done in the R environment version 3.4.4. ANOVA and posthoc Tukey's HSD tests were used to compare prawn isotope values between each sample area. Standard elliptical areas corrected for small sample sizes (SEA c ) were calculated for each sample area and used to estimate isotopic niche overlap between groups [56]. Each SEA c encompassed 40% of the data in each sample area. Corrected ellipses were fitted with maximum likelihood estimations. Overlap between each SEA c was measured as a percent (%). Standard elliptical area and overlap were calculated using the package SIAR [57] and SIBER [56]. SEA analysis was used in conjunction with the ANOVA and Tukey's HSD tests to identify isotopically distinct regions of Moreton Bay.
Discriminant analysis. Linear discriminant analysis was used to isolate the isotopes and elements (together referred to as tracers) that could distinguish the eight sampling areas and correctly assign prawn composites to known sample areas. Prior to linear discriminant analysis (LDA), Z-scores were calculated for each tracer at each site within each sample area. Z-scores were calculated by subtracting the mean M. plebejus Moreton Bay value for each tracer from the sample site value, and dividing the output by the standard deviation of the mean [51] XÞ=s " X ). This approach converts all values to a similar standard deviation scale and creates more reliable comparisons between tracers. Mean, standard deviations, and Z-scores of measured values for isotopes for each sample site are listed in S2 Table. Measured and Z-score  values for major and trace elements are listed in S3 and S4 Tables, respectively. Tracers were removed from the LDA model using a step-wise subtraction process that eliminated tracers that exhibited high collinearity or low importance in sample area discrimination. Collinearity was assessed using the tolerance statistic (1-R 2 , where R 2 is produced from a regression of one predictor variable against the remaining predictors). Tolerance values ranged from 0 to 1, where values close to 0 indicated variables were strongly correlated. The relative importance of a tracer in discriminating between sites was assessed using an F-to-remove statistic. Tracers with small F-to-remove statistics were deemed the least useful in discriminating between sites and were removed from the analysis. "Leave one out" cross validation was used to calculate assignment error rates and assess the accuracy of each subsequent version of the model. LDA models were fitted and assessed using the klaR [58] package. The results of LDA were used to identify important biogeochemical regions and tracers in Moreton Bay.

Results
A total of 249 prawns were collected from 25 sites inside (20) and outside (5) Moreton Bay. Only nine prawns were collected from Site 16-CS. Prawn length ranged from 47-156 mm (mean = 90.7 mm), and wet mass ranged from 0.48-25.6 g (mean = 5.7 g). There were significant differences in prawn length between areas (F (7, 251) = 12.43, P < 0.05; Table 1). Prawns from Deception Bay and the Southern Bay were significantly smaller compared to all other areas. Prawns from Waterloo Bay were also smaller than other sites, but were only significantly smaller than prawns collected from the Eastern Bay. The largest individuals were caught Offshore and in the Eastern Bay. Overall, prawns increased in size with distance from the mainland. This pattern was expected as western areas of Moreton Bay are key nursery grounds for M. plebejus, and as prawns grow, most move offshore through the northern passage of the bay [43,44].

Spatial isotope analysis
δ 13 C and δ 15 N analyses consistently demonstrated that prawn isotope values varied significantly between coastal areas. Circle plots and isoscapes indicated that δ 15 N values gradually decreased with increased distance from shore, while δ 13 C gradually increased (Fig 2A and 2B). In contrast, δ 34 S values were highly variable between sites with no consistent offshore trend ( Fig 2C), although mean prawn values were highest in Deception Bay. Offshore areas had significantly lower δ 15 N (F (1, 247) = 75.67, P < 0.05) and significantly higher δ 13 C than inshore areas (F (1, 247) = 36.59, P < 0.05). There was no significant difference in δ 34 S between inshore and offshore sites (F (1, 247) = 1.11, P > 0.05). There were also significant differences in prawn δ 13 C (F (7, 241) = 13.79, P < 0.05) and δ 15 N (F (7, 241) = 35.67, P < 0.05) between sample areas within Moreton Bay (Table 2). There were no significant differences in prawn δ 34 S between any sample areas (F (7, 241) = 1.65, P > 0.05). Prawns from nearshore areas, i.e. Bramble Bay, Deception Bay, and Waterloo Bay, generally had higher δ 15 N and lower δ 13 C compared to most other areas. Bramble Bay prawns had significantly higher δ 15 N values compared to all other areas ( Table 2; Fig 3). Bramble Bay prawns also had significantly lower δ 13 C values than all other areas with the exception of Deception Bay, and exhibited low isotopic overlap with all sample areas (0-10%). Waterloo Bay and Deception Bay prawns exhibited moderate isotopic overlap (13-41%), where prawns from Deception Bay had lower δ 13 C and higher δ 15 N than prawns from Waterloo Bay. However values were not significantly different between these two bays.  Central Bay and Central South prawns had mean δ 13 C and δ 15 N values that transitioned between nearshore and offshore sites, but individual prawns from these areas exhibited a wide range of δ 13 C and δ 15 N values. Central Bay and Central South prawns also exhibited a nearly identical range of isotope values (58-100% overlap). Central Bay and Central South prawns also exhibited high isotopic overlap with all other groups except for Bramble Bay and Offshore prawns (23-78%). Central and Central South prawn isotope values were significantly different from Bramble Bay and Offshore prawns. Central Bay prawns also had significantly higher δ 13 C than prawns collected from Deception Bay.
The Eastern Bay, Southern Bay, and Offshore sites had significantly lower δ 15 N and higher δ 13 C compared to nearshore sites and exhibited low isotopic overlap with nearshore areas (0-32%). Thus Offshore, Eastern Bay and the Southern Bay sites could be distinguished from nearshore areas using δ 13 C and δ 15 N alone. However, these three areas exhibited moderate to high isotopic overlap (24-79%) and were not significantly different. Eastern Bay and Southern prawns were also not significantly different from Central Bay or Central South prawns. δ 15 N and δ 34 S analysis indicated δ 34 S values were highly variable between sites within the same sample area, and δ 34 S poorly distinguished different areas of Moreton Bay (Fig 4). Overall, isotope analysis indicated there were three broad biochemical regions within Moreton Bay: (1) the three nearshore areas; (2) the two central areas; and (3) Eastern Bay, Southern Bay, and Offshore areas.

Major and trace element analysis
The elements Al, Fe, Mn, Pb, Cu, Ce, La, and Y had higher concentrations inside Moreton Bay compared to offshore sites, based on circle plots and tracerscapes (Fig 5). Manganese was generally elevated in all nearshore sites. Major elements Fe and Al, rare earth elements La, Ce, and the trace element Y (together referred to as Rare Earths plus Yttrium, REY), had correlated distribution patterns. Southern Bay and Deception Bay prawns had the highest concentrations of each of these elements. Prawn Pb values were highest in the Central and Southern Bay. Cu concentrations were highest in nearshore areas and decreased with distance from shore. Co, Ni, and V concentrations were also high within the bay, most notably in Waterloo Bay, however these trace elements also had high concentrations at some offshore sites. In contrast to most other elements, prawn Cd concentrations were higher in offshore than inshore areas. Cd concentrations were low and consistent in inshore areas, while offshore prawn Cd concentrations were high. Prawn As values also increased further from the mainland and were highest in prawns collected from the Eastern Bay and Offshore. However prawns collected from Bramble Bay, Waterloo Bay, and Site 21-SB within the Southern Bay area also had moderately high As values.

Linear discriminant analysis
Based on initial LDA, no combination of tracers could distinguish composites from Central and Central South areas. Central and Central South areas were therefore combined into one larger area (hereafter referred to as Central Bay). Initial analysis revealed strong correlations between Al, Fe, Ce, La, and Y. To prevent autocorrelation in the LDA, all but Ce were removed from the analysis. F-to-remove statistics indicated δ 34 S reduced the discriminatory power of the model and was removed. The subsequent removal of Ni, V, Cd, Mn, and Cu produced the simplest and most well-resolved model with 100% assignment accuracy (Fig 6). Arsenic, δ 13 C, Ce, and Co were the most important tracers in LD1, while δ 15 N was the most important tracer in LD2 (Table 3). Lead made smaller contributions to each LD function. Using combined results from each analysis, we identified seven distinct biochemical regions within Moreton Bay. Each area was defined by a unique "fingerprint" or combination of elements: (1) Bramble Bay composites were defined by the highest δ 15 N, low δ 13 C, and high Cu values; (2) Deception Bay composites had low δ 13 C, high δ 15 N, high Cu, low As, and distinctly high levels of REY; (3) Central Bay composites were defined by increased δ 13 C, decreased δ 15 N, high Pb, but relatively lower levels of other trace metals; (4) Eastern Bay composites had low δ 15 N, high δ 13 C, high As, and low levels of most other trace metals; (5) Offshore composites were distinguished by the lowest δ 15 N, high δ 13 C, high Cd, and relatively high As compared to some other areas; (6) the Southern Bay composites had high δ 13 C, low δ 15 N, relatively high Pb, the highest levels of REY, and low As and Cd; finally, (7) Waterloo Bay composites had intermediate δ 15 N and δ 13 C values compared to nearshore and offshore sites, but was primarily defined by high levels specific trace elements, most notably Co.
The inconsistent or patchy distribution of Ni, V, Co and As suggested these tracers may not be as reliable as tracers that exhibited more consistent and continuous distribution patterns, or those that were correlated with available sediment data. Therefore, additional LDA models were developed where Ni, V, Co and As were pre-emptively removed in a step-wise fashion and the remaining tracers were used to construct a more conservative assessment of potential biogeochemical regions. This approach helped to ensure that the supervised LDA technique was as objective as possible.
The pre-emptive exclusion of V and Ni had no effect on the final LDA model or model accuracy. This was not surprising as the original LDA indicated these tracers had low F-toremove values. However, the subsequent removal of Co changed the LD functions and reduced assignment accuracy to 85% (Table 4). Without Co, the LDA model could not distinguish Central Bay and Waterloo Bay composites. The subsequent removal of As resulted in another unique set of LD functions (Table 5), but maintained 85% assignment accuracy. However, without As it was more difficult to distinguish composites from Deception Bay and Southern sites. This would suggest that the removal of Co and As had a substantial effect on overall LDA accuracy. When Ni, V, Co, and As were excluded from analysis, the final LDA model indicated there were five to six conservative bioregions in Moreton Bay; (1) Bramble Bay, (2) Eastern Bay (3) Offshore sites, (4) Central Bay and Waterloo Bay, which could not be distinguished, (5) the Southern Bay and (6) Deception Bay, however the two latter regions were more difficult to distinguish compared to the original LDA model. Comparisons between different LDA models indicates that while Co and As were important for complete and accurate site discrimination, models that only included the more reliable tracers were still able to differentiate the most prominent biogeochemical regions in Moreton Bay.

Discussion
Our study identified significant spatial variation in a range of tracers in a dynamic, urban coastal environment. Most tracers displayed a gradual but clear west to east trend away from the coast, and the strongest biochemical distinction was between prawns collected from more oceanic (Offshore and Eastern Bay) and nearshore environments. Almost all of the tracers in this study could be used to monitor this transition. Based on these results, the shift from inshore to offshore environments likely strongly influenced the biochemical values of prawns in Moreton Bay. This trend is consistent with the well-established coastal influences typically used to define nearshore ecosystems. However, inshore areas and embayments were also readily distinguished by variation in trace elements. Stable isotopes were used to identify broader biochemical patterns and a few specific areas (e.g. Bramble Bay), while less commonly used tracers, such as REY, identified smaller biochemical regions. Discriminant analysis also revealed distinct biochemical "fingerprints" that could be used to define each area. Previous assessments of Moreton Bay have divided the bay into four regions, based on water quality: Western Bay, Central Bay, Eastern Bay and Southern Bay [39]. Our study extended past evaluations and produced a similar but slightly more differentiated view of the bay, with three regions along the western shore and two regions in the central bay. For example, Bramble Bay prawns had high levels of potential pollutants such as Cu and δ 15 N (which can indicate high levels of wastewater input [59]), and low δ 13 C (which suggests low levels of seagrass cover [11]). This combination of tracers suggests Bramble Bay is best defined by a "polluted nearshore" fingerprint. In contrast, the low levels of some urban pollutants and high δ 13 C in the Southern Bay prawns suggests this area is best defined as a relatively "pristine seagrass"  [11]. "Oceanic" prawns were readily identified by a combination of low Cu, Pb, and δ 15 N, but high Cd. Different regions of the bay were ultimately identified using a variety of disparate biochemical tracers that in-turn were indicative of variable land-based or oceanic inputs and local processes. These findings demonstrate that using stable isotopes and trace elements with a spatially explicit, multivariate approach can substantially increase the accuracy and precision of coastal spatial analysis. The use of M. plebejus as a bio-indicator was supported by the general consistency between prawn tracer values and abiotic values from previous studies in Moreton Bay. For example, Costanzo et al. [60,61] found that wastewater treatment plants were the primary source of elevated nitrogen loading and δ 15 N in Moreton Bay. δ 15 N was highest near coastal river mouths and adjacent bays, then decreased with distance from shore. Sediments within Bramble Bay and other nearshore areas also have relatively high δ 15 N content compared with other areas (Coates-Marnane pers. comm. 2016). This is consistent with spatial trends in prawn isotope values. Heavy metals that are often used as indicators of anthropogenic pollution were highest in prawns collected along the western shore and within the bay. Higher levels of these tracers in nearshore prawns is consistent with current water quality assessments and sediment patterns in the bay [37, 39]. The spatial regularity and agreement of δ 13 C, δ 15 N, Cu, Pb, and Cd also suggests these tracers are robust indicators of biochemical patterns in Moreton Bay and that M. plebejus is an effective bio-indicator at a regional level. However, it is important to note prawn migration from inshore nurseries to central bay areas had the potential to affect tracer analysis. The presence of migrant prawns in the Central Bay sites from biochemically distinct nursery areas may have increased tracer overlap between sample areas resulting in inaccurate tracer values and trends. Preliminary analysis of prawn isotope values indicated that one to three individuals from each of the Central Bay sites were recent migrants from nearshore nurseries. Therefore Central Bay composites included a few migrant individuals. However, given the low number of migrants at each site, composite values are likely accurate representations of local values and bay-wide trends. Using a less mobile species, such as sessile suspensionfeeder, would have been one way to eliminate this potential complication; however preliminary assessments of the bay suggested it would have been difficult to collect the same sessile suspension-feeding species from all areas of the bay. Because M. plebejus feed on plants, algae, animals and are found throughout Moreton Bay we believed M. plebejus tracer values would be a holistic representation of the different tracer sources found in both the sediment and in the water column and that values could be readily compared between sample sites.
M. plebejus size may also have affected our benthic biochemical maps of Moreton Bay. Previous studies have indicated crustacean size and growth rate can affect tracer uptake and accumulation, although size-based trends are often inconsistent between environments, tracers, and species [62][63][64]. Juvenile and sub-adult M. plebejus growth patterns and condition are relatively consistent between different regions of Moreton Bay (Munroe et al. unpublished data), suggesting that the effect of growth on spatial variation in M. plebejus tracer uptake are small relative to the effects of location and tracer bioavailability. Size-based variation in diet may also have affected tracer values between prawns from different sample areas [65]. Unfortunately, it is difficult to discriminate between size and location-based effects (i.e. food availability) on diet because M. plebejus ontogeny and habitat use are tightly linked. There is currently no research available analysing changes in M. plebejus diet with size or location. Given the broad scale of our analysis, the diverse diet of M. plebejus, and the similarity in M. plebejus and sediment tracer value trends throughout the bay, it is more likely prawn tracer values were primarily driven by location and bioavailability, rather than M. plebejus size. Nonetheless, effects of size on M. plebejus diet is an area of research that requires further study. As is detailed in the following sections, in-depth consideration of each tracer can lead to a better understanding of the potential effects of M. plebejus as a bio-indicator, the unique environmental influences that define each region, and help determine how different tracers can be applied to future coastal research.

Stable isotopes
SEA c overlap and discriminant analysis demonstrated δ 13 C and δ 15 N were highly effective and significant tracers. Spatial patterns and biochemical groups derived using δ 13 C and δ 15 N were likely the most reliable because isotope integration is relatively conservative and predicable, and values were measured at the individual level. The most distinctive isotopic pattern was the gradual change in δ 13 C and δ 15 N from the coastline to the offshore area. Relatively low δ 13 C in nearshore areas is often attributed to terrestrial organic and dissolved inorganic carbon (DIC) input into estuaries and coastal basins from coastal rivers. This material is more depleted in 13 C than marine or offshore sources [17], and can be indirectly incorporated into the diet of local consumers [66,67]. However, the Brisbane River, which is the largest coastal river in the bay, has low flow rates and high water residency time [36], suggesting that freshwater input from coastal rivers is typically not large enough to inject high volumes of terrestrial organic matter throughout individual bays. Therefore, with the exception of large flood events, direct terrestrial organic input is likely not the primary cause of the depleted 13 C that was observed in nearshore prawns. However, respiration in marine sediment by bacteria also produces depleted DIC. When shallow water sediment is disturbed and resuspended, the depleted DIC is released from the sediment and re-enters nearshore food chains, thus indirectly lowering the δ 13 C values of nearshore prawns [17].
Prawn δ 15 N values also helped identify areas with high anthropogenic input. Natural terrestrial sources of nitrogen, such as soil, typically have lower δ 15 N than human sources, such as treated wastewater from sewage [17,60]. As a result, δ 15 N is considered an extremely reliable metric with which to measure wastewater input in coastal environments [59]. Prawns with the highest δ 15 N values were collected near coastal river mouths and within adjacent bays. Prawn δ 15 N then decreased with increasing distance from shore. This pattern strongly suggests that the high δ 15 N in nearshore prawns was the result of wastewater input near the mouths of multiple coastal rivers. This is consistent with previous studies that showed elevated δ 15 N in the water column was due to wastewater input into Moreton Bay [60,61]. Prawn δ 15 N and δ 13 C may have also been influenced by seagrass habitat distribution. Low δ 15 N and high δ 13 C values are often associated with seagrass habitats [11]. The Southern Bay has the largest and heathiest seagrass beds in Moreton Bay [68]. Southern Bay and Eastern Bay prawns had higher δ 13 C and lower δ 15 N values compared to other nearshore prawns. This could indicate, similar to previous studies of benthic invertebrates in the area, that seagrass habitats and associated food resources are important in these regions [11,18]. δ 15 N analysis also indicated that benthic and pelagic habitats may respond differently to wastewater nitrogen input over time. Previous studies have shown that modern upgrades to wastewater treatment plants have substantially reduced the amount of excess nitrogen and δ 15 N in Moreton Bay [61,69]. However, both prawn and sediment δ 15 N patterns suggest that nearshore benthic habitats still contain wastewater N. There are several explanations for this apparent contradiction. Prawns may still be able to access sewage effluent that was buried before the upgrades were installed [69,70]. Periodic flooding of major coastal rivers can also resuspend and release contaminated sediment from rivers into adjacent coastal areas [61]. Consistent but low levels of wastewater input into Moreton Bay could also increase benthic δ 15 N values, while having no obvious effect on pelagic values. Whereas pelagic environments will quickly process and flush contaminated water out of the bay, contaminated sediment and particulate matter is more likely to settle and persist [71]. It is also possible that elevated nearshore δ 15 N reflects denitrification by benthic bacteria in areas enriched by N inputs from local rivers [72][73][74]. A combination of these processes, primarily controlled by historical and modern day wastewater input, is the most likely explanation for the observed δ 15 N benthic trend in Moreton Bay.
It was anticipated that δ 34 S would increase model accuracy and provide complementary analytical power [75]. Pelagic and offshore sources, such as seawater, typically have greater δ 34 S values compared to benthic food sources [76,77]. Therefore, it was expected that M. plebejus δ 34 S would exhibit a decreasing east-to-west trend, similar to δ 15 N. However, the mean prawn δ 34 S values did not exhibit any consistent pattern, and individual data indicated high within-site variance in all areas of the bay. While it is difficult to explain this result, δ 34 S cycling can be highly complex in estuarine and coastal environments and δ 34 S values can vary substantially between food sources [78]. For example, Oakes and Connolly [79] observed small scale δ 34 S variation in Moreton Bay seagrasses. If δ 34 S sources are not evenly distributed, this can result in dramatically different consumer δ 34 S values over a relatively small area. Each prawn trawl could have contained several resource patches with distinct δ 34 S. This may explain the inconsistent M. plebejus δ 34 S patterns observed in this study.
Ultimately, Moreton Bay isoscapes were able to delineate prawn populations that likely integrated greater amounts of wastewater derived N and nearshore C (Bramble Bay, Deception Bay, Waterloo Bay), populations that likely relied primarily on more offshore marine food sources (Offshore, Eastern Bay), and prawns that were transitioning between these two environments (Central Bay). Although prawns from Southern Bay sites likely incorporated seagrass food resources, isotope analysis could not accurately separate Southern Bay and offshore areas. It is important to note that most geographically distinct sample areas in this study could not be distinguished using stable isotopes. As previously discussed, distinct habitats can have similar isotopic values, which can make it almost impossible to distinguish different sources or environments without additional tracers. If prawns consume isotopically distinct food sources, this can also lead to high individual variability and isotopic overlap [80,81]. The presence of migrant prawns from isotopically distinct habitats can also increase overlap between distinct areas, although the effect of migration was likely relatively low based on the results of individual isotope analysis. Regardless of the ultimate cause(s), isotopic overlap between sample areas emphasises the importance of a multivariate approach to coastal biochemical mapping.

Major and trace elements
Major and trace elements, alongside stable isotopes, substantially improved separation between sampling areas. Although trace element distributions were only reported for one season, we had confidence in most tracerscapes because results broadly matched trace element patterns in sediments [37]. The M. plebejus distribution patterns of several major elements (Al, Mn, and Fe) and trace elements (Cu, Pb, and REY) strongly suggested these elements were derived from terrestrial sources. This is consistent with previous sediment analyses that showed local coastal rivers, most notably the Brisbane River, were the primary entry points for terrestrial soils and trace elements into Moreton Bay [37,82]. High inshore Mn, Fe and Al concentrations in M. plebejus likely reflects the higher Mn, Fe and Al found in nearshore, terrestrially-derived silts and clays [37]. Similarly, higher inshore prawn trace element values were likely derived from the higher concentrations of trace elements found in inshore, muddy sediments. Thus, there was broad agreement between prawn and sediment trace element patterns between inshore and offshore sites.
Within the bay, elements exhibited unique inshore distribution patterns. For example, high prawn Pb values were concentrated in the Central and Southern Bay, while high Cu values were dispersed along the mainland coast. This variation was critical to the high resolution of each LDA. The distribution of each tracer can be linked to different sources and processes, including different physiological processing specific to M. plebejus. Therefore, individual consideration of each tracer can help identify the distinctive chemical and environmental conditions that distinguish different areas of Moreton Bay. Evaluation of trace element distribution is partially limited by a poor understanding of tracer sources and uptake by marine biota. Trace element analysis of individuals would also have increased the analytical power of this study. However, exploratory results in the present study will help identify tracers for individual analysis of coastal bioindicators. In future, trace element distributions in different potential indicators, including sediment and species at different trophic levels, should be combined for a more complete and powerful understanding of element cycling in the bay.
Lead. High Pb values helped to distinguish Central and Southern Bay M. plebejus and was important in the final LDA model. The correspondence between sediment and prawn Pb patterns indicates sediment was likely the primary Pb source for M. plebejus in Morton Bay, and that M. plebejus were good bioindicators of benthic Pb distribution. As Pb is a toxic, persistent, immobile element with low solubility, fine sediments are often a significant Pb sink in aquatic habitats [83,84]. The Brisbane and Logan River are the most significant contributors of anthropogenic, Pb-rich sediment to Moreton Bay [37,82]. Pb values in both M. plebejus and sediment in the Central and Southern Bays can likely be attributed to urban run-off and industrial effluent transported via these rivers [26,37,82]. However, prawn Pb values from the Southern Bay were higher than expected given the low Pb levels found in Southern Bay sediment. This discrepancy may be due to differences in tracer bioavailability and prawn diet in a seagrass ecosystem as compared to the muddy benthic habitats which dominate other parts of the bay. It is also possible that elevated Pb in Southern Bay M. plebejus was originally accumulated closer to the mouth of the Logan River. Future fine-scale sampling could resolve these possibilities.
Copper. Nearshore prawns (i.e. Bramble, Deception, and Waterloo Bay) had higher Cu concentrations than prawns from all other sample areas. This trend was consistent with previous estimates of Cu values in Moreton Bay sediment [37]; however, elevated Cu was also detected in most inshore prawns. Cu is a common agricultural and industrial pollutant [85][86][87], and most rivers in Moreton Bay drain catchments with substantial agricultural and urban development. Cu is also a common pollutant near shipyards and marinas as it is a primary component of antifouling paint [88][89][90]. However, unlike Pb, Cu is a highly mobile, essential element. Elevated Cu in most inshore prawns likely reflects increased Cu exposure from a broad range of more urban and industrial sources in more mobile dissolved and particulate forms. Crustaceans can accumulate and retain more Cu than other biota exposed to similar concentrations of Cu [91]. Therefore, inshore prawns likely also had elevated Cu values because they retained high Cu concentrations as they moved away from more enriched nearshore habitats.
Arsenic. Prawn As values increased further from the mainland, and helped separate Eastern Bay composites in the full LD model. This trend is consistent with previous studies that found oceanic prawns in Australia have higher As concentrations than inshore populations (Burford, unpublished data), and prawns imported from other countries [92]. Based on these combined results, there may be a correlation between prawn As and marine habitat exposure in Australian waters. Although some forms of As are toxic [82,83,93,94], As is also a naturally ubiquitous element often found in high concentrations in seawater and marine primary producers, such as phytoplankton and algae [91,92,[95][96][97][98]. Thus, high As in offshore M. plebejus may indicate increased marine and planktonic food web contributions to its diet. This supposition is supported by similar trends in Moreton Bay prawn δ 13 C. Overall, bulk aquatic animal As values could be used to indicate increased marine or phytoplankton contributions to animal diet, as compared to contributions from estuarine or seagrass habitats which have much lower As values [97,99,100]. However, As trends in Moreton Bay were inconsistent between nearshore areas, and will likely fluctuate with seasonal changes in freshwater flow and phytoplankton productivity [95,98,101]. Thus, As trends need to be validated with replicate samples over time before it can be included in coastal biogeochemical assessments.
Cadmium. Prawn Cd distribution was particularly distinct as inshore prawns had substantially lower concentrations compared to offshore prawns. Although high Cd concentrations in aquatic biota can be associated with human based pollution, our findings suggest elevated prawn Cd was not derived from local anthropogenic sources [83]. High offshore Cd concentrations are consistent with previous studies which have shown north Australian oceanic prawns caught above 30˚N have higher Cd concentrations than imported or south Australian prawns [92,102]. At a national scale, geographic location is the largest determinant of prawn Cd concentrations [102]. However, our results suggest that at a local scale, Cd in northern prawns may be controlled by marine habitat exposure. The mechanism(s) behind this phenomena is not well understood at this time. Potential explanations may include increased bioavailability of Cd in oligotrophic waters, or increased Cd uptake in offshore habitats as a result of the competitive relationship for uptake routes with essential trace elements (Zn, Cu, Mn) [103][104][105][106]. Previous work has also indicated that bivalves, mussels and oysters integrate Cd from pelagic (phytoplankton) rather than benthic (sediment) sources [91,107]. Thus as M. plebejus transition to an offshore diet (characterized by a greater contribution from pelagic food webs), Cd levels may increase. More research on Cd uptake and distribution is required to better comprehend this local and national trend.
Rare earth elements and yttrium. The REYs (Ce, La, and Y), and the correlated major elements Al and Fe, differentiated inshore and offshore sites. REY are typically supplied to the coast through terrestrial run-off from coastal rivers [108,109] and are consistently higher in freshwater, estuarine, and nearshore environments [110]. However, prawn REY patterns in Moreton Bay were more complex than a simple inshore-offshore trend. Prawns from the Southern Bay and Deception Bay were distinguished from other inshore sites by high REY, Fe, and Al concentrations. Research on REY accumulation in marine animals is limited, but available data indicates that prawn REY concentrations were likely influenced by local differences in tracer availability and environmental conditions (i.e. pH, salinity, redox conditions). REY in freshwater systems are influenced by parent rock concentrations, lithology and weathering, water chemistry, and anthropogenic modification [13,108,111,112]. REY variation between coastal rivers directly influences REY variation between adjacent nearshore habitats [13]. This variation could be passed on to contiguous estuaries and reflected in local benthic species. However, research also indicates that animal REY uptake is nonlinear and can be affected by local conditions, such as changes in pH or alkalinity [110,113]. Therefore, prawn REY values are not necessarily an accurate representation of abiotic REY distribution.
The correlations between Fe, Al, and REY suggest prawn accumulation of these elements may be related. In marine habitats, REY preferentially bind to Fe-oxyhydroxide colloids found in organic, muddy, claylike sediments, and as a result Fe and REY sediment values are usually correlated [108,109,114,115]. Prawn REY concentrations in Moreton Bay cannot be compared to REY in sediment because sediment REY data are not available for this site. Nonetheless, spatial correlations between prawn Fe, Al, and REY values indicate that prawn REY accumulation may mimic expected correlations in sediment. Similar correlations have been identified in freshwater [116], estuarine [Burford et al. unpublished data] and coastal crustaceans [63]. Our results are also consistent with past work which suggested Fe-oxyhydroxides may be the main forms for REY uptake in amphipods [113]. However, while prawn Al, Fe and REY values had similar distribution patterns, Fe prawn patterns did not correspond with mud and Fe distribution in Moreton Bay sediment, suggesting that REY, Fe, and Al integration is more complex than simple association with muddy sediments and Fe-hydroxides. Despite these unknowns, our results have shown REY elements are useful additions to coastal spatial analysis, and including REY in future studies may prove highly beneficial.
Cobalt, nickel, and vanadium. Prawn Co, Ni, and V exhibited incongruent distribution patterns between sample sites, making it difficult to determine the environmental controls that effected tracer distribution. Ni and V were ultimately excluded from all the discriminate models, potentially due to their inconsistent distribution patterns; however Co was important and helped to distinguish Waterloo Bay from Central Bay samples. Although elevated levels of these tracers are usually associated with urban input [117,118], inconsistent distribution patterns suggest higher concentrations were not exclusively derived from mainland sources. Moreover, elevated prawn Ni concentrations were not associated with elevated Ni concentrations in sediment [37], which suggests prawn Ni was not exclusively controlled by sediment exposure. The biological role of Ni and V in marine animals remains unclear, and the ability of crustaceans to regulate Ni concentrations appears species specific [24, [119][120][121][122][123][124]. Marine animals can integrate Ni and V from solution [120,125,126] and their diet [122,126,127], so prawn Ni and V concentrations may be a combination of various sources, resulting in less predictable accumulation and distribution patterns. In contrast, Co is an essential metabolic element [128]. Prawns may accumulate Co to satisfy their metabolic requirements [128,129]. However, more research on Co uptake and distribution in marine environments is needed before the relevance of Co in this study and its potential applications can be more clearly discussed. Co, V and Ni results should be interpreted with caution because trends were defined by high values at only a few disconnected sites. It is possible these trends were distorted by outliers within the composite sample. Therefore, while discriminant analysis ultimately separated all of the original sample areas, the five to six biogeochemical regions defined by δ 13 C, δ 15 N, REY, Cu, Pb, and Cd, are probably the most robust and reliable groups. Any explanation about Co, Ni, and V distribution should be carefully considered within this context.

Conclusions
Our study has helped address key knowledge gaps linking tracer dispersal patterns and how tracers reflect sources and processes by identifying different biochemical regions in an urbanised coastal environment using a generalist bioindicator, M. plebejus. This species was an effective bioindicator of δ 13 C, δ 15 N and some trace elements, such as Pb and Cu. However, discrepancies between some sediment and prawn trace element distribution patterns (e.g. Ni) indicated that local environmental conditions and prawn physiology and diet also had a measurable effect on prawn values. Therefore, juveniles were useful indicators of environmental and chemical conditions at a regional level, but the results of some tracers should be interpreted with caution. The final analysis showed it was possible to distinguish seven of the eight original sample areas within Moreton Bay. These "fingerprints" were able to help distinguish the dominant influences, i.e. urban, catchment and offshore, that affected different sample areas. The majority of tracers showed clear transitions between nearshore and offshore environments. Inshore prawn δ 15 N, Cu, and Pb values indicated prawns were exposed to higher levels of different land-based pollution. Changes in δ 13 C, Cd, and As appeared to indicate a transition from terrestrial to marine resources. However, less commonly used tracers, such as REY and Co, were required to distinguish prawns from specific embayments. Greater study of the distribution and biological uptake of less commonly used biochemical tracers (i.e. REY), could help increase our understanding of coastal regions and species. In future, continued monitoring of Moreton Bay using the strategies presented in this study would allow researchers and managers to observe and track environmental changes in relation to local conditions and development. Ultimately, tracerscapes can lead to more accurate assessments of how changes in local processes and catchment inputs may affect coastal estuaries and basins.
Supporting information S1