Translator Disclaimer
1 November 2013 Macroinfauna Dynamics and Sediment Parameters of a Subtropical Estuarine Lake—Coombabah Lake (Southern Moreton Bay, Australia)
Author Affiliations +

Dunn, R.J.K.; Lemckert, C.J.; Teasdale, P.R., and Welsh, D.T., 2013. Macroinfauna dynamics and sediment parameters of a subtropical estuarine lake—Coombabah Lake (Southern Moreton Bay, Australia).

The distribution, composition, density, and biomass of benthic macrofauna within estuarine environments typically exhibit significant variations attributable to heterogeneity in and interactions between physical, biological, and chemical processes. The spatial and temporal dynamics of benthic macroinfauna assemblages and physicochemical sediment parameters within the intertidal mudflats of a subtropical estuarine lake (Coombabah Lake, Southern Moreton Bay) were studied at four sites from August 2006 to April 2007. No significant seasonal changes were observed at any site for all physical sediment parameters. The northern sample sites were characterised by fine- to medium-grained to moderately to poorly sorted sediments and the southern sample sites by fine-grained to moderately well to well-sorted. A total of 1029 individuals representing species from three orders, including deposit feeding and filter feeding macroinfaunal groups, were collected. The highest combined species densities occurred in the fine-grained southern sites, with the greatest combined species density occurring at Site 4 during winter. Amphipods (Victoriopisa australiensis) and polychaete worms (Simplisetia aequisetis) dominated the lake-wide faunal community with V. australiensis, representing 49% of the total retrieved macroinfauna. Significant correlations between mean macroinfauna densities, biomassDW, sediment parameters, and seasonal maximum monthly temperatures were identified during the study. Seasonal trends in combined site densities were observed at each of the lake sites, with the highest combined density occurring during winter. Spatial and temporal variations might also be partially explained by the predation pressures of fish and migratory wading birds within the lake, with the seasonal presence of migratory wading birds coinciding with the minimum observed macroinfauna densities at each sample site.


Despite ever-changing physical and chemical conditions, estuaries are characterised by high primary productivity (Alongi, 1998; Heip et al., 1995; Kennish, 2002) and concentrations of organic matter (Cifuentes, 1991; Dunn et al., 2008). Consequently, estuarine sediments often support high abundances and biomasses of benthic macrofauna (Moreira, Aldea, and Troncoso, 2010; Snelgrove, 1999), which provide vital food for many crustacean, fish, and shorebird species (Wolff, 1987; Ysebaert et al., 2005). As a result, intertidal habitats are of high conservation value (Fujii, 2007; Kennish, 2002) and often function as nursery grounds for important recreational and commercial fish species (Gray, McElligott, and Chick, 1996; Ysebaert et al., 2005).

The distribution, composition, density, and biomass of benthic macrofauna within estuarine environments typically exhibit significant variations attributable to the heterogeneity in and interactions between physical, biological, and chemical processes. To date, extensive efforts have been made to investigate macrofauna community dynamics and the role of influencing variables, such as hydrodynamic conditions, water depth, sediment composition, temperature, biogeochemical conditions, productivity and primary producer community structure, water column conditions, predator-prey relationships, and anthropogenic disturbances (e.g., Arvanitidis et al., 1999; Jones, Watson-Russell, and Murray, 1986; Kanaya and Kikuchi, 2008; McKindsey, Archambault, and Simard, 2012; Rossi, Castelli, and Lardicci, 2006; Ysebaert et al., 2002, 2003). Furthermore, investigations regarding temporal variations of benthic macrofauna are also of great interest (e.g., Cheng, 1995; Cheng and Chang, 1999; Jones, 1987; Kanaya, Suzuki, and Kikuchi, 2011; Morrisey et al., 1992; Sánchez-Moyano, García-Asencio, and García-Gómez, 2010; Ysebaert and Herman, 2002).

Because of their importance to estuarine food web dynamics and biogeochemical processes (Kristensen, 2000; Welsh, 2003), benthic macroinfauna are regarded as a key feature in estuarine monitoring programmes (Ysebaert and Herman, 2002). Through their feeding, burrow construction, bioturbation, and burrow irrigation activities, infauna can influence rates of organic matter inputs to the sediment, the vertical distribution of sedimentary organic matter, rates and pathways of organic matter mineralisation, and the fluxes of the regenerated dissolved nutrients back to the overlying water (Dunn et al., 2009, 2012; Jordan et al., 2009; Welsh, 2003). The extent of these influences depends on the functional groups of animals present and their abundance, biomass density, population dynamics, and individual size (Kristensen, 2000; Welsh, 2003). Therefore, investigations regarding macroinfauna community structure and dynamics are important when formulating models and quantitative predictions about the functioning of soft-bottom marine systems.

The present study investigated the spatial and temporal dynamics of benthic macroinfauna assemblages along with physicochemical sediment parameters within the intertidal mudflats of a subtropical estuarine lake of conservation significance. This work is intended to provide an initial assessment of the benthic macroinfauna assemblages in conjunction with physicochemical sediment and environmental parameters in an attempt to identify relationships and spatial and temporal variations within Coombabah Lake.


Site Description

