Eigenrays in 3D heterogeneous anisotropic media,Part II: Dynamics

This paper is the second in a sequel of two papers and dedicated to the computation of paraxial rays and dynamic characteristics along the stationary rays obtained in the first paper. We start by formulating the linear, second‐order, Jacobi dynamic ray tracing equation. We then apply a similar finite‐element solver, as used for the kinematic ray tracing, to compute the dynamic characteristics between the source and any point along the ray. The dynamic characteristics in our study include the relative geometric spreading and the phase correction due to caustics (i.e. the amplitude and the phase of the asymptotic form of the Green's function for waves propagating in 3D heterogeneous general anisotropic elastic media). The basic solution of the Jacobi equation is a shift vector of a paraxial ray in the plane normal to the ray direction at each point along the central ray. A general paraxial ray is defined by a linear combination of up to four basic vector solutions, each corresponds to specific initial conditions related to the ray coordinates at the source. We define the four basic solutions with two pairs of initial condition sets: point–source and plane‐wave. For the proposed point–source ray coordinates and initial conditions, we derive the ray Jacobian and relate it to the relative geometric spreading for general anisotropy. Finally, we introduce a new dynamic parameter, similar to the endpoint complexity factor, presented in the first paper, used to define the measure of complexity of the propagated wave/ray phenomena. The new weighted propagation complexity accounts for the normalized relative geometric spreading not only at the receiver point, but along the whole stationary ray path. We propose a criterion based on this parameter as a qualifying factor associated with the given ray solution. To demonstrate the implementation of the proposed method, we use several isotropic and anisotropic benchmark models. For all the examples, we first compute the stationary ray paths, and then compute the geometric spreading and analyse these trajectories for possible caustics. Our primary aim is to emphasize the advantages, transparency and simplicity of the proposed approach.     

Use of a High-Resolution Sodar to Study Surface-layer Turbulence at Night

Measurements in the atmospheric surface layer are generally made with point sensors located in the first few tens of metres. In most cases, however, these measurements are not representative of the whole surface layer. Standard Doppler sodars allow a continuous display of the turbulent thermal structure and wind profiles in the boundary layer up to 1000 m, with a few points, if any, in the surface layer. To overcome these limitations a new sodar configuration is proposed that allows for a higher resolution in the surface layer. Because of its capabilities (echo recording starting at 2 m, echo intensity vertical resolution of approximately 2 m, temporal resolution of 1 s) this sodar is called the surface-layer mini-sodar (SLM-sodar). Features and capabilities of the SLM-sodar are described and compared with the sodar. The comparison of the thermal vertical structure given by the SLM-sodar and the sodar provides evidence that, in most cases, the surface layer presents a level of complexity comparable to that of the entire boundary layer. Considering its high vertical resolution, the SLM-sodar is a promising system for the study of the nocturnal surface layer. The nocturnal SLM-sodar measurements have shown that, depending on wind speed, the structure of the surface layer may change substantially within a short time period. At night, when the wind speed is greater than 3 m s⁻¹, mechanical mixing destroys the wavy structure present in the nocturnal layer. Sonic anemometer measurements have shown that, in such cases, also the sensible heat flux varies with height, reaching a peak in correspondence with the wind speed peak. Under these conditions the assumption of horizontal homogeneity of the surface layer and the choice of the averaging time need to be carefully treated.     

Determining the main atmospheric factor on ozone concentrations

This paper presents techniques for studying the influence of the atmospheric variables for the explanation of ozone concentrations. A methodology is described for separating the different time scale components in time series of ozone, namely, the global-term component which describes the long and seasonal variations, and the synoptic-scale component. We prove that the decomposition of the time series of ozone and other atmospheric variables is essential for the analysis. The results indicate that solar radiation is the only main factor between all the atmospheric variables for the global and the synoptic-scale component of the ozone time series. Most of the total variance of ozone can be explained by the atmospheric variables, which makes possible the precise estimation of human effect on ozone. In case the solar radiation data are not available, lagged temperature can be used as a supplement of solar radiation.     

Geo–ecology and management of sensitive montane landscapes

Montane (alpine) areas are generally of high value for nature conservation. Such environments and the habitats they support are dynamic and often fragile. They are vulnerable to disturbance from a range of human activities and are responsive to climate changes over short and long timescales. Biodiversity and conservation values are closely linked to geological history, geomorphological processes and soils, and it is crucial that management systems are based on understanding these links.
There are many similarities between the Cairngorm Mountains (Scotland), the Giant Mountains (Czech Republic) and Abisko Mountains (Sweden) in terms of geology, geomorphology, ecology, links with biodiversity and high conservation importance. Comparable pressures and management issues involve, to varying degrees, a history of human use and impacts from deforestation, pasturing, grazing, recreation and atmospheric pollution. Landscape change therefore involves a complex interplay between natural and anthropogenic factors. Managing such change requires better understanding of the geo–ecological processes involved and the factors that determine landscape sensitivity. This is illustrated through a simple framework and examples from the three areas. Comparison of landscape sensitivity between similar montane areas, but in different geographic locations and climatic environments, should allow more informed management planning and a precautionary approach in advance of further changes in human activity and from predicted global warming scenarios.     

Current estimates of freshwater flux through Arctic and subarctic seas

As the world warms, the expectation is that the freshwater outflows from the Arctic Ocean to the North Atlantic will strengthen and may act to suppress the rate of the climatically-important Atlantic meridional overturning circulation. Hitherto, however, we have lacked the system of measurements required to estimate the totality of the freshwater flux through subarctic seas. Though observations remain patchy and rudimentary in places, we piece-together the results from recent large-scale observational programmes together with associated modelling, to establish preliminary maps of the rates and pathways of freshwater flux through subarctic seas. These fluxes are calculated according to two reference salinities, S = 34.8 to conform with the majority of estimates reported in the literature, and S = 35.2, the salinity of the inflowing Atlantic water, to calculate the freshwater balance of the ‘Arctic Mediterranean’. We find that 148 mSv of freshwater enters the Nordic Seas across its northern boundary. There it is supplemented by around 54 mSv of freshwater from Baltic runoff, Norwegian runoff, P − E and Greenland ice melt, so that the total freshwater contribution to the Nordic Seas from all sources is 202 mSv. Of this, around 51 mSv of freshwater is estimated to pass south to the deep Atlantic in the dense water overflows leaving an assumed balance of 151 mSv to leave the Nordic Seas in the upper water export through Denmark Strait. The corresponding estimate for the freshwater outflow west of Greenland is 103 mSv relative to 35.2 so that the total freshwater flux reaching the North Atlantic through subarctic seas is around 300 mSv.     

Finite element modeling of net panels using a consistent net element

A consistent finite element is proposed to model the hydrodynamic response of net panels to environmental loading. This equivalent net element is constructed to reproduce the drag, buoyancy, inertial and elastic forces exerted on the netting by current and waves. To evaluate the accuracy of the proposed finite element modeling, numerical predictions have been compared with the experimental observations and (simplified) analytical results of other authors. This new modeling technique has been applied to evaluate the performance of a tension leg fish cage in the open ocean environment.     

Effects of the Anticyclonic Eddies on Water Masses, Chemical Parameters and Chlorophyll Distributions in the Oyashio Current Region

Data from the R/V Mirai cruise (May–June 2000) have been examined to discover how mesoscale processes associated with eddy dynamics direct affect the water masses, the distributions and the vertical fluxes of the dissolved oxygen, nutrients and dissolved inorganic carbon in the western subarctic Pacific. Using maps of the temperature, salinity, dissolved oxygen, nutrients, chlorophyll and sea-air pCO₂ difference we show that the boundaries of the anticyclone eddies in the study region were composed of high productivity coastal Oyashio water. The coastal waters were wrapped around the anticyclone eddies (thus creating a high productivity belt) and intruded inside of them. Using SeaWifs data we demonstrate that temporal variations in the position and the strength of anticyclone eddies advected the Kuril island coastal high productivity waters to the pelagic part, resulting in temporal variations of the chlorophyll in the Oyashio region. Computed vertical fluxes of the dissolved oxygen (DO), inorganic carbon (DIC) and silicate show that the anticylonic eddies in the Kuroshio-Oyashio Zone are characterized by enhanced vertical fluxes of the DO and DIC between the upper (σ_θ = 26.7–27.0) and lower (σ_θ = 27.1–27.5) intermediate layer, probably due to the intrusions of the Oyashio waters into the eddies. This revised version was published online in August 2006 with corrections to the Cover Date.     

Modelling hydrographic changes in the Labrador sea over the past five decades

Inter-annual to inter-decadal changes of hydrographic structure and circulation in the subpolar North Atlantic are studied using a coarse resolution ocean circulation model. The study covers 1949 through 2001, inclusive. A “time-mean state nudging” method is applied to assimilate the observed hydrographic climatology into the model. The method significantly reduces model biases in the long-term mean distribution of temperature and salinity, which commonly exist in coarse-resolution ocean models. By reducing the time-mean biases we also significantly improve the model’s representation of inter-annual to inter-decadal variations. In the central Labrador Sea, the model broadly reproduces the heat and salt variations of the Labrador Sea Water (LSW) as revealed by hydrographic observations. Model sensitivity experiments confirm that the low-frequency hydrographic changes in the central Labrador Sea are closely related to changes in the intensity and depth of deep convection. Changes in surface heat flux associated with the winter North Atlantic Oscillation (NAO) index play a major role in driving the changes in T–S and sea surface height (SSH). Changes in wind stress play a secondary role in driving these changes but are important in driving the changes in the depth-integrated circulation. The total changes in both SSH and depth-integrated circulation are almost a linear combination of the separate influences of variable buoyancy and momentum fluxes.     

Geochronology of the Larderello geothermal field: new data and the “closure temperature” issue

The Larderello geothermal field is generally accepted to have been produced by a granite intrusion at 4–9 km depth. Hydrothermal parageneses and fluid inclusions always formed at temperatures greater than or equal to the current ones, which implies that the field has always undergone a roughly monotonic cooling history (fluctuations < 40 K) since intrusion of the granite at 4 Ma. The heat required to maintain the thermal anomaly over such a long period is supplied by a seismically anomalous body of 32000 km³ rooted in the mantle. Borehole minerals from Larderello are thus a unique well-calibrated natural example of thermally induced Ar and Sr loss under geological conditions and time spans. The observations (biotites retain Ar above 450°C) agree well with other, albeit less precise, geological determinations, but contrast with laboratory determinations of diffusivity from the literature. We therefore performed a hydrothermal experiment on two Larderello biotites and derived a diffusivity D _Lab(370°C)=5.3·10^-18 cm²s^-1, in agreement with published estimates of diffusivity in annite. From D _Lab and the rejuvenation of the K/Ar ages we calculate maximum survival times at the present in-hole temperatures. They trend smoothly over almost two orders of magnitude from 352 ka to 5.3 ka, anticorrelating with depth: laboratory diffusivities are inconsistent not only with geological facts, but also among themselves. From the geologically constrained lifetime of the thermal anomaly we derive a diffusivity D _G(370°C)=3.81·10²¹ cm²s^-1, 3±1 orders of magnitude lower than D _Lab. The cause of these discrepancies must be sought among various laboratory artefacts: overstepping of a critical temperature T ^*; enhanced diffusivities in wet experiments; presence of fast pathway (dislocation and pipe) diffusion, and of dissolution/reprecipitation reactions, which we imaged by scanning electron microscopy. These phenomena are minor in geological settings: in the absence of mineral transformation reactions, complete or near-complete resetting is achieved only by volume diffusion. Therefore, laboratory determinations will necessarily result in apparent diffusivities that are too high compared to those actually effecting the resetting of natural geochronometers.This word is dedicated to the memory of Aldo Valori (1958–1991)     

Management of karst aquifers in Serbia for water supply

Karst aquifers represent one of the main sources of water supply in Serbia. Wide distributions of karstic areas, abundant reserves, and excellent quality of karst groundwater have been the reasons for its extensive use in water supply systems throughout the country. In Serbia, 70 karstic sources have been tapped for centralized water supply, with the assessed minimal capacity exceeding 4.5 m³/s. Most of the large cities of Eastern and Western Serbia currently use karst groundwater for water supply; however, due to unstable flow regime when only natural springflow is tapped, numerous problems arise during the recession period (summer–autumn). During the last two decades, after favorable conditions had been confirmed through hydrogeological survey and feasibility studies, several newly constructed systems for artificial control of karst aquifers (mostly in Eastern Serbia) resulted in significantly improved water supply.