Heat and moisture budgets from airborne measurements and high-resolution model simulations

High-resolution simulations with a mesoscale model are performed to estimate heat and moisture budgets of a well-mixed boundary layer. The model budgets are validated against energy budgets obtained from airborne measurements over heterogeneous terrain in Western Germany. Time rate of change, vertical divergence, and horizontal advection for an atmospheric column of air are estimated. Results show that the time trend of specific humidity exhibits some deficiencies, while the potential temperature trend is matched accurately. Furthermore, the simulated turbulent surface fluxes of sensible and latent heat are comparable to the measured fluxes, leading to similar values of the vertical divergence. The analysis of different horizontal model resolutions exhibits improved surface fluxes with increased resolution, a fact attributed to a reduced aggregation effect. Scale-interaction effects could be identified: while time trends and advection are strongly influenced by mesoscale forcing, the turbulent surface fluxes are mainly controlled by microscale processes.

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

Balance model for equatorial long waves

Geophysical fluid models often support both fast and slow motions. As the dynamics are often dominated by the slow motions, it is desirable to filter out the fast motions by constructing balance models. An example is the quasi geostrophic (QG) model, which is used widely in meteorology and oceanography for theoretical studies, in addition to practical applications such as model initialization and data assimilation. Although the QG model works quite well in the mid-latitudes, its usefulness diminishes as one approaches the equator. Thus far, attempts to derive similar balance models for the tropics have not been entirely successful as the models generally filter out Kelvin waves, which contribute significantly to tropical low-frequency variability. There is much theoretical interest in the dynamics of planetary-scale Kelvin waves, especially for atmospheric and oceanic data assimilation where observations are generally only of the mass field and thus do not constrain the wind field without some kind of diagnostic balance relation. As a result, estimates of Kelvin wave amplitudes can be poor. Our goal is to find a balance model that includes Kelvin waves for planetary-scale motions. Using asymptotic methods, we derive a balance model for the weakly nonlinear equatorial shallow-water equations. Specifically we adopt the ‘slaving’ method proposed by Warn et al. (Q. J. R. Meteorol. Soc., vol. 121, 1995, pp. 723–739), which avoids secular terms in the expansion and thus can in principle be carried out to any order. Different from previous approaches, our expansion is based on a long-wave scaling and the slow dynamics is described using the height field instead of potential vorticity. The leading-order model is equivalent to the truncated long-wave model considered previously (e.g. Heckley & Gill, Q. J. R. Meteorol. Soc., vol. 110, 1984, pp. 203–217), which retains Kelvin waves in addition to equatorial Rossby waves. Our method allows for the derivation of higher-order models which significantly improve the representation of Rossby waves in the isotropic limit. In addition, the ‘slaving’ method is applicable even when the weakly nonlinear assumption is relaxed, and the resulting nonlinear model encompasses the weakly nonlinear model. We also demonstrate that the method can be applied to more realistic stratified models, such as the Boussinesq model.

Contribution of land use changes to near-surface air temperatures during recent summer extreme heat events in the Phoenix Metropolitan Area

The impact of 1973–2005 land use–land cover (LULC) changes on near-surface air temperatures during four recent summer extreme heat events (EHEs) are investigated for the arid Phoenix, Arizona, metropolitan area using the Weather Research and Forecasting Model (WRF) in conjunction with the Noah Urban Canopy Model. WRF simulations were carried out for each EHE using LULC for the years 1973, 1985, 1998, and 2005. Comparison of measured near-surface air temperatures and wind speeds for 18 surface stations in the region show a good agreement between observed and simulated data for all simulation periods. The results indicate consistent significant contributions of urban development and accompanying LULC changes to extreme temperatures for the four EHEs. Simulations suggest new urban developments caused an intensification and expansion of the area experiencing extreme temperatures but mainly influenced nighttime temperatures with an increase of up to 10 K. Nighttime temperatures in the existing urban core showed changes of up to 2 K with the ongoing LULC changes. Daytime temperatures were not significantly affected where urban development replaced desert land (increase by 1 K); however, maximum temperatures increased by 2–4 K when irrigated agricultural land was converted to suburban development. According to the model simulations, urban landscaping irrigation contributed to cooling by 0.5–1 K in maximum daytime as well as minimum nighttime 2-m air temperatures in most parts of the urban region. Furthermore, urban development led to a reduction of the already relatively weak nighttime winds and therefore a reduction in advection of cooler air into the city.

Multi-scale sensible heat fluxes in the urban environment from large aperture scintillometry and eddy covariance

Sensible heat fluxes (QH) are determined using scintillometry and eddy covariance over a suburban area. Two large aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 km and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2-4 km2, whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5-10 km2. These large-scale heat fluxes are compared with local-scale eddy covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly QH is related to the meteorological conditions. At shorter time scales the response of QH to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation (Q*) coincide with greater differences between heat fluxes. For clear days QH lags Q*, thus the ratio of QH to Q* increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.

