The cross‐channel flow at the entrance of Lancaster Sound

Abstract

A westward current flows along the northern side of Lancaster Sound and an eastward current flows along the southern side. A cross‐channel flow is commonly observed to link them near the eastern entrance of Lancaster Sound; this flow is modelled assuming inviscid flow and conservation of potential vorticity. It is shown that the westward decrease in depth is sufficient to cause a cross‐channel flow that couples the inflow to the outflow. The modelled cross‐channel flow takes place at a distance inside the entrance that is less than that observed for the surface current. Obviously stratification reduces the coupling of the surface current to the bathymetry. A more realistic result is obtained with the barotropic model if the bottom slope is halved.

An inviscid mean barotropic flow out of the channel is also modelled and found to be concentrated on the southern bank in order to conserve potential vorticity. It seems that barotropic instability and friction would limit the narrowing of the flow.     

Boundary‐layer influence on the development of a mesoscale rainstorm

Abstract

A major convective rainstorm slowly developed in the largely barotropic warm sector of an open‐wave cyclone. Meso‐beta‐scale boundary‐layer conditions played a significant role in the location of storm development, and provided favourable thermodynamic gradients and persistently convergent flow in preferred areas. Minimal storm movement and a strong association between the surface divergence field and cloud evolution permitted the estimation of that fraction of the moisture fed into the storm from the boundary layer that was returned to the ground as rain.     

Impressions of the turbulence variability in a weakly stratified,flat-bottom deep-sea ‘boundary layer’

The character of turbulent overturns in a weakly stratified deep-sea is investigated in some detail using 144 high-resolution temperature sensors at 0.7 m intervals, starting 5 m above the bottom. A 9-day, 1 Hz sampled record from the 912 m depth flat-bottom (<0.5% bottom-slope) mooring site in the central-north Alboran Sea (W-Mediterranean) demonstrates an overall conservative temperature range of only 0.05 °C, a typical mean buoyancy period as large as 3 h and a 1 Hz-profile-vertically-averaged turbulence dissipation rate maximum of only 10⁻⁸ m² s⁻³. Nonetheless, this ‘boundary layer’ varies in height between <6 and >104 m above the bottom and is thus not homogeneous throughout; the temperature variations are seldom quiescent and are generally turbulent in appearance, well exceeding noise levels. The turbulence character is associated with small-scale internal waves; examples are found of both shear- and convection-driven turbulence; particular association, although not phase-locked, is found between turbulence variations and tidal rather than with inertial motions; the mean buoyancy frequency of a few times the inertial frequency implies the importance of ‘slantwise convection’ in the direction of the earth rotational vector rather than in the direction of gravity. Such convection is observed both in near-homogeneous and weakly stratified form.     

Sea surface temperature variability in the shelf‐slope region of the Northwest Atlantic

Abstract

Sea surface temperature (SST) variability in the shelf‐slope region of the northwest Atlantic is described and then explained in terms of latent and sensible heat exchange with the atmosphere. The basic data are primarily engine‐intake temperature measurements made by merchant ships over the period 1946–80. The data have been grouped by month and area and an empirical orthogonal function analysis has been performed to determine the dominant modes of variation. The first two modes account for 44% of the total variance. The first mode corresponds to in‐phase changes of SST from the Grand Banks to Mid‐Atlantic Bight; the second mode corresponds to opposite changes of SST on the Grand Banks and Mid‐Atlantic Bight. The time‐dependent amplitudes of these large‐scale modes have pronounced low‐frequency components; the associated changes in SST are typically 3°C. It is also shown that winter anomalies last longer than summer anomalies; their typical decay scales are 6 and 3 months, respectively.

The onshore component of geostrophic wind is significantly correlated with the amplitude of the first mode in winter. We note the strong land‐sea contrast of temperature and humidity in this region during winter and explain the wind‐SST correlation in terms of latent and sensible heat exchanges. The second mode (i.e. the difference in SST between the Grand Banks and Mid‐Atlantic Bight) also appears to be related to changes in atmospheric circulation during the winter. A stochastic model for mixed layer temperature is finally used to model the SST autocorrelation functions. Following Ruiz de Elvira and Lemke (1982), it includes a seasonally‐varying feedback coefficient. The model successfully reproduces the extended persistence of winter anomalies with physically realistic parameter values but it cannot account for the summer reinforcement of winter anomalies on the Scotian Shelf. We speculate that this is due to the occasional entrainment of water, cooled the previous winter, into the shallow summer mixed layer.     

