Processes controlling surface, bottom and lateral melt of Arctic sea ice in a state of the art sea ice model

We present a modelling study of processes controlling the summer melt of the Arctic sea ice cover. We perform a sensitivity study and focus our interest on the thermodynamics at the ice–atmosphere and ice–ocean interfaces. We use the Los Alamos community sea ice model CICE, and additionally implement and test three new parametrization schemes: (i) a prognostic mixed layer; (ii) a three equation boundary condition for the salt and heat flux at the ice–ocean interface; and (iii) a new lateral melt parametrization. Recent additions to the CICE model are also tested, including explicit melt ponds, a form drag parametrization and a halodynamic brine drainage scheme. The various sea ice parametrizations tested in this sensitivity study introduce a wide spread in the simulated sea ice characteristics. For each simulation, the total melt is decomposed into its surface, bottom and lateral melt components to assess the processes driving melt and how this varies regionally and temporally. Because this study quantifies the relative importance of several processes in driving the summer melt of sea ice, this work can serve as a guide for future research priorities.

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating.

Surface Urban Energy and Water Balance Scheme (SUEWS): development and evaluation at two UK sites

The Surface Urban Energy and Water Balance Scheme (SUEWS) is evaluated at two locations in the UK: a dense urban site in the centre of London and a residential suburban site in Swindon. Eddy covariance observations of the turbulent fluxes are used to assess model performance over a twoyear period (2011-2013). The distinct characteristics of the sites mean their surface energy exchanges differ considerably. The model suggests the largest differences can be attributed to surface cover (notably the proportion of vegetated versus impervious area) and the additional energy supplied by human activities. SUEWS performs better in summer than winter, and better at the suburban site than the dense urban site. One reason for this is the bias towards suburban summer field campaigns in observational data used to parameterise this (and other) model(s). The suitability of model parameters (such as albedo, energy use and water use) for the UK sites is considered and, where appropriate, alternative values are suggested. An alternative parameterisation for the surface conductance is implemented, which permits greater soil moisture deficits before evaporation is restricted at non-irrigated sites. Accounting for seasonal variation in the estimation of storage heat flux is necessary to obtain realistic wintertime fluxes.

Critical flux in ultrafiltration of skimmed milk

he best operating conditions, using the critical flux concept during ultrafiltration of skimmed milk, were evaluated for tubular membranes. It was found that irreversible fouling was greatly reduced by operating at or below the critical flux, but was not totally eliminated. The critical flux of skimmed milk was found to be the weak form. The critical flux at cross flow velocity 3.4 in s(-1) for MWCO 200 kDa membrane was 56.9 kg m(-2) h(-1) while for MWCO 25 kDa membranes it was 45 kg m(2) h(-1) suggesting that membrane pore size influenced the flux. The critical flux increased with increasing wall shear stress and decreased with increasing protein concentration. Empirical equations, for predicting the critical flux (J(crit)) for skimmed milk with a protein concentration (c(b)) in the range 3-7% w/w and wall shear stress (tau(w)) in the range 7-60 Pa for MWCO 200 kDa and 25 kDa membranes were J(crit) = 5.1 (tau(w)/c(b)) and J(crit) = 4.0 (tau(w)/c(b)) respectively. In general, the rejections of protein and lactose at the critical flux were not affected by protein concentration, wall shear stress and membrane used, and they were similar to those found when operating at the limiting flux.

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.

Heat-inactivated peroxidases and the role of calcium abstraction as a cause of their enhanced lipid oxidation activity: potential effects on the flavour quality of heat-processed vegetables

