Abstract
A large proportion of estuaries along microtidal wave-dominated coastlines worldwide have entrances that intermittently close to the ocean when tidal currents and fluvial discharge are insufficient to erode sediment delivered onshore by waves. In this study, these systems are termed “intermittently open/closed estuaries” (IOCE) in order to include all estuaries which intermittently close to the ocean. IOCE do not fit neatly into existing generalized estuary classification models and have been traditionally recognized as a single estuary type that constitute a rare subset of wave-dominated estuaries. In this study, 111 estuaries in Victoria, Australia, are used to develop a classification model that delineates between different IOCE types. This was undertaken using historic aerial imagery and quantification of the estuary channel width, catchment area, lagoon dimensions and tidal prism derived from remotely sensed data. Field surveying of entrance morphology was undertaken for a subsample of 35 IOCE characteristic of each section of the coast and which had detailed entrance condition records. Using this subset, IOCE were classified into three distinct types using multiple methods of statistical delineation (non-metric multidimensional scaling, hierarchical cluster analysis and distribution analysis). These three types are: (1) Type A, the largest IOCE which both close and open infrequently but for the longest durations; (2) Type B, medium sized IOCE which open and close several times per year for weekly to monthly durations; and (3) Type C (tidal creeks), the smallest IOCE located specifically in high rainfall, mountainous catchments and which exist in a predominantly open state. The three types of IOCE showed an order of magnitude difference in entrance closure duration as controlled by variations in the catchment area, tidal prism volume, dimensions of the estuarine lagoon and the entrance channel at the mouth. The classification is also applicable to wave-dominated coastlines internationally where IOCE are present.
I Introduction
Estuaries are broadly classified based on the relative dominance of marine and fluvial processes (Dalrymple et al., 1992; Hume and Herdendorf, 1988; Kench, 1999; Kennedy, 2011; Woodroffe et al., 1989). In microtidal regions, wave-dominated estuaries are most common and occur in three forms: (i) drowned river valleys, (ii) barrier estuaries, and (iii) saline coastal lakes (Allen and Posamentier, 1994; Boyd et al., 1992; Roy, 1984). The specific morphology of each type further varies based on the degree of sedimentary infill. Youthful estuaries have limited infill and often deep basins, whereas mature estuaries have little accommodation space in the central basin and tend to be shallow (Heap et al., 2004; Roy, 1984; Saintilan, 2004). The estuary accommodation space refers to the available space for potential sediment accumulation (Slagle et al., 2006). As these systems infill, the entrance often progressively shallows and may close as the accommodation space becomes occupied with sediment. For this reason, entrance condition is often used as a proxy for defining estuarine type (Cooper, 2001; Roy, 1984).
In certain estuarine systems, the entrance is not permanently open to the sea, but oscillates between an open and a closed state (Perissinotto et al., 2010a; Potter et al., 2010; Ranasinghe et al., 1999). Estuaries which intermittently close to the ocean are common on microtidal, wave-dominated coastlines such as southwest and southeast Australia, where they are termed “intermittently closed/open lakes and lagoons” (ICOLLs) (Dye and Barros, 2005; Gale et al., 2006; Haines et al., 2006; Morris and Turner, 2010), South Africa, where they are referred to as “temporarily open/closed estuaries” (TOCE) (Perissinotto et al., 2010b; Whitfield, 1992), New Zealand (Hart, 2009; Kirk and Lauder, 2000), Mexico (Ruiz-Luna and Berlanga-Robles, 2003), south-east Asia (Piazza et al., 2010), and parts of North and South America (Behrens et al., 2013; Hanes et al., 2011; Kraus et al., 2002; Oliveira and Kjerfve, 1993). In this study, we refer to these systems as “intermittently open/closed estuaries” (IOCE), a universally applicable term that covers all estuaries which intermittently close to the ocean. Despite their ubiquity, IOCE are traditionally thought to be a transitional stage of permanently open estuaries on the trajectory to forming a fluvial delta and therefore do not fit neatly into generalized classification models (Cooper, 2001; Roy 1984, Roy et al., 2001).
General estuary classifications based on geomorphic processes are less common than those focusing on the physio-chemical, ecological and management attributes of estuaries. Within existing geomorphic classifications, IOCE are often excluded or recognized merely as a homogenous estuary type that are a rare subset of wave-dominated systems (Boyd et al., 1992; Dalrymple et al., 1992; Hume and Herdendorf, 1988), and there have been few attempts to differentiate between any sub-types of IOCE that might exist. At a continental scale, attempts at differentiating between IOCE sub-types tend to be conceptual classifications only (Cooper, 2001; Jacobs et al., 2011; Kjerfve, 1994). For South African systems, IOCE have been conceptually classified into two sub-types based on the frequency of their connection with the ocean. Intermittently open estuaries (IOE) are closed for > 50% of the time, on average, whereas intermittently closed estuaries (ICE) are closed for < 50% of the time (Bell et al., 2001; Cowley et al., 2001; James et al., 2007; Whitfield and Bate, 2007). In Australia, IOCE are becoming increasingly recognized as an estuary type in their own right (Harris and Heap 2003; Harris et al., 2002), although classification predominantly focuses on estuary ecology, or their degree of modification (Bucher and Saenger, 1991; Digby et al., 1999; Edgar et al., 2000; Murray et al., 2006). Along the coast of Victoria, Australia, which forms the study region for this project, IOCE are common and show considerable diversity in form and entrance condition (Kench, 1999). Site specific and broad geomorphic descriptions of a small subset of IOCE in Victoria have been used to develop classification models (Arundel et al., 2008; Barton and Sherwood, 2004; Barton et al., 2008; Mondon et al., 2003). However, these classifications again focus on the ecology and local physio-chemical variations in these systems. What is needed is a classification of IOCE based on geomorphic processes that differentiates between sub-types of these estuaries from easily measurable variables of marine and fluvial parameters.
Estuary entrance condition is linked to a suite of physical factors including: (a) tidal velocity and tidal prism (Escoffier, 1940; Hinwood et al., 2012; Jarrett, 1976; O’Brien, 1931); (b) the wave climate and annual littoral drift (Bruun, 1986; Ranasinghe et al., 1999); (c) river discharge (Cooper, 1994); (d) entrance orientation (Arundel et al., 2008; Barton and Sherwood, 2004), and (e) basin accommodation space (Roy, 1984; Roy et al., 2001). For example, estuaries with larger tidal prisms, in combination with a wider and deeper entrance channel, are likely to remain open for long periods due to the more efficient scour of channel and berm sediments (Bruun, 1986; Hume and Herdendorf, 1993). Despite these broad boundary conditions being identified, there has yet to be defined a specific combination of variables that controls the diversity in entrance closure regimes of IOCE. Identifying how these catchment-scale variables act as geomorphic controls on entrance condition is important for informing future management decisions regarding IOCE.
This study aims to use catchment-scale geomorphic variables to statistically delineate IOCE and identify potential subsets within IOCE. The microtidal coast of Victoria, Australia, containing over 100 estuaries, is the principle site of investigation. Remotely sensed and field data, in combination with records of estuary entrance condition, are used to identify the geomorphic variables that are key controls on the diversity in entrance functioning. A classification model of IOCE types is developed using 35 IOCE in Victoria, Australia, with detailed records of entrance condition, and then applied to all estuaries in Victoria. Critically, as IOCE are currently managed as a single estuarine type, it is important to determine whether they are simply an end member state or a diverse form of estuary with a range of sub-types. The geomorphic based classification model presented in this study provides an important tool for predicting not only the typical frequency and duration of entrance closures, but also which estuaries are at the most at risk of declining water quality, circulation and ecological condition based on their typical closure regime. It is also an important tool that is applicable on other wave-dominated coastlines internationally where IOCE are present.
II Methods
This study maps the morphology of all of the 111 estuaries along the coastline of Victoria. Several catchment-scale geomorphic variables were measured at each site using remotely sensed data, including: entrance channel width at the mouth, estuary perimeter, estuary surface water area, catchment area, estuary length, beach width and tidal prism (Table 1). A subset of 35 IOCE with detailed entrance condition data were selected for statistical analysis. For these 35 estuaries, the beach berm elevation was surveyed in the field. The characteristics of these 35 estuaries were then used to develop a classification of IOCE based on a statistical analysis of the catchment-scale controls of the diversity in morphology and entrance function. This model was then applied to all 111 estuaries in Victoria.
Definition of parameters and methodology of measurement for data used to classify estuaries in Victoria. For all estuaries in Victoria, the entrance channel width, estuary perimeter, estuary length, beach width, estuary surface water area, catchment area and tidal prism were measured using remote sensing methods. Mean annual rainfall was also included. For a subset of 35 IOCE with entrance condition data, the entrance channel and berm elevation were also measured by field surveying.
*Lake Tyers: http://www.laketyersbeach.net.au/openings.html
Using aerial photography from 1946–2010 from the Victorian State Library, supplemented with digital globes (GoogleEarth and DigiGlobe), the entrance condition of all estuaries in Victoria was firstly classified as being either (a) permanently open or (b) intermittently closed. An estuary was classified as being open when there was a free and ongoing exchange of sediment and water between the estuarine lagoon and the ocean, and closed when a continuous subaerial berm was present at the entrance. Any estuary that was observed to have been both closed and opened at any time in the past 70 years was classified as an IOCE, as per the temporal extent of aerial imagery. Google Earth imagery was available from 1984–2016 for all estuaries in Victoria, and had an average time interval between images of six months. The total amount of historic Google Earth images available for all estuaries was 15 on average (± 8 images).
For each estuary in Victoria, measurements of the entrance channel width at the mouth, estuary perimeter, estuary surface water area, catchment area, estuary length, beach width and tidal prism were undertaken using airborne light detection and ranging (LiDAR) data, which was collected in 2007–2009 by the Victorian State Government (Table 1; Figure 1). The LiDAR survey was conducted using a LADS Mk II system coupled with a GEC-Marconi FIN3110 inertial motion sensing system and a dual frequency kinematic global positioning system (kGPS). This dataset was processed to produce a seamless terrestrial–marine 2.50 m resolution mosaic from elevations of +10 m to depths of −25 m (Quadros and Rigby, 2010) and a 1 m resolution dataset for the upper catchment with > 4 points per m2. Estuary tidal prisms were calculated as per Barber (2003) with the spring tidal range taken from the Bureau of Meteorology (BOM). Obtaining field measurements of tidal current velocity and cross-sectional area for all sites state-wide was outside the scope of this study. The river flow per tidal cycle was added to the tidal prism for permanently open estuaries with large rivers (as per Ranasinghe et al., 2013), however for IOCE, due to the absence of flow gauging stations on 75% of sites, the tidal prism provides a minimum estimate as river flow was not added. Mean annual rainfall was taken from BOM with all stations having > 40 years of data and being within 10 km of each estuary. Mean annual rainfall was then analyzed relative to catchment area in order to provide an indication of the likely inflow at each estuary due to many sites not having flow gauging data.

