Self-gravity in curved mesh elements

The local character of self-gravity along with the number of spatial dimensions are critical issues when computing the potential and forces inside massive systems like stars and disks. This appears from the discretisation scale where each cell of the numerical grid is a self-interacting body in itself. There is apparently no closed-form expression yet giving the potential of a three-dimensional homogeneous cylindrical or spherical cell, in contrast with the Cartesian case. By using Green’s theorem, we show that the potential integral for such polar-type 3D sectors—initially, a volume integral with singular kernel—can be converted into a regular line-integral running over the lateral contour, thereby generalising a formula already known under axial symmetry. It therefore is a step towards the obtention of another potential/density pair. The new kernel is a finite function of the cell’s shape (with the simplest form in cylindrical geometry), and mixes incomplete elliptic integrals, inverse trigonometric and hyperbolic functions. The contour integral is easy to compute; it is valid in the whole physical space, exterior and interior to the sector itself and works in fact for a wide variety of shapes of astrophysical interest (e.g. sectors of tori or flared discs). This result is suited to easily providing reference solutions, and to reconstructing potential and forces in inhomogeneous systems by superposition. The contour integrals for the 3 components of the acceleration vector are explicitely given.     

Temperature effect on the photodissociation rates in the atmospheres of Mars and Venus

The available solar flux at a given altitude in the atmospheres of Mars and Venus is attenuated mainly by CO₂ (molecular absorption and Rayleigh scattering) with an extra contribution due to SO₂ on Venus. The dissociation cross section of CO₂ depends on temperature. At temperatures appropriate for these atmospheres (～250°K), the cross sections are about 15% lower than those at room conditions (Y.L. Yung and W.B. De More, 1982, Icarus, 51, 199). It is shown that this temperature effect cannot be neglected in the evaluation of photolysis rates. Calculations of the photodissociation coefficients of CO₂, SO₂, HCl, and H₂O are presented. For example, at the surface of Mars, the coefficient of H₂O is nearly multiplied by a factor of 10!     

Impact of spatial resolution on simulated surface water mass transformations in the Atlantic

We analyze the water mass transformation in coarse (1°) and high (1/6°) resolution ocean simulations with the identical configuration of the CLIPPER model and interannual ERA15 forcing function. Climatological characteristics of surface water mass transformation in the two experiments are quite different. The high resolution experiment exhibits a stronger surface transformation in equatorial and tropical regions, in the Gulf Stream area and in the location of the formation of Subtropical Mode Water (STMW), associated with high levels of eddy kinetic energy. The coarse resolution experiment shows a better representation of the transformation rates corresponding to the densest subpolar mode waters and Labrador Sea Water (LSW). This is explained by the differences in lateral mixing procedures between high and coarse resolution experiments. The high resolution 1/6° run is eddy-resolving only in the tropics and mid-latitudes. In these areas eddies are found to enhance the process of water mass transformation compared to the isopycnal diffusion used to parameterized the eddies in the 1° model. Despite its 1/6° resolution, the high resolution model does not adequately represent eddies in the subpolar gyre and Labrador Sea. In these areas the high resolution model fails to correctly simulate water mass transformation because the lateral mixing (provided through the bi-harmonic sub-gridscale parameterization) of newly ventilated waters with surrounding waters is not efficient enough. In contrast in the coarse 1° resolution model, the strong lateral mixing and the unrealistically broad boundary currents imposed by the high diffusivity required for numerical stability mixes newly formed LSW waters with the warmer and saltier waters of the rim current. Finally, it results in a more effective representation of the surface water mass transformation in high latitudes in the 1° model. A possible impact of the increased lateral diffusion in high resolution experiment on the representation of re-stratification in the Labrador Sea was studied in sensitivity experiments with different lateral diffusion coefficients compared to the regional eddy-resolving 1/15° simulation in the subpolar North Atlantic. If the eddies are not resolved in subpolar latitudes (as in the case of 1/6° model), the GM90 parameterization with the coefficient close to 800 m² s⁻¹ provides the closest agreement with the solution of eddy-resolving 1/15° model.     

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.     

