The Parameterisation of Scalar Transfer over Rough Ice

We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z _s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z ₀). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z ₀ below a threshold value of approximately 10^?3 m, the model performs well and in accord with earlier studies. However, with growing hummock size, z ₀ increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z ₀ = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z _s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z ₀  >  10^?3 m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z _s /z ₀)  =  1.5  ? 0.2 ln (Re _*)  ? 0.11(ln (Re _*))².     

Volatile emissions and gas geochemistry of Hot Spring Basin,Yellowstone National Park,USA

We characterize and quantify volatile emissions at Hot Spring Basin (HSB), a large acid-sulfate region that lies just outside the northeastern edge of the 640 ka Yellowstone Caldera. Relative to other thermal areas in Yellowstone, HSB gases are rich in He and H₂, and mildly enriched in CH₄ and H₂S. Gas compositions are consistent with boiling directly off a deep geothermal liquid at depth as it migrates toward the surface. This fluid, and the gases evolved from it, carries geochemical signatures of magmatic volatiles and water–rock reactions with multiple crustal sources, including limestones or quartz-rich sediments with low K/U (or ^40?Ar/^4?He). Variations in gas chemistry across the region reflect reservoir heterogeneity and variable degrees of boiling. Gas-geothermometer temperatures approach 300 °C and suggest that the reservoir feeding HSB is one of the hottest at Yellowstone. Diffuse CO₂ flux in the western basin of HSB, as measured by accumulation-chamber methods, is similar in magnitude to other acid-sulfate areas of Yellowstone and is well correlated to shallow soil temperatures. The extrapolation of diffuse CO₂ fluxes across all the thermal/altered area suggests that 410 ± 140 t d^− 1 CO₂ are emitted at HSB (vent emissions not included). Diffuse fluxes of H₂S were measured in Yellowstone for the first time and likely exceed 2.4 t d^− 1 at HSB. Comparing estimates of the total estimated diffuse H₂S emission to the amount of sulfur as SO₄²⁻ in streams indicates ~ 50% of the original H₂S in the gas emission is lost into shallow groundwater, precipitated as native sulfur, or vented through fumaroles. We estimate the heat output of HSB as ~ 140–370 MW using CO₂ as a tracer for steam condensate, but not including the contribution from fumaroles and hydrothermal vents. Overall, the diffuse heat and volatile fluxes of HSB are as great as some active volcanoes, but they are a small fraction (1–3% for CO₂, 2–8% for heat) of that estimated for the entire Yellowstone system.     

Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke     

Mathematical Analysis of a Model of River Channel Formation

The study of overland flow of water over an erodible sediment leads to a coupled model describing the evolution of the topographic elevation and the depth of the overland water film. The spatially uniform solution of this model is unstable, and this instability corresponds to the formation of rills, which in reality then grow and coalesce to form large-scale river channels. In this paper we consider the deduction and mathematical analysis of a deterministic model describing river channel formation and the evolution of its depth. The model involves a degenerate nonlinear parabolic equation (satisfied on the interior of the support of the solution) with a super-linear source term and a prescribed constant mass. We propose here a global formulation of the problem (formulated in the whole space, beyond the support of the solution) which allows us to show the existence of a solution and leads to a suitable numerical scheme for its approximation. A particular novelty of the model is that the evolving channel self-determines its own width, without the need to pose any extra conditions at the channel margin.     

The influence of alongshore and cross‐shore wave energy flux on large‐ and small‐scale coastal erosion patterns

A strong low‐pressure system traveled along the Japanese main island Honshu in October 2006. High waves and storm surge attacked the Kashima Coast resulting in huge erosion over the area. Airborne laser data measured in October 2005 and November 2006 were analyzed to estimate cross‐sectional changes within the subaerial zone. The results of the alongshore distribution of the changes of cross‐sectional area indicate that the amount of erosion of the 38 km‐long northern and 15 km‐long southern parts decreased toward the south in each part and that the amount of erosion was smaller in protected areas with artificial headlands than in unprotected areas. The local alongshore variation of the erosion and accretion patterns showed wavy fluctuations of several hundreds of meters. The total amounts of the estimated eroded volume of the subaerial zone over the northern and southern parts were 620 000 m³ and 600 000 m³, respectively. The Simulating Waves Nearshore (SWAN) wave model was applied to estimate wave conditions along the coast during the storm. The computational results were verified, and then the alongshore distribution of wave energies, expressed as the alongshore and cross‐shore components of the wave energy flux, was compared with the alongshore distribution of cross‐sectional change. The results show that the distribution of energy flux explains the distribution of erosion well: The alongshore variability in the cross‐shore energy flux is responsible for the large‐scale variability in erosion, and shorter‐scale variability is due to gradients in the alongshore energy fluxes, especially for the areas without coastal works. Copyright © 2011 John Wiley & Sons, Ltd.     

