Effects of the large-scale atmospheric circulation on the onset and strength of urban heat islands: a case study

Air temperature was monitored at 13 sites across the urban perimeter of a Brazilian midsize city in winter 2011. In this study, we show that the urban heat island (UHI) develops only at night and under certain weather conditions, and its intensity depends not only on the site's land cover but also on the meteorological setting. The urban heat island intensity was largest (6.6 °C) under lingering high-pressure conditions, milder (3.0 °C) under cold anticyclones and almost vanished (1.0 °C) during the passage of cold fronts. The cooling rates were calculated to monitor the growth and decay of the UHI over each specific synoptic setting. Over four contiguous days under the effect of a lingering high-pressure event, we observed that the onset of cooling was always at about 2 h before sunset. The reference site attained mean cooling rate of ?2.6 °C h^?1 at sunset, whilst the maximum urban rate was ?1.2 °C h^?1. Under a 3-day cold anticyclone episode, cooling also started about 2 h before sunset, and the difference between maximum rural (?2.0 °C h^?1) and urban (?1.0 °C h^?1) cooling rates diminished. Under cold-front conditions, the cooling rate was homogeneous for all sites and swang about zero throughout the day. The air temperature has a memory effect under lingering high-pressure conditions which intensified the UHI, in addition to the larger heat storage in the urban area. Cold anticyclone conditions promoted the development of the UHI; however, the cold air pool and relatively light winds smoothed out its intensity. Under the influence of cold fronts, the urban fabric had little effect on the city's air temperature field, and the UHI was imperceptible.     

Permeability-porosity Relationships in Rocks Subjected to Various Evolution Processes

— It is well known that there is no “universal” permeability-porosity relationship valid in all porous media. However, the evolution of permeability and porosity in rocks can be constrained provided that the processes changing the pore space are known. In this paper, we review observations of the relationship between permeability and porosity during rock evolution and interpret them in terms of creation/destruction of effectively and non-effectively conducting pore space. We focus on laboratory processes, namely, plastic compaction of aggregates, elastic-brittle deformation of granular rocks, dilatant and thermal microcracking of dense rocks, chemically driven processes, as a way to approach naturally occurring geological processes. In particular, the chemically driven processes and their corresponding evolution permeability-porosity relationships are discussed in relation to sedimentary rocks diagenesis.     

Development and sensitivity analysis of a model for assessing stratification and safety of Lake Nyos during artificial degassing

To prevent the recurrence of a disastrous eruption of carbon dioxide (CO₂) from Lake Nyos, a degassing plan has been set up for the lake. Since there are concerns that the degassing of the lake may reduce the stability of the density stratification, there is an urgent need for a simulation tool to predict the evolution of the lake stratification in different scenarios. This paper describes the development of a numerical model to predict the CO₂ and dissolved solids concentrations, and the temperature structure as well as the stability of the water column of Lake Nyos. The model is tested with profiles of CO₂ concentrations and temperature taken in the years 1986 to 1996. It reproduces well the general mixing patterns observed in the lake. However, the intensity of the mixing tends to be overestimated in the epilimnion and underestimated in the monimolimnion. The overestimation of the mixing depth in the epilimnion is caused either by the parameterization of the k-epsilon model, or by the uncertainty in the calculation of the surface heat fluxes. The simulated mixing depth is highly sensitive to the surface heat fluxes, and errors in the mixing depth propagate from one year to the following. A precise simulation of the mixolimnion deepening therefore requires high accuracy in the meteorological forcing and the parameterization of the heat fluxes. Neither the meteorological data nor the formulae for the calculation of the heat fluxes are available with the necessary precision. Consequently, it will be indispensable to consider different forcing scenarios in the safety analysis in order to obtain robust boundary conditions for safe degassing. The input of temperature and CO₂ to the lake bottom can be adequately simulated for the years 1986 to 1996 with a constant sublacustrine source of 18 l s^–1 with a CO₂ concentration of 0.395 mol l^–1 and a temperature of 26 °C. The results of this study indicate that the model needs to be calibrated with more detailed field data before using it for its final purpose: the prediction of the stability and the safety of Lake Nyos during the degassing process.Responsible Editor: Hans Burchard     

