A verification of an analytical formula for estimating the height of the stable internal boundary layer

A stable thermal internal boundary layer (IBL) develops when warm air is advected from warmer land upstream to a cooler sea downstream. It is shown that the analytical model for estimating the height (h) of this stable IBL as formulated by Garratt (1987) is verified. It is also demonstrated that a simpler equation, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGObGaeyisIS% RaaGymaiaaiAdacaWGybWaaWbaaSqabeaadaWcgaqaaiaaigdaaeaa% caaIYaaaaaaaaaa!390B!\[h \approx 16X^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} \] (where h is in meters and X, the fetch downwind, is in kilometers), is useful operationally as a first approximation.     

大气垂直运动对雷达估测降水的影响

运用雷达反射率因子（Z）与降水强度（R）之间的关系定量估测降水,降水云体中的大气垂直运动（w_a）不可忽视。PARSIVEL激光雨滴谱仪（简称PARSIVEL）在获取雨滴粒径分布的同时可以从测量的雨滴下降速度分布中提取w_a,用于分析PARSIVEL高度上的大气垂直运动对雷达Z估测降水强度影响。使用2014年5-6月华南季风降水观测试验期间广东阳江5次层状云、6次对流性降水过程中PARSIVEL数据分析大气垂直运动对定量估测降水影响,w_a对层状云和对流云降水强度影响范围分别为-0.18~1.05 mm·h-1和-5.44~24.81 mm·h-1,相对影响值分别为-13.61%~13.99%,-38.59%~25.92%。静止大气条件下,雷达估算降水Z-R关系式中系数A,b引起的层状云和对流云降水估测偏差平均分别为10.9%和25.5%。真实大气中雷达估测降水的偏差平均情况是层状云降水由于w_a的对消作用降低为9.2%,对流云降水则增加到51.2%。对流性降水中大气垂直运动对雷达估测降水的影响较大。     

The 2001 Mesoscale Convective Systems over Iberia and the Balearic Islands

Summary This paper characterizes Mesoscale Convective Systems (MCSs) during 2001 over Iberia and the Balearic Islands and their meteorological settings. Enhanced infrared Meteosat imagery has been used to detect their occurrence over the Western Mediterranean region between June and December 2001 according to satellite-defined criteria based on the MCS physical characteristics. Twelve MCSs have been identified. The results show that the occurrence of 2001 MCSs is limited to the August–October period, with September being the most active period. They tend to develop during the late afternoon or early night, with preferred eastern Iberian coast locations and eastward migrations. A cloud shield area of 50.000 km² is rarely exceeded. When our results are compared with previous studies, it is possible to assert that though 2001 MCS activity was moderate, the convective season was substantially less prolonged than usual, with shorter MCS life cycles and higher average speeds. The average MCS precipitation rate was 3.3 mm·h⁻¹ but a wide range of values varying from scarce precipitation to intense events of 130 mm·24 h⁻¹ (6 September) were collected. The results suggest that, during 2001, MCS rainfall was the principal source of precipitation in the Mediterranean region during the convective season, but its impact varied according to the location. Synoptic analysis based on NCEP/NCAR reanalysis show that several common precursors could be identified over the Western Mediterranean Sea when the 2001 MCSs occurred: a low-level tongue of moist air and precipitable water (PW) exceeding 25 mm through the southern portion of the Western Mediterranean area, low-level zonal warm advection over 2 °C·24 h⁻¹ towards eastern Iberia, a modest 1000–850 hPa equivalent potential temperature (θ_e) difference over 20 °C located close to the eastern Iberian coast, a mid level trough (sometimes a cut-off low) over Northern Africa or Southern Spain and high levels geostrophic vorticity advection exceeding 12·10⁻¹⁰ s⁻² over eastern Iberia and Northern Africa. Finally, the results suggest that synoptic, orographic and a warm-air advection were the most relevant forcing mechanisms during 2001.     

Structure of the marine atmospheric boundary layer during two cold air outbreaks of varying intensities: GALE 86

