Identification of micro-scale anthropogenic CO2, heat and moisture sources – processing eddy covariance fluxes for a dense urban environment

Anthropogenic emissions of heat and exhaust gases play an important role in the atmospheric boundary layer, altering air quality, greenhouse gas concentrations and the transport of heat and moisture at various scales. This is particularly evident in urban areas where emission sources are integrated in the highly heterogeneous urban canopy layer and directly linked to human activities which exhibit significant temporal variability. It is common practice to use eddy covariance observations to estimate turbulent surface fluxes of latent heat, sensible heat and carbon dioxide, which can be attributed to a local scale source area. This study provides a method to assess the influence of micro-scale anthropogenic emissions on heat, moisture and carbon dioxide exchange in a highly urbanized environment for two sites in central London, UK. A new algorithm for the Identification of Micro-scale Anthropogenic Sources (IMAS) is presented, with two aims. Firstly, IMAS filters out the influence of micro-scale emissions and allows for the analysis of the turbulent fluxes representative of the local scale source area. Secondly, it is used to give a first order estimate of anthropogenic heat flux and carbon dioxide flux representative of the building scale. The algorithm is evaluated using directional and temporal analysis. The algorithm is then used at a second site which was not incorporated in its development. The spatial and temporal local scale patterns, as well as micro-scale fluxes, appear physically reasonable and can be incorporated in the analysis of long-term eddy covariance measurements at the sites in central London. In addition to the new IMAS-technique, further steps in quality control and quality assurance used for the flux processing are presented. The methods and results have implications for urban flux measurements in dense urbanised settings with significant sources of heat and greenhouse gases.

Controls on sensible heat and latent energy fluxes from a short-hydroperiod Florida Everglades marsh