A new calcium-aluminate from a refractory inclusion in the Leoville carbonaceous chondrite

The first meteoritic occurrence of CaAl₄O₇ is described from a Ca-Al-rich inclusion (CAI) in the Leoville carbonaceous chondrite. This CAI consists mainly of gehlenitic melilite, spinel, perovskite, and hibonite. CaAl₄O₇ is a minor component and occurs within melilite preferentially in portions rich in perovskite.The CAI is enveloped by a succession of three rims (from inside out): (a) hibonite+melilite+spinel+perovskite, (b) diopside, and (c) olivine.On the basis of mineral associations found and from the presence of moderately volatile elements (Fe and Cr) we conclude that the CaAl₄O₇-bearing CAI from Leoville is of residual nature. CaAl₄O₇ is apparently stable in the very Mg- and Si-poor environment of this CAI and is probably of igneous origin.The rims are interpreted as products of partial evaporation (rim (a)) and associated re-condensation (rims (b) and (c)).     

Turbulent heat fluxes by profile and inertial dissipation methods: analysis of the atmospheric surface layer from shipboard measurements during the SOFIA/ASTEX and SEMAPHORE experiments

Heat flux estimates obtained using the inertial dissipation method, and the profile method applied to radiosonde soundings, are assessed with emphasis on the parameterization of the roughness lengths for temperature and specific humidity. Results from the inertial dissipation method show a decrease of the temperature and humidity roughness lengths for increasing neutral wind speed, in agreement with previous studies. The sensible heat flux estimates were obtained using the temperature estimated from the speed of sound determined by a sonic anemometer. This method seems very attractive for estimating heat fluxes over the ocean. However allowance must be made in the inertial dissipation method for non-neutral stratification. The SOFIA/ASTEX and SEMAPHORE results show that, in unstable stratification, a term due to the transport terms in the turbulent kinetic energy budget, has to be included in order to determine the friction velocity with better accuracy. Using the profile method with radiosonde data, the roughness length values showed large scatter. A reliable estimate of the temperature roughness length could not be obtained. The humidity roughness length values were compatible with those found using the inertial dissipation method.     

Benthic primary producers--a neglected environmental problem in Mediterranean maricultures?

Marine fish farming is increasing rapidly in the Mediterranean and in contrast to the Atlantic the coastal zone in the Mediterranean is characterized by clear waters with high transparency. This allows benthic primary producers such as the slow-growing seagrass Posidonia oceanica to grow at large depths at locations suitable for fish farming and generating a conflict between the conservation of these meadows and the growth of aquaculture operations in the Mediterranean. In this paper we review the current knowledge on environmental interactions between fish farming and benthic primary producers with particular focus on P. oceanica, as this seagrass is a key component along Mediterranean coasts. The recovery times of P. oceanica are very long, in the order of centuries, and losses of this species are thus considered to be irreversible at managerial time scales.     

Wind-forced circulation model and water exchanges through the channel in the Bay of Toulon