Example measurements of: beach width, estuary entrance channel width, berm elevation, estuary perimeter, length and surface water area using Moggs Creek, Victoria. (a) LiDAR and (b) aerial imagery view.
The entrances of a subset of 35 estuaries were mapped in detail and surveyed in the field between 2012 and 2015 using a Sokkia S5X total station and Trimble R6 Real Time Kinematic (RTK) GPS. These were estuaries where entrance condition monitoring was undertaken by Estuary Watch Victoria (spanning six years on average, range 2–11 years). The entrance monitoring data included a record of the mouth status (open or closed) with supporting photographs for each time the site was visited and was supplemented with field observations during 2012–2015 and any available local databases (Table 1). The closure and opening duration of each site was obtained from taking the duration (in days) of each individual opening or closure event. Durations were totaled over the whole monitoring record for each site and then calculated as mean values. Opening and closure frequency were also calculated as an average over the extent of the monitoring period. The closure duration was recognized as the full period in which the estuary was separated from the ocean by a subaerial berm with the frequency being each individual instance of closure/opening. It is recognized that the infrequency (i.e. monthly vs weekly monitoring) of entrance monitoring at some sites may provide a minimum estimation of closure and opening frequency.
Statistical analysis of all geomorphic variables and entrance condition data was undertaken in IBM SPSS Statistics 23 and MATLAB R2015b. For the subset of 35 IOCE with records of entrance condition, all variables were tested for normality using the Kolmogorov-Smirnov coefficient and were below the threshold for normal distribution. Therefore, non-parametric methods were used to assess the relationships between variables and to classify IOCE into sub-types. A Pearson correlation analysis was performed to identify the correlations between catchment-scale variables and their influence on estuary entrance condition using the entrance condition data obtained from Estuary Watch Victoria (Table 1). Non-metric multidimensional scaling (nMDS) and hierarchal cluster analysis were undertaken using methods of Everitt and Dunn (1991); Clarke and Warwick (1994) and Saintilan (2004), to identify possible groupings between IOCE. Hierarchal cluster analysis was performed following nMDS to confirm the groupings identified using nMDS and to further identify any outliers. Following the groupings of IOCE which became apparent during nMDS and hierarchal cluster analysis, the distribution of all catchment-scale variables and entrance condition data (Table 1) were presented as boxplots for each group with descriptive statistics calculated within each group.
Multivariable linear regression using a stepwise backwards selection criteria was used to produce an equation to predict the likely percentage of time (PTO) an IOCE will exist in an open state annually. Data was log transformed to ensure a normal distribution and all measured geomorphic parameters (Table 1) for the subset of 35 IOCE were input into the multivariable linear regression to identify the strongest predictors of the independent variable (percentage of time open annually).
III Regional setting
The coast of Victoria, Australia, is 1700 km long and is divided into estuaries to the west and east of Wilson’s Promontory (Figure 2). On the western coast, mean annual significant wave height is 2.49 m with a mean annual maximum wave height of 4.13 m, mean annual wave period of 7.65 s and direction of 212.80° at Port Campbell (Hughes and Heap, 2010; Water Technology, 2004). The spring tidal range at Lorne, in the central western coast, is 1.60 m, while in the far western coast, the spring tidal range is 0.90 m at Portland. In west Victoria, south and southwesterly winds are dominant for most of the year; however, during spring and summer, winds from the south and southeast increase in frequency. This results in a predominantly easterly longshore drift for about nine months annually with a westerly occurring drift over summer (Barton and Sherwood, 2004). In eastern Victoria, (Figure 2), the mean annual significant wave height is 1.03 m with a mean maximum wave height of 1.73 m and period of 8–10 s (Short, 1996). The dominant direction of wave approach is from the south and southeast (Water Technology, 2004) with some influence from the northeast associated with east coast lows (25% of annual swell frequency) (Bird, 1993; Short, 1996). This results in a dominant longshore drift towards the northeast however the direction of drifting is reversed when north-easterly winds dominate (Bird, 1961; Short, 2010). The spring tidal range at Lakes Entrance, in the central eastern coast is 0.95 m. Wave heights peak in June to August along both coastlines.