Little is known of energy balance in low latitude wetlands where there is a year-round growing season and a climate best defined by wet and dry seasons. The Florida Everglades is a highly managed and extensive subtropical wetland that exerts a substantial influence on the hydrology and climate of the south Florida region. However, the effects of seasonality and active water management on energy balance in the Everglades ecosystem are poorly understood. An eddy covariance and micrometeorological tower was established in a short-hydroperiod Everglades marsh to examine the dominant environmental controls on sensible heat (H) and latent energy (LE) fluxes, as well as the effects of seasonality on these parameters. Seasonality differentially affected H and LE fluxes in this marsh, such that H was principally dominant in the dry season and LE was strongly dominant in the wet season. The Bowen ratio was high for much of the dry season (1.5–2.4), but relatively low (H and LE fluxes across nearly all seasons and years (). However, the 2009 dry season LE data were not consistent with this relationship () because of low seasonal variation in LE following a prolonged end to the previous wet season. In addition to net radiation, H and LE fluxes were significantly related to soil volumetric water content (VWC), water depth, air temperature, and occasionally vapor pressure deficit. Given that VWC and water depth were determined in part by water management decisions, it is clear that human actions have the ability to influence the mode of energy dissipation from this ecosystem. Impending modifications to water management under the Comprehensive Everglades Restoration Plan may shift the dominant turbulent flux from this ecosystem further toward LE, and this change will likely affect local hydrology and climate.

Seasonal changes of heat and moisture fluxes from the Caspian Sea surface and its effect on synoptic systems with heavy rainfalls in southern shores of the sea during 2005-2010

Caspian Sea with its unique characteristics is a significant source to supply required heat and moisture for passing weather systems over the north of Iran. Investigation of heat and moisture fluxes in the region and their effects on these systems that could lead to floods and major financial and human losses is essential in weather forecasting. Nowadays by improvement of numerical weather and climate prediction models and the increasing need to more accurate forecasting of heavy rainfall, the evaluation and verification of these models has been become much more important. In this study we have used the WRF model as a research-practical one with many valuable characteristics and flexibilities. In this research, the effects of heat and moisture fluxes of Caspian Sea on the synoptic and dynamical structure of 20 selective systems associated with heavy rainfall in the southern shores of Caspian Sea are investigated. These systems are selected based on the rainfall data gathered by three local stations named: Rasht, Babolsar and Gorgan in different seasons during a five-year period (2005-2010) with maximum amount of rainfall through the 24 hours of a day. In addition to synoptic analyses of these systems, the WRF model with and without surface flues was run using the two nested grids with the horizontal resolutions of 12 and 36 km. The results show that there are good consistencies between the predicted distribution of rainfall field, time of beginning and end of rainfall by the model and the observations. But the model underestimates the amounts of rainfall and the maximum difference with the observation is about 69%. Also, no significant changes in the results are seen when the domain and the resolution of computations are changed. The other noticeable point is that the systems are severely weakened by removing heat and moisture fluxes and thereby the amounts of large scale rainfall are decreased up to 77% and the convective rainfalls tend to zero.

An analysis of MONTBLEX data on heat and momentum flux at Jodhpur

Parameterization of sensible heat and momentum fluxes as inferred from an analysis of tower observations archived during MONTBLEX-90 at Jodhpur is proposed, both in terms of standard exchange coefficients C-H and C-D respectively and also according to free convection scaling. Both coefficients increase rapidly at low winds (the latter more strongly) and with increasing instability. All the sensible heat flux data at Jodhpur (wind speed at 10m <(U)over bar (10)>, < 8ms(-1)) also obey free convection scaling, with the flux proportional to the '4/3' power of an appropriate temperature difference such as that between 1 and 30 m. Furthermore, for <(U)over bar (10)> < 4 ms(-1) the momentum flux displays a linear dependence on wind speed.

On the relationship between finger width, velocity, and fluxes in thermohaline convection

Double-diffusive finger convection occurs in many natural processes.The theories for double-diffusive phenomena that exist at present consider systems with linear stratification in temperature and salinity. The double-diffusive systems with step change in salinity and temperature are, however, not amenable to simple stability analysis. Hence factors that control the width of the finger, velocity, and fluxes in systems that have step change in temperature and salinity have not been understood so far. In this paper we provide new physical insight regarding factors that influence finger convection in two-layer double-diffusive system through two-dimensional numerical simulations. Simulations have been carried out for density stability ratios (R-rho) from 1.5 to 10. For each density stability ratio, the thermal Rayleigh number (Ra-T) has been systematically varied from 7x10(3) to 7x10(8). Results from these simulations show how finger width, velocity, and flux ratios in finger convection are interrelated and the influence of governing parameters such as density stability ratio and the thermal Rayleigh number. The width of the incipient fingers at the time of onset of instability has been shown to vary as Ra-T-1/3. Velocity in the finger varies as Ra(T)1/3/R-rho. Results from simulation agree with the scale analysis presented in the paper. Our results demonstrate that wide fingers have lower velocities and flux ratios compared to those in narrow fingers. This result contradicts present notions about the relation between finger width and flux ratio. A counterflow heat-exchanger analogy is used in understanding the dependence of flux ratio on finger width and velocity.

Canopy processes, fluxes and microclimate in a pine forest

Interaction between forests and the atmosphere occurs by radiative and turbulent transport. The fluxes of energy and mass between surface and the atmosphere directly influence the properties of the lower atmosphere and in longer time scales the global climate. Boreal forest ecosystems are central in the global climate system, and its responses to human activities, because they are significant sources and sinks of greenhouse gases and of aerosol particles. The aim of the present work was to improve our understanding on the existing interplay between biologically active canopy, microenvironment and turbulent flow and quantify. In specific, the aim was to quantify the contribution of different canopy layers to whole forest fluxes. For this purpose, long-term micrometeorological and ecological measurements made in a Scots pine (Pinus sylvestris) forest at SMEAR II research station in Southern Finland were used. The properties of turbulent flow are strongly modified by the interaction between the canopy elements: momentum is efficiently absorbed in the upper layers of the canopy, mean wind speed and turbulence intensities decrease rapidly towards the forest floor and power spectra is modulated by spectral short-cut . In the relative open forest, diabatic stability above the canopy explained much of the changes in velocity statistics within the canopy except in strongly stable stratification. Large eddies, ranging from tens to hundred meters in size, were responsible for the major fraction of turbulent transport between a forest and the atmosphere. Because of this, the eddy-covariance (EC) method proved to be successful for measuring energy and mass exchange inside a forest canopy with exception of strongly stable conditions. Vertical variations of within canopy microclimate, light attenuation in particular, affect strongly the assimilation and transpiration rates. According to model simulations, assimilation rate decreases with height more rapidly than stomatal conductance (gs) and transpiration and, consequently, the vertical source-sink distributions for carbon dioxide (CO2) and water vapor (H2O) diverge. Upscaling from a shoot scale to canopy scale was found to be sensitive to chosen stomatal control description. The upscaled canopy level CO2 fluxes can vary as much as 15 % and H2O fluxes 30 % even if the gs models are calibrated against same leaf-level dataset. A pine forest has distinct overstory and understory layers, which both contribute significantly to canopy scale fluxes. The forest floor vegetation and soil accounted between 18 and 25 % of evapotranspiration and between 10 and 20 % of sensible heat exchange. Forest floor was also an important deposition surface for aerosol particles; between 10 and 35 % of dry deposition of particles within size range 10 30 nm occurred there. Because of the northern latitudes, seasonal cycle of climatic factors strongly influence the surface fluxes. Besides the seasonal constraints, partitioning of available energy to sensible and latent heat depends, through stomatal control, on the physiological state of the vegetation. In spring, available energy is consumed mainly as sensible heat and latent heat flux peaked about two months later, in July August. On the other hand, annual evapotranspiration remains rather stable over range of environmental conditions and thus any increase of accumulated radiation affects primarily the sensible heat exchange. Finally, autumn temperature had strong effect on ecosystem respiration but its influence on photosynthetic CO2 uptake was restricted by low radiation levels. Therefore, the projected autumn warming in the coming decades will presumably reduce the positive effects of earlier spring recovery in terms of carbon uptake potential of boreal forests.

Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

Turbulent heat flux variation over the Monsoon-Trough region during MONTBLEX-90

Tower data collected during the Monsoon-Trough Boundary Layer Experiment (MONTBLEX-90) have been analysed to understand the observed structure of the surface layer over an arid region (Jodhpur) and a moist region (Kharagpur) during active and weak phases of the 1990 southwest monsoon. Turbulent heat and momentum fluxes are estimated by the eddy correlation method using sonic data. The turbulent momentum flux at both Jodhpur and Kharagpur was larger when the winds were stronger, reaching a maximum of the order of 0.5 N m(-2) on 5 and 6 August when a low pressure system was located over the region. The heat flux at Jodhpur is high during weak monsoon days, the maximum being 450 W m(-2), whereas during active days the flux never exceeds 200 W m(-2). At Kharagpur, the flux does not vary significantly between active and weak monsoon days, the maximum in either phase being 160 W m(-2) At Jodhpur, there is significant contrast in the near-surface air temperature, being higher during weak monsoon days as compared to active days. Cloud cover did not vary significantly in both the regions. The turbulent heat flux variation at both the sites appears to be correlated mainly with soil mixture, and less sensitive to cloud cover.

Transient ventilation dynamics induced by heat sources of unequal strength

We examine theoretically the transient displacement flow and density stratification that develops within a ventilated box after two localized floor-level heat sources of unequal strengths are activated. The heat input is represented by two non-interacting turbulent axisymmetric plumes of constant buoyancy fluxes B1 and B2 > B1. The box connects to an unbounded quiescent external environment of uniform density via openings at the top and base. A theoretical model is developed to predict the time evolution of the dimensionless depths λj and mean buoyancies δj of the 'intermediate' (j = 1) and 'top' (j = 2) layers leading to steady state. The flow behaviour is classified in terms of a stratification parameter S, a dimensionless measure of the relative forcing strengths of the two buoyant layers that drive the flow. We find that dδ1/dτ α 1/λ1 and dδ2/dτ α 1/λ2, where τ is a dimensionless time. When S 1, the intermediate layer is shallow (small λ1), whereas the top layer is relatively deep (large λ2) and, in this limit, δ1 and δ2 evolve on two characteristically different time scales. This produces a time lag and gives rise to a 'thermal overshoot', during which δ1 exceeds its steady value and attains a maximum during the transients; a flow feature we refer to, in the context of a ventilated room, as 'localized overheating'. For a given source strength ratio ψ = B1/B2, we show that thermal overshoots are realized for dimensionless opening areas A < Aoh and are strongly dependent on the time history of the flow. We establish the region of {A, ψ} space where rapid development of δ1 results in δ1 > δ2, giving rise to a bulk overturning of the buoyant layers. Finally, some implications of these results, specifically to the ventilation of a room, are discussed. © Cambridge University Press 2013.

N2O fluxes from the native and grazed semi-arid steppes and their driving factors in Inner Mongolia, China

This paper presents results of 2 years (from January 2005 to December 2006) of measurement of N2O fluxes from the native and grazed Leymus chinensis (LC) steppes in Inner Mongolia, China using the static opaque chamber method. The measurement was at a frequency of twice per month in the growing season and once per month in the non-growing season. In addition, the possible effect of water-heat factors on N2O fluxes was statistically analyzed. The results indicated that there were distinct seasonal patterns in N2O fluxes with large fluxes in spring, summer, and autumn but negative fluxes in winter. The annual net emission of N2O ranging from 0.24 to 0.30 kg N2O-N ha(-1) and from 0.06 to 0.26 kg N2O-N ha(-1) from the native and grazed LC steppe, respectively. Grazing activities suppressed N2O production. In the growing season, soil moisture was the primary driving factor of N2O fluxes. The high seasonal variation of N2O fluxes was regulated by the distribution of effective rainfall, rather than precipitation intensity. Air temperature or soil temperature at 0, 5, and 10 cm depth was the most restricting factor of N2O fluxes in the non-growing season.

Air–sea fluxes of oxygenated volatile organic compounds across the Atlantic Ocean

We present air-sea fluxes of oxygenated volatile organics compounds (OVOCs) quantified by eddy covariance (EC) during the Atlantic Meridional Transect cruise in 2012. Measurements of acetone, acetaldehyde, and methanol in air as well as in water were made in several different oceanic provinces and over a wide range of wind speeds (1-18 m s(-1)). The ocean appears to be a net sink for acetone in the higher latitudes of the North Atlantic but a source in the subtropics. In the South Atlantic, seawater acetone was near saturation relative to the atmosphere, resulting in essentially zero net flux. For acetaldehyde, the two-layer model predicts a small oceanic emission, which was not well resolved by the EC method. Chemical enhancement of air-sea acetaldehyde exchange due to aqueous hydration appears to be minor. The deposition velocity of methanol correlates linearly with the transfer velocity of sensible heat, confirming predominant airside control. We examine the relationships between the OVOC concentrations in air as well as in water, and quantify the gross emission and deposition fluxes of these gases.

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.

The cumulus-capped boundary layer. II: Interface fluxes

This paper considers the relationship between the mean temperature and humidity profiles and the fluxes of heat and moisture at cloud base and the base of the inversion in the cumulus-capped boundary layer. The relationships derived are based on an approximate form of the scalar-flux budget and the scaling properties of the turbulent kinetic energy (TKE) budget. The scalar-flux budget gives a relationship between the change in the virtual potential temperature across either the cloud base transition zone or the inversion and the flux at the base of the layer. The scaling properties of the TKE budget lead to a relationship between the heat and moisture fluxes and the mean subsaturation through the liquid-water flux. The 'jump relation' for the virtual potential temperature at cloud base shows the close connection between the cumulus mass flux in the cumulus-capped boundary layer and the entrainment velocity in the dry-convective boundary layer. Gravity waves are shown to be an important feature of the inversion.

Dependence of Indian monsoon rainfall on moisture fluxes across the Arabian Sea and the impact of coupled model sea surface temperature biases

The Arabian Sea is an important moisture source for Indian monsoon rainfall. The skill of climate models in simulating the monsoon and its variability varies widely, while Arabian Sea cold sea surface temperature (SST) biases are common in coupled models and may therefore influence the monsoon and its sensitivity to climate change. We examine the relationship between monsoon rainfall, moisture fluxes and Arabian Sea SST in observations and climate model simulations. Observational analysis shows strong monsoons depend on moisture fluxes across the Arabian Sea, however detecting consistent signals with contemporaneous summer SST anomalies is complicated in the observed system by air/sea coupling and large-scale induced variability such as the El Niño-Southern Oscillation feeding back onto the monsoon through development of the Somali Jet. Comparison of HadGEM3 coupled and atmosphere-only configurations suggests coupled model cold SST biases significantly reduce monsoon rainfall. Idealised atmosphere-only experiments show that the weakened monsoon can be mainly attributed to systematic Arabian Sea cold SST biases during summer and their impact on the monsoon-moisture relationship. The impact of large cold SST biases on atmospheric moisture content over the Arabian Sea, and also the subsequent reduced latent heat release over India, dominates over any enhancement in the land-sea temperature gradient and results in changes to the mean state. We hypothesize that a cold base state will result in underestimation of the impact of larger projected Arabian Sea SST changes in future climate, suggesting that Arabian Sea biases should be a clear target for model development.