Environmental fluid mechanics of minimum energy loss weirs: hydrodynamics and self-aeration at Chinchilla MEL weir during the November–December 2021 flood event

During the recent decades, a number of overflow embankment protection systems were implemented. One design is the Minimum Energy Loss (MEL) weir, developed to pass large flood events with minimum energy loss and low erosion. Several MEL weirs have successfully operated for decades in Australian catchments affected by heavy tropical and sub-tropical rainfalls with very flat gradients. Their historical performances are discussed and confirmed the design capability to pass large floods with small afflux and very small energy loss and little environmental impact including no significant erosion at the abutments. During a major flood, visual and quantitative observations were undertaken at the Chinchilla MEL weir. On the smooth converging chute, the observations showed that the inception of self‐aeration was a three‐dimensional process with a gradual change in free‐surface roughness, from a smooth glassy free-surface to a very-rough choppy surface. An optical technique (OF) was implemented to derive the contour maps of surface velocities based upon video movies taken from a sturdy tripod. The self-aerated flow velocity data were close to the backwater equation in the self‐aerated region, while highlighting regions of high- and low-velocities across the chute. The OF data showed large streamwise surface velocity fluctuations in the aerated flow region, consistent with the broad literature on self-aerated flow measurements. The ratio of transverse to streamwise surface turbulence intensity indicated a strong anisotropy of the free-surface turbulence. Minimum energy loss (MEL) weirs were developed for river catchment with heavy rainfall and very flat gradient. Field observations were performed at the Chinchilla MEL weir spillway in late 2021 during a major flood. Surface velocity data showed large transverse difference in longitudinal surface velocities across the chute Minimum energy loss (MEL) weirs were developed for river catchment with heavy rainfall and very flat gradient. Field observations were performed at the Chinchilla MEL weir spillway in late 2021 during a major flood. Surface velocity data showed large transverse difference in longitudinal surface velocities across the chute


Introduction
During the recent decades, a number of overflow embankment protection systems were developed and implemented [18,27]. These include the earth dam spillway with precast concrete blocks [28,48], some Roller Compacted Concrete (RCC) protection of the downstream slope of embankment dams [13,44] and the Minimum Energy Loss (MEL) weir [43]. The MEL weir structures are designed to pass large flood events with minimum energy loss and low erosion, and several MEL weirs have successfully operated for decades in Australian catchments affected by heavy tropical and sub-tropical rainfalls with very flat gradients, i.e. S o ~ 0.1% [43,55].
In most overflow spillways, the upstream flow region is brown and non-aerated, and a strong air-water mixture develops downstream of the onset region of air entrainment [30,37]. At the upstream end of the spillway, the flow is accelerated in a relatively smooth manner [29,36]. A turbulent boundary layer develops along the invert from the upstream end [32,34]. On a steep overflow spillway, the inception of free-surface aeration typically occurs when the outer edge of the bottom boundary layer starts to interfere with the water surface [31,36,50]. The interactions between boundary layer turbulence and free-surface can be explosive accompanied by strong projections of air-water packets [16,21,57]. In contrast, on a flat chute, the onset region exhibits a free surface with the gradual development of small ripple and surface instabilities, and the water surface tends to present a choppy appearance when self-aeration starts to occur [11,12,46]. Levi [40,41] and Anwar [1] argued that the inception conditions then combine the effects of longitudinal vortices' breakdown and of gravity waves produced by boundary irregularities and roughness. Recently, Toro et al. [54] and Chanson [19,21] demonstrated that the two theories are not exclusive using two different approaches. That is, threedimensional Computational Fluid Dynamics (CFD) computations [54] and detailed prototype observations [19,21] respectively. The inception region of free-surface aeration is illustrated on two prototype spillways in Fig. 1, using high-shutter speed (1/8000 s) photographs of the onset region in relatively steep and flat overflow chutes, i.e. Figure 1A, B respectively. The location of both structures is shown in Fig. 2.
To date, the literature on prototype spillway operation, including minimum energy loss (MEL) weirs, is very limited, with negligible quantitative field data. Herein, a series of field observations was conducted in a large protype MEL weir during a major flood event in late 2021. The contribution first aims to describe the process that led to the MEL weir design concept and to review the performance of MEL weirs taking into account recent observations of the effects of major flood events on the MEL weir at Chinchilla. The latter is complemented by an unique series of field observations including detailed quantitative measurements performed in terms of the surface velocity field in the highspeed, high-Reynolds-number overflow. Some key features of the inception region of freesurface aeration were described, showing the onset of air entrapment being a complicated process, while the surface velocity field differed from the assumed one-dimensional flow (1D) pattern.