Application of a coupled ice‐ocean model to the Labrador Sea

Abstract

The steady, coupled ice‐ocean circulation model of Willmott and Mysak (1989) for a meridional channel is applied to the Labrador Sea for the winter season. The model consists of a thermodynamic reduced‐gravity ocean combined with a variable thickness ice cover that is in thermal equilibrium. Upon specifying the forcing fields of surface air temperature, wind stress and water temperature along the open southern boundary, the winter climatological ice‐edge position, ice thickness, ocean circulation and temperature fields are determined in the channel domain. The sensitivity of the results to the various model parameters is examined. In particular, the optimum heat exchange coefficients for the interfaces of air‐water, ice‐water and air‐ice are found.

The model ice‐edge position compares favourably with the 50% winter climatological ice concentration isoline obtained from an analysis of 32 years (1953–84) of sea‐ice concentration data. The simulations of the ocean temperature and ice thickness are also quite realistic according to the observed records available. The model is also applied to two specific winters (1981 and 1983) during which anomalous sea‐ice and weather conditions prevailed in the Labrador Sea.     

Large‐eddy simulation of small‐scale surface effects on the convective boundary‐layer structure

Abstract

The effects of small‐scale surface inhomogeneities on the turbulence structure in the convective boundary layer are investigated using a high‐resolution large‐eddy simulation model. Surface heat flux variations are sinusoidal and two‐dimensional, dividing the total domain into a checkerboard‐like pattern of surface hot spots with a 500‐m wavelength in the x and y directions, or 1/4 of the domain size. The selected wind speeds were 1 and 4 m s^‐l, respectively. As a comparison, a simulation of the turbulence structure was performed over a homogeneous surface.

When the wind speed is light, surface heat flux variations influence the horizontally averaged turbulence statistics, including the higher moments despite the small characteristic length of the surface perturbation. Stronger mean wind speeds weaken the effects of inhomogeneous surface conditions on the turbulence structure in the convective boundary layer.

Results from conditional sampling show that when the mean wind speed is small, weak mean circulations occur, with updraft branches above the high heat flux regions and down‐draft branches above the low heat flux regions. The inhomogeneous surface induces significant differences in the turbulence statistics between the high and low heat flux regions. However, the effect of the surface perturbations weaken rapidly when the mean wind speed increases. This research has implications in the explanation of the large‐scale variability commonly encountered in aircraft observations of atmospheric turbulence.     

Prairie agroclimate boundary‐layer model: A simulation of the atmosphere/crop‐soil interface

Abstract

This paper describes a 1‐D agroclimatic model of the atmosphere/crop‐soil interface. Vertical profiles of wind, potential temperature and water vapour are constructed twice daily for the overnight‐low and maximum temperature times by combining 1200 and 0000 UTC upper‐air standard‐level grid‐point data with climatological observations. The vertical structure of the atmospheric boundary layer has a surface constant‐flux layer that is usually topped by a mixed layer by day but not at night. The crop‐soil boundary layer consists of a shallow top‐zone and a growing root‐zone. Vegetation cover and root depth depend upon crop type and phenological stage. Water‐balance accounting tracks the moisture contents of both the top‐ and root‐zones. Evapotranspiration or the vertical flux of water vapour in the atmospheric boundary layer is tied to the evolution of the crop‐soil boundary layer.

The model was calibrated using field data from the Regional Evaporation Study's primary site in an agricultural area of central Saskatchewan. The evolution of 1991's wheat‐soil boundary layer from the crop's heading to ripe stages was then successfully simulated at two additional sites in the same geographical area.     

Estimate of the steric contribution to global sea level rise from a comparison of the WOCE one‐time survey with 2006–2008 Argo observations

Abstract