Map of the coast of Victoria with the location of 35 case study IOCE marked. Mean spring tidal range and the dominant swell direction are shown. Major towns are indicated by dots and with names in italics.
The coastal climate of Victoria is temperate with median annual precipitation ranging between 900–1500 mm and peak rainfall occurring in the mountainous regions, particularly in the Otway Ranges (> 1400 mm/yr). Rainfall is highest in late winter to early spring and lowest in summer, and the river flow of catchments draining the coast also exhibit high inter-annual variability (Peel et al., 2001). The current coastline is tectonically stable (Murray-Wallace and Belperio, 1991) and the Holocene sea level is considered to be primarily driven by eustatic processes, with the Holocene highstand being ∼2.30 m above present sea level (Bryant, 1992).
IV Results
4.1 Geomorphic diversity of estuaries in Victoria
There are 111 estuaries in Victoria of which 103 (93%) are IOCE, having been both closed and open to the ocean at least once between 1946–2016. IOCE in Victoria typically have a shallow (< 2 m deep) lagoon impounded by a subaerial berm (Figure 3). The catchments behind IOCE are on average, < 1000 km2, have entrance channels < 500 m wide and small backing estuarine dimensions (e.g. lengths < 10 km, perimeters < 10 km, surface water areas < 10 km2; Table 2). The small lagoon size leads to a mean tidal prism for these estuaries of 5.01 × 106 m3 (Table 2).

