One-dimensional Mixed Layer Model Research Articles

Sensitivity of the Short-Range Precipitation Forecast in South China Region to Sea Surface Temperature Variations

Ocean variability plays an essential role in the climate system at different time scales through air–sea interactions. Recent studies have addressed the importance of the ocean mixed layer in cooling feedback to tropical cyclones (TCs). However, using constant sea surface temperature (SST) in short-range weather forecasts remains common, especially in high-resolution regional models. This study investigates the role of subsurface ocean mixing in the short-range forecast of non-TC extreme rainfall with the Weather Research and Forecast (WRF) model. In the simulations of 26 heavy rainfall cases, we found that using a one-dimensional mixed layer model leads to a 15% enhancement (reduction) of rainfall maximum in six (two) cases compared to using constant SST. When the initial depth of the mixed layer model is perturbed by the amount of daily variability, 13 cases exhibit larger than 15% increases or decreases. A detailed analysis of one case suggests that the radiative process dominates the overall response of SST. The warming and moistening of boundary layer air cause significant strengthening of updrafts in convection. Although the SST change in most cases due to varying mixed layer model setups is less than 0.5 K, convective motions in some cases are surprisingly sensitive to small changes.

Salinity response to atmospheric forcing of the Terra Nova Bay polynya, Antarctica

AbstractThe density and salinity of High Salinity Shelf Water, a key component of Antarctic Bottom Water emanating from the Ross Sea, are intensified by brine rejection induced by ice formation within the Terra Nova Bay (TNB) polynya. Ocean mooring data from 2007, meteorological observations from automatic weather stations and a satellite-derived history of the opening of TNB polynya delineate variability in water column salinity linked to atmospheric forcing, with a period on the order of 10 days. Lagged correlation analysis indicates that on average salinity response lags the polynya opening by 2 days and the wind forcing by 5 days. We find stronger correlations of salinity with the wind during March through May and with the polynya open-water fraction during June through October, with decreasing lags in the latter period. A one-dimensional mixed-layer model incorporating thermodynamic ice formation captures the oscillations in salinity. A process study shows that the variability in the polynya open-water fraction governs the final salinity attained by the model as well as the salinity cycling. Variability in surface heat fluxes modulates that effect. Our work suggests that there is a more complex relationship between salinity, the polynya open-water fraction, and atmospheric forcing than previously suggested.

Taiwan Earth System Model Version 1: description and evaluation of mean state

Abstract. The Taiwan Earth System Model (TaiESM) version 1 is developed based on Community Earth System Model version 1.2.2 of National Center for Atmospheric Research. Several innovative physical and chemical parameterizations, including trigger functions for deep convection, cloud macrophysics, aerosol, and three-dimensional radiation–topography interaction, as well as a one-dimensional mixed-layer model optional for the atmosphere component, are incorporated. The precipitation variability, such as diurnal cycle and propagation of convection systems, is improved in TaiESM. TaiESM demonstrates good model stability in the 500-year preindustrial simulation in terms of the net flux at the top of the model, surface temperatures, and sea ice concentration. In the historical simulation, although the warming before 1935 is weak, TaiESM captures the increasing trend of temperature after 1950 well. The current climatology of TaiESM during 1979–2005 is evaluated by observational and reanalysis datasets. Cloud amounts are too large in TaiESM, but their cloud forcing is only slightly weaker than observational data. The mean bias of the sea surface temperature is almost 0, whereas the surface air temperatures over land and sea ice regions exhibit cold biases. The overall performance of TaiESM is above average among models in Coupled Model Intercomparison Project phase 5, particularly in that the bias of precipitation is smallest. However, several common discrepancies shared by most models still exist, such as the double Intertropical Convergence Zone bias in precipitation and warm bias over the Southern Ocean.

Diurnal Sea surface temperature response to tropical cyclone Dahlia in the Eastern tropical Indian Ocean in 2017 revealed by the Bailong buoy

