Remote sensing of upwelling off Australia ' s north-east coast Mochamad

Cross-shelf processes drive the exchange of water between the continental shelf and western boundary currents, 10 leading to the import and export of heat, freshwater, sediments, nutrients, plankton, fish larvae, and other properties. Upwelling is an important process which modulates those exchanges. It regulates primary productivity, which in turn promotes higher trophic levels and fisheries. In this paper, we investigate upwelling events in the East Australian Current (EAC) intensification zone off Southeast Queensland through the analysis of remotely-sensed Chlorophyll-a (Chl-a) and Sea Surface Temperature (SST) as well as wind and ocean reanalysis products. A particular focus is on identifying the likely mechanisms that drive 15 upwelling events during the austral autumn to winter which are evident from cold SST and enhanced Chl-a concentrations. Four complementary Upwelling Indices (UIs) are derived. Chl-a (UIChla) and SST (UISST) based indices characterize the oceanic response to upwelling, while indices based on wind (UIw) and current (UIc) data capture the forcing of upwelling. The spatial and temporal variability of all UIs is examined over the continental shelf. It reveals distinct seasonal patterns. For the northern region, UIs identify the well-known Southeast Fraser Island Upwelling System. It prevails during the austral spring 20 to early summer and is driven by currentand upwelling favourable wind. In contrast, upwelling is enhanced over the southern shelf during austral autumn to winter. About 70% of all UISST and UIChla identified upwelling events occur during this period. A case study is presented that provides observational evidence for the existence of a shelf-break upwelling. Simultaneous downwelling favourable wind stress and upwelling favourable current-driven bottom stress establish a flow convergence in the bottom boundary layer (BBL). These convergent BBL flows force upwelling of cold and nutrient-rich slope waters as 25 evident from negative SST anomaly and enhanced Chl-a in austral autumn to winter. It is evident from these results that the shelf region is characterised by two distinct seasonally reoccurring upwelling regimes.


Introduction
Cross-shelf processes play an important role in the exchange and mixing of coastal waters over the continental shelf and the offshore waters of adjacent boundary currents.It leads to the export and import of properties such as sediments, nutrients, 30 carbon species and fish larvae (e.g.Brink, 2016).Coastal upwelling is a physical process of particular importance since it is Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License.associated with exchanges both vertically and horizontally.It regulates primary productivity and drives the major fisheries of the Eastern Boundary Upwelling Systems (EBUS) due to persistent winds (e.g.Wang et al., 2015;Kämpf and Chapman, 2016).
Upwelling also occurs in Western Boundary Currents (WBCs) but appears to be more sporadic, weaker in intensity and of 5 moderate spatial extent.It results from a combination of current-and wind-driven boundary layer stresses as found for the East Australian Current (EAC, Rossi et al., 2014;Brieva et al., 2015), the Brazil Current (Aguiar et al., 2014), the Agulhas Current (Leber et al., 2017) and the South Atlantic Patagonian shelf (Carranza et al., 2017).Interactions between along-shore wind stress oscillations and shelf-break fronts commonly observed in WBCs can have a significant influence on upwelling through Ekman dynamics (e.g.Siedlecki et al., 2011).Gibbs (2000) and Gibbs et al. (2000) found that concurent actions of downwelling 10 favourable winds and the EAC bottom stress along the shelf-break can establish flow convergence in the Bottom Boundary Layer (BBL) inducing vertical motion and upwelling.This provides a mechanism for the transient but substantial nutrient supply into the euphotic zone resulting in enhanced chlorophyll-a concentrations (e.g.Everett et al., 2014).This study aims to explain upwelling events off the Southeast Queensland coast, Australia (Fig. 1), during the austral autumn to winter and to identify the mechanisms that would explain remotely observed elevated Chlorophyll-a (Chl-a) concentrations and cold surface 15 waters.
The EAC is the WBC of the South Pacific Ocean and forms in the southern Coral Sea at about 15 o S. It intensifies further south into a narrow swift southward flow along the shelf-break extending to about 31 o S, where it separates from the shelf and turns eastward into the Tasman Sea (Ridgway and Dunn, 2003).The EAC exhibits seasonal variability in its transport and 20 hydrographic properties.Its flow is more intense and closer to the coast during summer (Ridgway and Godfrey, 1997).The EAC shelf encroachment drives shoreward cross-shelf transport of cold bottom slope water in the BBL which in turn preconditions upwelling (Schaeffer et al., 2013;Schaeffer et al., 2014;Schaeffer and Roughan, 2015).
The preconditioning of upwelling by the EAC and the simultaneous occurrence of upwelling favourable wind stress is known 25 to lift cold and nutrient-rich slope water into the euphotic zone (Rossi et al., 2014).This initiates and sustains Chl-a blooms (e.g.Everett et al., 2014) as well as drives cross-shelf phytoplankton composition and distribution (Armbrecht et al., 2015).Rossi et al. (2014) document the latitudinal difference in the frequency of upwelling events along the east Australian shelf.A maximum of eight days per month of current-driven favourable conditions was found between 25 o -32 o S. Wind-driven upwelling favourable conditions occurred on about 3-6 days per month with a maximum observed during austral spring and 30 summer (Rossi et al. 2014).
Most studies of upwelling along the southeast coast of Australia focus on the EAC separation zone located to the south of about 31 o S (e.g.Gibbs et al., 1998;Gibbs et al., 2000;Roughan andMiddleton, 2002, 2004;Schaeffer et al., 2013;Wood et al., Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License. 2016).Previous studies of elevated Chl-a, nutrients and Sea Surface Temperature anomalies (SSTa) and associated coastal upwelling off southeast Queensland by Middleton et al. (1994) and Brieva et al. (2015) focus on the austral spring to summer.
In this paper, we investigate the causes of remotely-sensed seasonally reoccurring mid-shelf elevated Chl-a and negative SSTa during austral autumn to winter.Furthermore, Brieva et al.'s (2015) results are based only on the analysis of remotely sensed Chl-a.Negative SSTa, indicative of upwelled cooler subsurface waters, was not investigated further and this paper expands on 5 that analysis.Rossi et al. (2014) provided a broader view of sporadic upwelling along the eastern Australian shelf without focussing in detail on the shelf of southeast Queensland.In addition, complex bathymetry, convoluted coastline, and a series of islands along 25°-27° S advocate for the need to re-evaluate and refine Rossi et al.'s (2014) findings in this region.
In this study, we apply a combination of ocean and climate datasets to investigate the causes of upwelling and to understand 10 the spatial and temporal evolution of Chl-a and SSTa in a less well studied region of the EAC.Wind-and EAC boundary layer stresses are the two main physical processes that drive upwelling along the coastal shelf of Australia (e.g.Schaeffer et al. 2013, Rossi et al. 2014, Brieva et al. 2015).Therefore, we utilize wind and current induced BBL stresses to compute upwelling indices and to quantify the characteristics of upwelling events based on these two drivers.Due to the temporal resolution of the data used in this study, we focus on long-lived (≥ 8 days) upwelling events.We analyse and compare upwelling indices 15 derived from satellite remotely-sensed Chl-a, SST, scatterometer wind, and ocean reanalysis to describe upwelling events in southeast Queensland shelf waters.The findings contribute to an improved understanding of the physical processes that maintain the prominent ecological and marine biodiversity which characterises the shelf waters of the region (e.g. Ward et al., 2003;Neira and Keane, 2008;Young et al., 2011;Dambacher et al., 2012).