A range of estuaries in Victoria, Australia, showing differences in morphology. (a) Barwon River; (b) Andersons Inlet; (c) Wingan Inlet; (d) Tamboon Inlet; (e) Fitzroy River; (f) Darby River. Photos courtesy of Neville Rosengren.
Descriptive statistics for geomorphic parameters of all estuaries in Victoria (n = 111), all permanently open estuaries (non-modified entrances) (n = 4) and IOCE (n = 103). Parameters measured using remotely sensed data with mean annual rainfall from BOM (2016). Data parameters as per Table 1 and include: (CW) entrance channel width (m), (CA) catchment area (km2), (EP) perimeter (km), (EL) estuary length (km), (SWA) estuary surface water area (km2), (BW) beach width (m), mean annual rainfall (MAR), (TP) tidal prism (x 106 m3).
On the west coast, there are 46 IOCE, the largest being Curdies Inlet with a channel width of 314 m and a catchment area of 1110 km2. The estuary has a large central basin (12 km2) which is up to 2 m deep. In contrast, the smallest IOCE is Moggs Creek (channel width of 34 m, surface water area < 1 km2). This estuary is found at the base of a steep catchment of 14 km2, which forms in the Otways Ranges and exits onto the beach via a channel < 0.25 m deep. The Fitzroy River could be considered to have a typical morphology for a west coast IOCE (Figure 3(e)). It has an entrance channel width of 160 m, a catchment area of 260 km2 and an estuarine lagoon of 1.20 km2 which is, on average, 1.25 m deep.
On the east coast, 57 IOCE are present with a similar range of dimensions to those found on the west coast (Table 2). The largest IOCE is Sydenham Inlet with a channel width of 600 m, perimeter of 22 km and a surface water area of 19 km2. The estuary has a large central basin (22 km2) which is up to 3 m deep. The smallest IOCE in east Victoria is Davis Creek, which has an entrance channel width of 35 m, a catchment area of 15 km2 and a surface water area < 1 km2. There are numerous small IOCE (40% of all east Victorian IOCE) around the granitic Wilsons Promontory, occupying steep, small catchments with areas < 30 km2, channel widths < 50 m and surface water areas < 1 km2 (e.g. Figure 3(f)).
Of the eight permanently open systems present on the open coast of Victoria, four estuaries had unmodified entrances (e.g. Barwon River, Andersons Inlet, Shallow Inlet, Mallacoota Inlet) with the other half being kept open through human intervention (e.g. entrance training walls at Moyne River). Permanently open estuaries tend to have large catchments (> 1000 km2) with entrance channel widths (> 450 m) and basin dimensions all above the mean values for estuaries in Victoria (Table 2). For example, the Barwon River in west Victoria (Figure 3(a)) has a catchment area of 8590 km2, a length of 19 km, a depth in the central basin of > 3 m and an entrance channel 480 m wide. Similarly, Andersons Inlet in east Victoria (Figure 3(b)) has a length of 17.65 km, a channel width of 2 km, catchment area of 1600 km2 and a surface water area of 13 km2 (Table 2). Permanently open estuaries in Victoria without entrance modification had tidal prisms ranging between 23.40 and 29.80 × 106 m3 (Table 2) and tended to have mouths which were orientated away from the dominant wave direction (Table 3).
Entrance configuration and entrance condition of a subsample of 35 IOCE and four permanently open estuaries without entrance modification relative to the dominant direction of wave approach on both the east and west coasts of Victoria. PTO refers to the percentage of time open annually.
4.2 Controls on entrance condition and estuary morphology
A continuum of catchment, channel and estuary basin sizes is evident across the 35 field-surveyed IOCE in Victoria showing that the estuary dimensions of these IOCE are positively scaled to catchment area and channel width (Figure 4; Table 4; Table 5). Using a Pearson correlation analysis, the relationship between the catchment-scale variables measured for the 35 field-surveyed IOCE were tested for correlations to determine relationships between the dimensions of the estuary basin, channel and catchment and their influence on entrance functioning. The highest positive correlations occur between tidal prism and surface water area (p = 0.981, γ = 0.01 level), estuary length and perimeter (p = 0.923, γ = 0.01 level), estuary perimeter and channel width (p = 0.921, γ = 0.01 level), estuary perimeter and surface water area (p = 0.906, γ = 0.01 level), perimeter and tidal prism (p = 0.902, γ = 0.01 level), and catchment area and tidal prism (p = 0.901, γ = 0.01 level) (Figure 4; Table 5). Berm height also showed a positive relationship with all other geomorphic variables (p = 0.425– 0.555, γ = 0.01 level) (Table 5).

Linear relationships of parameters influencing morphology and entrance condition of IOCE in Victoria: (a) catchment area and surface water area; (b) catchment area and channel width; (c) perimeter and channel width; (d) estuary length and channel width; (e) berm elevation and channel width; (f) channel width and tidal prism; (g) channel width and percentage of time open annually; (h) mean annual closure frequency and channel width; (i) mean annual rainfall: catchment area (proxy for fluvial inflow scaled per catchment) and channel width. A sample size of 35 case study IOCE is presented as these sites had detailed records of entrance condition. The boxplot distributions on the x and y axes show the distribution of parameters.
Descriptive statistics for geomorphic and entrance parameters of the 35 case study IOCE sites in Victoria. Data is derived from sites which were visited and surveyed as part of field work, and with entrance condition monitoring data. Entrance condition monitoring recorded by Estuary Watch Victoria online database system and field observations.
Pearson correlation analysis of geomorphic parameters and their relation to IOCE morphology and entrance condition. Parameters as per Table 1 and data for 35 case study IOCE with detailed entrance condition data (n = 35 for each parameter).
*Correlation is significant at the 0.05 level (2-tailed).
**Correlation is significant at the 0.01 level (2-tailed).
4.3 Dynamics of estuary mouths
As channel width at the estuary entrance has a consistently strong positive relationship with all other geomorphic variables, this was used to illustrate the influence of all other scaled geomorphic parameters in relation to entrance functioning (Figure 4(a)–(i)). The 35 field-surveyed IOCE closed between 1 and 10 times annually (mean 2.70 ± 2.40) (Table 4) with the lowest closure frequency occurring on smaller estuaries (channel widths < 75 m; catchments < 75 km2) (e.g. Kennett River, Skenes Creek and Wye River). The mean closure frequency showed a negative relationship with all geomorphic parameters aside from berm height (p = 0.571, γ = 0.01 level) (Table 5). The total percentage of time open per year was highly variable ranging between 18% and 98% (mean 61% ± 31%) with the highest percentage of time open occurring at Wild Dog Creek and Cumberland River, both with small channel, basin and catchment dimensions (Figure 4 (g)). The lowest proportion of time open was at Curdies Inlet, the largest IOCE in western Victoria (18%). Surprisingly, the percentage of time open annually showed a strong negative relationship with all geomorphic variables (Table 5). The mean closure duration was also variable and ranged between one and 372 days (mean of 64 days). Closure duration showed a strong positive relationship with all geomorphic variables (p = 0.768–0.634, γ = 0.01) (Table 5). The longest closure duration occurred at larger IOCE such as Curdies Inlet, while IOCE with the smallest entrance channels and dimensions have the shortest closure durations, that persisted for days on average. Additionally, the mean duration of each entrance opening was 184.90 ± 126.60 days with the maximum occurring at both the larger and smaller end member IOCE (Table 4).
All surveyed estuary mouths had a broad U-shaped cross section. At the mouth, the base of the widest channels (> 200 m) was < 0.50 m below MSL (Figure 5a). In long profile, the berm crest elevation was +0.5–1 m above MSL with a landward decrease in elevation to below MSL upon entering the lagoon (Figure 5b). In channels between 75 and 200 m wide, the channel bed depth at the mouth was generally above MSL. Berms were steep rising to > +1 m above MSL (Figure 5b). The base of the narrowest channels (< 75 m) were generally the deepest, incising to < −1 m below MSL where the thalweg was present (Figure 5c). In long profile, a distinct berm was not evident in channels of < 75 m size (Figure 5c). The long profile morphology may, however, be more a function of the entrance conditions at the time of surveying because all IOCE with channel widths < 75 m were open and surveyed in their dominant entrance state. At the time of surveying, all other IOCE were closed, so closed profiles are used for comparison.

