Frequency of quasi-geostrophic modes on hexagonal grids

Summary Finite-difference analysis of Rossby modes has been performed for two staggered hexagonal grids. The solutions are compared with those obtained in analytical case and for rectangular grids. The result for one of the selected hexagonal grids better fits to the analytical solution then the results for the other considered grids. The obtained results may contribute to better understanding of the appropriateness of hexagonal grids in atmospheric and oceanographic modeling and numerical computations.With 6 Figures     

Rječina River sediments (Croatia): from captured spring to polluted prodelta

Surface sediments (fraction <63 μm) from the source to the mouth of the Rječina, short (18.3 km) karst allogenic river in Croatia, which is an important source of drinking water, were studied to investigate their mineral (by XRD) and chemical (by ICP-MS) composition to check possible anthropogenic influence at the lower course due to paper industry and mills, and in the prodelta area from untreated municipal sewage and the large harbor of Rijeka town. In all analyzed sediment samples and in the sandstone source, rock quartz is a major mineral, while feldspar and mica group minerals are less abundant. Chlorite is a minor or trace mineral in all samples. Calcite and dolomite are abundant in the river prodelta, reflecting changes in bedrock lithology from flysch to carbonates. In river sediments, Fe is the most abundantly analyzed element, while Ca is the most abundant in prodelta sediments. Concentrations of Al, Mn, Ni, Cr, Co, La and Nd decrease downstream, while Mg, S, Na, B, Pb, Zn, As, Sn, U, Mo, Hg and Ag have relatively higher concentration in prodelta sediments. The results are compared with sediments of other rivers in the area: Raša, Rižana and Dragonja, as well as with those of the Rosandra Creek (Italy). Sediments in the Raša River showed similar behavior as those in the Rječina, as the highest concentration of metals was found in the restricted upper part of the estuary, characterized by rapid deposition of clay particles and terrestrial sedimentary organic matter. The comparison also showed that the most contaminated were the sediments from the Rižana, followed by those from the Rječina and Rosandra Creek, which had similar results. Among the studied elements, As was present in all sediment samples at concentrations >6 ppm that might have the lowest toxic effects. At the lower Rječina and in prodelta sediments, Pb was also present at slightly elevated concentrations (>31 ppm) that could cause such effects. Concentrations of Zn in the prodelta correspond to those occurring in moderately polluted sediments (90–200 ppm). In the prodelta sediments, Hg is slightly below toxicity threshold (1 ppm), while Ag is present at toxicity threshold (0.5 ppm) or close to it. Rječina River could act as a good illustrative example for behavior of toxic metals in allogenic karstic rivers, in which accumulation of anthropogenically introduced pollutants usually occurs in their estuaries, as a result of transport and deposition of fine particles.     

Hydrogeochemical stratification of the unconfined Samobor aquifer (Zagreb,Croatia)

In the western part of the city of Zagreb, Croatia, hydrogeochemical and isotopic investigations of the Samobor aquifer were carried out with the aim of determining the differences in hydrogeochemical characteristics at increasing aquifer depths. The aquifer comprises 40-m thick gravelly–sandy deposits, with lenses and interlayers of silt and clay. The analyses have proven that with increasing aquifer depth, there are decreases in groundwater temperature and the values of electrical conductivity and increases in the sodium, iron and manganese concentrations. The δ¹³C distribution shows an evident increase in biogenic carbon concentrations with increasing aquifer depth. The measured specific ¹⁴C activities showed that the deeper part of the aquifer is characterized by slow water exchange, while the shallower part is influenced by current recharge, although the pumping wells located on the well-field downstream penetrate the aquifer fully. A direct exchange of water from the Sava River and groundwater occurs in the near vicinity of the river. This exchange weakens further away, while the difference in hydrogeochemical characteristics between the Sava River water and groundwater increases.     

A note on the performance of the multiply-upstream semi-Lagrangian advection schemes for one-dimensional nonlinear momentum conservation equation

Summary Two-time-level multiply-upstream semi-Lagrangian schemes were examined in the case of the self-advecting, one-dimensional nonlinear momentum conservation equation. The shock formation process was analyzed. It is pointed out that the shocks cannot be created in the truncated systems satisfying the Pudykiewicz, Benoit and Staniforth criterion.The numerical integrations were restricted to 12 h. It was shown that, at least in the sub-CFL range, increased complexity of the scheme can compensate reduced horizontal resolution. A considerable sensitivity of the schemes with respect to the time step was detected. In the super-CFL mode, several windows on various time scales were found within which the Pudykiewicz, Benoit and Staniforth criterion was satisfied. The time step of 1.44 times the maximum time step allowed by the CFL criterion was used in the semi-Lagrangian runs.The super-CFL, semi-Lagrangian solutions were diverging progressively from the sub-CFL ones as the forecasts advanced. This was also reflected in the energy spectra.Unacceptably large energy losses were encountered in the super-CFL, semi-Lagrangian runs. Most of these losses could be explained by the reduced mean wind speed, i.e., the amplitude of the zero wavenumber wave. At the same time, the energy content in the shorter waves increased. In a more complex model, such a situation would resemble a loss of zonal, and an increase of transient eddy kinetic energy.A trajectory error measure was defined as the maximum absolute value of the distance between the actual arriving point of the particle originating at the estimated departure point, and the grid point assumed to be the arrival point in the semi-Lagrangian procedure. In contrast to the sub-CFL regime, this measure could reach a considerable fraction of the grid distance in the computations with the super-CFL time steps.In the physical system considered, the trajectories are determined only by the velocities at the departure points. With the semi-Lagrangian schemes the distances traveled by the particles are estimated on the basis of the velocities at the points downstream with respect to the departure points. Thus, unless the solution is smooth (in space and time) on the scales of the extrapolation distances/times, the upstream extrapolation does not promise the convergence of the solution.With 16 Figures     

Improvement of Vertical Diffusion Analytic Schemes Under Stable Atmospheric Conditions

Based on gradient transport theory or K-theory, turbulent transport in the atmosphere has long been parameterized using the eddy diffusivity. Due to its simplicity, this approach has often been applied in many numerical models but rarely tested with observations. Here, the widely used O’Brien cubic polynomial approach has been validated together with an exponential approach against eddy diffusivity profiles determined from measurements and from large-eddy simulation data in stable conditions. Verification is completed by analyzing the variability effects on pollutant concentrations of two different vertical diffusion (K(z)) schemes incorporated in an atmospheric chemical model. It is shown that the analytical, exponential solution agrees better with observations than the O’Brien profile and should be used henceforth in practical applications.     

Magnetic field measurements in India during the last two centuries

Summary The direction of the geomagnetic field has changed substantially less in India than in Europe since the end of the 18th century.Whereas the European geomagnetic time series indicate a tendency of secular variation from which a remarkable westward drift of the non-axial field can be deduced, in India the field direction seems to have remained more or less fixed over the last two centuries. The drift of the non-axial field seems to be very small, or even eastward.