Aircraft, surface, upper air and satellite measurements have been used to observe the evolution and growth of the convective Marine Atmospheric Boundary Layer (MABL) offshore of North Carolina in close proximity to the Gulf Stream, during the intense cold air outbreak of 28 January 1986 and the moderate event of 12 February 1986, as part of the Genesis of Atlantic Lows Experiment (GALE). Air mass modification processes, driven primarily by the ocean-atmosphere exchanges of surface turbulent sensible and latent heat fluxes, caused the overlying air mass to warm and moisten as it advected over the warmer waters of the eastern United States continental shelf. Maximum observed near-surface total heat fluxes were 1045 and 811 W·m^–2 over the core of the Gulf Stream, for 28 January and 12 February 1986, respectively. The observed changes in the overlying air mass occurred almost instantaneously as the ambient flow traversed different underlying SST conditions.The turbulent structure showed a buoyancy-dominated MABL below approximately 0.8z/h. However, shear was also observed to be an important production term above 0.8z/h and below 0.1z/h for the 28 January 1986 event. Dissipation of turbulent kinetic energy was the dominant destruction term in the budgets, but vertical transport of energy was a strong contributor below 0.5z/h, above which this term became a source of turbulent energy. Additionally, the normalized standard deviations of the horizontal velocity components showed a near-equal contribution to the turbulence, while the vertical velocity components displayed the characteristic mid-layer maximum profile observed for a convective, well-mixed boundary layer.     

Comparisons of xylem sap flow and water vapour flux at the stand level and derivation of canopy conductance for Scots pine

Summary Simultaneous measurements of xylem sap flow and water vapour flux over a Scots pine (Pinus sylvestris) forest (Hartheim, Germany), were carried out during the Hartheim Experiment (HartX), an intensive observation campaign of the international programme REKLIP. Sap flow was measured every 30 min using both radial constant heating (Granier, 1985) and two types of Cermak sap flowmeters installed on 24 trees selected to cover a wide range of the diameter classes of the stand (min 8 cm; max 17.5 cm). Available energy was high during the observation period (5.5 to 6.9 mm.day^–1), and daily cumulated sap flow on a ground area basis varied between 2.0 and 2.7 mm day^–1 depending on climate conditions. Maximum hourly values of sap flow reached 0.33 mm h^–1, i.e., 230 W m^–2.Comparisons of sap flow with water vapour flux as measured with two OPEC (One Propeller Eddy Correlation, University of Arizona) systems showed a time lag between the two methods, sap flow lagging about 90 min behind vapour flux. After taking into account this time lag in the sap flow data set, a good agreement was found between both methods: sap flow = 0.745^* vapour flux,r ² = 0.86. The difference between the two estimates was due to understory transpiration.Canopy conductance (g _c ) was calculated from sap flow measurements using the reverse form of Penman-Monteith equation and climatic data measured 4 m above the canopy. Variations ofg _c were well correlated (r ² = 0.85) with global radiation (R) and vapour pressure deficit (vpd). The quantitative expression forg _c =f (R, vpd) was very similar to that previously found with maritime pine (Pinus pinaster) in the forest of Les Landes, South Western France.With 6 Figures     

在不同粗糙表面上龙卷涡旋特性的实验室模拟实验研究

本工作是在日本京都大学,防灾研究所的龙卷模拟实验室,用流体动力学模似实验方法研究三种不同粗糙表面对龙卷涡旋的影响。实验中,用热线风速表及微型风向标,测定了五个不同高度上的水平风速和风向分布,揭示出在不同粗糙表面上,龙卷涡旋的风速和风向分布特征。通过不断改变模拟装置的涡旋比S,发现当S值超过某一数值后(S>1.5),在不同粗糙表面上,龙卷涡旋的涡核尺寸将趋于一致。实验结果可进一步解释自然界中的观测事实,并对某些理论研究提供实验依据。     

The stably stratified internal boundary layer for steady and diurnally varying offshore flow