Field surveys of IOCE entrance morphologies (cross section and long profile) for typical examples of IOCE in Victoria. (a) Examples of estuaries with channel widths > 200 m (Group 1), (b) Examples estuaries with channel widths > 75–200 m (Group 2) (c) Examples of estuaries with channel widths < 75 m (Group 3). Cross sections taken at the seaward extent of the lagoon for estuaries closed at the time of surveying (a), (b) and across the berm for estuaries surveyed in an open state (c). Long profiles taken from swash zone, over the berm and landward into the lagoon.
4.4 Classification of IOCE in Victoria
Using a non-metric multidimensional scaling analysis (nMDS), three distinct groupings of IOCE are apparent (Figure 6). The 2D stress factor for nMDS was 0.01 and a value of < 0.1 indicates a significant relationship between variables (Clarke and Warwick, 1994). The first group includes the larger estuaries such as Curdies Inlet and Lake Tyers, which have entrance channel widths > 200 m, catchment area > 500 km2 and tidal prism > 10 × 106 m3 (Figure 6, Group 1). The entrances of this group are open less frequently, but for longer periods of time. In the second group, which includes estuaries of mid-size range (75–200 m channel width, 75–1000 km2 catchment area and tidal prism 1.50–10 × 106 m3), the entrance condition is highly variable. Entrances are open for 19–90% of the time annually with 1–10 openings per year. Examples in Victoria include Merriman Creek and Aire River (Figure 6, Group 2). Within this group, the nMDS plot shows that the Gellibrand River could be considered an outlier. This is due to the larger basin and catchment dimensions (i.e. catchment > 1000 km2, perimeter > 10 km and length > 10 km) but lower surface water area as attributed to the potentially higher level of sedimentary infill present. The final group is most tightly clustered and includes the smallest IOCE with entrance channel widths < 75 m, catchments < 75 km2 and tidal prisms < 1.5 × 106 m3 (Figure 6, Group 3).

nMDS plot showing distinct groups of IOCE in terms of channel width, catchment area, perimeter, length and surface water area and entrance functioning. For nMDS analysis, 35 IOCE were used with detailed records of entrance condition data. Legend shows IOCE site. Group 1 (largest sized) have channel widths > 200 m (black dots), Group 2 are > 75–200 m channel widths (grey dots), Group 3 (smallest size) have channel widths < 75 m (open dots).
Hierarchical cluster analysis produces similar results to nMDS (Figure 7). The larger estuaries of Sydenham, Wingan, Tamboon and Curdies Inlets and Lakes Tyers and Bunga are clustered together (Figure 7, Group 1). Group 1 estuaries can be further subdivided between Lakes Bunga and Tyers, Curdies Inlet, Wingan Inlet and Sydenham and Tamboon Inlets as based on the differences in basin and catchment size (i.e. Sydenham and Tamboon Inlets are most similar in surface water area and have the largest catchments > 1100 km2; Figure 7). As per nMDS analysis, the Gellibrand River is considered an outlier. Relative to all other IOCE, however, dimensions of all Group 1 estuaries were an order of magnitude larger than the smallest IOCE (Group 3) and > 5 km2 larger than the next highest value in Group 2. The next most distinct cluster (Group 2: groups B1, B2, B3) subdivide between Thurra River and Sydenham Inlet (Figure 7).

Dendrogram using average linkage (between groups) showing hierarchal analysis for case study sites in Victoria, Australia. Study sites falling within each class are listed. The dendrogram is created using the 35 IOCE with detailed records of entrance condition data as per nMDS analysis.
Group 2 represents the mid-sized IOCE which have channel widths 75–200 m and the most variable entrance condition. Within Group 2 there is a secondary divergence, which is likely attributed to the varying degrees of basin accommodation space associated with these IOCE, ranging between the two end member types. Group B3 (Group 3) represents the smaller IOCE with channel widths < 75 m (Figure 7). Compared to nMDS results, the dendrogram shows more overlap between groups. This is explained by the split between Mueller River (B2) and Shipwreck Creek (B3) being predominantly based on estuary and catchment size, as opposed to recognizing the different entrance conditions and topographic setting that may otherwise differentiate between these estuaries. For example, Cumberland, Wye, Kennett and St George rivers are open for > 90% of the time annually and have limited to no central basin accommodation space, but have catchments ranging between 20 and 36 km2. Barnham River and Aireys Inlet, also included in Group B3, are both extensively infilled and open for < 60% of time annually; however, because they have catchments 79–81 km2, they too classify as Group 3 IOCE. As a result, Group B3 is included in Group 3 in the nMDS plot, the smaller end member of all IOCE but with some overlap as hierarchal analysis does not account for infill and topographic setting, and appears to prioritize spatial scale over the dominant entrance condition.
The distribution of variables within these three groupings of IOCE further confirms the nMDS and hierarchal cluster results with distinct differences between the distribution and median values for each group (Figure 8a–l). It is evident that Group 1 consistently show the largest channel, catchment and basin dimensions, Group 2 occupy the mid-range of scale but have the most variable range of values, and Group 3 are the smallest, with the most tightly distributed values on average (Figure 8a–k). When considering the impact of these groupings on entrance condition, Group 1 close to the ocean less often; however, entrance closure and opening both persist for months at a time (Figure 8i–l). Group 2 are highly variable in terms of both opening and closure frequency and duration, and Group 3 are open for the majority of the time, with closures only persisting for days (Figure 8i–l).