Presentation
The Minimum Energy Loss (MEL) weir design approach was developed in response to the failure to achieve a satisfactory design of a "conventional" weir for a site on Sandy Creek near Clermont, Queensland (QLD), Northern Autralia. Through physical model tests of many conventional weirs, late Professor Gordon Reinecke ("Mac") McKay  progressively refined the weir design. In essence, his strategy had as its main features: (i) arrange the level of the crest in several steps, rising progressively on each side from the central section in order to reduce the difference between the flow distribution over the weir and that in the natural stream and thereby reduce the potential for bank erosion; (ii) reduce as much as possible the afflux across the weir when the upstream flood level overtopped the banks in order to reduce the potential scour at the abutments: (iii) design adequate scour protection for the embankments to prevent scour there. An example of such a design is the weir on Wararba Creek near Caboolture, Queensland, Eastern Australia (Fig. 3). However, eventually it proved impossible to develop a satisfactory design based on these principles for the site on Sandy Creek, despite a large number of attempts. The competing requirements for this weir and for many others constructed in streams with erodible banks are: (1) maximum in-stream storage, (2) protection against erosion of the stream banks at the abutments of the weir and downstream and (3) no increase in the frequency of outof-bank flooding. The conventional type of weir caused excessive afflux and over-bank flow and massive erosion of the alluvial banks at the abutments and downstream, even for flows considerably less than the design flow. When Apelt [4] was invited to observe the performance of the last of these unsatisfactory model tests, it was clear that the weir was causing too much constriction of the flow with large energy losses. He noted that the flow over the weir was analogous to that past an orifice plate in a pipe and suggested an approach analogous to a Venturi tube. The idea was to minimise the energy loss across the structure by using a more streamlined geometry everywhere and a gradual expansion downstream to recover as much as possible of the kinetic energy of the flow at the weir crest. A large number of trials finally produced a satisfactory design for the site on Sandy Creek, Clermont and the first MEL weir was built there in 1962 (Fig. 4) [4].
The essential principles of the approach, described as the Minimum Energy Loss (MEL) weir design, are simple: (i) the crest is made long enough to pass the bank-full flow at critical conditions with the simplifying assumption there is no energy loss and therefore no change in the upstream water level; (ii) the weir crest in plan is a circular arc of the length required by the first principle and concave downstream in order to converge the flow horizontally towards the centre of the stream after it has passed the crest; (iii) the faces of the weir must have relatively flat slopes to reduce energy losses by avoiding rapid lateral convergence horizontally and expansion of the flow in the vertical plane as it passes over the weir; the downstream face has a much flatter slope than does the upstream face. In order to satisfy the first two principles, the upper levels of the banks immediately upstream from the weir crest are excavated sufficiently to provide for efficient flow normal to the crest, uniformly distributed along its full length. Figure 5, based on the design drawings of the Chinchilla MEL weir, illustrates these principles.