Tropical Cyclone (TC) Dahlia occurred adjacent to over the equatorial southeastern Indian Ocean during the period 26 November – 3 December 2017 and was observed by the Bailong buoy, which provides in situ observations of high-frequency variations in the upper ocean environment. The diurnal sea surface temperature (dSST) variabilities during different stages of the passage of TC Dahlia are studied. The dSST variability is rather weak during the TC passing stage in contrast to the strong ranges before (0.35 °C) and after (0.57 °C) the TC. Before the influence of TC Dahlia, the dSST presented significant regular variability with a peak in the afternoon and minimum value in the morning, which is similar to the even larger range that occurred after TC Dahlia. During the passage of TC Dahlia, dSST decreased dramatically, and a uniform variation was presented due to the absence of strong heat fluxes and stirring and upwelling induced by strong winds. Further analysis through a one-dimensional mixed layer model (Price-Weller-Pinkel, PWP) indicated that the dominant elements responsible for the different dSST variations during distinct stages of TC Dahlia were shortwave radiation and surface wind, which strongly impacted the dSST evolution during TC Dahlia. The asymmetrical wind strength was responsible for the asymmetry of dSST variation.

Southern Ocean Seasonal Restratification Delayed by Submesoscale Wind–Front Interactions

AbstractOcean stratification and the vertical extent of the mixed layer influence the rate at which the ocean and atmosphere exchange properties. This process has direct impacts for anthropogenic heat and carbon uptake in the Southern Ocean. Submesoscale instabilities that evolve over space (1–10 km) and time (from hours to days) scales directly influence mixed layer variability and are ubiquitous in the Southern Ocean. Mixed layer eddies contribute to mixed layer restratification, while down-front winds, enhanced by strong synoptic storms, can erode stratification by a cross-frontal Ekman buoyancy flux. This study investigates the role of these submesoscale processes on the subseasonal and interannual variability of the mixed layer stratification using four years of high-resolution glider data in the Southern Ocean. An increase of stratification from winter to summer occurs due to a seasonal warming of the mixed layer. However, we observe transient decreases in stratification lasting from days to weeks, which can arrest the seasonal restratification by up to two months after surface heat flux becomes positive. This leads to interannual differences in the timing of seasonal restratification by up to 36 days. Parameterizing the Ekman buoyancy flux in a one-dimensional mixed layer model reduces the magnitude of stratification compared to when the model is run using heat and freshwater fluxes alone. Importantly, the reduced stratification occurs during the spring restratification period, thereby holding important implications for mixed layer dynamics in climate models as well as physical–biological coupling in the Southern Ocean.

Water Mass Transformation in the Greenland Sea during the Period 1986–2016

AbstractHydrographic measurements from ships, autonomous profiling floats, and instrumented seals over the period 1986–2016 are used to examine the temporal variability in open-ocean convection in the Greenland Sea during winter. This process replenishes the deep ocean with oxygen and is central to maintaining its thermohaline properties. The deepest and densest mixed layers in the Greenland Sea were located within its cyclonic gyre and exhibited large interannual variability. Beginning in winter 1994, a transition to deeper (>500 m) mixed layers took place. This resulted in the formation of a new, less dense class of intermediate water that has since become the main product of convection in the Greenland Sea. In the preceding winters, convection was limited to <300-m depth, despite strong atmospheric forcing. Sensitivity studies, performed with a one-dimensional mixed layer model, suggest that the deeper convection was primarily the result of reduced water-column stability. While anomalously fresh conditions that increased the stability of the upper part of the water column had previously inhibited convection, the transition to deeper mixed layers was associated with increased near-surface salinities. Our analysis further suggests that the volume of the new class of intermediate water has expanded in line with generally increased depths of convection over the past 10–15 years. The mean export of this water mass from the Greenland Sea gyre from 1994 to present was estimated to be 0.9 ± 0.7 Sv (1 Sv ≡ 106m3s−1), although rates in excess of 1.5 Sv occurred in summers following winters with deep convection.