A critical revision of the estimation of the latent heat flux from two-wavelength scintillometry

Simultaneous scintillometer measurements at multiple wavelengths (pairing visible or infrared with millimetre or radio waves) have the potential to provide estimates of path-averaged surface fluxes of sensible and latent heat. Traditionally, the equations to deduce fluxes from measurements of the refractive index structure parameter at the two wavelengths have been formulated in terms of absolute humidity. Here, it is shown that formulation in terms of specific humidity has several advantages. Specific humidity satisfies the requirement for a conserved variable in similarity theory and inherently accounts for density effects misapportioned through the use of absolute humidity. The validity and interpretation of both formulations are assessed and the analogy with open-path infrared gas analyser density corrections is discussed. Original derivations using absolute humidity to represent the influence of water vapour are shown to misrepresent the latent heat flux. The errors in the flux, which depend on the Bowen ratio (larger for drier conditions), may be of the order of 10%. The sensible heat flux is shown to remain unchanged. It is also verified that use of a single scintillometer at optical wavelengths is essentially unaffected by these new formulations. Where it may not be possible to reprocess two-wavelength results, a density correction to the latent heat flux is proposed for scintillometry, which can be applied retrospectively to reduce the error.

Freshwater and heat transports from global ocean synthesis

[1] An eddy-permitting ¼° global ocean reanalysis based on the Operational Met Office FOAM data assimilation system has been run for 1989–2010 forced by ERA-Interim meteorology. Freshwater and heat transports are compared with published estimates globally and in each basin, with special focus on the Atlantic. The meridional transports agree with observations within errors at most locations, but where eddies are active the transports by the mean flow are nearly always in better agreement than the total transports. Eddy transports are down gradient and are enhanced relative to a free run. They may oppose or reinforce mean transports and provide 40–50% of the total transport near midlatitude fronts, where eddies with time scales <1 month provide up to 15%. Basin-scale freshwater convergences are calculated with the Arctic/Atlantic, Indian, and Pacific oceans north of 32°S, all implying net evaporation of 0.33 ± 0.04 Sv, 0.65 ± 0.07 Sv, and 0.09 ± 0.04 Sv, respectively, within the uncertainty of observations in the Atlantic and Pacific. The Indian is more evaporative and the Southern Ocean has more precipitation (1.07 Sv). Air-sea fluxes are modified by assimilation influencing turbulent heat fluxes and evaporation. Generally, surface and assimilation fluxes together match the meridional transports, indicating that the reanalysis is close to a steady state. Atlantic overturning and gyre transports are assessed with overturning freshwater transports southward at all latitudes. At 26°N eddy transports are negligible, overturning transport is 0.67 ± 0.19 Sv southward and gyre transport is 0.44 ± 0.17 Sv northward, with divergence between 26°N and the Bering Strait of 0.13 ± 0.23 Sv over 2004–2010.

Meteorology, air quality, and health in London: the ClearfLo project

The ClearfLo project provides integrated measurements of the meteorology, composition and particulate loading of London's urban atmosphere to improve predictive capability for air quality. Air quality and heat are strong health drivers and their accurate assessment and forecast are important in densely populated urban areas. However, the sources and processes leading to high concentrations of main pollutants such as ozone, nitrogen dioxide, and fine and coarse particulate matter in complex urban areas are not fully understood, limiting our ability to forecast air quality accurately. This paper introduces the ClearfLo project's interdisciplinary approach to investigate the processes leading to poor air quality and elevated temperatures. Within ClearfLo (www.clearflo.ac.uk), a large multi-institutional project funded by the UK Natural Environment Research Council (NERC), integrated measurements of meteorology, gaseous and particulate composition/loading within London's atmosphere were undertaken to understand the processes underlying poor air quality. Long-term measurement infrastructure installed at multiple levels (street and elevated), and at urban background, kerbside and rural locations were complemented with high-resolution numerical atmospheric simulations . Combining these (measurement/modeling) enhances understanding of seasonal variations in meteorology and composition together with the controlling processes. Two intensive observation periods (winter 2012 and summer Olympics 2012) focus upon the vertical structure and evolution of the urban boundary layer, chemical controls on nitrogen dioxide and ozone production, in particular the role of volatile organic compounds, and processes controlling the evolution, size, distribution and composition of particulate matter. The paper shows that mixing heights are deeper over London than in the rural surroundings and the seasonality of the urban boundary layer evolution controls when concentrations peak. The composition also reflects the seasonality of sources such as domestic burning and biogenic emissions.

An optimal inverse method using Doppler lidar measurements to estimate the surface sensible heat flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average −63%) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.

Using Distributed Temperature Sensing to monitor field scale dynamics of ground surface temperature and related substrate heat flux