Discussion
The MEL weir design presents a number of advantages. The weir passes the bank-full flow with little afflux. The term "afflux" is used here to describe the increase in the actual water level upstream from the weir above the undisturbed water level that would have existed in the absence of the weir (see discussion Sect. 2.3). The simple design approach implies that no energy loss occurs as the flow approaches and passes over the weir crest. Clearly, this is an oversimplification and, inevitably, some energy loss occurs though it is much smaller than that for a conventional weir. Consequently, the bank-full flow with the weir in place will be slightly smaller than for the undisturbed stream. As a consequence of the small energy loss, the crest level of the MEL weir can be much higher than would be feasible for a conventional weir for the same design flow and the water storage capacity is much increased. The horizontal convergence of the flow downstream from the crest results in a flow distribution much closer to that of the natural stream than is achievable with a conventional weir and the potential for bank erosion is much reduced.
There are however a number of disadvantages. Because of the flatter cross-sectional shape, the construction of a MEL weir involves a much larger volume of material than is needed for a conventional concrete gravity weir. Typically, the approach has been to construct the weir as a compacted earth embankment covered everywhere by concrete slabs with provision for drainage to relieve any water pressure that develops in the embankment. Inevitably, there is risk of damage should a significant flow occur during construction. The MEL weir on Sandy Creek at Clermont was built in 1962. When it was nearly completed a minor flood overtopped it and caused considerable damage. The weir was completed at the next attempt. The construction of the MEL weir on the Condamine River at Chinchilla began in late 1972 but was interrupted by a large flood in October 1972 that caused extensive damage when the weir was nearly completed. Unfortunately, when the renewed construction was almost finished, another flood in July 1973 caused extensive damage. Eventually the weir was completed successfully in December 1973 [55].
Whether an MEL weir is appropriate for a particular site will depend on the design criteria and on the relative costs-that must include some insurance against possible damage during construction if there is a risk of flooding.

Historical performance of MEL weirs
The weir at Clermont has been in operation for about 60 years, that at Chinchilla for nearly 50 years and a review of their performance is instructive. The significant criteria in such review are whether there have been issues concerning erosion and the hydraulic performance.
The MEL weir on Sandy Creek, Clermont was designed for a design head of 2.9 m. McKay [43] reported that the MEL weir "has been overflowed many times, but the size of the maximum flow is not known. No scour or erosion has taken place downstream". In Apelt [2], J.D. Turnbull reported that "no further problems have been reported from the weir" since completion of construction. The weir is still in use as an emergency water supply [14].
The Chinchilla minimum energy loss weir (QLD, Australia) (Figs. 1B, 5) is located in the Western Downs, along the Condamine River to provide both town and irrigation water supply [49]. The longitudinal river bed slope is S o = 0.215% in average between Brigalow and Chinchilla weir. Completed in 1973, the weir provides irrigation water, and it is listed as a "large dam". With a catchment area of 19,192 km 2 , the weir is a 14 m high earth fill embankment with a 410 m long dam crest including abutments. The overflow spillway (Fig. 1B) has a design capacity is of 850 m 3 /s corresponding to bank full at a design head 1.83 m. The spillway system consists of a broad crest that is 214 m long, followed by a smooth converging chute with a 1 V:5H slope. There is no stilling basin. After completion, grass was planted in the top soil and protective mesh that had been placed on the subsidiary earth embankments on each bank. Before this had time to grow, the weir was subjected to a major out-of-bank flood of order 1,130 m 3 /s, one of the highest on record, and that overtopped the embankments [55]. At the peak of the flood, it was difficult to see where the weir was, because it was completely drowned out. As the flood receded, it was wonderful to see the structure emerge from the flood waters relatively unscathed. The top soil and grass planted over the plastic mesh and gravel on the top of these embankments were washed away for the most part, but the mesh 1 3 held and the damage was minor. The soil and grass were replaced and the grass became well established before the next flood [2]. Since the weir's completion, the Condamine River at Chinchilla has carried a large number of floods in that period, many of them being major out-of-bank floods. Between 1973 and 2022, the weir was overtopped by a number of large flood events, including events larger than the design flow ( Table 1). The weir operated safely and properly, and inspections after the flood showed no damage. The first author has inspected the weir several times and found no evidence of erosion. Most recently he visited soon after a major out-of-bank flood in November/December 2021. The flood level had risen to at least 2.5 m above the right bank at the carpark, and no erosion had occurred.
The Chinchilla MEL weir was designed to give no afflux at bank full flow, 850 m 3 /s, when the head over the weir is only 1.83 m, [55]. Although the design of the weir itself was based on the assumption that no energy losses occurred, "the level of the supplementary embankments was determined from theory and the model as the lowest head water level at which there would be an afflux of less than one foot (300 mm)" [55]. The term "afflux" was used here to describe the increase in the actual water level upstream from the weir above the undisturbed water level that would have existed in the absence of the weir. Implicit in this definition is the assumption of steady flow in a uniform channel. During the major flood in 1973 of order 1,130 m 3 /s, the measured afflux at the bank full flow was only 100 mm, i.e. the head over the crest was 1.93 m, only slightly larger than the simple design assumption [55]. It was not stated how this was measured-it would be a very difficult procedure during a major flood. Elsewhere, the cited authors stated that, for the same conditions, the measured afflux was approximately six inches (0.150 m) compared with the model estimate of over one foot (0.305 m).
With the uncertainty about the stated magnitudes of afflux based on observation during the flood, it is thought that the most likely magnitude is that measured on the model-"over" 300 mm. Photographs of the model [2,5] showed that this measurement was made just upstream from the weir. In contrast to the usage here, the usual meaning given to afflux is the difference between the observed headwater and tailwater elevations at a structure such as a weir [3,15]. This will be discussed further in the context of flood levels recorded during a major flood in November/December 2021. The historical performance of both MEL weirs discussed has shown that, even during large out-of-bank floods, they cause no significant erosion at their abutments or downstream-the regions that are most at risk and require extensive protection for conventional weirs. There is evidence that they do cause some small increase in flood levels upstream during what would have been a bank-full flood, but it is difficult to assess the magnitude of this. Overall, the MEL weirs have much less environmental impact than conventional weirs.