Box plot distribution of geomorphic parameters and measurements of entrance opening and closure regime for the three sub-groups of IOCE: Type A, B and C as corresponding with the groups identified in nMDS and statistical analysis. (a) Variation in channel width; (b) catchment area; (c) perimeter; (d) length; (e) surface water area; (f) tidal prism; (g) berm elevation; (h) mean annual rainfall; (i) mean closure duration; (j) percentage of time open annually; (k) mean closure duration; (l) mean opening duration. Data is from 35 case study IOCE with detailed entrance condition records.
4.5 Predicting entrance condition in IOCE
Multivariable linear regression was used to produce an equation to predict the likely percentage of time an IOCE will exist in an open state annually prior to any human modification to the catchment or entrance morphology (equation (1)). Channel width (CW), catchment area (CA) and surface water area (SWA) were found to be the strongest drivers of the annual percentage of time open (PTO) in IOCE. Multivariate analysis had a standard error of +0.19, degrees of freedom (df) value of 30 and a R Square value of 0.61 (i.e. 61% of variation in the response variable: PTO is predicted using these parameters). PTO can therefore be accounted for using the three predictor values (CW, CA and SWA). As a result, equation (1) is used to model PTO in IOCE in Victoria where CW is channel width, CA is catchment area and SWA is surface water area
Channel width, as a predictive parameter represents the capacity of the entrance channel to facilitate outflow and exchange with the ocean. Surface water area relative to catchment area is effectively a proxy for the amount of basin accommodation space along with representing the size of the backing lagoon as a whole. Catchment area provides an approximation of how much fluvial input into the estuary basin occurs in terms of runoff, discharge and groundwater flow.
V Discussion
5.1 Classification model of IOCE in Victoria
93% of estuaries in Victoria close intermittently; therefore, this is a fundamental property of Victorian estuaries. Of the small number of permanently open systems, most remain open due to human intervention through features designed to counteract infilling of the entrance channel. The larger tidal prism and entrance channels at those sites which remain permanently open without entrance modification reflects the well-known AP empirical stability relationship where a larger entrance cross-sectional area (A) and tidal prism (P) corresponds with a more stable channel (Hume and Herdendorf, 1993; Jarrett, 1976; O’Brien, 1931). As wave energy along the Victorian coast is persistently high due to the influence of uninhibited swell propagation from the Southern Ocean, littoral drift would also be high by global standards (Bruun, 1986; Davis, 1989). For example, median annual wave power exceeds 30 kW/m in west Victoria with a mean annual Hs > 2 m (Hughes and Heap, 2010; Water Technology, 2004). While these sites may experience similar wave conditions to IOCE, the larger tidal prisms mean that a higher proportion of sediment can be removed from the entrance channel due to a larger volume of tidal exchange and higher ebb-tidal velocities (Bruun, 1978; Bruun and Gerritsen, 1960). The larger catchments also add more river flow, thus strengthening the ebb-tidal prism (Table 2). In Victoria, these sites also tend to have entrances which do not directly face the dominant direction of wave approach, meaning local topographic features may shelter these sites from the direct influence of waves (Table 3).
Within IOCE in Victoria, the continuum of entrance channel width, basin dimensions (perimeter, basin surface water area, estuary length), tidal prism and catchment size can be used as a basis for classifying IOCE (Table 6; Figure 9). Statistical analysis indicates that three clear sub-types of IOCE can be identified in Victoria: Type A – large IOCE (corresponding to Group 1); Type B – medium sized IOCE (Group 2); and Type C – small, tidal creeks (Group 3) (Figure 6; Figure 7; Figure 9). As the variability in these geomorphic parameters in turn controls entrance functioning, each sub-type of IOCE exhibits a different characteristic opening and closure regime. For example, there is at least one order of magnitude difference between mean closure duration in the three types of IOCE.
Parameters evident in Victoria for conceptual classification of IOCE groups. The top line for each group indicates the thresholds for the classification for estuaries in Victoria. Parameters are taken from 35 case study IOCE sites and are as defined in Table 1 including: entrance channel width (CW), catchment area (CA), estuary surface water area (SWA), estuary perimeter (km) (EP), estuary length (EL), tidal prism (TP), mean annual rainfall (MAR), mean closure duration (MCD) and percentage of time open annually (PTO).

Geomorphic based classification model of wave-dominated microtidal estuaries in Victoria.
5.2 Type A – large IOCE
Type A IOCE exhibit the largest estuary dimensions in terms of channel width (> 200 m), estuary water surface area (> 5 km2), length (> 10 km), perimeter (> 12 km), catchment size (> 500 km2) and tidal prism (> 10 × 106 m3) (Table 6; Figure 8; Figure 10). This increase in catchment size, basin water surface area and tidal prism (Figure 8) is consistent with past literature (Haines et al., 2006; Mondon et al., 2003; West et al., 1985). Type A estuaries form on low-gradient coastal plains and have a basin morphology which tends to be circular, with a large central basin and ample accommodation space. They have a distinct beach berm at an elevation (when closed) of +0.50–1.50 m MSL (Figure 5a).