It is well known from observations by altimetric satellites (predominantly TOPEX/Poseidon and Jason‐1) that global sea level is rising. What is less well known is exactly how the observed sea level rise is partitioned between a steric contribution (sea level rising because of changes in ambient temperature and salinity) and a contribution arising from the addition of new water mass to the oceans. Strictly speaking, such a separation is not possible because of the non‐linearity in the equation of state for sea water, but in practice the non‐linearities are sufficiently small to allow this separation as a very good first approximation.

A careful comparison of the World Ocean Circulation Experiment (WOCE) one‐time survey with recent observations by the Argo array indicate a steric component to sea level rise of 2.2 mm y^–1 between the early 1990s and 2006 to 2008. This is a significantly larger rise rate than previously estimated and, along with recent estimates of melt rate from ice sheets, is in much closer agreement with the total rise rate as reported by altimetric satellites, 3.2 ± 0.4 mm y^–1 over this period.     

On the efficiency of horizontal diffusion and numerical filtering in an Arakawa‐type model

Abstract

In a series of 5‐day forecasts with a 3‐layer/2.8° hemispheric model, horizontal diffusion schemes of the second and fourth degree are compared with a numericalfilter technique. The results, which are discussed mainly in terms of spectral energetics in zonal wavenumber space, indicate that fourth‐degree diffusion is more scale selective than second‐degree and equivalent to filtering. The seventh‐order filter applied only intermittently to the prognostic variables is superior to fourth‐degree diffusion from the viewpoint of computational economy. Excessive dissipation of the long waves may inhibit the production of eddy kinetic energy from eddy available potential energy.     

Surface‐absorbed and top‐of‐atmosphere radiation fluxes for the Mackenzie River basin from satellite observations and a regional climate model and an evaluation of the model

Abstract

Both the earth‐reflected shortwave and outgoing longwave radiation (OLR) fluxes at the top of the atmosphere (TOA) as well as surface‐absorbed solar fluxes from Canadian Regional Climate Model (CRCM) simulations of the Mackenzie River Basin for the period March 2000 to September 2003 are compared with the radiation fluxes deduced from satellite observations. The differences between the model and satellite solar fluxes at the TOA and at the surface, which are used in this paper to evaluate the CRCM performance, have opposite biases under clear skies and overcast conditions, suggesting that the surface albedo is underestimated while cloud albedo is overestimated. The slightly larger differences between the model and satellite fluxes at the surface compared to those at the TOA indicate the existence of a small positive atmospheric absorption bias in the model. The persistent overestimation of TOA reflected solar fluxes and underestimation of the surface‐absorbed solar fluxes by the CRCM under all sky conditions are consistent with the overestimation of cloud fraction by the CRCM. This results in a larger shortwave cloud radiative forcing (CRF) both at the TOA and at the surface in the CRCM simulation. The OLR from the CRCM agrees well with the satellite observations except for persistent negative biases during the winter months under all sky conditions. Under clear skies, the OLR is slightly underestimated by the CRCM during the winter months and overestimated in the other months. Under overcast conditions the OLR is underestimated by the CRCM, suggesting an underestimation of cloud‐top temperature by the CRCM. There is an improvement in differences between model and satellite fluxes compared to previously reported results largely because of changes to the treatment of the surface in the model.     

Lightning‐flash locations related to the precipitation pattern of the storm

Abstract

Triangulation from paired sferics (100 ± 50 kHz) received at stations 57 km apart gave 42 map locations of lightning flashes belonging to a thunderstorm well‐defined on weather radar. This storm was one of many in a long line along the leading edge of a band of light stratiform precipitation (snow above 3 km). The band was 9 km high and 60 km wide, and the storms were embedded in it. The convective region of the storm studied had a north‐south extent of 20 km, the east‐west extent averaging two‐thirds as much. Cumulus heights progressed from 3 km on the west to 7 ± 1 km on the east. Two‐thirds of the lightning flashes were in the convective region. The remaining third were in 160 km² of the light stratiform precipitation immediately west of it.     

Application of the delta‐Eddington method to the absorption of solar radiation in the atmosphere

Abstract