Cationic swede and anionic turnip peroxidases were partially purified by ion-exchange and gel-filtration chromatography, respectively. Heat treatment of these enzymes and of a commercial high purity horseradish peroxidase (HRP) caused a loss of enzyme activity and a corresponding increase in linoleic acid hydroperoxide formation activity. The hydroperoxide levels in model systems increased only in the early stages of the oxidation reaction and then declined as degradation became more significant. The presence of a dialysed blend of cooked swede markedly lowered the hydroperoxide level formed. Analysis of volatile compounds formed showed that hexanal predominated in a buffer system and in a blend of cooked turnip. In dialysed blends of cooked swede, hexanol was the primary volatile compound generated. After inactivation under mild conditions in the presence of EDTA, the peroxidases showed hydroperoxide formation activity and patterns of volatile compounds from linoleic acid that were similar to those found on heat-inactivation. This suggested that calcium abstraction from the peroxidases was critical for the enhancement of lipid oxidation activity. (C) 2002 Elsevier Science Ltd. All rights reserved.

The climate response to CO2 of the Hadley Centre coupled AOGCM with and without flux adjustment

Coupled atmosphere‐ocean general circulation models have a tendency to drift away from a realistic climatology. The modelled climate response to an increase of CO2 concentration may be incorrect if the simulation of the current climate has significant errors, so in many models, including ours, the drift is counteracted by applying prescribed fluxes of heat and fresh water at the ocean‐atmosphere interface in addition to the calculated surface exchanges. Since the additional fluxes do not have a physical basis, the use of this technique of “flux adjustment” itself introduces some uncertainty in the simulated response to increased CO2. We find that the global‐average temperature response of our model to CO2 increasing at 1% per year is about 30% less without flux adjustment than with flux adjustment. The geographical patterns of the response are similar, indicating that flux adjustment is not causing any gross distortion. The reduced size of the response is due to more effective vertical transport of heat into the ocean, and a somewhat smaller climate sensitivity. Although the response in both cases lies within the generally accepted range for the climate sensitivity, systematic uncertainties of this size are clearly undesirable, and the best strategy for future development is to improve the climate model in order to reduce the need for flux adjustment.

Investigation of flux in flat-plate modules for membrane distillation

A theoretical model for predicting the behaviour of membrane distillation by incorporating mass and heat transfer equations has been used to find permeate fluxes, and has been validated experimentally. The model accurately predicts mass and heat transfer. The main work studied the effect of module design using a flat-plate module in laminar flow conditions. Areas of investigation included the use of channels across the membrane surface, decreasing the available membrane surface area, and widening the inlet and outlet channels. The work showed that widening the channels increased the flux. Increased flux was also obtained by the use of channels on the permeate side, though not on the feed side.

Characterization of energy flux partitioning in urban environments: links with surface seasonal properties

A better understanding of links between the properties of the urban environment and the exchange to the atmosphere is central to a wide range of applications. The numerous measurements of surface energy balance data in urban areas enable intercomparison of observed fluxes from distinct environments. This study analyzes a large database in two new ways. First, instead of normalizing fluxes using net all-wave radiation only the incoming radiative fluxes are used, to remove the surface attributes from the denominator. Second, because data are now available year-round, indices are developed to characterize the fraction of the surface (built; vegetation) actively engaged in energy exchanges. These account for shading patterns within city streets and seasonal changes in vegetation phenology; their impact on the partitioning of the incoming radiation is analyzed. Data from 19 sites in North America, Europe, Africa, and Asia (including 6-yr-long observation campaigns) are used to derive generalized surface–flux relations. The midday-period outgoing radiative fraction decreases with an increasing total active surface index, the stored energy fraction increases with an active built index, and the latent heat fraction increases with an active vegetated index. Parameterizations of these energy exchange ratios as a function of the surface indices [i.e., the Flux Ratio–Active Index Surface Exchange (FRAISE) scheme] are developed. These are used to define four urban zones that characterize energy partitioning on the basis of their active surface indices. An independent evaluation of FRAISE, using three additional sites from the Basel Urban Boundary Layer Experiment (BUBBLE), yields accurate predictions of the midday flux partitioning at each location.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Energy exchange in a dense urban environment Part II: impact of spatial heterogeneity of the surface