The Tyrrhenian Intermediate Water (TIW): Characterization and formation mechanisms

This work focuses on the Tyrrhenian Intermediate Water, a water mass present in the Tyrrhenian Sea (Western Mediterranean Sea), which approximately lies between 100 m and 200 m of depth, and is characterized by a local minimum of temperature. Here, for the first time, a thorough investigation of the Tyrrhenian Intermediate Water has been performed, based on the analysis of long time series of observations (XBT, Argo float, dissolved oxygen data) and on modeling results. It was found that this water is present in a large part of the Tyrrhenian basin, and is a persistent feature of the basin hydrology. It is formed by winter convection, with some differences between the southern part of the basin and the northern part, where strong winter cooling produces deeper convection and mixing. The process of formation of the Tyrrhenian Intermediate Water has been investigated in detail, through the analysis of the experimental datasets, of dedicated numerical experiments with a one-dimensional mixed layer model, and of the outputs of a three-dimensional ocean circulation model. Such analysis has excluded remote contributions from intermediate waters produced in the western portion of the Western Mediterranean basin. On the other hand, the presence of the Tyrrhenian Intermediate Water has been found to impact not only the Tyrrhenian Sea hydrology, but also that of the adjacent Liguro-Provençal basin, because this water mass outflows from the Corsica Strait, and is transported as far as the Gulf of Lion, in the core of the region where deep waters are formed. Recent changes in the properties of this intermediate water mass have also been highlighted, which result from the warming of the neighboring water layers.

Salinity Variability Associated with the Positive Indian Ocean Dipole and Its Impact on the Upper Ocean Temperature

Abstract Both surface and subsurface salinity variability associated with positive Indian Ocean dipole (pIOD) events and its impacts on the sea surface temperature (SST) evolution are investigated through analysis of observational/reanalysis data and sensitivity experiments with a one-dimensional mixed layer model. During the pIOD, negative (positive) sea surface salinity (SSS) anomalies appear in the central-eastern equatorial Indian Ocean (southeastern tropical Indian Ocean). In addition to these SSS anomalies, positive (negative) salinity anomalies are found near the pycnocline in the eastern equatorial Indian Ocean (southern tropical Indian Ocean). A salinity balance analysis shows that these subsurface salinity anomalies are mainly generated by zonal and vertical salt advection anomalies induced by anomalous currents associated with the pIOD. These salinity anomalies stabilize (destabilize) the upper ocean stratification in the central-eastern equatorial (southeastern tropical) Indian Ocean. By decomposing observed densities into contribution from temperature and salinity anomalies, it is shown that the contribution from anomalous salinity stratification is comparable to that from anomalous thermal stratification. Furthermore, impacts of these salinity anomalies on the SST evolution are quantified for the first time using a one-dimensional mixed layer model. Since enhanced salinity stratification in the central-eastern equatorial Indian Ocean suppresses vertical mixing, significant warming of about 0.3°–0.5°C occurs. On the other hand, stronger vertical mixing associated with reduced salinity stratification results in significant SST cooling of about 0.2°–0.5°C in the southeastern tropical Indian Ocean. These results suggest that variations in salinity may potentially play a crucial role in the evolution of the pIOD.

Global Estimates of Lateral Springtime Restratification