Presentation and instrumentation
The Chinchilla MEL weir (QLD, Australia) is located on the Condamine River, on the northern part of the Murray-Darling basin (Fig. 2). The catchment area is 19,192 km 2 . The weir experienced a large flood event in Nov./Dec. 2021 and field observations took place. The maximum tailwater elevation and afflux at the peak of the 2021 flood event are reported in Table 1 and compared with several documented major floods, including several major floods in 2022 ( Table 1).
The water discharge was deduced from the measured reservoir elevation, using the predicted discharge coefficient, based upon the drawings as constructed (Fig. 5) [33]. Visual, photographic and cinematographic records were undertaken from the right bank of the spillway and from the right river bank downstream of the spillway toe. Figures 1B  and 6 present hand-held photographs taken from several locations. The observations were conducted using three dSLR Pentax™ cameras with sensor resolutions between 12 and 24 Mpx, a digital camera Sony™ RC100VA and an iPhone XI. The dSLR cameras were equipped with full-frame prime lenses producing photographs and movies with negligible barrel distortion. The dSLR camera movies were recorded in high definition (1920 × 1080 px) at 30 fps and 60 fps.
The analyses of surface features were conducted using both high-shutter speed photographs and video movies. For all the cases, the tracking and measurements of the water surface features were conducted manually to guarantee the best quality control, owing to the complexity of the prototype flow's turbulent motion, characterised by very rapid and unpredictable changes with time and space. Further, a number of movies were taken from downstream and analysed using an optical flow (OF) technique, with the camera fixed on a sturdy professional-grade tripod. The study area covered by the non-intrusive optical measurements of surface velocity was restricted to the left spillway bay, because of optical access limitations. The Optical Flow (OF) is a set of tools, detecting the flow motion between consecutive frames based upon brightness intensity gradients [9,59]. In the current study, Farnebäck's [26] OF technique was applied to calculate the surface velocity field based upon movies recorded from a camera fixed on a sturdy tripod. The Farnebäck OF technique was applied with OF parameters previously validated in laboratory for surface velocity field [6] and used for the Hinze Dam prototype spillway data sets [20]. It is acknowledged that the OF calculations were not truly 'calibrated' owing to the significant number of intrinsic difficulties with field observations during flood events, including uncontrolled optical conditions [19,20]. Since the camera field of view was perpendicular to the spillway chute, the raw data included the vertical and horizontal surface velocity components recorded at various vertical elevations. The streamwise surface velocity was deduced from the vertical surface velocity component and invert slope based upon geometric considerations, assuming a two-dimensional flow.