A model to compute rapidly the absorption of solar radiation in the atmosphere is described. The model is based partially on the parameterization of Lacis and Hansen and also makes use of the delta‐Eddington method. In addition to absorption by ozone and water vapour, and scattering by air molecules and clouds, the mode1 includes absorption and scattering by aerosols. Good agreement is found in comparison with the Lacis and Hansen parameterization in the absence of clouds and aerosol. The present model represents an improvement in the treatment of scattering by clouds. Its main advantage though, is in its flexibility in allowing for interactions with the atmospheric aerosol.     

Near‐bottom currents and sediment transport on the inner Scotian Shelf: Sea‐floor response to winter storms during CASP

Abstract

As part of the Canadian Atlantic Storms Program (CASP), near‐bottom current velocity, pressure, light transmission (as a measure of suspended sediment concentration) and water temperature were recorded using a variety of instruments deployed in water depths of 20 to 37 m on the inner Scotian Shelf, during February and March 1986. Detailed mapping of a 12‐km² area encompassing the instrument mooring sites revealed a variety of bottom types. These include sand and gravel (both forming ripples at various scales), cobble‐boulder lags, and bedrock, resulting in bottom roughness estimates that vary widely (10^?4 m < k < 10° m) over short horizontal distances (of the order 10² m). The velocity data provided information on the near‐bottom current response to winter storms anda basis for computations of sediment load and transport rates. The near‐bottom mean flow showed distinct storm‐driven circulation patterns, with velocities roughly parallel to alongshore wind stress but opposed to shore‐normal wind. Wave‐induced oscillatory motions also showed marked increases during storms and frequently dominated the near‐bottom flow. Sediment load (depth‐integrated concentration) and transport were computed using a model in which the load is related to the excess normalized shear stress. The computed mean concentrations were compatible with the optical transmis someter data. These computations yielded estimates ranging up to 0.7 kg m^?2 for the mean sediment load and 443 kg m^?1 h^?1 for the net transport. Hindcast scour rates, of the order of 1 mm h^?1 under moderate storm conditions were generally compatible with depths of scour measured by divers.     

A global available potential energy‐kinetic energy budget in terms of the two‐dimensional wavenumber for the FGGE year

Abstract

A global vertically integrated available potential energy‐kinetic energy budget in terms of the two‐dimensional wavenumber is formulated using spherical harmonics. Results of the budget equations applied to the four mid‐season months of the FGGE year are given.     

Reply to comments on and addenda to “some properties and comparative performance of the semi‐lagrangian method of Robert in the solution of the advection‐diffusion equation”

Abstract

Droughts are major natural disasters for many parts of the world. Dry areas, where the precipitation pattern is markedly seasonal, or is otherwise highly variable, are the most susceptible. The Canadian Prairies, together with the U.S. Great Plains, are one such area. While immediate loss of life is seldom a feature of most droughts, malnutrition and even starvation do follow severe droughts in some parts of the world. In Canada, economic losses, particularly in the agricultural sector, may reach several hundred millions of dollars in a drought year, with major socio‐economic repercussions affecting the entire region. Environmental damages include soil degradation and erosion, vegetation damage, slough and lake deterioration and wildlife loss. Unlike most other natural disasters, drought onset is difficult to identify. Droughts develop slowly, and until human activity begins to be affected by an on‐going reduction of precipitation, their existence is unrecognized. Development and application of specific soil moisture and drought indices based on cumulative precipitation deficits have enhanced drought monitoring programs. These in turn provide guidance on the need for mitigative measures that can be initiated early in the course of a drought. Any improvement in the timely application of these measures is, however, strongly dependent on being able to determine the drought's course, extent and likely severity at the earliest stage possible. The identification of precursor conditions for past drought has raised the possibility that the likelihood of a drought occurring in a particular year or growing season might be predictable. Teleconnections between North American precipitation patterns and ENSO events and other surface boundary conditions in the North Pacific have been detected. Forecasting seasonal temperature and perhaps precipitation anomalies appears to be potentially feasible using a suitable merging of precursor parameters and modelling methodologies. Clearly, future research must focus both on those precursors that have been identified and on a search for possible new ones. Development of better forecasting methodology is also essential. Research activity to identify and evaluate new mitigative measures should also be increased to keep pace with the prospects of drought predictability.     