Examples of typical Type A (large) IOCE (a) Lake Tyers, east Victoria; (b) Curdies Inlet, west Victoria; (c) Schematic image of characteristic Type A IOCE morphology.
Type A IOCE experience few openings per year (1–2 per year) being open on average for 32% of the time annually; however, they are distinguished by long periods of both opening and closing (Table 6). When entrance openings occur, they persist for more than a month, while closures typically last for months to years at a time. For example, on average, Curdies Inlet was closed for just 18% of the time annually, but had a mean opening duration of 200 days and mean closure duration of 210 days. The monthly scale duration of entrance opening is due to the larger tidal prism and increased hydraulic efficiency through a larger entrance channel (Bruun, 1978; Gao and Collins, 1994; O’Brien, 1931). Evidence of this is the deposition of prominent flood tidal deltas indicating periods of sustained tidal exchange during past openings (Figure 10a). As the lagoon drains, and in the absence of sufficient fluvial flow to add to the ebb-tidal prism, the high-wave energy of the Victorian coastline is likely to close these larger IOCE faster than if they were located on coasts with a less energetic wave climate and lower drift rate. The frequency and duration of entrance opening and closure is therefore not only dictated by contemporary processes such as the relative balance of wave and fluvial energy, but also largely by the underlying basin topography and relative accommodation space. This is a fundamental scaling feature of IOCE in Victoria as, while the marine conditions are similar across the state, those estuaries with larger catchments and lagoons can maintain an open entrance for longer due to the larger tidal prism. Once isolated from the ocean and tidal influence, however, the larger basin volume takes longer to fill, until overtopping from the catchment side reinitiates opening.
In Australia, Type A IOCE are a counterpart of the large, immature ‘saline coastal lagoons’ found in New South Wales (NSW) (Roy, 1984; Roy et al., 2001 – Type IV 8). These sites have similar morphologies; for example, Smiths Lake (extensive tidal delta deposits, channel width > 200 m, surface water area 10 km2) and Coila Lake (surface water area 7 km2, closed 95% of time) (Roy et al., 2001). A fundamental difference between these sites and Victorian Type A IOCE is that in NSW these IOCE showed large basin and channel dimensions, but often with small catchments; for example, Smiths Lake (33 km2) and Coila Lake (52 km2). As a larger catchment area has a positive influence on the percentage of time open annually (equation (1)), Type A IOCE with small catchments would be expected to be closed for a higher proportion of the time compared to those with larger catchments, as their basin size is similar but there is less inflow to the estuary. This underpins the role of estuary basin accommodation space as a key control on opening frequency and duration. Internationally, Type A IOCE may also be considered a counterpart of the ‘Waituna’ coastal lakes found in New Zealand, such as Waituna Lagoon and Lakes Ellsemere and Forsyth, but forming on sandy beaches. Waituna have large surface water areas (> 5 km2), are typically closed for months to years at a time and are associated with low river flows relative to catchment size (Kirk and Lauder, 2000; NIWA, 2016; Schallenberg et al., 2010). The large ‘coastal lagoons’ of northern California are another international form of Type A IOCE, including Stone Lagoon and Lake Earl (Kraus et al., 2008). These sites have catchments > 500 km2, entrance channels > 200 m, tidal prisms > 5 × 106 m3 and close every 1–2 years for months to years at a time (Kraus et al., 2002, 2008). Consistent with Victoria, both the national and international counterparts of Type A IOCE exist in a predominantly closed entrance state, with closures persisting for months to years.
5.3 Type B – medium sized IOCE
Type B IOCE occupy the mid-spectrum in terms of size, and are smaller than Type A, with channel widths 75–200 m, tidal prisms 1.50–10 × 106 m3, catchment areas 75–1000 km2, estuary surface water areas 1–5 km2, lengths 2–10 km and perimeters 2–12 km (Table 6; Figure 11). The catchment and estuary dimensions of Type B IOCE showed more variance than other types of IOCE, likely due to the wider range of coastal and topographic settings that they occupy and their varying degrees of sedimentary infill. For example, the Gellibrand River has a catchment area > 1000 km2, a channel width of 120 m, a channelized (< 3 m depth) central basin with extensive infill, and a widespread fringing swamp (Barton and Sherwood, 2004). The Anglesea River, in contrast, has a catchment area of 125 km2 and an entrance channel of 90 m with a shallow (< 1.50 m) central basin. Throughout the central basin of Type B IOCE, width increases laterally, but then becomes narrow again at the head of the estuary. Type B IOCE exhibit a steep beach berm, but the range of berm heights is variable due to the variety of beach types forming at the entrance. For example, both Sherbrook River and Spring Creek in west Victoria occupy rocky embayments and have berm heights of +1.50–3.15 m above MSL. In contrast, the Thurra, Aire and Hopkins estuaries enter the open coast and have berm heights of +2.20–3.05 m above MSL. Berms increase in elevation rapidly over the first 20 m landward of the swash zone where after 20 m the channel bed elevation only decreases slightly (±0.30 m) to indicate further sediment storage behind the berm (Figure 5b).

Typical examples of Type B (medium sized) IOCE (a) Hopkins River, west Victoria; (b) Gellibrand River, west Victoria; (c) Schematic image of characteristic Type B IOCE morphology.
Type B IOCE entrances open at a weekly to monthly frequency (1–10 times/year) (Table 6). Entrances stay open for less time than Type A IOCE, typically for weeks to months. For example, the Aire River opens on average four times per year for a mean duration of 21 days, while the Anglesea River opens five times per year with a mean opening duration of 14 days. Entrance closures are more frequent than in Type A IOCE, but of a shorter duration as the lower volume of flood accommodation space requires less freshwater input before the berm is overtopped from the catchment side. Entrance openings also last for weeks to months as the smaller tidal prism means that marine depositional processes can dominate more rapidly in the absence of a rapid influx of river flow.
In Australia, Type B IOCE are a counterpart of a semi-mature “saline coastal lake” where extensive infill of the central basin is evident (Roy, 1984). Internationally, they are similar to the “perched, normally closed estuaries” defined by Cooper (2001), which have steep berms, an estuary bed maintained at or above MSL at the mouth, channels ∼100 m wide and openings that typically persist for weeks at a time. The Mhlanga, Mdloti, Groot Brak and East Kleinemonde rivers in South Africa are examples of estuaries which can be classified as Type B IOCE. These sites have catchments ranging between 43 and 550 km2, basin surface water areas of 0.50–3 km2 and entrance channels 100–200 m wide at the mouth (Cooper, 2001; Nunes and Adams, 2014; Riddin and Adams, 2008; Whitfield and Bate, 2007). These sites open between 1 - 6 times per year for a duration of weeks - days on average (Anandraj et al., 2007; Mundree et al., 2003; Wooldridge, 1994). In the USA, Type B IOCE are similar in morphology and entrance functioning to the larger “bar-built” estuaries of California (Kraus et al., 2008). Examples include the Russian River and Redwood Creek, which also have high steep berms (> +2 m MSL), similar catchment and estuary dimensions, and close several times annually for weeks to months at a time (Behrens et al., 2013).
5.4 Type C – small tidal creeks
Type C IOCE exhibit the smallest estuary dimensions with channel widths < 75 m, tidal prisms < 1 × 106 m3, surface water areas < 1 km2, lengths < 2 km, perimeters < 3 km and catchments < 75 km2 (Table 6; Figure 12). The typical entrance morphology of Type C IOCE shows a low berm (all < +0.50 m above MSL) and smaller cross-sectional areas at the mouth (Figure 5c). Type C IOCE exhibit a narrow and meandering channel up to 200 m upstream of the mouth which, upon flowing onto the beach, becomes shallow (< 0.30 m) and fan-shaped. Tidal terraces are evident along the lower reaches of the channel, which indicate the tidal flow regularly enlarging the channel under different tidal elevations (Figures 5c). There is a generally constant depth landward of the entrance, except for deeper sections where the channel is confined by areas of bedrock or boulders (Figure 5c). Type C IOCE occupy steep, mountainous catchments with high rainfall and where the topographic setting limits the size of the estuary itself with a smaller inherited river valley.