AbstractSubmesoscale frontal dynamics are thought to be of leading-order importance for stratifying the upper ocean by slumping horizontal density gradients to produce vertical stratification. Presented here is an investigation of submesoscale instabilities in the mixed layer—mixed layer eddies (MLEs)—as a potential mechanism of frontal slumping that stratifies the upper ocean during the transition from winter to spring, when wintertime forcings weaken but prior to the onset of net solar warming. Observations from the global Argo float program are compared to predictions from a one-dimensional mixed layer model to assess where in the world’s oceans lateral processes influence mixed layer evolution. The model underestimates spring stratification for ~75% ± 25% of the world’s oceans. Relationships between vertical and horizontal temperature and salinity gradients are used to suggest that in 30% ± 20% of the oceans this excess stratification can be attributed to the slumping of horizontal density fronts. Finally, 60% ± 10% of the frontal enhanced stratification is consistent with MLE theory, suggesting that MLEs may be responsible for enhanced stratification in 25% ± 15% of the world’s oceans. Enhanced stratification from frontal tilting occurs in regions of strong horizontal density gradients (e.g., midlatitude subtropical gyres), with a small fraction occurring in regions of deep mixed layers (e.g., high latitudes). Stratification driven by MLEs appears to constrain the coexistence of sharp lateral gradients and deep wintertime mixed layers, limiting mixed layer depths in regions of large lateral density gradients, with an estimated wintertime restratification flux of order 100 W m−2.

Turbulence structure in the upper ocean: a comparative study of observations and modeling

Observations of turbulent dissipation rates measured by two independent instruments are compared with numerical model runs to investigate the injection of turbulence generated by sea surface gravity waves. The near-surface observations are made by a moored autonomous instrument, fixed at approximately 8 m below the sea surface. The instrument is equipped with shear probes, a high-resolution pressure sensor, and an inertial motion package to measure time series of dissipation rate and nondirectional surface wave energy spectrum. A free-falling profiler is used additionally to collect vertical microstructure profiles in the upper ocean. For the model simulations, we use a one-dimensional mixed layer model based on a k–ε type second moment turbulence closure, which is modified to include the effects of wave breaking and Langmuir cells. The dissipation rates obtained using the modified k–ε model are elevated near the sea surface and in the upper water column, consistent with the measurements, mainly as a result of wave breaking at the surface, and energy drawn from wave field to the mean flow by Stokes drift. The agreement between observed and simulated turbulent quantities is fairly good, especially when the Stokes production is taken into account.

The effect of non-uniform depth on the baroclinic response to wind in elongated basins

We examined the effect of along-thalweg depth variability on the baroclinic response to wind in elongated narrow basins with a sharp thermocline. The effect of depth variability was examined by deriving a modal-based forced model with two density layers and applying the model to a symmetric curved-bottom basin (CB), an asymmetric wedge-shaped basin (with a sloping bottom towards a vertical wall, WB), and a flat-bottom basin (FB). The baroclinic responses of CB, WB, and FB to uniform wind were found to differ in time-scale, number and energy of excited modes, and temporal pattern and along-thalweg structure of baroclinic flow and thermocline deflection. For all bottom profiles that were examined, the fundamental mode was found to dominate the response to spatially-uniform wind. Compared to FB, the asymmetric depth variability in WB increased the number and energy of excited higher modes and localized the interface shear, while the symmetric deviation from flat bottom in CB caused the opposite effects. Linear deviation from uniform wind was found to feed energy into higher baroclinic modes for the symmetric CB, but was found to reduce the energy of higher baroclinic modes for WB when the deviation from uniform wind is comparable to the spatial-average magnitude. Our results can explain the observation of the second baroclinic mode and irregular wave patterns in some lakes and reservoirs. Further, our results suggest that one-dimensional vertical mixed-layer models provide better results for shear entrainment in curved-bottom basins than in wedge-shaped basins.

Mechanisms behind the Temporary Shutdown of Deep Convection in the Labrador Sea: Lessons from the Great Salinity Anomaly Years 1968–71

