Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.     

Biochemical characteristics of alga-bacterial mats and invertebrates from shallow-water hydrothermal fields of the West Pacific Ocean

In this study we assessed the biochemical adaptations of biotic specimens from the hydrothermal vent fields from different geographical zones. Algabacterial and bacterial mats and specimens of bivalves (9 species) and sea urchins (2 species) from different coastal shallow-water fields of hydrothermal activity along the West Rift zones of the Pacific Ocean were sampled to compare biochemical activities to volcanogenic chemical characteristics. The algabacterial and bacterial mat cells and subcellular fractions of bivalve gills and sea urchin guts were measured for 5-aminolevulinic acid (ALA), cytochromes b, (c + c₁), (a + a₃) and P450 (CYP), and activities of ALA synthetase and CYP-dependent aminopyrine N-demethylase (APND). The algabacterial mats from the Kuril Islands and New Zealand displayed comparable ALA and cytochrome levels and enzyme activities. Similar ALA synthesis, APND activity, and levels of cytochromes b and (c + c₁) were observed in bacterial mats from hydrothermal fields from Kuriles, New Britain Island and New Zealand. More active ALA synthesis, APND activity and enhanced content of cytochromes b and (c + c₁) and lower levels of CYP were found in bacterial mats than in algabacterial mats. Further, all bivalves influenced by hydrothermal discharges had elevated CYP levels and APND activity as well as a significant increase in ALA content and ALA synthesis (P < 0.05). Statistically sufficient (P < 0.05) alterations in ALA synthesis, APND activity and total CYP level were found in sea urchins under hydrothermal influence. Our findings contribute to the global assessment of hydrothermal effluents on biota and indicate that living conditions near hydrothermal vents accelerate metabolism of bivalves and sea urchins in examined hydrothermal zones in the northern and southern West Pacific Ocean.     

The Gunn–Peterson effect and the Lyman alpha forest

We show that spatial correlations in a stochastic large-scale velocity field in an otherwise smooth intergalactic medium (homogeneous comoving density) superposed on the general Hubble flow may cause a 'line-like' structure in QSO spectra similar to the population of unsaturated Lyα forest lines which usually are attributed to individual clouds with 10¹¹ ≲ N _{H i}  5 × 10¹³ cm⁻². Therefore there is no clear observational distinction between a diffuse intergalactic medium and discrete intergalactic clouds. It follows that the H  i density in the diffuse intergalactic medium might be substantially underestimated if it is determined from the observed intensity distribution near the apparent continuum in high-resolution spectra of QSOs. Our tentative estimate implies a diffuse neutral hydrogen opacity τ_GP ∼ 0.3 at z  ∼ 3 and a current baryon density Ω_IGM ≃ 0.08, assuming a Hubble constant H ₀ = 70 km s⁻¹ Mpc⁻¹.     

Scattering parameters of the lithosphere below the Massif Central, France, from teleseismic wavefield records

Teleseismic P -wave recordings are analysed in the frequency range 0.3–6  Hz to derive structural (statistical) parameters of the lithosphere underneath the French Massif Central. For this we analyse differences in frequency-dependent intensities of the mean wavefield and the fluctuation wavefield. It is possible to discriminate a weak fluctuation regime of the wavefield in the frequency range below 1  Hz and a strong fluctuation regime starting above 1  Hz and continuing to higher frequencies. The observed wavefield fluctuations in the frequency range 0.3–3  Hz can be explained by scattering of the teleseismic P wave front at elastic inhomogeneities in the lithosphere. A statistical distribution of the inhomogeneities is assumed and the concept of random media is applied. The lithospheric structure under the Massif Central can be described as a 70  km thick heterogeneous layer with velocity fluctuations of 3–7 per cent and correlation lengths of the heterogeneities of 1–16  km.     

Influence Of The Atmosphere On Earth Rotation: What New Can Be Learned From The Recent Atmospheric Angular Momentum Estimates?