Prototype flow conditions
In November 2021, more than 100 mm of rainfall was recorded in the Condamine River catchment area upstream of Chinchilla, mostly from mid-November. 1 The Chinchilla weir overflowed for more than one month, with the peak headwater observed on 5 December 2021 morning. Observations of overflow discharges are presented in Fig. 7A. For seven days, the spillway discharged a flow rate greater than the bank full design discharge, with a peak flow about 1,930 m 3 /s corresponding to 2.3 times the design discharge. Note that the inspection on 15 December 2021 confirmed that the abutments were overtopped by nearly 3 m of water at the peak of the flood, without damage.
The Chinchilla weir overflow was documented on 27 November 2021 and 15 December 2021. The first data set was obtained at the start of the flood event, while the second took place after the peak of the flood (Fig. 7). The flow conditions at the time of observations are detailed in Table 2. Note that the weir overflow could not be witnessed at the peak of the flood by the authors, because the access road become submerged when the tailwater elevation exceeds about 290-291 m AHD. The afflux, difference between recorded High Water Rating Level (HWRL) and Tailwater Rating Level (TWRL), was 3.23 m at the peak of the flood event (Fig. 7B), which was somewhat consistent with some photographic observation by a local on 5 December 2021 morning. The maximum tailwater elevation and afflux during the 2021 flood event is compared with several documented major floods, including several major floods in 2022 in Table 1.

Discussion on the afflux
The differences between the magnitudes of the afflux given in Table 1 and those reported in the preceding section on historical performance are due, at least partly, to the difference between the meaning of the term implicit in the discussion of the afflux reported to have been observed during the major flood in 1973 [55] and its meaning as specified for Table 1. Clearly, the assumption of steady flow in a uniform channel, implicit in the use of the term in Turnbull and McKay [55] is not valid for a time varying flood event in a real riverthe Condamine River is very nonuniform in crosssection. In Table 1, the headwater and tailwater elevations are those observed at the respective flood gauges maintained by the Australian Bureau of Meteorology. These are some distance from the weir and part of the afflux reported in Table 1 would be due to the flood gradient at the time. Further, it is not known whether all of the flow diverted from the river upstream from the weir during overbank stages returns to the main channel upstream from the tailwater gauge. Nevertheless,  it is clear that the increase in flood levels upstream due to the weir, while small, are larger than those implicit in its design, but it is not possible to estimate the magnitude of this difference from the available information.