A two-dimensional numerical mesoscale model is used to investigate the internal structure and growth of the stably stratified internal boundary layer (IBL) beneath warm, continental air flowing over a cooler sea. Two situations are studied — steady-state and diurnally varying offshore flow. In the steady-state case, vertical profiles of mean quantities and eddy diffusion coefficients (K) within the IBL show small, but significant, changes with increasing distance from the coast. The top of the IBL is well defined, with large vertical gradients within the layer and a maximum in the coast-normal wind component near the top. Well away from the coast, turbulence, identified by non-zero K, decreases to insignificant levels near the top of the IBL; the IBL itself is characterised by a critical value of the layer-flux Richardson number equal to 0.18. The overall behaviour of the mean profiles is similar to that found in the horizontally homogeneous stable boundary layer over land.A simple physical model is used to relate the depth of the layer h to several relevant physical parameters viz., x, the distance from the coast and U, the large-scale wind (both normal to the coastline) and g/, being the temperature difference between continental mixed-layer air and sea surface, is the mean potential temperature and g is the acceleration due to gravity. Excellent agreement with the numerical results is found, with h = 0.014x ^1/2 U (g/)^–1/2.In the diurnally varying case, the mean profiles within the IBL show only small differences from the steady-state case, although diurnal variations, particularly in the wind maximum, are evident within a few hundred kilometres of the coast. A mesoscale circulation normal to the coast, and superimposed upon the mean offshore flow, develops seawards of the coastline with maximum vertical velocities about sunset, of depth about 2 km and horizontal scale 500 km. The circulation is related to the advection, and subsequent decay, of daytime convective turbulence over the sea.     

Sensitivity analysis of the equations for a convective mixed layer

Jump or slab models are frequently used to calculate the depth of the convectively mixed layer and its potential temperature during the course of a clear day. Much attention has been paid theoretically to the parameterization of the budget for turbulent kinetic energy that is required in these models. However, for practical applications the sensitivity of the solutions of the model equations to variations in the entrainment formulation and in the initial and boundary conditions is also very important. We analyzed this sensitivity on the basis of an analytical solution for the model which uses the well-known constant heat flux ratio. The initial conditions for the mixed-layer height (h) and potential temperature ( _m ) quickly lose their influence. Only the initial temperature deficit is important. The mixed-layer temperature at noon on convective days is insensitive to the entrainment coefficient c. It is governed by the integral of the heat input and by the stable lapse rate. A change in c from 0.2 to 0.5 leads to a variation of 20% in h. This is not very much considering the accuracy in the determination of h from actual observations.     

Mean circulations of boundary-layer rolls in lake-effect snow storms

An observational study of the roll-average flow fields of boundary-layer rolls is presented. Data used for this purpose were collected by dual-Doppler radar and aircraft measurements taken over southern Lake Michigan during the 1983/1984 field operations of Project Lake Snow. The roll circulations agreed well with findings of past observational, numerical and theoretical studies, with cross-roll components roughly 10% of the convective internal boundary layer (CIBL)-mean wind speeds and weaker vertical components. Along-roll winds were systematically stronger in the rollupdraft regions than in the roll-downdraft regions, probably due to distortion of the along-roll wind profile by the rolls. Comparison of observed wind profiles to those required by roll formation mechanisms found by past numerical and theoretical studies suggested that the observed rolls were formed by the along-roll wind shear (Asai, 1970) or wind shear curvature (Kuettner, 1971) in the lowest 0.2Z _i, whereZ _i is the height of the top of the CIBL.     

Thermal effect of the sea breeze on the structure of the boundary layer and the heat budget over land

The daytime boundary-layer heating process and the air-land heat budget were investigated over the coastal sea-breeze region by means of observations over the Sendai plain in Japan during the summer. In this area, the onset of the sea breeze begins at the coast around 0900 LST, intruding about 35 km inland by late afternoon. The cold sea breeze creates a temperature difference of over 10°C between the coastal and inland areas in the afternoon. On the other hand, warm air advection due to the combination of the counter-sea breeze and land-to-sea synoptic wind occurs in the layer above the cold sea breeze in the coastal region. Owing to this local warm air advection, there is no significant difference in the daytime heating rate over the entire atmospheric boundary layer between the coastal and inland areas. The sensible heat flux from the land surface gradually decreases as distance from the coastline increases, being mainly attributed to the cold sea breeze. The daytime mean cold air advection due to the sea breeze is estimated asQ _adv ^local =–29 W m^–2 averaged over the sea breeze region (035 km from the coastline). This value is 17% of the surface sensible heat fluxH over the same region. The results of a two-dimensional numerical model show that the value ofQ _adv ^local /H is strongly affected by the upper-level synoptic wind direction. The absolute value ofQ _adv ^local /H becomes smaller when the synoptic wind has the opposite direction of the sea breeze. This condition occurred during the observations used in the present study.     