20
This paper is organized as follows: Section 2 describes the geographical and physical characteristics of the study site, the data used and methods applied to compute four upwelling indices from Chl-a, SST, wind and current data.Our results are presented in Section 3, where the spatial features of the seasonally reoccurring elevated Chl-a and negative SSTa distributions as well as their temporal variability are discussed in Sections 3.1 and 3.2, respectively.Section 3.3 examines a census of upwelling events using the four upwelling indices.This is followed by a characterisation of the upwelling mechanisms from exploring the 25 temporal climatological variation of spatially averaged Chl-a, SSTa, Ekman transport and bottom layer stress in Section 3.4.
Discussion and conclusions are presented in Section 4 with a particular focus on a case study that underpins the conceptual model of shelf-break front promoted upwelling off southeast Queensland.
2 Study Site, Data and Methodology

Study Site 30
We refer to the study region as the Southeast Queensland Marine Coastal Zone (SEQMCZ).It is defined as the coastal shelf ocean between 25 o -27 o S and it extends eastward to the 200 m isobath (Fig. 1).The on-shelf circulation is characterised by the Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License.
seasonally appearing and wind-driven Fraser Gyre (Azis Ismail et al., 2017), which influences cross-shelf exchanges.It possibly facilitates the offshore transport of cold shelf water, Chl-a and larval fish (Ward et al., 2003;Keane and Neira, 2008;Ribbe et al., 2018).The climate of the region is subtropical with maximum rainfall during the austral spring and summer (e.g.Ribbe, 2014).Here, the EAC approaches the shelf-break and intensifies with the strongest flow occurring during austral spring and summer and the region is often referred to as the intensification zone of the EAC (Ridgway and Godfrey, 1997).5 High Chl-a concentrations and cold surface water are frequently observed south of Fraser Island (Brieva et al., 2015;Wijffels et al., 2018).This is probably due to different mechanisms which include tidal flushing and export of high-nutrient coastal mangrove water (Middleton et al., 1994), a combination of wind-and current-driven upwelling (Rossi et al., 2014) and possibly high river discharge (Brieva et al., 2015).In addition, Brieva et al. (2015) identified two distinct quasi-climatological features 10 of high Chl-a concentrations: the Southeast Fraser Island Upwelling System (Fig. 2a) and the mid-shelf blooms (Fig. 2b).The former has been described in detail and its causes identified by Brieva et al. (2015).Mean Chl-a values of > 0.8 mg m -3 in the northern region of the SEQMCZ during austral spring and summer are well above the background value of 0.2 mg m -3 which characterises Coral Sea water advected within the EAC.The upwelling is driven by the Ekman BBL stress of the EAC encroaching the shelf with the northerly (southward direction) along-shelf wind stress playing a secondary forcing role (Brieva 15 et al., 2015).However, in the southern region of the shelf, Chl-a values are close to the background value during this period and instead enhanced during autumn and winter.
The upwelled nutrient-rich water off Fraser Island generates coastal phytoplankton blooms, which support commercial and recreational fisheries (Ward et al., 2003;Young et al., 2011).It establishes the Southeast Fraser Island Upwelling System 20 (Brieva et al., 2015) as a key ecological region of high primary productivity and an important fishery along Australia's eastern marine region (Dambacher et al., 2012;Evans et al., 2017).The austral autumn to winter mid-shelf bloom identified by Brieva et al. (2015) downstream of Fraser Island has Chl-a values well above the background level of 0.2 mg m -3 .The bloom extends far offshore and its eastern limit follows closely the 200 m isobath (Fig. 2b).The austral spring to summer and autumn to winter difference in the mean Chl-a distribution (Fig. 2c) emphasises the different timing of the Chl-a maxima in the northern 25 (25 o -26 o S) as opposed to the southern (26 o -27 o S) region of the SEQMCZ.Brieva et al. (2015) found that the mid-shelf blooms appear to be unrelated to the along-shelf wind or EAC-driven BBL stress, leaving unexplained that local biological hotspot.Furthermore, in-situ wind measurements near 30 o S by Wood et al. (2016) show that from austral autumn to winter, the wind is dominated by south-easterly direction which is downwelling favourable.Thus, this study aims to investigate and identify the likely drivers of the mid-shelf elevated autumn-winter Chl-a concentrations.30 Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License.