Basic flow observations and self-aeration
On both occasions (Table 2), the approach flow was very smooth, and this was confirmed by aerial footages on 21 December 2018 and on 3 December 2021. The reservoir flow converged smoothly towards the weir crest and the water surface was waveless and almost still. On both 27 November 2021 and 15 December 2021, the flow was critical at the spillway crest which acted as a broad-crested weir, with a ratio of crest length to upstream head over crest of 6.2 and 5.54 respectively. The smooth change on water surface elevation between the upstream reservoir and spillway crest was nicely documented with photographs (Fig. 6A). Immediately downstream of the broad-crest, the water surface was smooth and glassy (Figs. 1B, 6). On 27 November 2021, a few logs were trapped at the crest, and their wake was clearly seen at the surface of the non-aerated flow. Further downstream, the wake progressively disappeared in the self-aerated air-water flow region (Fig. 6B).
As the flow accelerated down the smooth converging chute, the glassy free-surface became rough and choppy, before becoming self-aerated, i.e. the inception region. The inception region showed a progressive transition of the surface roughness, from a smooth glassy surface to a rough, coarse-sand-paper-like appearance and further downstream to a very-rough choppy surface (Fig. 6). The self-aerated region presented a beige colour, evidence of a three-phase air-water sediment motion. Such an observation was well-documented visually on both 27 November and 15 December 2021, as well as during an earlier overflow on 8 November 1997 [53].
The average location x I of the inception region of self-aeration was recorded during both overflows, and the data are reported in Table 3 (column 5). The corresponding water depth d I was deduced for the backwater calculations (Table 3, column 6). The observations are compared to the correlations of Cain and Wood [10] for a smooth prismatic channel, assuming k s = 1 mm as used in the backwater calculations (see below). The comparative results (Table 3) showed a close agreement in terms of the water depth at inception, while the observed location data presented more scatter, possibly caused by the wide spread of the inception region, while the flow convergence was expected to induce a slower growth in bottom boundary layer compared to prismatic chutes. Physically, the converging chute geometry has two major effects on the free-surface aeration inception conditions. First, the chute convergence yields a greater water depth at a given vertical elevation for a given discharge than in absence of convergence, based upon conservation of mass since the free-stream velocity is unchanged, Second, the contraction of the chute affects the turbulent boundary layer development, as predicted by the von Karman momentum integral equation [17,42]. The former might have a more sizeable impact, leading to a further downstream selfaeration inception region location. The video movies further showed the occurrence of some form of two-dimensional roll waves in the self-aerated flow region. The phenomenon was clearly seen on 27 November 2021, and less visible on 15 December 2021. Although the same camera system (camera body and lens) was used on both occasions, the difference might be linked to differences in light conditions as well as slightly different water discharges. For completeness, the quasi-two-dimensional roll wave patterns was not visually seen in-situ, but it was clearly evidenced during the movie replay, as well as in the OF data. For both flow conditions, the instability criteria of Keulegan and Patterson [35] was fulfilled [25]. The Vedernikov number [24,47]  At the toe of the spillway chute, the high-velocity self-aerated flow impinged into the tailwater and a hydraulic jump formed. The jump was stationary and very stable, located on the underwater sloping chute on both days. The position of the jump was controlled by the tailwater conditions and the very-flat Condamine River channel bed slope downstream of Chinchilla. The jump roller was highly turbulent and fluctuating about its mean position. Both visual and video movie observations indicated strong three-dimensional large-scale vortices in the roller. The interactions of these large turbulent structures with the roller free-surface created large scars and surface features, with length-scales comparable to the roller height.