We present one of the first studies of the use of Distributed Temperature Sensing (DTS) along fibre-optic cables to purposely monitor spatial and temporal variations in ground surface temperature (GST) and soil temperature, and provide an estimate of the heat flux at the base of the canopy layer and in the soil. Our field site was at a groundwater-fed wet meadow in the Netherlands covered by a canopy layer (between 0-0.5 m thickness) consisting of grass and sedges. At this site, we ran a single cable across the surface in parallel 40 m sections spaced by 2 m, to create a 40×40 m monitoring field for GST. We also buried a short length (≈10 m) of cable to depth of 0.1±0.02 m to measure soil temperature. We monitored the temperature along the entire cable continuously over a two-day period and captured the diurnal course of GST, and how it was affected by rainfall and canopy structure. The diurnal GST range, as observed by the DTS system, varied between 20.94 and 35.08◦C; precipitation events acted to suppress the range of GST. The spatial distribution of GST correlated with canopy vegetation height during both day and night. Using estimates of thermal inertia, combined with a harmonic analysis of GST and soil temperature, substrate and soil-heat fluxes were determined. Our observations demonstrate how the use of DTS shows great promise in better characterising area-average substrate/soil heat flux, their spatiotemporal variability, and how this variability is affected by canopy structure. The DTS system is able to provide a much richer data set than could be obtained from point temperature sensors. Furthermore, substrate heat fluxes derived from GST measurements may be able to provide improved closure of the land surface energy balance in micrometeorological field studies. This will enhance our understanding of how hydrometeorological processes interact with near-surface heat fluxes.

Alternative body sites for heat stress measurement in milking cows under tropical conditions and their relationship to the thermal discomfort of the animals

This study was conducted to determine the relationship among temperatures measured at different anatomical sites of the animal body and their daily pattern as indicative of the thermal stress in lactating dairy cows under tropical conditions. Environmental dry bulb (DBT) and black globe (BGT) temperatures and relative humidity (RH) were recorded. Rectal temperature (RT), respiratory frequency (RF), body surface (BST), internal base of tail (TT), vulva (VT) and auricular temperatures (AT) were collected, from 37 Black and White Holstein cows at 0700, 1300 and 1800 hours. RT showed a moderately and positive correlations with all body temperatures, ranging from 0.59 with TT to 0.64 with BST. Correlations among AT, VT and TT with RF were very similar (from 0.63 to 0.64) and were greater than those observed for RF with RT (0.55) or with BST (0.54). RF and RT were positively correlated to TT (0.63 and 0.59, respectively), AT (r = 0.63 for both) and VT (r = 0.64 and 0.63, respectively). Positive and very high correlations were observed among AT, VT and TT (from 0.94 to 0.97) indicating good association of temperatures measured in these anatomical sites. Correlations of BST with AT and VT were positive and very similar (0.71 and 0.72, respectively) and lower with TT (0.66). The AT, TT, VT and BST presented similar patterns and follow the variations of DBT through the day. Temperatures measured at different anatomical sites of the animal body have the potential to be used as indicative of the thermal stress in lactating dairy cows.

An assessment of the latent and sensible heat flux on the simulated regional climate over Southwestern South Atlantic Ocean

A Regional Climate Model (RegCM3) 10-year (1990-1999) simulation over southwestern South Atlantic Ocean (SAO) is evaluated to assess the mean climatology and the simulation errors of turbulent fluxes over the sea. Moreover, the relationship between these fluxes and the rainfall over some cyclogenetic areas is also analyzed. The RegCM3 results are validated using some reanalyses datasets (ERA40, R2, GPCP and WHOI). The summer and winter spatial patterns of latent and sensible heat fluxes simulated by the RegCM3 are in agreement with the reanalyses (WHOI, R2 and ERA40). They show large latent heat fluxes exchange in the subtropical SAO and at higher latitudes in the warm waters of Brazil Current. In particular, the magnitude of RegCM3 latent heat fluxes is similar to the WHOI, which is probably related to two factors: (a) small specific humidity bias, and (b) the RegCM3 flux algorithm. In contrast, the RegCM3 presents large overestimation of sensible heat flux, though it simulates well their spatial pattern. This simulation error is associated with the RegCM3 underestimation of the 2-m air temperature. In southwestern SAO, in three known cyclogenetic areas, the reanalyses and the RegCM3 show the existence of different physical mechanisms that control the annual cycles of latent/sensible heating and rainfall. It is shown that over the eastern coast of Uruguay (35A degrees-43A degrees S) and the southeastern coast of Argentina (44A degrees-52A degrees S) the sea-air moisture and heat exchange play an important role to control the annual cycle of precipitation. This does not happen on the south/southeastern coast of Brazil.

Origin of Convectively Coupled Kelvin Waves over South America