Typical examples of Type C (small tidal creek) IOCE (a) Darby River, east Victoria (channel width on beach 15 m); (b) Wye River, west Victoria (Image from Estuary Watch (2016) photo point site 1; (c) Schematic image of characteristic Type C IOCE morphology.
At all study sites, Type C IOCE were open for > 90% of the year (Figure 8j). Entrance closures were also infrequent (1.18 per/ year on average) and of a short duration, persisting only over a daily or tidal scale (Figure 8g–i; Table 6). Although Type C IOCE have small basin dimensions and tidal prisms, which has previously indicated an estuary that is more prone to closure (Bruun, 1978; Bruun and Gerritsen, 1960; Jarrett, 1976), their steep, mountainous catchments with high runoff rates means that they are maintained open by more constant fluvial discharge. Type C IOCE in Victoria are located in areas with high rainfall relative to catchment size, meaning that the more constant runoff is likely to dominate over marine processes of deposition under normal conditions (Figure 8i). For example, Type C IOCE on Wilson’s Promontory (e.g. Darby River) and the Otways Ranges (e.g. the Wye and Kennet rivers) receive the highest annual rainfall state-wide of > 1400 mm (Table 6). In comparison, mean annual rainfall at Type A IOCE is between 800 and 1000 mm and at Type B IOCE is between 600 and 1050 mm (Table 6). Type C IOCE are also located seaward of coastal ranges, with a more impermeable substrate, meaning rainfall is channeled directly downslope with little opportunity to disperse through the water table. For example, the underlying granitic rock on Wilsons Promontory assists runoff from the mountainous topography (Hill, 1994). As rainfall and river flow is a main control on entrance opening, any change in precipitation or catchment characteristics will be rapidly reflected in the entrance. Many Type C IOCE in Victoria also tend to enter the ocean alongside rocky headlands, which may act to promote an area of deeper scour of the channel thalweg alongside the headland, thus extending the entrance opening duration under normal wave conditions (Figure 12a–c).
In Victoria, Type C IOCE are a counterpart of the “intermittently open small coastal creeks” described by Roy et al. (2001) occurring in NSW (Type IV 9 in their classification). A key distinction is that many small coastal creeks in NSW are infilled (i.e. mature), whereas Type C IOCE in Victoria are small on account of the topographic constraints of their confined valley setting, which limits central basin development. Internationally, Type C IOCE are a smaller counterpart of the “normally open, river-dominated” estuaries of South Africa in the sense that they are open to the ocean for the majority of the time due to the influence of river flow, lack a distinct flood tidal delta and are found in steep catchments (Cooper, 2001). Again, a key distinction between normally open, river-dominated estuaries and Type C IOCE is that, in Victoria, the confined valley setting is the primary factor limiting the estuary size as opposed to progressive sedimentary infill (Cooper, 1994; Cooper, 2001). In southern California, the “steep coastal drainage (W-C)” coastal lagoons are similar in morphology to Type C IOCE (Jacobs et al., 2011). These estuaries have small (< 75 km2), high gradient catchments and drain the face of coastal ranges, with IOCE of the Santa Lucia Range (e.g. Toro Creek) being good examples (Bell et al., 2011; Jacobs et al., 2011).
5.5 Regional distribution of estuary types
Across Victoria, permanently open estuaries are rare and the three types of IOCE identified in the classification model are relatively evenly distributed along both the east and west coasts (Figure 13a and b). Type A IOCE are less common due to their characteristic large basin area whereby these systems require an extensive past valley space to occupy. They are also found only on open coastlines. Type B IOCE occupy the most diverse topographic settings. They are present along sheltered embayments and on open coasts. Type C IOCE are concentrated in areas of high gradient catchment topography and areas with narrow coastal plains. In Victoria, there are two regions with a high proportion of Type C IOCE: the Otways coast and Wilson’s Promontory.

(a) Location of all estuaries in Victoria, Australia including IOCE and permanently open types (both from human intervention to the mouth and naturally occurring). (b) Distribution map of IOCE and estuary types when applying the classification model to all estuaries in Victoria, Australia.
VI Conclusion
In Victoria, IOCE are the dominant type of estuary and constitute 93% of all estuaries along the 1700 km long open coastline. IOCE in Victoria exhibit considerable morphological diversity in terms of spatial scale, topographic setting and entrance functioning, thus comprising a suite of estuaries in their own right. This study presents a geomorphic classification of IOCE types that is both locally and internationally relevant as it enables any IOCE to be classified based on easily obtainable spatial and field data to predict entrance functioning. The delineation between IOCE types is an important advance from previous models which recognize IOCE as a homogenous estuary type despite clear differences in the frequency and duration of entrance closures. IOCE in Victoria can be classified into three distinct sub-groups: Type A (large IOCE); Type B (medium sized IOCE); and Type C (small, tidal creeks). These three types of IOCE showed an order of magnitude difference in entrance closure duration. The channel width, catchment area and lagoon size proved to be key geomorphic controls on opening and closure frequency and duration whereby a larger channel area, estuarine basin and tidal prism resulted in less frequent openings but both openings and closures of a longer duration.
Identifying the typical frequency and duration of entrance closures occurring in each type of IOCE provides a valuable asset for informing management decisions such as predicting the potential risk of each estuary type to the adverse effects of prolonged closure. By being able to predict the closure frequency and duration of IOCE, this information can be used to develop inventories of sites which would be at a high risk of water quality issues associated with longer water residence times, decreased fish passage to the ocean and flooding of surrounding land. An important consideration is that as many IOCE internationally exist in increasingly modified catchment settings where dams and irrigation practices reduce river flow into the estuary, a decrease in the current opening frequency and duration can be expected to occur in the future, with further implications for IOCE ecology and management.
Footnotes
Acknowledgements
Neville Rosengren is thanked for providing photos of estuaries and Estuary Watch Victoria is thanked for records of entrance condition. We thank the two anonymous reviewers for helpful comments on an earlier draft of this paper.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The Victorian Department of Environment, Land, Water & Planning provided financial support for this project and the first author was supported by an Australian Postgraduate Award and a Faculty of Science Postgraduate Writing up Award.