A Parameterization of Bowen Ratio with Respect to Soil Moisture Availability

The Bowen ratio (B) is impacted by 5 environmental elements: soil moisture availability, m, the ratio of resist-ances between atmosphere and soil pores, ra/rd, atmospheric relative humidity, h, atmospheric stability, ΔT, and environment temperature. These impacts have been investigated over diverse surfaces, including bare soil, free water surface, and vegetation covered land, using an analytical approach. It was concluded that: (a) B is not a continuous function. The singularity exists at the condition αhcb=h, occurring preferably in the following conditions: weak turbulence, stable stratified stability, dry soil, and humid air, where hcb, defined by Eq.(11) is a critical variable. The existence of a singularity makes the dependence of B on the five variables very complicated. The value of B approaches being inversely proportional to m under the conditions m≥mfc (the soil capacity) and / or ra/rd→0. The proportional coefficient changes with season and latitude with relatively high values in winter and over the poles; (b) B is nearly independent of ra/rd during the day. The impact of m on B is much larger as compared to that of ra/rd on B, (c) when h increases, the absolute value of B also increases; (d) over bare soil, when the absolute surface net radiation increases, the absolute value of B will increase. The impact of RN on B is larger at night than during the day, and (e) over plant canopy, the singularity and the dependcies of B on m, ra , and h are modified as compared to that over bare soil. Also (i) during the daytime unstable condition, m exerts an even stronger impact on B, at night, however, B changes are weak in response to the change in m; (ii) the value of B is much more sensitive in response to the changes of turbulent intensity; (iii) the B response to the variation of h over a vegetation covered area is weaker; and (iv) the singularity exists at the condition hcp=h instead of αhcb=h as over bare soil, where hcp is defined by Eq.(49). The formulas derived over bare soil also hold the same when applied to free water bodies as long as they are visualized as a special soil in which the volumetric fraction of soil pore is equal to one and are fully filled with water. Finally, the above discussions, are used to briefly study the impact on the thermally induced mesoscale circulations.     

Measurements of the humidity structure function parameters,C q 2and C Tq,over the ocean

Atmospheric turbulence measurements, including temperature and humidity fluctuations, were made from the R/V Acania off the coast of California in June, 1979. The purpose of the experiment was to investigate the scaling properties of the humidity structure function parameter (C _q ²) and temperaturehumidity cospectrum structure parameter (C _Tq) in the marine surface layer. The bulk parameterization method was used to obtain Monin-Obukhov Similarity (MOS) scaling parameters u _*, T _*, q _*and L. Assuming a neutral stability humidity drag coefficient c _qn = 1.3 × 10^-3the dimensionless humidity structure function parameter C _q ²Z^2/3/q_* ²was found to be 18% lower than the corresponding temperature function obtained by Wyngaard et al. (1971). Furthermore, the measurements indicate that the temperature-humidity fluctuations are highly coherent well into the inertial subrange. The results have direct application to turbulent scattering of waves propagating in the atmosphere (particularly microwaves) and methods of estimating air-sea surface fluxes.     

Fluctuations of the horizontal wind under unstable conditions

The similarity relations for _u/u_* proposed by Panofskyet al. (1977) and Højstrup (1982) have been verified using eddy-correlation data collected during the EFEDA-experiment, conducted over the extensive plain of La Mancha (Spain), where vine plants form a primary crop. Also, the standard Monin-Obukhov relation is considered. It is found that the expressions by Panofskyet al. and Højstrup both yield almost identical results, and are better than the Monin-Obukhov expression. Also, _u measured with a cup anemometer obeys the similarity expressions well, provided that the effect of the variation of wind direction on propellor wind speed is accounted for. The relationship of Panofskyet al. works rather well even when the boundary-layer height scaleh is replaced by a fixed height,h _c. Best results were obtained forh _c=1800 m. This height scale is possibly associated with the horizontal variability in the surface sensible heat flux pattern.     

A simple model of drainage flow on a slope