Convectively coupled Kelvin waves over the South American continent are examined through the use of temporal and spatial filtering of reanalysis, satellite, and gridded rainfall data. They are most prominent from November to April, the season analyzed herein. The following two types of events are isolated: those that result from preexisting Kelvin waves over the eastern Pacific Ocean propagating into the continent, and those that apparently originate over Amazonia, forced by disturbances propagating equatorward from central and southern South America. The events with precursors in the Pacific are mainly upper-level disturbances, with almost no signal at the surface. Those events with precursors over South America, on the other hand, originate as upper-level synoptic wave trains that pass over the continent and resemble the ""cold surges`` documented by Garreaud and Wallace. As the wave train propagates over the Andes, it induces a southerly low-level wind that advects cold air to the north. Precipitation associated with a cold front reaches the equator a few days later and subsequently propagates eastward with the characteristics of a Kelvin wave. The structures of those waves originating over the Pacific are quite similar to those originating over South America as they propagate to eastern South America and into the Atlantic. South America Kelvin waves that originate over neither the Pacific nor the midlatitudes of South America can also be identified. In a composite sense, these form over the eastern slope of the Andes Mountains, close to the equator. There are also cases of cold surges that reach the equator yet do not form Kelvin waves. The interannual variability of the Pacific-originating events is related to sea surface temperatures in the central-eastern Pacific Ocean. When equatorial oceanic conditions are warm, there tends to be an increase in the number of disturbances that reach South America from the Pacific.

Anthropogenic heat in the city of So Paulo, Brazil

The main goal of this work is to describe the anthropogenic energy flux (Q (F)) in the city of So Paulo, Brazil. The hourly, monthly, and annual values of the anthropogenic energy flux are estimated using the inventory method, and the contributions of vehicular, stationary, and human metabolism sources from 2004 to 2007 are considered. The vehicular and stationary sources are evaluated using the primary consumption of energy based on fossil fuel, bio fuel, and electricity usage by the population. The diurnal evolution of the anthropogenic energy flux shows three relative maxima, with the largest maxima occurring early in the morning (similar to 19.9 Wm(-2)) and in the late afternoon (similar to 20.3 Wm(-2)). The relative maximum that occurs around noontime (similar to 19.6 Wm(-2)) reflects the diurnal pattern of vehicle traffic that seems to be specific to So Paulo. With respect to diurnal evolution, the energy flux released by vehicular sources (Q (FV)) contributes approximately 50% of the total anthropogenic energy flux. Stationary sources (Q (FS)) and human metabolism (Q (FM)) represent about 41% and 9% of the anthropogenic energy flux, respectively. For 2007, the monthly values of Q (FV), Q (FS), Q (FM), and Q (F) are, respectively, 16.8 +/- 0.25, 14.3 +/- 0.16, 3.5 +/- 0.03, and 34.6 +/- 0.41 MJ m(-2) month(-1). The seasonal evolution monthly values of Q (FV), Q (FS), Q (FM), and Q (F) show a relative minimum during the summer and winter vacations and a systematic and progressive increase associated with the seasonal evolution of the economic activity in So Paulo. The annual evolution of Q (F) indicates that the city of So Paulo released 355.2 MJ m(-2) year(-1) in 2004 and 415.5 MJ m(-2) year(-1) in 2007 in association with an annual rate of increase of 19.6 MJ m(-2) year(-1) (from 2004 to 2006) and 30.5 MJ m(-2) year(-1) (from 2006 to 2007). The anthropogenic energy flux corresponds to about 9% of the net radiation at the surface in the summer and 15% in the winter. The amplitude of seasonal variation of the maximum hourly value of the diurnal variation increases exponentially with latitude.

Wolf-Rayet optically thick winds with Alfven waves

The Wolf-Rayet (WR) stars are hot luminous objects which are suffering an extreme mass loss via a continuous stellar wind. The high values of mass loss rates and high terminal velocities of the WR stellar winds constitute a challenge to the theories of radiation driven winds. Several authors incorporated magnetic forces to the line driven mechanism in order to explain these characteristics of the wind. Observations indicate that the WR stellar winds may reach, at the photosphere, velocities of the order of the terminal values, which means that an important part of the wind acceleration occurs at the optically thick region. The aim of this study is to analyze a model in which the wind in a WR star begins to be accelerated in the optically thick part of the wind. We used as initial conditions stellar parameters taken from the literature and solved the energy, mass and momentum equations. We demonstrate that the acceleration only by radiative forces is prevented by the general behavior of the opacities. Combining radiative forces plus a flux of Alfven waves, we found in the simulations a fast drop in the wind density profile which strongly reduces the extension of the optically thick region and the wind becomes optically thin too close its base. The understanding how the WR wind initiate is still an open issue. (C) 2010 COSPAR. Published by Elsevier Ltd. All rights reserved.