Chl-a
MODerate resolution Imaging Spectroradiometer (MODIS) Chl-a estimates with a resolution of about 1 km were obtained from the Integrated Marine Observing System (IMOS, 2018), which is a National Collaborative Research Infrastructure (NCRIS) supported by the Australian Government.The data is available from the Australian Ocean Data Network (see URL: 5 https://portal.aodn.org.au/).Satellite ocean color derived Chl-a is used as a proxy for phytoplankton biomass to study the biological response to upwelling (Chen et al., 2012).Chl-a distributions are assessed to describe periods of high Chl-a concentrations as well as their spatial distribution (see Fig. 2).The 8-day composite Chl-a data from 2003 to 2016 is used to maximise the temporal resolution and minimise data gaps e.g.due to cloud coverage (Everett et al., 2014).The data was previously utilised by Brieva et al. (2015) who also discussed further data limitations.For example, in near coastal waters, 10 remotely sensed Chl-a is unreliable as bottom reflectance may cause errors.Therefore, the analysis and interpretation presented in this paper is restricted to water deeper than 40 m (see location of isobath in Fig. 1).

SST
Advanced Very High Resolution Radiometer (AVHRR) SST with a spatial resolution of 2 km is available from IMOS too (see URL: https://portal.aodn.org.au/).The cloud-free data comply with Group for High Resolution Sea Surface Temperature 15 (GHRSST) requirements and further technical details are provided in Griffin et al. (2017).The quality of the IMOS data product is similar to that of the daily SST analysed by Aguiar et al. (2014), Leber et al. (2017) and Jayaram and Kumar (2018) to detect upwelling off Brazil, northern Patagonia Shelf, southeast coast of South Africa and at southwest coast of India, respectively.The 6-day composite IMOS SST data is used in this study to identify long-lived upwelling events in the SEQMCZ during the 2003 to 2016 period.The data were previously used by Brieva et al. (2015) to identify the southeast Fraser Island 20 Upwelling System and its driving forces.

Scatterometer winds
An index from wind measurements is developed to characterise upwelling events and quantify associated magnitudes.Wind variability has been found to play a significant role in controlling phytoplankton growth (e.g.Kämpf, 2015).Daily data with a resolution of 25 km was extracted from two scatterometer satellite products for periods 1999 to 2009 (QuikSCAT) and 2007 25to 2016 (ASCAT) from www.remss.com.QuikSCAT wind speed and direction has accuracies similar to ASCAT data (Bentamy et al., 2012).

Ocean reanalysis
Daily zonal velocities from BRAN3p5 for the period 2003 to 2012 is used to quantify the EAC-driven upwelling variability.
BRAN3p5 is an ocean modelling data assimilation system (Oke et al., 2008;Oke et al., 2013)  general circulation off the east coast of Australia and other shelf regions (see e.g.Oke and Griffin, 2011;Wang et al., 2013;Everett et al., 2014;Brieva et al., 2015;Azis Ismail et al., 2017).Azis Ismail et al. (2017) utilized the model data to identify the seasonally appearing Fraser Gyre and found that BRAN3p5 data captured the observed time-varying regional coastal ocean circulation.

Methodologies
We develop four different upwelling indices to identify and quantify upwelling process.Chl-a and SST derived upwelling indices are used to examine the physical and biological response to upwelling since an increased Chl-a concentration and a negative SSTa is an indicator of upwelling (Section 2.3.1 and 2.3.2).Wind and BBL current stress are the two physical processes previously identified as forcing upwelling along the eastern shelf of Australia (e.g.Rossi et al., 2014;Schaeffer et 10 al., 2013).Here, we utilize wind-and EAC-derived upwelling indices to quantify the mechanisms of upwelling forced by along-shelf wind stress and EAC-driven BBL stress.The methods to obtain the indices are presented in Section 2.3.3 and 2.3.4 respectively.

Chl-a derived upwelling index (UIChla)
Turbidity of riverine input, re-suspended sediments, and bottom albedo can cause errors and result in overestimated Chl-a data 15 from remote sensing (e.g.Everett et al. 2014). Brieva et al. (2015) found that river runoff into the SEQMCZ appears to be limited to the near-shore inner-shelf zone of 10 km width coinciding with the approximate location of the 40 m isobath.Thus, all Chl-a estimates from locations with water depths of less than 40 m were excluded from analysis (Brieva et al. 2015).The same approach is used in this study.20 Brieva et al. (2015) applied a Chl-a background value of 0.2 mg m -3 for the SEQMCZ and determined the average length of a Chl-a event or bloom with about 7-8 days.The background value is similar to that found in the Coral Sea (Condie and Dunn, 2006) and to the east off Fraser Island (Dambacher et al., 2012).The same threshold value is used in this study.Data from two zonal cross sections located at 25.5 o S and 26.5 o S and extending eastward from the 40 m to the 200 m isobath (see Fig. 1) were selected to investigate the occurrence of upwelling favourable events that last at least 8 days.25

SST derived upwelling index (UISST)
Following the approach by Leber et al. (2017), we define upwelling favourable condition by examining the occurrence of "cold events".The method identifies the two mean SSTa threshold values of -0.57o C and -0.48 o C for the two previously identified representative zonal cross-sections and SSTa with lower values are indicative of upwelling.We also use Aguiar et al.'s (2014) method by substracting the observed SST from its centered 90-day mean to remove seasonality from SST.A spatial filtering Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License. is applied by excluding the near coastal pixels shallower than 40 m in the SST data to avoid the potential influence from bottom reflectance and turbidity.

Wind derived upwelling index (UIw)
The Ekman transport estimated from along-shelf wind stress is used to quantify upwelling (e.g.Cropper et al., 2014;Lamont et al., 2017).In this study, the wind upwelling index (UIw) is obtained from the Bakun Index modification (Bakun and Nelson, 5 1991).It is defined as the fraction of the Ekman transport that is perpendicular to the coast (Alvarez et al., 2011;Rossi et al., 2014).The method used to compute the surface wind stress is described by e.g.Azis Ismail et al. (2017).The UIw was computed according to Rossi et al. (2014): where Qx (Qy) is the zonal (meridional) Ekman transport (m 2 s -1 ); τy (τx) is the meridional (zonal) wind stress (N m -2 ); ρw is the mean density of sea water (1025 kg m -3 ); f is the Coriolis parameter (m s -2 ); α is the angle of the unitary vector pointing northward at the shoreline (degree).Positive and negative values of UIw correspond to upwelling and downwelling favourable conditions, respectively.

EAC derived upwelling index (UIc) 15
The upwelling index derived from EAC BBL stress (UIc) indicates whether the EAC drives strong onshore Ekman BBL flow along the shelf-break.This potentially uplifts cold and nutrient-rich slope water toward the shelf.Everett et al. (2014) use a threshold value of ≤ -0.08 N m -2 with smaller values indicating current driven upwelling favourable conditions.The same threshold value is applied in this study in order to identify current-driven upwelling favourable conditions.The  (N m -2 ) follows from: 20 where v b (u b ) is the meridional (zonal) velocity (m s -1 ) at the shelf-break and Cd is the drag coefficient (0.0025).The meridional velocity component was not rotated to the along-shelf direction which potentially underestimates the along-shelf bottom stress.
However, the difference in the bottom stress between the rotated and non-rotated meridional velocity in the SEQMCZ is only about 10 % (Brieva et al. 2015).25

Spatial Chl-a and SST features
The mean monthly Chl-a and SSTa climatologies capture the spatial structure of possible upwelling into the SEQMCZ (Fig. 2 and 3).During austral spring to early summer, increased Chl-a occupies predominantly the region just to the southeast of Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License.
Fraser Island with values exceeding 0.8 mg .m -3 (Fig. 2a).This Chl-a bloom characterises the Southeast Fraser Island Upwelling System (Brieva et al. 2015).Chl-a values are below the threshold value of 0.2 mg .m -3 in the southern region of the shelf.In contrast, enhanced Chl-a concentrations are found across all of the shelf from austral autumn to winter and approach the threshold value just to the east of the 200 m isobath (Fig. 2b).Two distinct larger-scale Chl-a patterns can be identified in Fig. 2c which show the difference between the two periods.The negative anomalies in the northern region (blue shading) 5 correspond to the Southeast Fraser Island Upwelling System, whereas positive anomalies in the southern region (red shading) correspond to the mid-shelf blooms.Similar to the Chl-a spatial structure from austral spring and early summer, we find negative SSTa to the southeast of Fraser Island.These anomalies are confined by the 40 m isobath and the shelf-break and extend southward to about 26.5 o S (Fig. 3a).In austral autumn and winter, SSTa gradients are weak and SSTa (colder than average) is negative across all of the southern region of the shelf (Fig. 3b). 10 Negative SSTa values (Fig. 3) coincide with enhanced Chl-a values (Fig. 2c).This is particularly evident during the austral spring and early summer period (Fig 3a,Fig 2c) and is an indication of seasonal upwelling.Positive SSTa values (Fig. 3a) characterise the pathway of warm EAC water that flows along the shelf-break.The presence the EAC appears to serve as an eastern barrier of the continental shelf Chl-a blooms.The climatological views reveal the signatures of two distinct possible 15 upwelling regimes that appear to occur within the SEQMCZ albeit at different seasons.

Temporal variability of Chl-a and SST
A longitude-time Hovmöller diagram of the estimated Chl-a concentrations (mg m -3 ) is presented in Fig. 4a for the northern region and in Fig 4b for the southern region of the SEQMCZ.The most noticeable feature is the inter-annual, seasonal and monthly increase of the Chl-a concentrations exceeding the background value of 0.2 mg m -3 .The Chl-a blooms extend to about 20 100 km offshore and mainly during austral autumn and winter.This is indicative of frequent enhancements of the phytoplankton biomass as measured via Chl-a across the shelf and throughout the year.The eastward extend of the Chl-a bloom provides evidence for the existence of a nutrient enrichment process such as upwelling that controls Chl-a values up to 100 km offshore.Alternatively, enhanced Chl-a values near the shelf break could also result from the offshore transport and dilution of a bloom that initially originated closer to the coast.During the austral spring to summer, it is evident that the 25 dispersion of Chl-a is restricted to the shelf.The presence of the EAC restricts the eastward extend of the Chl-a blooms.
In the northern region (Fig. 4a), higher Chl-a concentrations typically develop in austral autumn (April and May) and expand fully across the shelf during austral winter to spring each year.In austral summer (December to February), the Chl-a concentrations decrease and is detected offshore to about 40 to 50 km eastward.The initial development of the Chl-a blooms 30 across the shelf in the southern region each year is quite similar to that of the northern region (Fig. 4b).However, the blooms appear to be shorter in duration and occur primarily during austral autumn and winter.In contrast to the northern region, the Chl-a concentrations in the southern region decrease throughout austral spring to summer (November to February).The cumulative number of UIChla detected Chl-a events indicates that concentrations from austral winter to spring are higher than from summer to autumn in both regions (Fig. 6).The period with a higher number of cumulative events is longer in the northern than the southern region.There is an increase in the number of UIChla events (> 10) from May to November with a maximum in October in the northern region (Fig. 6a).The number of events decreases in the following months with a minimum 15 in February.In contrast, the period with higher numbers in cumulative UIChla events ( > 10) is shorter and lasts from May to September in the southern region (Fig. 6b).A maximum of 14 events is found from June to August and in the following months, the number of events decreases to a minimum in January.Thus, the detected pattern in the distribution of events differs distinctly between the northern and southern region.This is also reflected in the total number of events (Table 1).A total of 145 and 99 events with about 52% and 70% in the northern and southern region, respectively, occur during austral autumn and 20 winter.An average about 10 events occurs in the northern and 7 in the southern region per year (Table 1).

Evidence of upwelling from SSTa (index UISST)
The mean monthly distribution of the cumulative number of UISST identified events resemblances the UIChla derived pattern (Fig. 7).The number of events increases during the autumn to winter period in both regions with maxima of 9 (Fig. 7a) and 11 (Fig. 7b) in July, yet the overall duration of that period is shorter than that found from UIChla identified events.The total number 25 of cold events in each region is 42 (Table 1).This is about three cold events per year and compares to the five cold events identified by Leber et al. (2017) for the Agulhas current system to the southeast off South Africa.About 52% and 69% of these events occur during austral autumn to winter in the northern and southern region, respectively.The autumn to winter period of the southern region is characterised by a higher number of total events than the spring to summer period, while the northern region has about the same total number of events during both periods (Table 1).

Wind forced upwelling (index UIw)
Figure 8 shows the cumulative number of UIw upwelling and downwelling events identified by applying the methodology presented in Section 2.3.3.The UIw (upwelling favourable) events are more frequent from austral spring to summer than from austral autumn to winter in both the northern and southern region (Fig. 8a and 8c).In the northern region, between 7 and 9 events occur from September to December (Fig. 8a), while in the southern region, about 10 events per month occur during this 5 period with a maximum number of 16 events in November (Fig. 8c).During the following austral autumn to winter, the number of UIw (upwelling favourable) events decreases to less than 2 events in both regions.In contrast, the cumulative number of UIw detected downwelling favourable events increases from a minimum period between September and November to a period characterised by a maximum number of events that lasts until about July (Fig 8b and 8d).

10
The cumulative number of UIw detected events indicating downwelling favourable conditions is consistently much larger than the number of UIw detected events indicating upwelling favourable conditions.From January to July more than 10 events are found for each month (Fig. 8b and 8d).A maximum number of 16 events is observed during May and June for the northern and in March for the southern region of the shelf.Overall, the total number of UIw (upwelling favourable) events for the northern and southern region is 43 and 57 events, respectively (see Table 1).In contrast, the total number of UIw downwelling 15 favourable events is 132 in the northern region and 129 in the southern region.This is about two to three times the number of UIw upwelling favourable events.Despite the wind direction being mostly southerly and therefore downwelling favourable during austral autumn to winter, we still find positive UIw (upwelling favourable) events which account for 9 % and 7 % of the total in northern and southern region.However, negative UIw (downwelling favourable) events make up 62 % and 65 % of the total in the northern and southern region.On average, about 3 to 4 upwelling events (positive UIw) occur each year, whereas 20 downwelling events (negative UIw) occur about 9 times per year (Table 1).

Current forced upwelling (index UIc)
The cumulative number of current derived upwelling favourable events (UIc) is shown in Fig. 9 for the period 2003-2012.
Upwelling favourable events occur every month throughout the year in both regions.The events seem to occur more frequent in the northern than in the southern region.The two maxima of 8 and 10 upwelling favourable events occur in May and 25 September for the northern region, while in the southern region 10 and 7 upwelling events appear in May and April/October.In total, there are 79 (64) upwelling favourable events and about 44 % (48 %) of them take place in austral autumn to winter in the northern (southern) region of the shelf (see Table 1).On average, there are 5 to 6 current-driven upwelling favourable events each year over the period 2003-2012.The results shown in the preceding section reveal the existence of episodic and seasonal upwelling within the SEQMCZ.In this section, the link between the forcing mechanisms (wind-and current driven) and the oceanic response as evident from changes in Chl-a and SSTa is explored further.A climatological view of daily Chl-a yields additional evidence for the existence of different seasonal patterns in the northern and southern region (Fig. 10a).Chl-a increases to above the background level of 5 0.2 mg m -3 in both regions simultaneously from austral autumn to spring, but maximum values occur at different times (Fig. 10a).In the northern region, Chl-a exhibits a bimodal distribution with two maxima.The first maximum of about 0.42 mg m - 3 occurs in June and a second maximum of about 0.47 mg m -3 follows in August and September.In October, Chl-a values begin to decrease.In contrast, maximum Chl-a values in the southern region occur from April to July with a maximum of 0.42 mg m -3 in about mid-May.Increases and decreases in Chl-a appear to correspond to decreases and increases in SSTa (Fig. 10b).10 The Chl-a and SSTa correlation is negative with a coefficient of 0.82 and 0.42 (95% confidence level) for the northern and southern section, respectively.These significant correlations highlight that most of the Chl-a increases are associated with a negative SSTa, particularly for the northern region.The concurrent occurrence of negative SSTa and enhanced Chl-a is in agreement with the Ekman transport driven downwelling (upwelling) during the period of April to August (September to 15 November) with differences in the latitudinal locations.Strong and persistent negative Ekman transport (inshore direction) is frequently observed during austral autumn and winter with a maximum transport of -0.84 m 2 s -1 and -0.88 m 2 s -1 in the northern and southern region, respectively (Fig. 10c).The transport changes direction to mostly positive in austral spring to summer and reaches its peak of 0.46 and 0.67 m 2 s -1 for the northern and southern region, respectively (Fig. 10c).Note that negative (positive) Ekman transport generates downwelling (upwelling) favourable condition . 20 The correlation between Ekman transport and Chl-a is generally negative and weak (strong), with a correlation coefficient of about 0.25 (0.60) and 95% confidence level for the northern (southern) region.This suggests that Ekman transport influences the Chl-a response more in the southern than the northern region.However, the correlation between Ekman transport and SSTa for both regions is weak (correlation coefficient < 0.05).The correlation becomes only significant (correlation coefficient > 25 0.5) when the Ekman transport data are lagged by 33 and 112 days for the northern and southern region, respectively.The daily climatology of  was identified as upwelling favourable ( < -0.08 N m -2 ) all year and particularly in the northern region (Fig. 10d).Moreover,  lower than -0.2 N m -2 develop intermittently during December to March and May to September (Fig. 10d).In the south,  in lower than -0.2 N m -2 primarily in January and February, April to June and September and October.Nevertheless, we find that there is no significant correlation between  and Chl-a or SSTa.30 Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License.

Discussion and Conclusions
Results derived from the analysis of four upwelling indices suggest that upwelling events occur in the SEQMCZ throughout the year, but due to seasonally varying forcing mechanisms.The number of total events is small as compared to the permanent upwelling systems of the world ocean.Nevertheless, this upwelling is crucial in driving regional and seasonally high productivity and biodiversity in contrast to the otherwise very low background productivity of the oligotrophic waters of the 5 EAC intensification zone.
The UIChla and UISST events detected from austral spring to summer in the northern region of the shelf characterise the Southeast Fraser Island Upwelling System which is predominately current-driven (Brieva et al., 2015).In contrast, the austral autumn and winter UIChla and UISST events detected in both the northern and southern regions of the shelf are not associated with the 10 Southeast Fraser Island Upwelling System as proposed by Brieva et al. (2015).In addition, during austral autumn and winter, strong and persistent downwelling favourable southerly winds prevail in the region (Rossi et al., 2014;Wood et al., 2016;Azis Ismail et al., 2017) and current-driven uplift of cold and nutrient rich-water onto the inner shelf through the Ekman BBL may not necessarily reach the surface (Roughan andMiddleton, 2002, 2004;Schaeffer et al., 2014).

15
The analysis of the UIw is consistent with earlier postulations of predominantly downwelling favourable conditions along the east coast of Australia (e.g.Rossi et al., 2014;Brieva et al., 2015;Wood et al., 2016;Azis Ismail et al., 2017).Although the bottom stress of the EAC is at its weakest during austral autumn and winter (Brieva et al. 2015), a small number of UIc events is found during April to June in both the northern and southern regions (Fig. 9).This potentially intensifies upwelling within the BBL.44 % and 48 % of UIc events occur during austral autumn and winter in the northern and southern region, respectively.20 In addition, 62 % and 65 % of UIw downwelling events occur in the northern and southern region from austral autumn to winter (Table 1), when the typical numbers of UIChla and UISST events are at a maximum for both regions.
More than 50% of UIChla (52% and 70% for the northern and southern region, respectively) and UISST (52 % and 69% for the northern and southern, respectively) occur from austral autumn to winter.These results suggest that the simultaneous 25 downwelling favourable wind stress and persistent bottom stress may play an important role in driving Chl-a blooms and negative SSTa.Downwelling favourable wind stress coinciding with the EAC bottom stress are likely to establish convergence in the BBL that in turn creates a shelf break front (Fig. 11).The convergent flow could potentially bring cold and nutrient-rich slope waters close to or reach the surface.A similar mechanism was invoked to explain upwelling along the narrow shelf off Sydney, Australia (Gibbs, 2000;Gibbs et al., 2000).30 Downwelling favourable wind stress induces onshore flows within the surface Ekman layer (Figure 11).This is compensated SSH anomaly and the establishment of a near coast geostrophic northerly flow in the direction of the wind.This northerly flow during the autumn to winter was previously identified and is referred to as the Fraser Gyre (Azis Ismail et al., 2017).
Furthermore, the deeper offshore directed flow converges with an onshore flow that is associated with the EAC's bottom stress in the BBL over the shelf.This bottom flow convergence leads to a steepening of the frontal isopycnals and is likely to enhance upwelling in the vicinity of the shelf (Carranza et al., 2017;Gibbs et al., 2000).Numerical modelling studies and observations 5 for the EAC off Sydney, Australia, have shown that the interaction of downwelling favourable wind stress with the EAC can establish a flow convergence in the BBL (Gibbs et al. 1998;Gibbs, 2000;Gibbs et al., 2000).Carranza et al. (2017) presented observational evidence of an upward nutrient pumping mechanism along a western boundary current shelf-break.This upward nutrient transport has the potential to increase Chl-a over the shelf.Everett et al. (2014) 10 demonstrated that the combination of current-driven upwelling and wind-driven downwelling can provide favourable conditions for the formation of upwelling that leads to observed positive Chl-a anomalies between 32 o and 33.5 o S. Thus, the bottom flow convergence potentially formed under simultaneous downwelling favourable winds and the along-shelf EAC explains also the occurrence of elevated Chl-a and cold surface waters found during austral autumn to winter in both northern and southern regions of the SEQMCZ (see Fig. 2c and 3b).This situation is similar to those found in other WBC regimes 15 including the Patagonian shelf (Carranza et al., 2017) and the Middle Atlantic Bight (Siedlecki et al., 2011).
A case study is investigated further to elucidate that this process appears to operate in the SEQMCZ (Figure 12).A specific upwelling event is discussed to distinguish factors that force upwelling ( and Ekman transport) from factors that result from upwelling (enriched Chl-a, negative SSTa).An upwelling event from June 2007 is selected.During this period,  (< -0.08 N 20 m -2 ) promoting current-driven upwelling corrsponds to wind-driven downwelling and appears to result in a Chl-a bloom in southern region several days later.Fig. 12a-12c provide observational evidence of the shelf-break front promoted upwelling in the southern region from June 10-12, 2007.The upwelling event is characterised by a distinct band of enhanced Chl-a and is observed near the shelf-break at 25 about 26.5 o S in latitude.The Chl-a is more than 1.2 mg m -3 higher than those of the surrounding water.In addition, the filament of elevated Chl-a in the southern region is not connected with the high Chl-a observed in proximity of Fraser Island or the near-coast region.Thus, the high Chl-a values within the proximity of the shelf-break appear to originate locally (Fig 12a -12c).It is the shelf-break upwelling that carries nutrient-rich water from the BBL to the euphotic layer and increases primary productivity . 30 Values for  and Ekman transport support the mechanism of shelf-break front induced upwelling (Fig. 12d-12e).Before the on-set of the Chl-a bloom,  is consistently lower than -0.08 N m -2 from June 1-9 at 26. persistent downwelling favourable conditions.Upwelling favourable  and downwelling favourable surface Ekman transports are likely to induce a flow convergence in the BBL that leads to upwelling of nutrient-rich water into the surface layer, which is followed by an increase in Chl-a about a week later.
Given the seasonally reoccurring upwelling events, it is no surprise that the region is characterized by enhanced pelagic 5 productivity (Evans et al., 2017) and was previously identified as a key ecological site off Australia's east coast (Dambacher et al., 2012).Neira and Keane (2008) suggest a relationship between upwelling of cooler nutrient-rich water and an increase in planktonic biomass, including for example, the large abundances of blue mackerel eggs during austral winter and spring.
Frequent upwelling events of low intensity boost primary production, enhance productivity and establish the region as an important fishery (e.g.Young et al., 2011). 10 This study uses complementary upwelling indices derived from remotely-sensed and reanalysis data.In-situ observations from the study site do not exist.Thus, it is suggested that future work is required to investigate bio-physical processes and explore upwelling indices using in-situ measurements along with high-resolution modelling.In addition, some limitations to this study need to be acknowledged and taken into account in interpreting the results.The analysis presented here focuses largely on the 15 seasonal climatology and mean state of upwelling off southeast Queensland, however, we also identify much longer lasting negative SST anomalies persisting for about 2 months (see lower panel of Fig. 5) during the period of July to August each year.This is the period when cyclonic eddies frequently encroach onto the shelf (Everett et al., 2014;Ribbe et al., 2018).The role of cyclonic frontal eddies in raising Chl-a and cooling shelf water is not addressed in this study.Yet, a recent study by Ribbe et al. (2018) observed frontal cyclonic eddies off SEQMCZ during austral autumn and winter that are characterised by 20 elevated Chl-a and cold SST.This is indicative of coastal water entrainment and possibly eddy-driven upwelling.All those upwelling factors (wind, currents, frontal eddies) and processes (surface/bottom convergence/divergence) may occur concomitantly, superimposing their hydrographic signatures.Hence, future studies should integrate both numerical and observational approaches to better understand how all the processes interact and how the shelf waters ultimately respond.

Data availability 25
The long-term Chl-a and SST data are provided by the Integrated Marine Observing System (IMOS) and can be retrieved from the Australian Ocean Data Network at https://portal.aodn.org.au/.(N m -2 ) in Fig. 10d indicates the threshold value of  considered upwelling favourable (Everett et al., 2014).15            (N m -2 ) in Fig. 10d indicates the threshold value of  considered upwelling favourable (Everett et al., 2014).
Fig. 5.It demonstrates the variability of negative SSTa events with different intensities as well as their meridional extend.For 5 example, a negative SSTa of less than -3 o C is shown between 25.5 o and 26.5 o S lasting from about March to April 2012 (see lower panel of Fig 5).Seasonal and monthly negative anomalies, particularly during the austral winter to spring, are clearly identifiable from the time series shown in the lower panel of Fig. 5. Negative anomalies appear from about late July to October and usually disappear in January to February.
5 o S and even reached a minimum value of -0.4 N m -2 between 1-2 June 2007.At the same time, wind-driven Ekman surface layer transport is negative indicating Ocean Sci.Discuss., https://doi.org/10.5194/os-2018-142Manuscript under review for journal Ocean Sci. Discussion started: 9 January 2019 c Author(s) 2019.CC BY 4.0 License.

Figure captions 25 Figure 1 .
Figure captions 25Figure 1.The geographic location of the study site along the north-east coast of Australia.Indicated are key physical oceanographic features.40 m (dotted), 200 m (dashed) and 1000 m (solid) isobaths are indicated.

Figure 2 .
Figure 2. Mean monthly climatology of Chl-a (mg m -3 ) and for austral (a) spring to summer, (b) autumn to winter and (c) the difference of data between (a) and (b).White areas indicate Chl-a values of <0.2 mg m -3 which is the background value of EAC advected Coral Sea water.

Figure 3 .
Figure 3. Mean monthly climatology of SST anomaly ( o C) and for austral (a) spring to summer and (b) autumn to winter.

Figure 4 .
Figure 4. Chl-a (mg m -3 ) across sections at (a) 25.5 o S and (b) 26.5 o S, which represent the northern and southern regions 5 respectively.Distance 0 km is located at the 40 m isobath.

Figure 5 .
Figure 5. Negative SST ( o C) anomaly averaged zonally between the 40m and 200m isobath and for the period 1993 to 2016 (upper panel) and an enlarged period of 2003 to 2009 in the lower panel.

Figure 6 .
Figure 6.Cumulative number of UIChla events at (a) 25.5 o S and (b) 26.5 o S.

Figure 7 . 10 Figure 8 .
Figure 7. Cumulative number of UISST events at (a) 25.5 o S and (b) 26.5 o S. 10 Figure 8. Cumulative number of UIw events at (a, b) 25.5 o S and (c, d) 26.5 o S.

Figure 9 .
Figure 9. Cumulative number of UIc events at (a) 25.5 o S and (b) 26.5 o S.

Figure 11 .
Figure 11.Schematic cross-shelf diagram of flow convergence that is formed under simultaneous onshore flow due to EAC's bottom stress and offshore flow due to downwelling favourable wind stress in the BBL.The flow convergence establishes a shelf break front which promotes uplift of nutrients.

Figure 12 .
Figure 12.Chl-a (mg m -3 ) for June (a) 10, (b) 11, and (c) 12, 2007.Cloud coverage limits the number of days availalbe for analysis.Daily time series of (d)  (N m -2 ) and (e) Ekman transport (m 2 s -1 ) at 25.5 o S (blue line) and 26.5 o S (red line) from 20 25 May to 30 June 2007.The dotted black line at -0.08 (N m -2 ) in Fig.12dindicates the threshold value of  considered as upwelling favourable(Everett et al., 2014).

Figure 1 .
Figure 1.The geographic location of the study site along the north-east coast of Australia.Indicated are key physical oceanographic features.40 m (dotted), 200 m (dashed) and 1000 m (solid) isobaths are indicated.

Figure 2 .
Figure 2. Mean monthly climatology of Chl-a (mg m -3 ) for the SEQMCZ during 2003-2016 and for austral (a) spring to summer, (b) autumn to winter and (c) the difference of data between (a) and (b).White areas indicate Chl-a values of <0.2 mg m -3 which is the background value of East Australian Currents (EAC) advected Coral Sea water.

Figure 3 .
Figure 3. Mean monthly climatology of SST anomaly ( o C) and for austral (a) spring to summer and (b) autumn to winter.

Figure 4 .
Figure 4. Chl-a (mg m -3 ) across sections at (a) 25.5 o S and (b) 26.5 o S, which represent the northern and southern regions respectively.Distance 0 km is located at the 40 m isobath.

Figure 5 .
Figure 5. Negative SST ( o C) anomaly averaged zonally between the 40m and 200m isobath and for the period 1993 to 2016 (upper panel) and an enlarged period of 2003 to 2009 in the lower panel.

Figure 6 .
Figure 6.Cumulative number of UIChla events at (a) 25.5 o S and (b) 26.5 o S.

Figure 7 .
Figure 7. Cumulative number of UISST events at (a) 25.5 o S and (b) 26.5 o S.

Figure 8 .
Figure 8. Cumulative number of UIw events at (a, b) 25.5 o S and (c, d) 26.5 o S.

Figure 9 .
Figure 9. Cumulative number of UIc events at (a) 25.5 o S and (b) 26.5 o S.

Figure 11 .
Figure 11.Schematic cross-shelf diagram of flow convergence that is formed under simultaneous onshore flow due to EAC's 5 bottom stress and offshore flow due to downwelling favourable wind stress in the BBL.The flow convergence establishes a shelf break front which promotes uplift of nutrients.

Figure 12 .
Figure 12.Chl-a (mg m -3 ) for June (a) 10, (b) 11, and (c) 12, 2007.Cloud coverage limits the number of days availalbe for analysis.Daily time series of (d)  (N m -2 ) and (e) Ekman transport (m 2 s -1 ) at 25.5 o S (blue line) and 26.5 o S (red line) from 25 May to 30 June 2007.The dotted black line at -0.08 (N m -2 ) in Fig.12dindicates the threshold value of  considered as upwelling favourable(Everett et al., 2014).5