Abstract From 1969 to 1971 convection in the Labrador Sea shut down, thus interrupting the formation of the intermediate/dense water masses. The shutdown has been attributed to the surface freshening induced by the Great Salinity Anomaly (GSA), a freshwater anomaly in the subpolar North Atlantic. The abrupt resumption of convection in 1972, in contrast, is attributed to the extreme atmospheric forcing of that winter. Here oceanic and atmospheric data collected in the Labrador Sea at Ocean Weather Station Bravo and a one-dimensional mixed layer model are used to examine the causes of the shutdown and resumption of convection in detail. These results highlight the tight coupling of the ocean and atmosphere in convection regions and the need to resolve both components to correctly represent convective processes in the ocean. They are also relevant to present-day conditions given the increased ice melt in the Arctic Ocean and from the Greenland Ice Sheet. The analysis herein shows that the shutdown was initiated by the GSA-induced freshening as well as the mild 1968/69 winter. After the shutdown had begun, however, the continuing lateral freshwater flux as well as two positive feedbacks [both associated with the sea surface temperature (SST) decrease due to lack of convective mixing with warmer subsurface water] further inhibited convection. First, the SST decrease reduced the heat flux to the atmosphere by reducing the air–sea temperature gradient. Second, it further reduced the surface buoyancy loss by reducing the thermal expansion coefficient of the surface water. In 1972 convection resumed because of both the extreme atmospheric forcing and advection of saltier waters into the convection region.

Water-mass transformation in the shelf seas

The rate at which water masses are transformed from one density class to another is assessed in the shelf seas using the Walin (1982) framework. For a tidal-mixing front, the transformation is estimated using air-sea density fluxes and the diapycnal mixing diagnosed from a series of one-dimensional mixed layer models running across the shelf. These transformation rates diagnosed from the air-sea fluxes and diapycnal mixing agree with volume changes dagnosed directly from the model. The transformation from air-sea fluxes reaches a maximum amplitude typically twice that provided by diffusive mixing. This framework is extended to estimate the rate at which nutrients are converted in nutrient space including the effects of biological consumption. The transformation in density and nutrient space are broadly related in the spring when there is a relatively tight relationship between density and nutrient concentrations. For a shelf-break front in the Celtic Seas, the transformation is estimated from a combination of observed air-sea fluxes, remotely-sensed sea surface temperature and ship-based measurements of density and turbulent mixing. The transformation is controlled by diapycnal mixing along dense surfaces and by air-sea fluxes for lighter surfaces, and each contribution reaches comparable magnitudes over a six-month period

Estimating Air–Sea Heat Fluxes in Semienclosed Basins: The Case of the Gulf of Elat (Aqaba)

Abstract Meteorological and oceanographic data collected at the head of the Gulf of Elat were used to compute the air–sea heat flux components and the heat storage in the water column, which are in turn used to estimate the heat balance of this semienclosed basin. The solar radiation was measured directly, whereas the longwave (LW) cooling and the turbulent heat fluxes (latent, LH; sensible, SH) were computed from commonly used bulk formulas. Nine formulas for LW and four formulas for LH + SH were tested for a total of 36 possible combinations. Independent estimates for the bounds on the advective heat flux through the straits and results from a one-dimensional mixed layer model provided criteria to help identify the best choice of bulk formulas for the gulf. It was concluded that the LW formula of Bignami together with the turbulent flux formulas of Kondo provide the best estimate of the heat balance of the gulf. Based on this, the annual mean evaporation is 1.6–1.8 m yr−1, with a minimum of 1 m yr−1 in (the long) summer and a maximum of 3–4 m yr−1 in (the short) winter. The increase in evaporation rate during the winter results from the instability of the atmosphere at that time when the sea surface temperature exceeds the air temperature; in the summer, when the air temperature is much higher than the sea surface temperature, evaporation nearly stops due to the atmospheric stability. This estimated evaporation rate for the gulf, which is similar for all four of the LH formulas considered, is significantly smaller than values commonly quoted in the literature. Finally, in contrast to previous studies, it is found that the advective heat flux from the Straits of Tiran is large and significant in spring, reaching an estimated value of over 125 W m−2, but its annually averaged value is only about 35–40 W m−2.

The influence of surface layer salinity on wintertime convection in Wilkinson Basin, Gulf of Maine