The interaction between the atmosphere and the underlying solid mantle is oneof the most important sources of changes in all three components of theEarth's rotation vector on different time scales. In this paper the NCEP/NCARreanalysis time series of four times daily atmospheric effective angularmomentum (EAM) estimates is used to investigate some selected aspects of theatmospheric influence on Earth rotation. Emphasis is placed on thecontroversial features which were difficult or impossible to study using theoperational EAM data, such as excitation of the free oscillations in polarmotion, the Chandler wobble (CW) and the free core nutation (FCN), or theroles of diurnal and semidiurnal atmospheric tides and atmospheric normalmodes in the rotational dynamics of the Earth.     

Post-collisional melting of crustal sources: constraints from geochronology,petrology and Sr,Nd isotope geochemistry of the Variscan Sichevita and Poniasca granitoid plutons (South Carpathians,Romania)

The Sichevita and Poniasca plutons belong to an alignment of granites cutting across the metamorphic basement of the Getic Nappe in the South Carpathians. The present work provides SHRIMP age data for the zircon population from a Poniasca biotite diorite and geochemical analyses (major and trace elements, Sr–Nd isotopes) of representative rock types from the two intrusions grading from biotite diorite to biotite K-feldspar porphyritic monzogranite. U–Pb zircon data yielded 311 ± 2 Ma for the intrusion of the biotite diorite. Granites are mostly high-K leucogranites, and biotite diorites are magnesian, and calcic to calc-alkaline. Sr, and Nd isotope and trace element data (REE, Th, Ta, Cr, Ba and Rb) permit distinguishing five different groups of rocks corresponding to several magma batches: the Poniasca biotite diorite (P₁) shows a clear crustal character while the Poniasca granite (P₂) is more juvenile. Conversely, Sichevita biotite diorite (S₁), and a granite (S₂*) are more juvenile than the other Sichevita granites (S₂). Geochemical modelling of major elements and REE suggests that fractional crystallization can account for variations within P₁ and S₁ groups. Dehydration melting of a number of protoliths may be the source of these magma batches. The Variscan basement, a subduction accretion wedge, could correspond to such a heterogeneous source. The intrusion of the Sichevita–Poniasca plutons took place in the final stages of the Variscan orogeny, as is the case for a series of European granites around 310 Ma ago, especially in Bulgaria and in Iberia, no Alleghenian granitoids (late Carboniferous—early Permian times) being known in the Getic nappe. The geodynamical environment of Sichevita–Poniasca was typically post-collisional of the Variscan orogenic phase.     

Distribution and changes of active layer thickness (ALT) and soil temperature (TTOP) in the source area of the Yellow River using the GIPL model

Active layer thickness (ALT) is critical to the understanding of the surface energy balance, hydrological cycles, plant growth, and cold region engineering projects in permafrost regions. The temperature at the bottom of the active layer, a boundary layer between the equilibrium thermal state (in permafrost below) and transient thermal state (in the atmosphere and surface canopies above), is an important parameter to reflect the existence and thermal stability of permafrost. In this study, the Geophysical Institute Permafrost Model (GIPL) was used to model the spatial distribution of and changes in ALT and soil temperature in the Source Area of the Yellow River (SAYR), where continuous, discontinuous, and sporadic permafrost coexists with seasonally frozen ground. Monthly air temperatures downscaled from the CRU TS3.0 datasets, monthly snow depth derived from the passive microwave remote-sensing data SMMR and SSM/I, and vegetation patterns and soil properties at scale of 1:1000000 were used as input data after modified with GIS techniques. The model validation was carried out carefully with ALT in the SAYR has significantly increased from 1.8 m in 1980 to 2.4 m in 2006 at an average rate of 2.2 cm yr^?1. The mean annual temperature at the bottom of the active layer, or temperature at the top of permafrost (TTOP) rose substantially from ?1.1°C in 1980 to ?0.6°C in 2006 at an average rate of 0.018°C yr^?1. The increasing rate of the ALT and TTOP has accelerated since 2000. Regional warming and degradation of permafrost has also occurred, and the changes in the areal extent of regions with a sub-zero TTOP shrank from 2.4×10⁴ to 2.2×10⁴ km² at an average rate of 74 km² yr^?1. Changes of ALT and temperature have adversely affected the environmental stability in the SAYR.