A hydrodynamic model of the Bay of Toulon has been developed for use as a post-accident radionuclide dispersion simulation tool. Located in a Mediterranean urban area, the Bay of Toulon is separated into two basins by a 1.4-km long seawall. The Little Bay is semi-enclosed and connected to the Large Bay by a fairway channel. This channel is the site of significant water mass exchange as a result of both wind-driven currents and bathymetry. It is therefore a focal point for marine contamination. As part of the model calibration and validation process, the first step consisted of studying the water mass exchange between the two basins. An Acoustic Doppler Current Profiler was moored in the channel for 1 year. The present study analyses in situ data to determine the current intensity and direction, and also to better understand the vertical current profile, which is highly correlated with meteorological forcing. Comparisons of model-generated and measured data are presented, and various atmospheric forcing datasets are used to enhance computed results. It appears that accurate meteorological forcing data is needed to enhance the accuracy of the hydrodynamic model. This channel is an important location for water mass renewal in the Bay of Toulon, and model results are used to quantify these exchanges. The mean calculated annual water exchange time is approximately 3.4 days. However, this duration is strongly wind dependent and shortens during windy winter months. It ranges from 1.5 days during strong wind periods to 7.5 days during calm weather. Residence time values calculated through tracer dispersion modelling after release at the back of the Little Bay are found to be comparable to the mean exchange time values, especially for windy conditions.     

An Experimental High-Resolution Forecast System During the Vancouver 2010 Winter Olympic and Paralympic Games

Environment Canada ran an experimental numerical weather prediction (NWP) system during the Vancouver 2010 Winter Olympic and Paralympic Games, consisting of nested high-resolution (down to 1-km horizontal grid-spacing) configurations of the GEM–LAM model, with improved geophysical fields, cloud microphysics and radiative transfer schemes, and several new diagnostic products such as density of falling snow, visibility, and peak wind gust strength. The performance of this experimental NWP system has been evaluated in these winter conditions over complex terrain using the enhanced mesoscale observing network in place during the Olympics. As compared to the forecasts from the operational regional 15-km GEM model, objective verification generally indicated significant added value of the higher-resolution models for near-surface meteorological variables (wind speed, air temperature, and dewpoint temperature) with the 1-km model providing the best forecast accuracy. Appreciable errors were noted in all models for the forecasts of wind direction and humidity near the surface. Subjective assessment of several cases also indicated that the experimental Olympic system was skillful at forecasting meteorological phenomena at high-resolution, both spatially and temporally, and provided enhanced guidance to the Olympic forecasters in terms of better timing of precipitation phase change, squall line passage, wind flow channeling, and visibility reduction due to fog and snow.     

Modeling the Phase Transition Associated with Melting Snow in a 1D Kinematic Framework: Sensitivity to the Microphysics

A simple 1D kinematic cloud model coupled to a two-moment bulk microphysics scheme is used to perform quasi-idealized simulations of snow, with a prescribed upper boundary snow field based on observed radar reflectivity and temperature, falling into a low-level melting layer. The model realistically simulates the formation of a nearly isothermal layer below the melting level, the surface precipitation rate, and the phase transition from liquid to solid, consistent with observations for this case. A series of test runs is performed to examine the sensitivity of modeling the timing and duration of the phase transition period to details of specific parameterization aspects related to snow in the microphysics scheme. The sensitivity tests include varying the number of prognostic moments, the mass–diameter relation, the fall velocity–diameter relation, the treatment of aggregation, and the lower limit for the slope of the size distribution. It is shown that the simulated transition period, for such a case with the initial melting level being close to the surface, can be quite sensitive to model parameters specified within realistic ranges and/or ranges within our physical understanding.     

Eurocode 8-compatible synthetic time-series as input to dynamic analysis

Nonlinear dynamic analysis of existing or planned structures often requires the use of accelerograms that match a target design spectrum. Here, our main concern is to generate a set of motions with a good level of fit to the Eurocode 8 design spectra for France. Synthetic time series are generated by means of a non-stationary stochastic method. To calibrate the input parameters in the stochastic approach, we select a reference set of accelerograms for a Eurocode 8 type B site category from the PEER Ground-Motion Database, which are then adjusted to the target spectrum through wavelet addition. Then, we compute nonlinear seismic responses of a soil column, including pore pressure effects, and brittle and ductile structures to the stochastic time-series, the natural accelerograms and time-series generated using stationary stochastic approaches. The results of these calculations reveal considerable variability in response despite the similarities in terms of spectral acceleration.