The centre of cities, characterised by spatial and temporal complexity, are challenging environments for micrometeorological research. This paper considers the impact of sensor location and heterogeneity of the urban surface on flux observations in the dense city centre of London, UK. Data gathered at two sites in close vicinity, but with different measurement heights, were analysed to investigate the influence of source area characteristics on long-term radiation and turbulent heat fluxes. Combining consideration of diffuse radiation and effects of specular reflections, the non-Lambertian urban surface is found to impact the measurements of surface albedo. Comparisons of observations from the two sites reveal that turbulent heat fluxes are similar under some flow conditions. However, they mostly observe processes at different scales due to their differing measurement heights, highlighting the critical impact of siting sensors in urban areas. A detailed source area analysis is presented to investigate the surface controls influencing the energy exchanges at the different scales

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

Thermospheric Control of the Auroral Source of O+Ions for the Magnetosphere

Linear theory, model ion-density profiles and MSIS neutral thermospheric predictions are used to investigate the stability of the auroral, topside ionosphere to oxygen cyclotron waves: variations of the critical height, above which the plasma is unstable, with field-aligned current, thermal ion density and exospheric temperature are considered. In addition, probabilities are assessed that interactions with neutral atomic gases prevent O+ ions from escaping into the magnetosphere after they have been transversely accelerated by these waves. The two studies are combined to give a rough estimate of the total O+ escape flux as a function of the field-aligned current density for an assumed rise in the perpendicular ion temperature. Charge exchange with neutral oxygen, not hydrogen, is shown to be the principle limitation to the escape of O+ ions, which occurs when the waves are driven unstable down to low altitudes. It is found that the largest observed field-aligned current densities can heat a maximum of about 5×1014 O+ ions m−2 to a threshold above which they are subsequently able to escape into the magnetosphere in the following 500s. Averaged over this period, this would constitute a flux of 1012 m−2 s−1 and in steady-state the peak outflow would then be limited to about 1013 m−2 s−1 by frictional drag on thermal O+ at lower altitudes. Maximum escape is at low plasma density unless the O+ scale height is very large. The outflow decreases with decreasing field-aligned current density and, to a lesser extent, with increasing exospheric temperature. Upward flowing ion events are evaluated as a source of O+ ions for the magnetosphere and as an explanation of the observed solar cycle variation of ring current O+ abundance.

Observational constraints on atmospheric and oceanic cross-equatorial heat transports: revisiting the precipitation asymmetry problem in climate models

Satellite based top-of-atmosphere (TOA) and surface radiation budget observations are combined with mass corrected vertically integrated atmospheric energy divergence and tendency from reanalysis to infer the regional distribution of the TOA, atmospheric and surface energy budget terms over the globe. Hemispheric contrasts in the energy budget terms are used to determine the radiative and combined sensible and latent heat contributions to the cross-equatorial heat transports in the atmosphere (AHT_EQ) and ocean (OHT_EQ). The contrast in net atmospheric radiation implies an AHT_EQ from the northern hemisphere (NH) to the southern hemisphere (SH) (0.75 PW), while the hemispheric difference in sensible and latent heat implies an AHT_EQ in the opposite direction (0.51 PW), resulting in a net NH to SH AHT_EQ (0.24 PW). At the surface, the hemispheric contrast in the radiative component (0.95 PW) dominates, implying a 0.44 PW SH to NH OHT_EQ. Coupled model intercomparison project phase 5 (CMIP5) models with excessive net downward surface radiation and surface-to-atmosphere sensible and latent heat transport in the SH relative to the NH exhibit anomalous northward AHT_EQ and overestimate SH tropical precipitation. The hemispheric bias in net surface radiative flux is due to too much longwave surface radiative cooling in the NH tropics in both clear and all-sky conditions and excessive shortwave surface radiation in the SH subtropics and extratropics due to an underestimation in reflection by clouds.