Observations suggest that the interannual variability of wintertime convection in Wilkinson Basin (WB) during 1977–2005 is related to the variation in the surface layer salinity in the western Gulf of Maine (wGOM). When the winter convection is stronger (weaker), the resulting deep-layer temperatures in the wGOM are colder (warmer), and the likelihood of deep winter mixing is greater when the wGOM salinities are high. This hypothesis was tested using a one-dimensional mixed layer model to simulate the water column structure over the cooling period. Comparisons were made between the convection potential for the range of observed late-fall salinity values. Model results indicated that the mixed layer depth could deepen by as much as 50% when the surface layer salinities in the wGOM were high. Sensitivity tests revealed that the surface layer salinity variability, and its influence on the water column density structure, was a significant factor, along with local climate variability, in determining the observed time series trend in springtime bottom layer temperatures. Density distributions in the wGOM in high salinity years suggested that winter-cooled water from the adjacent shallow coastal areas could become sufficiently dense to cascade into the deeper layers of WB and enhanced the convective cooling. The interannual variation in wintertime convection and in the resulting density structure has important implications for biological processes. The associated variation in deeper layer temperatures has been shown to influence the spatial distribution, timing of spawning, and growth rates of many demersal and mid-water fishery populations.

Winter Mixed Layer Development in the Central Irminger Sea: The Effect of Strong, Intermittent Wind Events

Abstract The impact of the Greenland tip jet on the wintertime mixed layer of the southwest Irminger Sea is investigated using in situ moored profiler data and a variety of atmospheric datasets. The mixed layer was observed to reach 400 m in the spring of 2003 and 300 m in the spring of 2004. Both of these winters were mild and characterized by a low North Atlantic Oscillation (NAO) index. A typical tip jet event is associated with a low pressure system that is advected by upper-level steering currents into the region east of Cape Farewell and interacts with the high topography of southern Greenland. Heat flux time series for the mooring site were constructed that include the enhancing influence of the tip jet events. This was used to force a one-dimensional mixed layer model, which was able to reproduce the observed envelope of mixed layer deepening in both winters. The deeper mixed layer of the first winter was largely due to a higher number of robust tip jet events, which in turn was caused by the steering currents focusing more storms adjacent to southern Greenland. Application of the mixed layer model to the winter of 1994–95, a period characterized by a high-NAO index, resulted in convection exceeding 1700 m. This prediction is consistent with hydrographic data collected in summer 1995, supporting the notion that deep convection can occur in the Irminger Sea during strong winters.

Observations and models of the energy flux from the wind to mixed-layer inertial currents

Observations of wind stress, upper-ocean currents, and stratification from a mooring at 10∘N in the Pacific are used along with two different one-dimensional mixed-layer models to investigate the horizontal kinetic energy (KE) balance for inertial motions in the mixed-layer. Results from four other sites spanning 25∘N to 60∘N in the Atlantic also are presented. Determining the work done by the wind on the ocean during strong resonant forcing events is of particular interest. The damped-slab mixed-layer model used in previous studies is appealing for this purpose, since it is possible to estimate the wind work using only the wind stress, prescribed mixed-layer depth, and an empirical damping constant. However, it is found that this model, which does not allow for interaction between the mixed-layer and the stratified water column below, does not reproduce the observed KE balance for strong forcing events and systematically over-estimates the wind work. A slightly more sophisticated model, which includes a shear instability mechanism to facilitate momentum transfer between the mixed-layer and a transition layer below, reproduces the observed KE balance and the cumulative wind work much more accurately. Mixed-layer inertial currents unrelated to wind forcing (e.g., due to upward propagating internal waves) are found in the observations, but are not accounted for in either model. When wind forcing is relatively weak, these background currents may be of comparable amplitude to wind-forced currents, and can have a significant impact on the inertial response.

Impact of the Diurnal Cycle of Solar Radiation on Intraseasonal SST Variability in the Western Equatorial Pacific