Study on the processing method of nighttime CO_2 eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the "true value" of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u., we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u.c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.     

Temperature and heat flux spectra in the turbulent buoyancy subrange

The physical nature of motions with scales intermediate between approximately isotropic turbulence and quasi-linear internal gravity waves is not understood at the present time. Such motions play an important role in the energetics of small scales processes, both in the ocean and in the atmosphere, and in vertical transport of heat and constituents. This scale range is currently interpreted either as a saturated gravity waves field or as a buoyancy range of turbulence.We first discuss some distinctive predictions of the classical (Lumley, Phillips) buoyancy range theory, recently improved (Weinstock, Dalaudier and Sidi) to describe potential energy associated with temperature fluctuations. This theory predicts the existence of a spectral gap in the temperature spectra and of an upward mass flux (downward buoyancy and heat fluxes), strongly increasing towards large scales. These predictions are contrasted with an alternate theory, assuming energetically insignificant buoyancy flux, proposed by Holloway.Then we present experimental evidences of such characteristic features obtained in the lower stratosphere with an instrumented balloon. Spectra of temperature, vertical velocity, and cospectra of both, obtained in homogeneous, weakly turbulent regions, are compared with theoretical predictions. These results are strongly consistent with the improved classical buoyancy range theory and support the existence of a significant downward heat flux in the buoyancy range.The theoretical implications of the understanding of this scale range are discussed. Many experimental evidences consistently show the need for an anisotropic theory of the buoyancy range of turbulence.     

Derivation Of Refractive Index And Temperature Gradients From Optical Scintillometry To Correct Atmospherically Induced Errors For Highly Precise Geodetic Measurements

Refraction effects of optical beams are generally caused by an inhomogeneous propagation medium and are a major source of systematic errors in the precise optical determination of angles and distances in the atmospheric surface layer. In this contribution a method for deriving vertical temperature and refractive index gradients from optical scintillation is presented. Knowledge of these gradients is required for the compensation of atmospherically induced errors for highly precise terrestrial geodetic measurements, like direct transfer and levelling. The advantage of the present optical method is, that temperature and refractive index gradients can be derived as line-averaged values over the propagation path, which is not possible by meteorological point measurements. Field observations have been carried out with a displaced-beam scintillometer over flat terrain and under different atmospheric conditions in order to verify this method. The experiments show, that this method allows to derive accurate correction values for precise terrestrial geodetic measurements.     

Characterization of seismic hazard and structural response by energy flux

Seismic safety of structures depends on the structure's ability to absorb the seismic energy that is transmitted from ground to structure. One parameter that can be used to characterize seismic energy is the energy flux. Energy flux is defined as the amount of energy transmitted per unit time through a cross-section of a medium, and is equal to kinetic energy multiplied by the propagation velocity of seismic waves. The peak or the integral of energy flux can be used to characterize ground motions. By definition, energy flux automatically accounts for site amplification. Energy flux in a structure can be studied by formulating the problem as a wave propagation problem. For buildings founded on layered soil media and subjected to vertically incident plane shear waves, energy flux equations are derived by modeling the building as an extension of the layered soil medium, and considering each story as another layer. The propagation of energy flux in the layers is described in terms of the upgoing and downgoing energy flux in each layer, and the energy reflection and transmission coefficients at each interface. The formulation results in a pair of simple finite-difference equations for each layer, which can be solved recursively starting from the bedrock. The upgoing and downgoing energy flux in the layers allows calculation of the energy demand and energy dissipation in each layer. The methodology is applicable to linear, as well as nonlinear structures.