The concept of a cold air ‘Parcel’ is introduced for describing the bulk properties of drainage flow. By means of a model based on the momentum and sensible heat transports under calm conditions, the thickness h and velocity u of the Parcel are derived in simple forms. It is shown that h and u correspond to the inversion height and maximum velocity of actual drainage flow. The governing parameters for h and u are the length and vertical drop of the slope, potential temperature difference between the ambient atmosphere and the Parcel, aerodynamic condition of the slope surface expressed by the mean bulk coefficients, and ambient stability. The mean bulk coefficients depend on the roughness lengths for the velocity and potential temperature profiles and are decreasing functions of the slope length. The Parcel Model agrees qualitatively with Manins and Sawford's (1979) model under neutral ambient stratification. But agreement is not so good under stable conditions. The thickness and velocity of drainage flow predicted by the Parcel Model agree with observations on slopes several tens of meters to several hundred kilometers long.     

Transfer coefficients of sensible heat on a snowmelt surface

Summary Field observations were carried out in order to determine the transfer coefficeient of sensible heat flux above a melting snow surface at the Moshiri experimental site. The coefficient is calculated as the ratio of the sensible heat flux determined by the eddy correlation method using a sonic anemometer to the product of the wind speed and the temperature difference between the air and snow surface. The sensible heat fluxes are also compared with the result of the precise heat balance observations. The nondimensional transfer coefficienth shows a good correlation with the atmospheric stability,Ri. The value ofh=2.3×10^–3 is obtained in the range of 0<Ri<0.1, however, it is smaller and scattered under stronger stable condition (Ri>0.1). The dimensional transfer coefficent for sensible heat flux is calculated, and a linear relationship is obtained as a function of the logarithm of atmospheric stability.With 11 Figures     

Bulk transfer coefficient over a snow surface

The drag coefficient C _D and the bulk transfer coefficient for sensible heat C _H over a flat snow surface were determined experimentally. Theoretical considerations reveal that C _D depends on the friction velocity u _* as well as on the geometrical roughness h of the snow surface. It is found that C _D increases with increasing u _* and/or h. The dependency of C _H on u _* and h is so small that it is possible to consider C _H as a constant for practical purposes: C _{H, 1} = 2.0 × 10^–3 for a reference height of 1 m. The bulk transfer coefficient for water vapor is estimated at C _{E, 1} = 2.1 × 10^–3 for a reference height of 1 m.     

A variable K planetary boundary-layer model

The steady-state, homogeneous and barotropic equations of motion within the planetary boundary layer are solved with the assumption that the coefficient of eddy viscosity varies as K(Z) = K _O(1–Z/h) ^p , where h is the height of the bounday layer and p is a parameter which depends on atmospheric stability. The solutions compare favourably with observed velocity profiles based on the Wangara data.     

Mean and turbulent structure of a baroclinic marine boundary layer during the 28 January 1986 cold-air outbreak (GALE 86)

Aircraft (NCAR Electra), ship (R/V Cape Hatteras), buoy (NCSU Buoy 2) and satellite (NOAA-7 and 9) measurements have been used to observe the structure of the Marine Boundary Layer (MBL) offshore of Wilmington, North Carolina, during the intense cold-air outbreak of 28 January, 1986, as part of the Genesis of Atlantic Lows Experiment (GALE). Air mass modification processes, driven primarily by the surface turbulent latent and sensible heat fluxes, caused the overlying air mass to warm and moisten as it advected over the warmer waters of the eastern United States continental shelf. Maximum observed total (latent + sensible) heat flux was 1045W/m² (at a height of 49 m) over the core of the Gulf Stream. Heat flux values decreased both east and west of this region, primarily in response to changes in the air-sea temperature difference.MBL height increased steadily in the offshore direction in response to increasing convection. The turbulent structure showed a buoyancy-dominated MBL between 0.1 z/h and 0.8 z/h; whereas shear was important above and below this level, vertical transport of kinetic energy (KE) was dominant as a source term only above 0.8 z/h. The normalized turbulent kinetic energy (TKE) budgets observed at different offshore locations showed general agreement at different flight levels. Thus the findings support the validity of the similarity relations under intense convective conditions.     

Simulating multi-decadal variability of Caspian Sea level changes using regional climate model outputs

