A systematic algorithm for obtaining SODAR wind profiles

Summary Current methods of obtaining wind profiles from SODARs generally employ ad hoc strategies for Doppler spectrum peak selection, profile smoothing, and consistency checks. This paper addresses the problem of systematically obtaining vertical profiles of wind vectors from a SODAR. A new approach is to regard every spectrum frequency as a possible estimator of the radial wind component, but to weight these Doppler estimates according to the spectral power at each frequency. This allows the peak detection, smoothing, and assimilation of a priori information to be combined into a constrained linear inversion methodology. The inversion process also allows an estimation of signal-to-noise ratio, degrees of freedom, and information content versus height. Examples are given of profiles obtained using this new inverse method. The relationship to the physical limitations of the SODAR instrument and the atmosphere is discussed.     

A Simple Parameterisation for Flux Footprint Predictions

Flux footprint functions estimate the location and relative importance of passive scalar sources influencing flux measurements at a given receptor height. These footprint estimates strongly vary in size, depending on receptor height, atmospheric stability, and surface roughness. Reliable footprint calculations from, e.g., Lagrangian stochastic models or large-eddy simulations are computationally expensive and cannot readily be computed for long-term observational programs. To facilitate more accessible footprint estimates, a scaling procedure is introduced for flux footprint functions over a range of stratifications from convective to stable, and receptor heights ranging from near the surface to the middle of the boundary layer. It is shown that, when applying this scaling procedure, footprint estimates collapse to an ensemble of similar curves. A simple parameterisation for the scaled footprint estimates is presented. This parameterisation accounts for the influence of the roughness length on the footprint and allows for a quick but precise algebraic footprint estimation.     

Numerical Studies with a Regional Atmospheric Climate Model Based on Changes in the Roughness Length for Momentum and Heat Over Antarctica

A regional atmospheric climate model is used toexamine the effect of changes in the roughnesslengths of momentum (_z0m) and heat (_z0h)on the structure of the lower atmosphere and on thesurface energy fluxes over Antarctica. Fourexperiments were carried out in which _z0mand/or _z0h were altered with respect to acontrol experiment. The changes consisted of (1) alowering of _z0m from a field aggregated froma vegetation map with an orographic correction basedon the European Centre for Medium-Range WeatherForecasts _z0m field, to a constant value of10^-3 m; and (2) a lowering of _z0h from a valueequal to _z0m to a constant value of 10^-3 mor a value dependent on the wind speed via a surfacerenewal model. A reduction of _z0m results in theexpected increase in near-surface wind speed. It alsoresults in an increase in the depth of the layer in whichsouth-easterly near-surface winds prevail, and in adecrease in the strength of the large-scale flow overthe continent, in particular in summer. In theescarpment region a decrease of _z0m is foundto result in too high wind speeds. Surface temperatureson average decrease while atmospheric temperaturesincrease, resulting in an increase of near-surfacestatic stability. Changes in roughness lengths donot significantly change the temperature profiles.The surface fluxes, on average found reduced, aremodelled best by using the _z0h based on thesurface renewal method.     

Tree-Ring Based Reconstructions of Precipitation for the Southern Canadian Cordillera

This paper reconstructs precipitation variability in the southern Canadian Cordillera over the past 3–400 years using dendroclimatologicaltechniques. Fifty-three total ring-width (RW) chronologies, 28 earlywood (EW) and 28 latewood (LW) chronologies were developed from open-grown, low-elevation stands of Pseudotsuga menziesii (Douglas-fir) and Pinusponderosa (ponderosa pine) across the southern Canadian Cordillera. RW, EW and LW chronologies from both species were used to develop 13 annual (prior July to current June) precipitation reconstructions across the region. The reconstructions range in length from 165 to 688 years, pass verification tests and capture 39–64% of the variancein the instrumental record. Coincident, prolonged intervals of dry conditions are estimated for the years: 1717–1732, 1839–1859, 1917–1941 and1968–1979. Shorter dry intervals are identified between 1581–1586, 1626–1630,1641–1653 1701–1708, 1756–1761, 1768–1772, 1793–1800,1868–1875, 1889–1897 and 1985–1989. The historic drought of the 1920–1930s was the longest but not the most intense across this area in the last 300 years. Wet conditions occur in the majority of reconstructions for the years: 1689–1700, 1750–55,1778–1789, 1800–1830, 1880–1890, 1898–1916 and 1942–1960. Thesedata, in conjunction with data from adjacent areas, are used to provide the first maps of decadal precipitation anomalies for the region between 1700 and 1990.     

A Least-squares Approach to Depth Determination from Numerical Horizontal Self-potential Gradients

We have developed a least-squares minimization approach to depth determination of a buried ore deposit from numerical horizontal gradients obtained from self-potential (SP) data using filters of successive window lengths (graticule spacings). The problem of depth determination from SP gradients has been transformed into the problem of finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for vertical and horizontal cylinders and spheres. Procedures are also formulated to estimate the electrical dipole moment and the polarization angle. The method is applied to synthetic data with and without random noise. Finally, the validity of the method is tested on two field examples. In both cases, the depth obtained is found to be in a very good agreement with that obtained from drilling information.     

Effective elastic properties of randomly fractured soils: 3D numerical experiments

This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results.     

Managed water levels and the expansion of emergent vegetation along a lakeshore

Water-level fluctuations may be used to promote the expansion of emergent vegetation along lakeshores. We present the case of the lake Volkerak-Zoommeer in the Netherlands, a fresh-water lake created in 1987 after the enclosure of an estuary. Using an experimental area in which the water level could be manipulated, it was shown that partial summer drawdown of the shoreline created suitable conditions for germination and growth of tall emergent species (in particular Phragmites australis). Plant survival and growth depended on subsequent water-level fluctuations and grazing by waterbirds. Based on the experiment and empirical data, a model was developed to predict the effects of the water-level regime on potential reed bed development. The model was applied for four hydrological scenarios that have been considered for the water-level management of the lake.