Presentation
Based upon the video movies of the left chute centreline, taken from downstream, the surface velocity field was measured in the plane parallel to the invert using the optical flow (OF) technique. Both the streamwise velocity V s and transverse velocity V t components were extracted. First, the present data provided physically meaningless results in the upstream non-aerated flow region. This was most likely caused by the shiny surface glare, seen in Figs. 1B and 6, upstream of the inception region. Second, physically meaningful surface velocity results were obtained in the self-aerated flow region and are presented below. This finding differed from the OF observations at the Hinze Dam spillway [19,20] and are discussed later.
For the left spillway section, the time-averaged centreline OF surface velocities were compared to the ideal fluid flow theory and to the backwater equation. On an uncontrolled chute spillway, the flow is accelerated by the gravity force component in the flow direction. The ideal velocity V max at a sampling location on the chute is derived from the Bernoulli principle: where H 1 is the upstream total head above the sampling point, θ is the channel slope and d is the water depth [15], pp. 407 & 486). For a steady non-uniform real fluid flow, the differential form of the energy equation is called the backwater calculations developed assuming that the flow is gradually varied. For a rectangular channel of non-constant slope and width, the backwater equation yields [15], pp. 377-381): with S f the friction slope, and B the free-surface width. At Chinchilla MEL weir, the channel cross-section is rectangular and the chute has a constant slope (θ = 11.31°) with a reducing width. Then, the backwater equation becomes: with cosθ = 0.9806, sinθ = 0.196, and B/ x = − 2.1805 at the Chinchilla MEL weir.
The surface velocity data showed a relatively close agreement between the timeaveraged streamwise OF surface velocity and the backwater-calculated velocity downstream of the onset region, both qualitatively and quantitatively (Fig. 8). Figure 8 presents some typical comparisons in terms of streamwise surface velocity component V s , with each graph corresponding to one movie. In Fig. 8, the ideal fluid flow estimates are shown with a solid green curve and the backwater calculations are the dashed red curve. For completeness, the streamwise surface velocity standard deviation v s ' and transverse surface velocity standard deviation v t ' are added in Fig. 8. In Fig. 8, the vertical axis is the vertical elevation in m AHD, where AHD stands for Australian Height Datum, and the horizontal axis is the surface velocity data in a dimensional form.
The time-averaged OF surface velocity data presented a number of key features. Firstly, the time-averaged longitudinal velocity component increased with decreasing vertical elevation, as it would be physically expected (Fig. 8). Secondly, the longitudinal velocity data showed some poor agreement with ideal fluid velocity estimates in the developing non-aerated flow region. For z o > 294 m AHD, the poor results were caused by the water surface reflections, seen in Figs. 1A and 6, preventing the OF technique to track the flow motion based upon the brightness intensity gradient. Thirdly, the time-averaged OF velocity data agreed relatively well with the backwater calculations developed for the converging chute (see above) and valid in the fully-developed self-aerated flow region. The comparison reflects well the adequacy of the backwater calculations assuming a gradually-varied chute flow. While the method was herein applied to the self-aerated flow region, drag reduction caused by self-aeration was not expected to occur for this flat slope (1 V:5H), because of the small values of mean void fraction down the chute slope, with nil to negligible effect on the friction factor, even if quasi uniform air-water flow would be attained (e.g. as per [12]. At the downstream end of the chute, i.e. for z o < 288 m AHD, a lack of match with backwater estimate might be caused by optical artifacts caused by spray and splashing above the hydraulic jump roller at the chute toe. Indeed, visual observations by the first author showed drops and splash reaching heights of more than 0.5-1 m above the hydraulic jump roller surface. Fourthly, the OF data output quality was closely linked to the quality of the original movies. Practical considerations encompassed the camera position, ideally placed perpendicularly to the spillway chute, a sturdy camera tripod, and the camera body and lens equipment, with all the present data set acquired with professional-grade full-frame prime lenses, following earlier experiences [6,19,21]. More, a higher movie definition delivered improved outputs with greater spatial resolutions, i.e. less millimetres per pixel. It is worth to mention that both OF field data sets were filmed under good weather conditions, hence decent lighting conditions.

Surface velocity field in the self-aerated flow region
The In the contour maps, the tailwater level is shown as well. As seen in Figs. 8, 9 and 10, the results were unreliable in the upstream non-aerated flow. Further, the surface velocity data were not reliable below elevation 288 m AHD, because of the splashing above the hydraulic jump roller. In the self-aerated flow region, the surface velocity data implied some regions of high-velocity and others of low-velocity. That is, longitudinal "canyons" of faster flowing water, with streaks of lower surface velocities in between, in the self-aerated flow region downstream of the inception region. The same data sets highlighted regions of high-and low-streamwise velocity fluctuations.
The observations were similar to recent field observations at the Hinze Dam Stage 3 [19,21] and in a relatively large-size spillway model [6]. The findings are important in terms of the design and operation of energy dissipators, e.g. hydraulic jump stilling basins and ski jumps. Indeed, these regions of high velocity are associated with concentrations in kinetic energy that must be dissipated safely in the stilling structure.
Further, the weighted averages of the velocity standard deviations are presented in Table 4, together with the number of analysed frames, the frame rate and the spatial  resolutions. Noteworthy, the data are based upon a large number of analysed frames in the current study (Table 4, 4th column). The streamwise velocity data presented large surface turbulent intensity Tu s = v s '/V s in the range 1.7 to 2 ( Table 4, 6th column). The large values were likely the resulting combination of (a) relatively large velocity fluctuations initiated by surface wave motion and of (b) true turbulence. The data were quantitatively of the same magnitude as a series of field measurements at the Hinze Dam spillway [20] and as self-aerated flow data in large-size laboratory facilities [6,23]. In contrast, the transverse velocity data exhibited much lower transverse surface turbulent intensity Tu t = v t '/V s within 0.11-0.17, likely induced by the smooth convergence of the flow. The ratio Tu t /Tu s was between 0.06 and 0.1 (Table 4, 8th column). Such a result implied some very strong anisotropy of the surface turbulence. Noteworthy, the present prototype observations were close to reported values of the ratio v t '/v s ' in rough turbulent boundary layers [51] and on the Hinze Dam prototype spillway [21].

Discussion
The present study results may be compared to earlier results in free-surface investigations in prototype and laboratory spillways, although most were conducted with steep stepped spillways ( Table 5). All the observations showed that the onset of self-aeration extended over a surface plane with some transient pattern, constantly fluctuating about a mean location. On the smooth relatively flat spillway of Chinchilla weir, the longitudinal transition from non-aerated to air-water flow region was somewhat gradual, while it was sudden and violent on steep stepped chutes. On steep stepped spillways, the onset was driven by the interactions of the developing boundary layer turbulence with the free-surface [10,21,45,58]. At Chinchilla MEL weir, it is believed that the inception of self-aeration resulted from the combined interactions between surface instabilities and developing boundary layer. The surface velocity data showed marked differences in the applicability of the OF calculation. Physically meaningful results were achieved in terms of the self-aerated surface velocity field in the large-size stepped spillway model and Chinchilla MEL weir spillway (Table 5). In both cases, the motion of the "grainy" air-water surface roughness was successfully captured by the OF method. Several studies discuss the influence of surface roughness on the inception of self-aeration (e.g. [7,8,56], while the interactions between self-aeration and surface waves were shown to have relatively little influence on the void fraction and air entrainment rate [53]. On prototype stepped spillways (Table 5), the air-water region exhibited an intensely bright white colour which prevented any detection of structure advection, with the camera equipment used in these studies, although it is conceivable that IR camera systems might be more successful.
Finally, it is acknowledged that the OF data were obtained with relatively low frame rates, compared to recent laboratory studies [6,9,59]. Indeed, the application of OF technique differs substantially between laboratory and field [19,20], and its implementation for prototype flow tends to be closer to large-sale particle image and infrared quantitative image velocimetry applications (e.g. [38,39,52]). There are however some key differences between laboratory and field applications of OF, linked to the differences in light intensity often limited indoors, surface feature sizes typically smaller in laboratory, and location of camera often restricted in the field. While ideally a faster frame rate would be desirable,  [6] other practical considerations tend to be more important and relevant to ensure successful surface velocity data in prototype spillway [20].

Conclusion
The Minimum Energy Loss (MEL) weir design was developed in Queensland, Australia specifically for the river catchments affected by heavy tropical and sub-tropical rainfalls where streams have very flat gradients and erodible banks. A number of MEL weirs have successfully operated for many decades and the operations of two of them were discussed. Inspections during and after major flood events indicated a reliable operation associated with minimum maintenance. The historical performances of the MEL weirs have confirmed the design capability to pass large floods with small afflux and very small energy loss and that they cause no significant erosion at their abutments or downstream-the regions that are most at risk and require extensive protection for conventional weirs.
During the November-December 2021 flood, visual and quantitative observations were undertaken at the Chinchilla MEL weir. The spillway chute was observed and documented during two days, at the start and towards the end of a major flood event peaking on 5 December 2021. The approach flow was very smooth and streamlined before it reached the broad crest, and then continued down a converging smooth chute. Downstream of the crest, the chute flow presented a very smooth and glassy free-surface. On 17 November 2021, the presence of large logs, blocked on the crest, had no major impact on the chute flow. Further downstream on the smooth converging chute, the observations showed that the inception of self-aeration was a three-dimensional process with a 'progressive' change in free-surface roughness, from a smooth glassy free-surface, to a coarse-sand-paper-like appearance and then a very-rough choppy surface. Both the upstream non-aerated and selfaerated flow regions exhibited a brown colour indicating the large amount of suspended sediment materials.
An optical technique was implemented to derive the contour maps of surface velocities based upon video movies taken from a sturdy tripod. The present data sets showed results close to the backwater equation in the self-aerated region, while highlighting large transverse difference in longitudinal surface velocities across the chute. The OF data presented large streamwise surface velocity fluctuations in the aerated flow region, with Tu s ~ 150-200%, consistent with other field data sets and some literature on self-aerated flow measurements using dual-tip phase detection probe in large-size laboratory models. The ratio of transverse to streamwise surface turbulence intensity Tu t / Tu s was between 0.06 and 0.1, indicating a strong anisotropy of the free-surface turbulence.