Abstract The mechanism by which the diurnal cycle of solar radiation modulates intraseasonal SST variability in the western Pacific warm pool is investigated using a one-dimensional mixed layer model. SSTs in the model experiments forced with hourly surface fluxes during the calm–sunny phase of intraseasonal oscillation are significantly warmer than those with daily mean surface fluxes. The difference in two experiments is explained by upper-ocean mixing processes during nighttime. Surface warming during daytime creates a shallow diurnal warm layer near the surface (0–3 m), which can be easily eroded by surface cooling during nighttime. Further cooling, however, requires a substantial amount of energy because deeper waters need to be entrained into the mixed layer. Since the shallow diurnal layer is not formed in the experiment with daily mean surface fluxes, the SST for the hourly forcing case is warmer most of the time due to the diurnally varying solar radiation. Sensitivity of the intraseasonal SST variation to the penetrative component of solar radiation is examined, showing that the diurnal cycle plays an important role in the sensitivity. Solar radiation absorbed in the upper few meters significantly influences intraseasonal SST variations through changes in amplitude of diurnal SST variation.

A modelling study of ice‐spectrum modes in deep frontal clouds

AbstractResults from a one‐dimensional mixed‐layer model with explicit microphysics and a discrete ice‐spectrum were compared to in situ data provided by the Met Office C‐130 aircraft. Bimodal ice‐spectra were observed. The first mode was found at diameters near the lower limit of the instrumentation resolution (25 µm) at a roughly constant concentration. The second mode was found initially at diameters exceeding 200 µm but it shifted to larger sizes and smaller peak concentrations as depths below the cloud top increased.The objectives of the paper are threefold: validation of the model with real cloud data; an analysis of parameters controlling the rate of aggregation, and to show that a secondary ice‐production process is needed to maintain the first mode and produce a second mode which is sensible.A reasonable first mode could only be maintained by producing a constant supply of small ice particles at a rate of at least 10 m−3s−1. Primary nucleation was unable to supply this rate of ice‐particle production at temperatures above −20°C. In the model, it was achieved by invoking a secondary ice‐particle production process, such as the break‐up of evaporating ice crystals. The rate of production required is in broad agreement with laboratory investigations of secondary ice‐particle production during evaporation.The relative importance of the processes of aggregation, deposition and sedimentation to the production of the bimodal spectra is also studied. While it is possible to grow ice particles to large sizes by deposition alone, the number concentration observed in these clouds produces so much competition for water vapour that deposition growth is very limited. The aggregation process is the strongest mechanism for the production of the larger ice‐particles, which then precipitate. Consequently, its strength controls the rate of small ice‐particle production needed, and so is mainly responsible for the second mode. Copyright © Royal Meteorological Society, 2003. P. R. Field's contribution is Crown copyright.

A Parameterization of the Roughness Length for the Air-Sea Interface in Free Convection

The response of the upper ocean to the parameterization of the roughness length z0 on the air side of the air-sea interface is studied using a one-dimensional mixed-layer model. In particular, it is shown that in the free convection limit when both the wind speed and the friction velocity approach zero, the familiar Charnock formula for the momentum roughness, which relies solely on wind gen- eration, can be modified to account for contributions arising from the thermally generated turbulence. Therefore, a new parameterization is proposed for the momentum roughness length which extends the Charnock formula down to zero friction velocity. The value of a parameter which enters in the new formulation is determined by making use of exsisting free convection surface flux parameterizations. The effect of the new parameterization on the model performance is tested using data from the ocean weathership station Papa (OWS P), and data from the Long-Term Upper-Ocean Study (LOTUS) experiment. Simulations were carried out using a recently developed one-dimensional, second-order, turbulence closure scheme over diurnal as well as seasonal time scales. The findings suggest that the new momentum roughness parameterization improves the overall agreement between the observed and simulated sea-surface temperature (SST).