Coombabah Lake (27°54′33″ S, 153°21′07″ E) is a subtropical, semiurbanised estuarine lake, which is part of Coombabah Creek, located in SE Queensland (Gold Coast, Australia) (Figure 1). The lake covers an area of approximately 2 km2; however, despite its modest dimensions, the lake and intertidal surrounds are recognised as a wetland of significant importance under the Convention on Wetlands as part of the Moreton Bay Ramsar wetland ( The site is also an important migratory bird habitat and is listed under the China–Australia Migratory Bird Agreement (1974) and Japan–Australia Migratory Bird Agreement (1986). Furthermore, the lake is a protected fish habitat area (Fisheries Act 1994) and provides a nursery for commercially and recreationally important fish, prawn and crab species. The lake is open to the Gold Coast Broadwater (a barrier island lagoon), through Coombabah Creek, and is subject to a mixed semidiurnal tidal regime with salinity ranging from ∼10 to 33 (Ali, Lemckert, and Dunn, 2010). The tidal and hydrodynamic conditions characteristics of the lake are described by Knight et al. (2008), Ali, Zhang, and Lemckert (2009), and Ali, Lemckert, and Dunn (2010). Annual rainfall for the region is approximately 1400 mm, with intense rainfall events common throughout the summer period. The minimum and maximum mean daily temperatures range from 9.2 to 28.5°C during July and January, respectively.

With the exception of shallow channels, Coombabah Lake is characterised by a relatively flat bathymetry, with water depths of 0 to ∼1 m at low water. During periods of low water, large portions of the lake sediments become exposed, which are characterised by mostly fine sediments away from the major drainage channels (Dunn et al., 2008).

The lake supports a mangrove dominated fringing flora, with Casuarina sp. and Melaleuca sp. on areas above the high-tide mark. Grey mangroves (Avicennia marina) form dense stands in low-lying areas connected to the lake in addition to Ceriops tagal, Aegiceras corniculatum, Bruguiera gymnorrhiza, and Rhizophora stylosa found in fringing mangrove stands. The lake edge–mangrove interface is ∼20 km in length, with mangroves extending landward between ∼5 and 700 m (commonly ∼400 m).

Sample Collection

Sediment samples were collected from four sampling sites within Coombabah Lake to determine macroinfauna abundance and biogeochemical properties, representing four regions of different surface sediment particle size distributions, organic matter sources, and hydrological features of the lake (see Dunn et al., 2012) (Figure 1). Six undisturbed sediment cores were collected during low water exposure from each of the four sites during consecutive seasons (winter [August] and spring [November] 2006; summer [January] and autumn [April] 2007) using sediment cores (plexiglass; 20 cm internal diameter × 33 cm length) for the determination of macroinfauna populations at each site during each season (i.e. total n = 96). Following collection, the sediment cores were transported to a constant temperature laboratory before being placed within a 220-L holding tank containing aerated site water, maintained under seasonal temperatures. Three additional sediment cores (PVC; 50 mm internal diameter × 400 mm length) were also collected at each of the four sites during each of the four seasons (i.e. total n = 48) for the determination of physical and biogeochemical parameters. Following collection these cores were immediately sliced into six depth horizons (0–1, 1–2, 2–4, 4–6, 6–10, and 10–15 cm) for determination of density, porosity, grain size distribution, organic matter content (LOI550), bioavailable ammonium ( i1551-5036-29-6a-156-ilm01.gif; porewater + exchangeable i1551-5036-29-6a-156-ilm12.gif), and bioavailable phosphate ( i1551-5036-29-6a-156-ilm18.gif; porewater + exchangeable i1551-5036-29-6a-156-ilm19.gif) concentrations. Subsamples of each depth horizon for the determination of i1551-5036-29-6a-156-ilm20.gif and i1551-5036-29-6a-156-ilm21.gif were immediately transferred to sample tubes containing 1 M KCl and 1 M MgCl2, respectively, and stored at <4°C in the field. All samples were returned to the laboratory within 2 hours of collection and stored at −20°C until analysis.

Sediment Analyses

Sediment wet-bulk density and water content were determined according to methods outlined by Percival and Lindsay (1997). Organic matter content (LOI550) was determined and calculated according to Heiri, Lotter, and Lemcke (2001). The dry sieving technique was used to determine the contribution of very coarse sand (1000–2000 μm), coarse sand (500-1000 μm), medium coarse sand (250–500 μm), fine sand (125–250 μm), very fine sand (63–125 μm), and mud (<63 μm) grain size fractions. The method of Folk (1968) was used to quantify the sorting characteristics of sediments. Grain size of sediments was recorded using the standard phi unit (Φ). A cumulative frequency plot of grain size (Φ) was used to measure the mean (Mz) and sorting characteristic (σI) graphically given by

where Mz is the graphical mean, σI is the sorting value (Inclusive Graphic Standard Deviation), and Φx is the grain size at the x% value on the cumulative frequency curve (Folk, 1968).

Descriptive terms of Folk (1968) were used to characterise Mz as calculated from Equation 1, including very coarse grained (−1–0.0 Φ), coarse grained (0.0–1.0 Φ), medium grained (1.0–2.0 Φ), fine grained (2.0–3.0 Φ), very fine grained (3.0–4.0 Φ), and coarse silt (4.0–5.0 Φ).

The sorting scale of Folk (1968) is used to characterise the sorting value σI as calculated from Equation 2: very well sorted (σI < 0.35), well sorted (σI = 0.35–0.5), moderately well sorted (σI = 0.5–0.7), moderately sorted (σI = 0.7–1.0), poorly sorted (σI = 1.0–2.0), very poorly sorted (σI = 2.0–4.0), and extremely poorly sorted (σI > 4.0).

Prior to nutrient analyses all laboratory glass- and plastic-ware were cleaned by soaking in 10% (v/v) HCl (>48 h) then rinsed three times with deionised water (Milli-Q; 18 MΩ cm). All reagents were of analytical grade purity, and all water used during dilutions was deionised water. Sediment i1551-5036-29-6a-156-ilm22.gif and i1551-5036-29-6a-156-ilm23.gif concentrations were determined by extracting aliquots of the homogenised sediment with 1 M KCl and 1 M MgCl2, respectively, by shaking. Following extraction, samples were centrifuged and the supernatant filtered (GF/F membrane, 47-mm internal diameter, Millipore) before spectrophotometric analysis at 640 nm using the phenate method for i1551-5036-29-6a-156-ilm24.gif concentrations and as molybdite reactive phosphorous (690 nm) for i1551-5036-29-6a-156-ilm02.gif concentrations (APHA, 1998).

Benthic Macroinfauna Analyses

Sediment cores placed in holding tanks were removed and the entire sediment core (20 cm internal diameter × 30 cm length) was sieved (250-μm stainless steel mesh) for the collection and identification of benthic macroinfauna species. Evaluation of the depth occurrences of the burrowing macroinfauna showed that more than 90% of the benthic infauna was concentrated in the top 0–10 cm interval in all cores sampled. Recovered macroinfauna were rinsed with freshwater to remove any adhering sediment or detritus and were preserved in 70% ethanol. Specimen identification (species level) and counts were performed under a low powered microscope. Dry weight (biomassDW) was determined after drying at 80°C for 48 h and reported as g m−2. Shannon diversity index (H′) was calculated according to Lee (1999).

Environmental Data

Regional air temperature (°C) and rainfall (mm) data obtained from the Australian Bureau of Meteorology automated Gold Coast weather station (Gold Coast Seaway, station number: 040764) were used as a complementary environmental characterisation of the study area. Previous hydrological measurements within the lake (data not shown) indicate that the approximate duration of sediment exposure per tidal cycle at Sites 1, 2, 3, and 4 was approximately 4, 3, 3, and 2 hours, respectively.

Statistical Analyses

Sediment parameter values for each sample depth horizon across the 15-cm depth profile are reported as mean values ± standard deviation (SD) from the triplicate cores for each site during each seasonal sampling event. Spatial and temporal variations (between site and/or season) were assessed using analysis of variance (ANOVAs), including Tukey's HSD or Games-Howell post hoc analysis. Prior to the analysis, the homogeneity of variance was evaluated using the Levene test, and, when violated, nonparametric analyses were performed. Pearson correlations (two-tailed; α = 0.05) were used to assess and identify relationships between physicochemical sediment conditions between sites, physicochemical sediment conditions (sample depth: 0–15 cm), environmental parameters, and benthic macroinfauna communities sampled (sample depth: 0–30 cm). Analyses were performed using SPSS Windows (SPSS, Inc., version 11.5.0). The infaunal assemblages and sediment characteristics at the sample sites were analysed using the PRIMER software package from the Plymouth Marine Laboratory (Clarke and Warwick, 2001).


Physicochemical Sediment Parameters

Daily temperature and rainfall data for the sampled period is shown in Figure 2. During the seasonal sampling events, regional temperatures ranged from 8°C (August) to 34.5°C (January). The maximum monthly rainfall occurred during November 2006 (122 mm).

No significant seasonal differences were observed at any site for all physical sediment parameters. Annual mean values of the sediment parameters with depth at each sampling site are shown in Table 1. Mean sediment wet-bulk densities ranged between 1.5 ± 0.1 (Site 4) to 2.0 ± 0.1 g cm−3 (Site 1) but were relatively consistent at each site with increasing depth. No significant difference was observed between sites. Water content throughout the sampled sediment profiles at all sites ranged from 20.4 to 69.6% and generally decreased with increasing sediment depth (Table 1). Water content significantly differed between all sample sites (p < 0.001, α = 0.05), with the exception of Site 2 compared to Site 3.

Mean grain size and sorting of sediments are presented in Table 1. The mean grain size ranged between very fine grained (Site 4, 0–1 cm and 1–2 cm depths) and coarse grained (Site 2, 10–15 cm depth). Northern sediments located toward the marine entrance of the lake (Sites 1 and 2) were typically described as fine- to medium-grained to moderately to poorly sorted sediments. In comparison, the southern sample sites are typically described as fine-grained to well to moderately well sorted. Sediment at the northern Sites 1 and 2 was significantly different in mean grain size compared to the southern Sites 3 and 4. For example, mean grain sizes throughout the sample profiles at Site 1 were significantly greater than Site 4 (p = 0.002, α = 0.05), and Site 2 was significantly greater than Site 3 (p = 0.004, α = 0.05) and Site 4 (p = 0.002, α = 0.05). At a depth of ∼10 cm, a patchy horizon of mollusc shells was present in the sediments at Sites 1 and 2.

Lake sediments were characterised by high organic matter contents with mean LOI550 values ranging from 1.21 ± 0.45% (Site 1, 1–2 cm depth, summer) to 9.57 ± 6.53% (Site 3, 10–15 cm depth, winter) (Figure 3). LOI550 values within the southern sites (3 and 4) were greater than the northern sites (1 and 2). Specifically, LOI550 values at Site 3 was significantly greater than Sites 1 (p < 0.001, α = 0.05) and 2 (p = 0.001, α = 0.05) values. Additionally Site 4 values were significantly greater then Site 1 (p <0.001, α = 0.05) and Site 2 (p = 0.003, α = 0.05) values. Maximum LOI550 contents within Site 3 and 4 sediments occurred during winter and within Site 1 and 2 sediments during autumn and spring, respectively.

i1551-5036-29-6a-156-ilm03.gif concentrations ranged between 1.8 ± 0.4 to 21.0 ± 8.1 nmol g−1 dry wt (Figure 4) and showed no significant difference between sample sites. Furthermore, no significant seasonal differences were observed. Additionally, i1551-5036-29-6a-156-ilm04.gif concentrations ranged from 13.2 ± 3.6 to 524.0 ± 26.0 nmol g−1 dry wt (Figure 5) and significantly changed between sites and seasons (p < 0.001, α = 0.05) during the study. Site 4 concentrations were significantly greater than Site 1 (p < 0.001, α = 0.05) and Site 2 (p = 0.001, α = 0.05), respectively. Maximal i1551-5036-29-6a-156-ilm05.gif concentrations at all sites throughout the lake occurred during the summer and autumn sampling periods.

Macroinfauna Dynamics

Over the course of the sampling programme, 1029 macroinfauna individuals were collected and identified, including six species from three orders representing deposit-feeding and filter-feeding macroinfaunal groups. Macrofauna abundance demonstrated significant spatial and temporal variations with the greatest combined seasonal abundances reported at the fine-grained Site 4 (p < 0.001, α = 0.05) and the greatest abundances across the lake recorded during winter (<0.001, α = 0.05). Seasonal trends in total site densities were observed at each of the four sites with temporal changes in macroinfauna densities characterised by a trending decline from winter through to autumn (Figure 6). The combined mean site densities were significantly correlated (negatively) with seasonal maximum monthly temperatures (r = −0.625, p = 0.010, n = 16). Mean combined densities ranged between 74.3 ± 74.5 ind. m−2 (Site 1, spring) and 955.5 ± 624.8 ind. m−2 (Site 4, winter; Figure 6). Victoriopisa australiensis (amphipods) followed by Simplisetia aequisetis (polychaete worms) were the most numerically dominant species collected (Figure 6). S. aequisetis was the greatest contributor to the burrowing faunal assemblages at Sites 1, 2, and 3, accounting for 60%, 87%, and 64% of the total individuals collected, respectively (Figure 6). V. australiensis accounted for 92% of all the total collected burrowing macrofauna at Site 4. The contribution of V. australiensis to mean densities decreased from the southern fine grained sites (Sites 3 and 4) in comparison to the northern fine to medium grained sites (Sites 1 and 2). Overall, V. australiensis represents 49% of the total retrieved macroinfauna across the four sites and four seasons sampled. No single species was recorded simultaneously at all four sites during any of seasonal sampling periods; however, V. australiensis, S. aequisetis, and H. cordiformis were recorded at all four sites during one or more seasons (Table 2). S. aequisetis demonstrated a seasonal pattern, with increased densities typically occurring during winter and spring.

Mean site biomassDW ranged between 0.16 ± 0.15 g m−2 (Site 1, summer) and 10.44 ± 11.69 (Site 1, winter; Figure 6). Combined biomassDW at Sites 1 and 2 were 91.5 g m−2 and 89.0 g m−2, respectively, compared to 18.0 g m−2 and 32.6 g m−2 for Sites 3 and 4, respectively. V. australiensis contributed greatest to Sites 3 and 4 biomassDW. The greatest contribution to biomassDW by S. aequisetis occurred at Site 3 (50%). Representing just 20% and 3% of the mean biomass at Sites 1 and 2, H. Cordiformis accounted for 61% and 82% of the recorded biomassDW at these sites, respectively.

Site 1 demonstrated the greatest seasonal Shannon diversity index values with increased values during winter (1.82) and autumn (1.74). Seasonal lake-wide diversity values ranged between 1.12 (spring) and 1.48 (autumn) with an annual mean value of 1.3 ± 0.16.

Relationships between Sediment Parameters and Macroinfauna Dynamics

Correlations between depth-averaged physicochemical sediment parameters are shown in Table 3. Significant correlations were observed between sediment parameters and macroinfauna mean density and biomassDW. The combined density of macroinfauna was significantly correlated with seasonal depth-averaged sediment density (r = −0.636, p = 0.008, n = 16), porosity (r = 0.570, p = 0.21, n = 16), and organic matter content (LOI550; r = 0.755, p = 0.001, n = 16). Mean densities of V. australiensis significantly correlated with mean i1551-5036-29-6a-156-ilm06.gif (0–4 cm; r = 0.505, p = 0.046, n = 16), mean i1551-5036-29-6a-156-ilm07.gif (r = 0.714, p = 0.002, n = 16), and LOI550 (r = 0.682, p = 0.004, n = 16). Additionally, mean density S. Aequisetis also significantly correlated with i1551-5036-29-6a-156-ilm08.gif (0–4 cm; r = −0.499, p = 0.049, n = 16) (Table 4). Furthermore, V. australiensis biomassDw significantly correlated i1551-5036-29-6a-156-ilm09.gif (0–4 cm; r = −0.508, p = 0.045, n = 16), mean i1551-5036-29-6a-156-ilm10.gif (r = 0.684, p = 0.003, n = 16), and LOI550 (r = 0.552, p = 0.027, n = 16). Mean biomassDw of H. cordiformis also significantly correlated with LOI550 values (r = −0.505, p = 0.046, n = 16, Table 4).


The northern sample sites (Sites 1 and 2) were characterised by fine to medium grained to moderately to poorly sorted sediments and the southern sample sites (Sites 3 and 4) by fine grained to moderately well to well-sorted sediments, as a result of the hydrodynamic regime and settling behaviour of sediment particles (Dronkers, 1986; van Leussen, 1999). The tidal regime of the lake and Coombabah Creek indicates pumping of sand-dominated sediments from the downstream marine entrance landward into the northern areas of the lake (Ali, Zhang, and Lemckert, 2009; Dunn et al., 2007a).

Nutrient concentrations measured during this study were comparable to surface sediment concentrations measured previously within the lake (Dunn et al., 2007b). Spatial and temporal differences in sedimentary inorganic nutrient concentrations throughout the lake are attributable to variations in sediment grain size, organic matter remineralisation, and external nutrient inputs. Sediment organic matter contents were greater within the southern lake fine-grained to moderately well to well-sorted sediments compared to the northern fine to medium grained to moderately to poorly sorted sediments. Organic matter content within the lake corresponded to the well-established and commonly observed correlation of increased organic matter with reducing sediment grain size (e.g., Dankers and Beukema, 1983; Dunn et al., 2008; Schrijvers, Van Gansbeke, and Vincx, 1995).

During a previous stable isotope and lipid biomarker study, Dunn et al. (2008) identified a greater contribution of autochthonous and labile organic matter to the sedimentary organic matter pool in the northern (marine entrance) sediments compared to the more allochthonous sourced organic matter of the southern region of the lake. Dunn et al. (2008) also observed that lower C/N ratios within the lake corresponded to regions of increased concentrations of chlorophyll-a and phaeopigment concentrations, with significant negative correlations occurring between chlorophyll-a and phaeopigment concentrations and C/N ratios, reflecting the influence of the low C/N ratio microalgal biomass and bacterial communities. In addition to stimulating microbial metabolism, sediment oxygen demand and nutrient regeneration (Ferguson, Eyre, and Gay, 2003; Slomp et al., 1993), sedimentary organic matter is a major food source for benthic infauna (Cheng and Chang, 1999; Spilmont et al., 2009) despite its often refractory nature (Cheng, 1995). Variations in the delivery and bioavailability of organic matter and regenerated nutrients, therefore, might influence benthic infauna population dynamics (Cheng and Chang, 1999; Ponti, Casselli, and Abbiati, 2010; Sánchez-Moyano, García-Asencio, and García-Gómez, 2010). During this study, increases in the combined macroinfauna density were observed to significantly correlate (r = 0.755, p = 0.001, n = 16) with increased sediment organic matter content (food supply) within the surface sediments of Coombabah Lake, indicating that sediment column organic matter content was a factor influencing the combined macrofauna population dynamics of Coombabah Lake on a relatively small spatial scale (e.g., ∼500 m). The significant correlation was influenced by the increased contribution of V. australiensis to overall densities in the fine-grained southern sites characterised by increased contents of organic matter. Additionally, significant correlations between V. australiensis density and meani1551-5036-29-6a-156-ilm11.gif andi1551-5036-29-6a-156-ilm13.gif concentrations were also identified (Table 4). The findings of the study illustrate associations between V. australiensis densities, LOI550,i1551-5036-29-6a-156-ilm14.gif, andi1551-5036-29-6a-156-ilm15.gif concentrations within Coombabah Lake, which themselves share significant correlations (Table 3). Additionally, significant correlations were also observed between LOI550 content and mean biomassDW of V.australiensis (r = 0.552, p = 0.027, n = 16), again reflecting the importance of organic matter content–food supply to the deposit feeding species. Previously several studies have also established sediment characteristics as significant factors influencing the distribution of macroinfauna species in estuarine and coastal waters (e.g., Cheng and Chang, 1999; Kanaya and Kikuchi, 2008; Ysebaert et al., 2002, 2003).

Although the sample sites were paucispecific, the lake included species commonly observed in intertidal surface sediments along the coastline of eastern Australia (e.g., Hirst, 2004; Ruello, 1975). Recorded densities are similar to those previously reported within the intertidal lake and surrounding environments (Dunn et al., 2009; GHD, 2003) and those reported in other shallow intertidal sediments (Kanaya and Kikuchi, 2008; Kanaya, Suzuki, and Kikuchi, 2011; Mucha, Vasconcelos, and Bordalo, 2004).

Species present demonstrated temporal variations with low abundance during summer conditions and significantly greater abundance during winter conditions, with mean site densities significantly correlating (negatively) with seasonal maximum monthly temperatures (r = −0.625, p = 0.010, n = 16). In addition to abundances, seasonal shifts in community structures were also observed. This finding suggests that seasonal temperature variation plays a potentially key role in the overall temporal macrofauna population dynamics of Coombabah Lake. Temporal evolutions observed within the lake have also been reported in temperate (e.g., Ysebaert and Herman, 2002) and subtropical intertidal environments (e.g., Chang, 1995; Da Cunha Lana and Guiss, 1991; Fonseca and Netto, 2006). The temporal variations observed within Coombabah Lake, however, are in contrast to those recorded in some other studies within Australian intertidal systems that reported no significant seasonal trends (Dittman, 2002; Jones, Watson-Russell, and Murray, 1986; Rainer, 1981). In many cases, intra-annual trends in macrofauna dynamics are attributable to patterns of species recruitment, water column and sediment conditions, competition, and predation stresses (Moreira, Aldea, and Troncoso, 2010). These facts could help explain the lower values of density observed during summer in Coombabah Lake, with recruitment potentially occurring during autumn, together with an increase in food supply (organic matter) leading to maximum macroinfauna densities during winter. This could be especially true for Site 3 and Site 4 where maximum LOI550 values, representing a food supply, occurred during winter.

Predation is most often presumed to be an important organising force of soft-sediment community structures (Bell, 1980; Brey, 1991). Furthermore, a feature of soft sediments is that any epibenthic predator feeding on the infauna is almost bound to cause a significant disturbance to the substrate (Glassom and Branch, 1997), potentially further influencing infauna dynamics. Macroinfauna within Coombabah Lake provides food sources for commercially and recreationally important fish and migratory wading bird species. Predation pressures presumably include resident species, Sillago ciliate (whiting) and Platycephalus fuscus (dusky flathead), which would typically feed in sandier sediments (northern region), and Acanthopagrus australis (bream) and Mugil cephalus (mullet), which generally feed in more mud-dominated sediments (GHD, 2003). Intertidal regions and sand/mud flats are also important feeding grounds for shorebirds, which are attracted to prey such as infaunal polychaetes and amphipods (Evans et al., 1998). In addition to resident shorebirds, the lake often serves as an important bird-staging area along avian migratory routes as an important declared migratory bird habitat. Migratory wading birds arriving from Europe and Asia during October (spring) to May (autumn) tend to roost in close proximity to surrounding feeding grounds and move out to feed on the tidal flats, often in hundreds (Dunn, personal observation), during ebb tides. Increased predation pressures resulting from the presence of migratory wading birds within the lake coincide with the minimum observed macroinfauna densities within the lake during summer. Avian predation resulting in reduced densities of prey species in intertidal soft-sediment environments has been demonstrated by many studies (e.g., Evans et al., 1998; Hulscher, 1982; Quammen, 1981; Rosa et al., 2008; Zwarts and Blomert, 1992). This observation suggests that avian predation is potentially a key factor influencing the temporal macrofauna population dynamics of Coombabah Lake.


In conclusion, this study investigated physicochemical sediment characteristics and macroinfauna dynamics within the sediments of Coombabah Lake. Seasonal sampling events indicated that physical sediment parameters demonstrated spatial variations but little seasonal variability in comparison to nutrient and organic matter concentrations. Spatial and temporal variations of the deposit-feeding and filter-feeding macroinfauna populations were observed according to sediment condition and sample season. Variations in the physicochemical surface parameters within the subtropical lake provide macroinfauna with different habitats that resulted in the coexistence of infauna with different habitat selectivity within a relatively small scale (e.g., approximately 500 m). During the study the highest combined species densities occurred in the southern fine-grained to moderately well to well-sorted sediments (Sites 3 and 4). Variations of mean macroinfauna density and biomassDW within the lake significantly correlated to nutrient and organic matter content. Strong temporal trends in total site densities were observed at each of the lake sites with significantly greater abundances and biomass observed during the winter.

The temporal variation of lake macroinfauna observed during this study showed trends that are hypothesised to be related to changing biotic and abiotic factors, including meteorological conditions, organic matter supply, patterns of species recruitment, and predation stresses. Specifically, the presence of migratory wading birds (October to May) within Coombabah Lake coincided with the minimum observed macroinfauna densities. In addition to interannual variations, biotic and abiotic factors would presumably determine long-term (intraannual) variations of macroinfauna assemblages throughout Coombabah Lake.


The authors thank W. Bennett, D. Dunn, K. Dunn, T. Dunn, M. Jordan, and D. Robertson for fieldwork assistance. The authors would also like to thank Dr. Kylie Pitt for guidance on the application of statistical analyses and comments of the early manuscript drafts.


  1. A Ali C.J Lemckertand R.J.K Dunn 2010. Salt fluxes within a very shallow sub-tropical estuary. Journal of Coastal Research, 26(3), 436–443. Google Scholar

  2. A Ali H Zhangand C.J Lemckert 2009. Numerical study of the hydrodynamics of a very shallow estuarine system—Coombabah Lake, Gold Coast, Australia. In: C.P da Silva (ed.), Proceedings of the 2009 International Coastal Symposium, Portugal. Journal of Coastal Research, Special Issue No. 56,922–926. Google Scholar

  3. D.M Alongi 1998. Coastal Ecosystem Processes. Boca Raton: CRC, p. Google Scholar

  4. American Public Health Association (APHA), 1998. Part 4000 inorganic nonmetallic constituents. In: M.A.H Franson (ed.), Standard Methods for the Examination of Water and Wastewater, 20th edition Washington, DC: American Public Health Association. Google Scholar

  5. C Arvanitidis D Koutsoubas C Dounasand A Eleftheriou 1999. Annelid fauna of a Mediterranean lagoon (Gialova Lagoon, south-west Greece): community structure in a severely fluctuating environment. Journal of the Marine Biological Association of the UK, 79(5), 849–856. Google Scholar

  6. S.S Bell 1980. Meiofauna-macrofauna interactions in a high salt marsh habitat. Ecological Monographs, 50(4), 487–505. Google Scholar

  7. T Brey 1991. The relative importance of biological and physical disturbance: an example from intertidal and subtidal sandy bottom communities. Estuarine, Coastal and Shelf Science, 33(4), 339–360. Google Scholar

  8. I.-J Cheng 1995. The temporal changes in benthic abundances and sediment nutrients in a mudflat of the Chuwei Mangrove Forest, Taiwan. Hydrobiologia, 295(1–3), 221–230. Google Scholar

  9. I.-J Chengand P.-C Chang 1999. The relationship between surface macrofauna and sediment nutrients in a mudflat of the Chuwei mangrove forest, Taiwan. Bulletin of Marine Science, 65(14), 603–616. Google Scholar

  10. L.A Cifuentes 1991. Spatial and temporal variations in terrestrially-derived organic matter from sediments of the Delaware estuary. Estuaries, 14(4), 414–429. Google Scholar

  11. K.R Clarkeand R.M Warwick 2001. PRIMER v5: User manual. Plymouth, UK: PRIMER-E Limited. Google Scholar

  12. P Da Cunha Lana and C Guiss 1991. Influence of Spartina alterniflora on structure and temporal variability of macrobenthic associations in a tidal flat of Paranaguá Bay (southeastern Brazil). Marine Ecology Progress Series, 73, 231–244. Google Scholar

  13. N Dankersand J.J Beukema 1983. Distributional patterns of macrozoobenthic species in relation to some environmental factors. In: W.J Wolff (ed.), Ecology of the Wadden Sea, Vol. 1. Balkema, The Netherlands: Rotterdam, pp. 69–103. Google Scholar

  14. S Dittmann 2002. Benthic fauna in tropical tidal flats of Hinchinbrook Channel, NE Australia: diversity, abundance and their spatial and temporal variation. Wetlands Ecology and Management, 10(4), 323–333. Google Scholar

  15. J Dronkers 1986. Tidal asymmetry and estuarine morphology. Netherlands Journal of Sea Research, 41, 109–118. Google Scholar

  16. R.J.K Dunn A Ali C.J Lemckert P.R Teasdaleand D.T Welsh 2007a. Short-term variability of physico-chemical parameters and the estimated transport of filterable nutrients and chlorophyll-a in the urbanised Coombabah Lake and Coombabah Creek system, southern Moreton Bay, Australia. In: C.J Lemckert (ed.), Proceedings of the 2007 International Coastal Symposium, Australia. Journal of Coastal Research, Special Issue No. 50, 1062–1068. Google Scholar

  17. R.J.K Dunn C.J Lemckert P.R Teasdaleand D.T Welsh 2007b. Distribution of nutrients in surface and sub-surface sediments of Coombabah Lake, southern Moreton Bay (Australia). Marine Pollution Bulletin, 54(5), 606–614. Google Scholar

  18. R.J.K Dunn D.T Welsh M.A Jordan J.M Arthur C.J Lemckertand P.R Teasdale 2012. Interactive influences of the marine yabby (Trypaea australiensis) and mangrove (Avicennia marina) leaf litter on benthic metabolism and nitrogen cycling in sandy estuarine sediment. Hydrobiologia, 693(1), 117–129. Google Scholar

  19. R.J.K Dunn D.T Welsh M.A Jordan P.R Teasdaleand C.J Lemckert 2009. Influence of natural amphipod (Victoriopisa australiensis) (Chilton, 1923) population densities on benthic metabolism, nutrient fluxes, denitrification and DNRA in sub-tropical estuarine sediment. Hydrobiologia, 628(1), 95–105. Google Scholar

  20. R.J.K Dunn D.T Welsh M.A Jordan N.J Waltham C.J Lemckertand P.R Teasdale 2012. Benthic metabolism and nitrogen dynamics in a sub-tropical coastal lagoon: microphytobenthos stimulate nitrification and nitrate reduction through photosynthetic oxygen evolution. Estuarine, Coastal and Shelf Science, 113, 272–282. Google Scholar

  21. R.J.K Dunn D.T Welsh S.Y Lee C.J Lemckert P.R Teasdaleand T Meziane 2008. Investigating the distribution and sources of organic matter in surface sediment of Coombabah Lake (Australia) using elemental, isotopic and fatty acid biomarkers. Continental Shelf Research, 28, 2535–2549. Google Scholar

  22. P.R Evans R.M Ward M Boneand M Leakey 1998. Creation of temperate-climate intertidal mudflats: factors affecting colonization and use by benthic invertebrates and their bird predators. Marine Pollution Bulletin, 37(8), 8–12. Google Scholar

  23. A.J.P Ferguson B.D Eyreand J.M Gay 2003. Organic matter and benthic metabolism in euphotic sediments along shallow sub-tropical estuaries, northern New South Wales, Australia. Aquatic Microbial Ecology, 33, 137–154. Google Scholar

  24. R.L Folk 1968. Petrology of Sedimentary Rocks. Austin, Texas: Hemphills, p. Google Scholar

  25. G Fonsecaand S.A Netto 2006. Shallow sublittoral benthic communities of the Laguna Estuarine System, South Brazil. Brazil Journal of Oceanography, 54(1), 41–54. Google Scholar

  26. T Fujii 2007. Spatial patterns of benthic macroinfauna in relation to environmental variables in an intertidal habitat in the Humber estuary, UK: developing a tool for estuarine shoreline management. Estuarine, Coastal and Shelf Science, 75(1–2), 101–119. Google Scholar

  27. D Glassomand G.M Branch 1997. Impact of predation by greater flamingos Phoenicopterus ruber on the macrofauna of two southern African lagoons. Marine Ecology Progress Series, 149, 1–12. Google Scholar

  28. C.A Gray D.J McElligottand R.C Chick 1996. Intra- and inter-estuary differences in assemblages of fishes associated with shallow seagrass and bare sands. Marine and Freshwater Research, 47(5), 723–735. Google Scholar

  29. Gutteridge, Haskins, and Davey Pty. Ltd. (GHD), 2003. Coombabah Creek Environmental Inventory. Brisbane: GHD Pty Ltd. Report prepared for Gold Coast City Council, p. Google Scholar

  30. C.H.R Heip N.K Goosen P.M.J Herman J Kromkamp J.J Middelbergand K Soetaert 1995. Production and consumption of biological particles in temperate tidal estuaries. Oceanography and Marine Biology: Annual Review, 33, 1–149. Google Scholar

  31. O Heiri A.F Lotterand G Lemcke 2001. Loss on ignition as a method for estimating organic and carbonate content in sediments: reproducibility and comparability of results. Journal of Paleolimnology, 25(1), 101–110. Google Scholar

  32. A.J Hirst 2004. Broad-scale environmental gradients among estuarine benthic macrofaunal assemblages of south-eastern Australia: implications for monitoring estuaries. Marine and Freshwater Research, 55(1), 79–92. Google Scholar

  33. J.B Hulscher 1982. The oystercatcher Haematopus ostralegus as a predator of the bivalve Macoma balthica in the Dutch Wadden Sea. Ardea, 70, 89–152. Google Scholar

  34. A.R Jones 1987. Temporal patterns in the macrobenthic communities of the Hawkesbury Estuary, New South Wales. Journal Marine and Freshwater Research, 38(5), 607–624. Google Scholar

  35. A.R Jones C.J Watson-Russelland A Murray 1986. Spatial patterns in the macrobenthic communities of the Hawkesbury Estuary, New South Wales. Marine and Freshwater Research, 37(4), 521–543. Google Scholar

  36. M.A Jordan D.T Welsh R.J.K Dunnand P.R Teasdale 2009. Influence of Trypaea australiensis population density on benthic metabolism and nitrogen dynamics in sandy estuarine sediment: a mesocosm simulation. Journal of Sea Research, 61(3), 144–152. Google Scholar

  37. G Kanayaand E Kikuchi 2008. Spatial changes in a macrozoobenthic community along environmental gradients in a shallow brackish lagoon facing Sendai Bay, Japan. Estuarine, Coastal and Shelf Science, 78(4), 674–684. Google Scholar

  38. G Kanaya T Suzukiand E Kikuchi 2011. Spatio-temporal variations in macrozoobenthic assemblage structures in a river-affected lagoon (Idoura Lagoon, Sendai Bay, Japan): influences of freshwater inflow. Estuarine, Coastal and Shelf Science, 92(1), 169–179. Google Scholar

  39. M.J Kennish 2002. Environmental threats and environmental future of estuaries. Environmental Conservation, 29(1), 78–107. Google Scholar

  40. J.M Knight P.E.R Dale R.J.K Dunn G.J Broadbentand C.J Lemckert 2008. Patterns of tidal flushing within a mangrove forest: Lake Coombabah, South East Queensland Australia. Estuarine, Coastal and Shelf Science, 76(3), 580–593. Google Scholar

  41. E Kristensen 2000. Organic matter diagenesis at the oxic/anoxic interface in coastal marine sediments, with emphasis on the role of burrowing animals. Hydrobiologia, 426(1), 1–24. Google Scholar

  42. S.Y Lee 1999. The effect of mangrove leaf litter enrichment on macrobenthic colonization of defaunated sandy substrates. Estuarine, Coastal and Shelf Science, 49(5), 703–712. Google Scholar

  43. C.W McKindsey P Archambaultand N Simard 2012. Spatial variation of benthic infaunal communities in baie de Gaspé (eastern Canada)—influence of mussel aquaculture. Aquaculture, 356–357, 48–54. Google Scholar

  44. J Moreira C Aldeaand J.S Troncoso 2010. Temporal dynamics of gastropod fauna on subtidal sandy sediments of the Ensenada de Baiona (NW Iberian Peninsula). Helgoland Marine Research, 64(4), 311–320. Google Scholar

  45. D.J Morrisey A.J Underwood L Howittand J.S Stark 1992. Temporal variation in soft-sediment benthos. Journal of Experimental Marine Biology and Ecology, 164(2), 233–245. Google Scholar

  46. A.P Mucha M.T.S.D Vasconcelosand A.A Bordalo 2004. Vertical distribution of the macrobenthic community and its relationships to trace metals and natural sediment characteristics in the lower Douro estuary, Portugal. Estuarine, Coastal and Shelf Science, 59(4), 663–673. Google Scholar

  47. J.B Percivaland P.J Lindsay 1997. Measurement of physical properties of sediments. In: A Mudroch J.M Azcueand P Murdoch (eds.), Manual of Physico-Chemical Analysis of Aquatic Sediments. Boca Raton: CRC, pp. 7–46. Google Scholar

  48. M Ponti C Casselliand M Abbiati 2010. Anthropogenic disturbance and spatial heterogeneity of macrobenthic invertebrate assemblages in coastal lagoons: the study case of Pialassa Baiona (northern Adriatic Sea). Helgoland Marine Research, 65(1), 25–42. Google Scholar

  49. M.L Quammen 1981. Use of exclosures in studies of predation by shorebirds on intertidal mudflats. The Auk, 98(4), 812–817. Google Scholar

  50. S Rainer 1981. Temporal patterns in the structure of macrobenthic communities of an Australian estuary. Estuarine, Coastal and Shelf Science, 13(6), 597–620. Google Scholar

  51. S Rosa J.P Granadeiro C Vinagre S França H.N Cabraland J.M Palmeirim 2008. Impact of predation on the polychaete Hediste diversicolor in estuarine intertidal flats. Estuarine, Coastal and Shelf Science, 78(4), 655–664. Google Scholar

  52. F Rossi A Castelliand C Lardicci 2006. Distribution of macrobenthic assemblages along a marine gradient in Mediterranean eutrophic coastal lagoons. Marine Ecology, 27(1), 66–75. Google Scholar

  53. N.V Ruello 1975. Geographical distribution, growth and breeding migration of the eastern Australian king prawn Penaeus plebejus Hess. Australian Journal of Marine and Freshwater Research, 26(3), 343–354. Google Scholar

  54. J.E Sánchez-Moyano I García-Asencioand J.C García-Gómez 2010. Spatial and temporal variation of the benthic macrofauna in a grossly polluted estuary from southwestern Spain. Helgoland Marine Research, 64(3), 155–168. Google Scholar

  55. J Schrijvers D Van Gansbekeand M Vincx 1995. Macrobenthic infauna of mangroves and surrounding beaches at Gazi Bay, Kenya. Hydrobiologia, 306(1), 53–66. Google Scholar

  56. C.P Slomp W Van Raaphorst J.F.P Malschaert A Kokand A.J.J Sandee 1993. The effect of deposition of organic matter on phosphorus dynamics in experimental marine sediment systems. Hydrobiologia, 253(1–3), 83–98. Google Scholar

  57. P Snelgrove 1999. Getting to the bottom of marine biodiversity: sedimentary habitats–ocean bottoms are the most widespread habitat on Earth and support high biodiversity and key ecosystem services. Bioscience, 49(2), 129–138. Google Scholar

  58. N Spilmont T Meziane L Seurontand D.T Welsh 2009. Identification of the food sources of sympatric ghost shrimp (Trypaea australiensis) and soldier crab (Mictyris longicarpus) populations using a lipid biomarker, dual stable isotope approach. Austral Ecology, 34(8), 878–888. Google Scholar

  59. W van Leussen 1999. The variability of settling velocities of suspended fine grained sediment in the Ems estuary. Journal of Sea Research, 41(1–2), 109–118. Google Scholar

  60. D.T Welsh 2003. It's a dirty job but someone has to do it: the role of marine benthic macrofauna in organic matter turnover and nutrient recycling to the water column. Chemistry and Ecology, 19(5), 321–342. Google Scholar

  61. W.J Wolff 1987. Flora and macroinfauna of intertidal sediments. In: M.J Bakerand W.J Wolff (eds.), Biological Surveys of Estuaries and Coasts. Cambridge, UK: Cambridge University, pp. 81–105. Google Scholar

  62. T Ysebaertand P.M.J Herman 2002. Spatial and temporal variation in benthic macrofauna and relationships with environmental variables in an estuarine, intertidal soft-sediment environment. Marine Ecology Progress Series, 244, 105–124. Google Scholar

  63. T Ysebaert M Fettweis P Meireand M Sas 2005. Benthic variability in intertidal soft-sediments in the mesohaline part of the Schelde estuary. Hydrobiologia, 540(1–3), 197–216. Google Scholar

  64. T Ysebaert P.M.J Herman P Meire J Craeymeersch H Verbeekand C.H.R Heip 2003. Large-scale spatial pattern in estuaries: estuarine macrobenthic communities in the Schelde estuary, NW Europe. Estuarine, Coastal and Shelf Science, 57, 335–355. Google Scholar

  65. T Ysebaert P Meire P.M.J Hermanand H Verbeek 2002. Macrobenthic species response surfaces along estuarine gradients: prediction by logistic regression. Marine Ecology Progress Series, 225, 79–95. Google Scholar

  66. L Zwartsand A.M Blomert 1992. Why knot Calidris canutus take medium-sized Macoma balthica when six prey species are available. Marine Ecology Progress Series, 83, 113–128. Google Scholar

Figure 1.

Location of sampling sites within Coombabah Lake (a). Insert shows location of the study area in Australia (b) and southern Moreton Bay (c).


Figure 2.

Daily recordings of air temperature (minimum and maximum) and rainfall in the southern Moreton Bay region (Australian Bureau of Meteorology weather station, Gold Coast Seaway, station number: 040764).


Figure 3.

Seasonal depth profiles of sediment LOI550 contents at each sampling site within Coombabah Lake. Data are mean values, and error bars show the standard deviation of the mean rate (n = 3).


Figure 4.

Seasonal depth profiles of sediment i1551-5036-29-6a-156-ilm16.gif concentrations at each sampling site. Data are mean values, and error bars show the standard deviation of the mean rate (n = 3).


Figure 5.

Seasonal depth profiles of sediment i1551-5036-29-6a-156-ilm17.gif concentrations at each sampling site. Data are mean values, and error bars show the standard deviation of the mean rate (n = 3).


Figure 6.

Burrowing macroinfaunal dynamics, including (a) mean density, (b) species contributions to mean density, (c) total biomassDW, and (d) species contribution to total biomassDW.


Table 1.

Annual averaged depth profiles of sediment parameters at each site (mean values ± standard deviation; n = 12 [three replicates × four seasons]).


Table 2.

Density (ind. m−2) and biomassDW (g m−2) of Coombabah Lake macroinfauna species (mean values ± standard deviation or single recorded density where n = 1).


Table 3.

Pearson correlation coefficients between mean annual sediment parameter values for each sediment depth horizon (n = 24 [six sample horizons × four seasons]).


Table 4.

Species-specific Pearson correlation coefficients between mean density and depth-averaged (0–4 cm and 0–15 cm) mean sediment nutrient concentrations and mean LOI550 content and biomassDW and depth-averaged (0–4 cm and 0–15 cm) mean sediment nutrient concentrations and mean LOI550 content (n = 16 [four sites × four seasons]).

2013, the Coastal Education & Research Foundation (CERF)
Ryan J.K. Dunn, Charles J. Lemckert, Peter R. Teasdale, and David T. Welsh "Macroinfauna Dynamics and Sediment Parameters of a Subtropical Estuarine Lake—Coombabah Lake (Southern Moreton Bay, Australia)," Journal of Coastal Research 29(6a), (1 November 2013).
Received: 24 October 2012; Accepted: 26 February 2013; Published: 1 November 2013

Back to Top