From rifting to spreading in the eastern Gulf of Aden: a geophysical survey of a young oceanic basin from margin to margin

A geophysical survey in the eastern Gulf of Aden, between the Alula–Fartak (52°E) and the Socotra (55°E) transform faults, was carried out during the Encens–Sheba cruise. The conjugate margins of the Gulf are steep, narrow and asymmetric. Asymmetry of the rifting process is highlighted by the conjugate margins (horst and graben in the north and deep basin in the south). Two transfer fault zones separate the margins into three segments, whereas the present‐day Sheba Ridge is divided into two segments by a transform discontinuity. Therefore segmentation of the Sheba Ridge and that of the conjugate margins did coincide during the early stages of oceanic spreading. Extensive magma production is evidenced in the central part of the western segment. Anomaly 5d was identified in the northern and southern parts of the oceanic basin, thus confirming that seafloor spreading in this part of Gulf of Aden started at least 17.6 Ma ago.     

Effects of hypoxia on benthic macrofauna and bioturbation in the Estuary and Gulf of St. Lawrence,Canada

The bottom water in the >300 m deep Lower St. Lawrence Estuary (LSLE) is persistently hypoxic in contrast to the normoxic bottom waters in the Gulf of St. Lawrence (GSL). We photographed the seabed at 11 stations in the Estuary and Gulf of St. Lawrence (EGSL) during the summers 2006 and 2007 and analysed the images to identify bioturbation traces (lebensspuren) and benthic macrofauna. The objective was to identify the environmental variables that influence the density and diversity of benthic macrofauna and bioturbation traces, and the differences that exist among regions with high, medium and low oxygen levels in the bottom water. The bottom water oxygen concentration is the variable that best explains the densities of total-traces as well as surface-traces. However, the density of these traces was higher in hypoxic regions than in well-oxygenated regions. The higher density of traces in the hypoxic region of the LSLE is mainly due to the activities of the surface deposit feeder Ophiura sp., which occurs in large numbers in this region. Possible explanations explored are stress behaviour of the organisms in response to hypoxia and different benthic macrofauna community structures between the hypoxic regions of the LSLE and the normoxic regions of the GSL. In the former, surface deposit feeders and low-oxygen tolerant species dominate over suspension feeders and low-oxygen intolerant species.     

The effects of ionic strength and organic matter on virus inactivation at low temperatures: general likelihood uncertainty estimation (GLUE) as an alternative to least-squares parameter optimization for the fitting of virus inactivation models

This study examined how the inactivation of bacteriophage MS2 in water was affected by ionic strength (IS) and dissolved organic carbon (DOC) using static batch inactivation experiments at 4 °C conducted over a period of 2 months. Experimental conditions were characteristic of an operational managed aquifer recharge (MAR) scheme in Uppsala, Sweden. Experimental data were fit with constant and time-dependent inactivation models using two methods: (1) traditional linear and nonlinear least-squares techniques; and (2) a Monte-Carlo based parameter estimation technique called generalized likelihood uncertainty estimation (GLUE). The least-squares and GLUE methodologies gave very similar estimates of the model parameters and their uncertainty. This demonstrates that GLUE can be used as a viable alternative to traditional least-squares parameter estimation techniques for fitting of virus inactivation models. Results showed a slight increase in constant inactivation rates following an increase in the DOC concentrations, suggesting that the presence of organic carbon enhanced the inactivation of MS2. The experiment with a high IS and a low DOC was the only experiment which showed that MS2 inactivation may have been time-dependent. However, results from the GLUE methodology indicated that models of constant inactivation were able to describe all of the experiments. This suggested that inactivation time-series longer than 2 months were needed in order to provide concrete conclusions regarding the time-dependency of MS2 inactivation at 4 °C under these experimental conditions.