The Caspian Sea is the largest enclosed body of water on earth, covering approximately 4×10⁵ km² and sharing its coast with five countries (Iran, Azerbaijan, Kazakhstan, Russia and Turkmenistan). Because it has no outlet to the ocean the Caspian Sea level (CSL) has undergone rapid shifts in response to climatic forcings, and these have been devastating for the surrounding countries. In this paper we present the initial results of a modeling effort aimed at building a regional climate model for the Caspian Sea basin suitable to study the response of the CSL to interdecadal climate variability and anthropogenic climate change. Simulations are performed using the International Centre for Theoretical Physics (ICTP) regional climate model RegCM at a 50 km grid spacing for the period 1948–1990. During this period an abrupt shift occurred in the sea level after 1977, when the CSL rose about two meters until the early 1990s. Using a simple equation of hydrologic balance for the Caspian Sea basin to predict the CSL, we show that the model is able to reproduce the observed CSL changes at interannual to multidecadal scales. The correlation coefficient between the simulated and observed annual CSL changes is 0.91 and the model is able to reproduce the abrupt shift in CSL which occurred after 1977. Analysis of the climatologies before and after 1977 indicate that the CSL rise was mostly due to an increase in precipitation over the northern basin and a decrease in evaporation over the sea, primarily during the warm season. We plan to apply our model to the investigation of the response of the CSL to anthropogenic climate forcings.     

Dissipative effects on inertial flows over a sill

A systematic investigation of the effects of various parametrizations of dissipation, e.g. quadratic and linear frictional drag, harmonic lateral viscosity, and harmonic lateral diffusion on inertial flow over a sill and possible hydraulic control is presented. Rotation effects are ignored and the geometry is assumed to vary only slowly with downstream distance so that the flow may be considered one-dimensional. Results are given both for a single-active layer and for two-active layers with a rigid lid.If the parametrization is only a function of the dependent variables and not of their spatial derivatives, then it may be possible to hydraulically control the flow. A general expression is derived for the possible control point and the two gradients there, which are functions of the slope and possibly of flow rate. Specific energy is irreversibly removed from the flow and non-controlled as well as controlled flows can exhibit significant asymmetry in fluid depth over a sill. The upstream specific energy, and hence depth of the lower layer, of the controlled flow is greater than for an ideal fluid. Frictional effects modify the behaviour of long gravity waves, such that they are dispersive and damped with time. The system will only exhibit hydraulic control if these effects are small.For a viscous single layer of fluid, the gradient in surface elevation is always uniquely defined, so classically defined hydraulic control, as such, cannot exist. However, for values of non-dimensional lateral eddy viscosity coefficient, , where q is the flow rate, there is a narrow band of specific energies centred around that for the control solution in an ideal fluid, E_crit, for which the surface elevation, h is very asymmetric over the sill; the solutions resemble the inviscid, hydraulically controlled solutions. Outside this range, either the fluid depth tends to zero, or the surface elevation is almost uniform over the sill. A ‘control’-type solution exists which has the conjugate values of the inviscid equation up- and downstream of the sill, where the gradient in fluid depth, and hence the viscous term, is zero. For larger values of A_M, the band of specific energies is much wider, and the upstream specific energy of the ‘control’-type solution is much lower than that for an inviscid fluid. Long gravity waves are dispersive and damped with time. There is a short-wave cut-off, k² > h/(4A_M²), above which waves are stationary in the flow. Longer waves, k² h/(4A_M²), are critical if , as for an ideal fluid. If these waves can propagate significant distances, then any observed asymmetry in h will be due to inertial and not to viscous effects. The behaviour of unidirectional, two-layer flow is similar. The governing equation for viscous, two-layer exchange flow is singular, and typically excludes the ‘control’-type solutions found for unidirectional flows.Establishing the existence and behaviour of steady inertial flows in the presence of lateral diffusion between layers is more difficult. It significantly alters the single-layer solutions once the non-dimensional coefficient A_H is large, i.e. . The flow rate may become zero on the downslope as all the fluid diffuses into the inert, infinitely deep, overlaying layer. The fluid depth is maintained by reverse flow from downstream. In this case, the depth of the active layer tends to zero downstream for all values of specific energy. For two-layer flow, both unidirectional and exchange, the governing equation is such that the lower-layer flow rate and interfacial height return to their upstream values.Motivation for the study is provided by the increasingly fine spatial resolution achievable in large-scale numerical models of the ocean general circulation, and the question of whether they are capable of simulating some form of hydraulic control. Application to modelling oceanic flows over a sill is discussed.