A decision‐analytic study of the joint value of seasonal precipitation and temperature forecasts in a choice‐of‐crop problem

Abstract

A static decision‐analytic method is used to investigate the economic value of bivariate ‐ precipitation and temperature ‐ seasonal forecasts of the form currently issued by the U.S. National Weather Service. This method is applied to a corn versus spring wheat choice‐of‐crop decision‐making problem by considering a transect of four counties across the northwestern margin of the North American corn belt. Numerical results indicate that seasonal forecasts of current quality can be of appreciable value (≥$1/ha) for some locations when the optimal action chosen on the basis of climatological information is only marginally preferred to another action. Increases in forecast value follow from hypothetical increases in the quality of both the precipitation and temperature components of the forecasts in the spring wheat region, whereas forecast value increases primarily as a function of the quality of the precipitation forecasts alone in the corn belt region. The results are very sensitive to absolute and relative crop prices.     

Implementation of a large‐scale variable velocity river flow routing algorithm in the Canadian Regional Climate Model (CRCM)

Abstract

Implementation and validation of a flow routing scheme for the North American domain of the Canadian Regional Climate Model (CRCM) is described. A variable velocity flow routing algorithm is used to transport runoff from the land surface to the continental edges and provide freshwater flux forcing for the oceans. The flow routing scheme uses Manning's equation to estimate flow velocities for river channels whose cross‐sections are assumed to be rectangular. Discretization of major North American river basins and their flow directions are obtained at the polar stereographic resolution of the CRCM using 5‐minute global river flow direction data as a template. In the absence of observation‐based gridded estimates of runoff, model runoff estimates from a global simulation of the Variable Infiltration Capacity (VIC) hydrological model (forced with observationbased meteorological data) are used to validate the flow routing scheme. Model results show that the inclusion of flow routing improves the comparison with observation‐based streamflow estimates when compared to the unrouted runoff. Monthly comparison of simulated streamflow with observation‐based estimates, and basin‐wide averaged flow velocities, suggests that the flow routing scheme performs satisfactorily.     

The climatology of the Canadian Climate Centre general circulation model as obtained from a five‐year simulation

Abstract

The climatology of the Canadian Climate Centre atmospheric general circulation model is presented and compared with the observed climatology of the atmosphere. The model climatology is obtained from a simulation over five annual cycles and the results are presented in terms of averages for the four seasons.

The climatology of the model is discussed in terms of zonally and time averaged values of the primary atmospheric variables as well as in terms of the spatial distributions of the important surface parameters and of the rotational and divergent components of the tropospheric flow. Some measure of model variability is also presented.

The model is generally quite successful in reproducing the mean observed climatology of the atmosphere.     

The use of σ T as a temperature scale in the atmospheric surface layer

The standard deviation of temperature _T is proposed as a temperature scale and as a velocity scale to describe the behaviour of turbulent flows in the Atmospheric Surface Layer (ASL), instead of _* andu _* of the Monin—Obukhov similarity theory, and ofT _f andU _f used for free convection stability conditions. On the basis of experimental evidence reported in the literature, it is shown that _T T _f andv _* U _f in the free convection region, and _{T *} andv _* U _* in nearneutral and stable conditions. This implies that the proposed scales can be applied for all stabilities. Furthermore, a new length scale is proposed and its relation with Obukhov length is given. Also, a simple semi-empirical expression is presented with which _T andv _* can be evaluated in a rather simple way. Some examples of practical applications are given, e.g., a stability classification for unstable conditions.     

Correction to “the vertical distribution of climate forcings and feedbacks from the surface to top of atmosphere”

We discuss here a mistake in the analysis of Previdi and Liepert (Clim Dyn, 2011). In that article, the surface albedo radiative kernels were calculated incorrectly. We present in this brief comment the corrected albedo kernels. We then use these kernels to compute the surface albedo radiative feedback in climate model simulations driven by increasing carbon dioxide, as in Previdi and Liepert (Clim Dyn, 2011). We find that the use of the corrected albedo kernels does not change the conclusions of our earlier work.