Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

Numerical Studies of Fluid Flow in Microgravity Conditions for Confined Crystal Growth

We study the convective flow induced by residual accelerations in microgravity conditions for different geometric arrangements which are relevant to crystal growth experiments. We consider both constantand oscillating acceleration and focus mostly on the transient relaxation dynamics. Results are relevant to estimate impact of more realistic residual accelerations in crystal growth experiments.     

Microlensing by cosmic strings

We consider the signature and detectability of gravitational microlensing of distant quasars by cosmic strings. Because of the simple image configuration such events will have a characteristic lightcurve, in which a source would appear to brighten by exactly a factor of 2, before reverting to its original apparent brightness. We calculate the optical depth and event rate, and conclude that current predictions and limits on the total length of strings on the sky imply optical depths of ≲ 10⁻⁸ and event rates of fewer than one event per 10⁹ sources per year. Disregarding those predictions but replacing them with limits on the density of cosmic strings from the cosmic microwave background fluctuation spectrum, leaves only a small region of parameter space (in which the sky contains about 3 × 10⁵ strings with deficit angle of the order of 0.3 milli-seconds) for which a microlensing survey of exposure 10⁷ source years, spanning a 20–40-year period, might reveal the presence of cosmic strings.     

Nutrients and zooplankton composition and dynamics in relation to the hydrological pattern in a confined Mediterranean salt marsh (NE Iberian Peninsula)

The role of the hydrological regime in the nutrients and zooplankton composition and dynamics has been analysed in five lagoons of La Pletera salt marshes (NE Iberian Peninsula) during a complete hydrological cycle (2002–2003). Two of the lagoons have their origin in the old river mouths while the other three were recently created in the framework of a Life Restoration project. This fact has also allowed us to study the effect of the lagoon age on nutrient and zooplankton composition and dynamics. The salt marsh hydrology is determined by a prolonged period of confinement without water inputs, irregularly interrupted by sudden water inputs due to flooding events (sea storms or intense rainfalls). While the dynamics of oxidized nitrogen compounds in the lagoons depends on the water inputs variability within each hydrological cycle, the internal load of phosphorus, total nitrogen and organic matter is related more to the cumulative mechanisms during the confinement periods. Accumulation processes may be easily related to lagoon age, since old lagoons have higher content of nutrients and organic matter, suggesting that these lagoons progressively accumulate nutrients during the successive confinement events. This is the usual case for most Mediterranean salt marshes without an artificially manipulated water regime. The zooplankton community in La Pletera integrates the effects of both the hydrological regime and the lagoon age since the former determines the temporal pattern of the main zooplankton species and the latter explains differences in composition and structure between old and new lagoons.     

Quality assessment of GPS reprocessed terrestrial reference frame

The International GNSS Service (IGS) contributes to the construction of the International Terrestrial Reference Frame (ITRF) by submitting time series of station positions and Earth Rotation Parameters (ERP). For the first time, its submission to the ITRF2008 construction is based on a combination of entirely reprocessed GPS solutions delivered by 11 Analysis Centers (ACs). We analyze the IGS submission and four of the individual AC contributions in terms of the GNSS frame origin and scale, station position repeatability and time series seasonal variations. We show here that the GPS Terrestrial Reference Frame (TRF) origin is consistent with Satellite laser Ranging (SLR) at the centimeter level with a drift lower than 1 mm/year. Although the scale drift compared to Very Long baseline Interferometry (VLBI) and SLR mean scale is smaller than 0.4 mm/year, we think that it would be premature to use that information in the ITRF scale definition due to its strong dependence on the GPS satellite and ground antenna phase center variations. The new position time series also show a better repeatability compared to past IGS combined products and their annual variations are shown to be more consistent with loading models. The comparison of GPS station positions and velocities to those of VLBI via local ties in co-located sites demonstrates that the IGS reprocessed solution submitted to the ITRF2008 is more reliable and precise than any of the past submissions. However, we show that some of the remaining inconsistencies between GPS and VLBI positioning may be caused by uncalibrated GNSS radomes.     

Using Two-Point Set Statistics to Estimate the Diameter Distribution in Boolean Models with Circular Grains

This paper deals with the counting and measuring of overlapping circular objects in binary images, a problem that arises in the mineral processing industry when estimating the distribution of bubble diameters in flotation cells in order to predict the geometallurgical performance of the flotation process. By viewing the images as realizations of a stationary planar Boolean model with circular grains and by using two-point set statistics, it is possible to jointly estimate the average number of objects per unit area and their diameter distribution. When applied to real and simulated images, this method is shown to give accurate estimates, to be robust to the presence of noise contaminating the images and of moderate drifts in the number of objects per unit area, and to speed up the processing time with respect to currently used methods. Combined, these properties serve to make the method suitable for an on-line monitoring of the flotation process.     

Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.     

INTEGRAL/XMM views on the MeV source GRO J1411-64

The COMPTEL unidentified source GRO J 1411-64 was observed by INTEGRAL and XMM-Newton in 2005. The Circinus Galaxy is the only source detected within the 4σ location error of GRO J1411-64, but in here excluded as the possible counterpart. At soft X-rays, 22 reliable and statistically significant sources (likelihood >10) were extracted and analyzed from XMM-Newton data. Only one of these sources, XMMU J141255.6-635932, is spectrally compatible with GRO J1411-64 although the fact the soft X-ray observations do not cover the full extent of the COMPTEL source position uncertainty make an association hard to quantify and thus risky. At the best location of the source, detections at hard X-rays show only upper limits, which, together with MeV results obtained by COMPTEL suggest the existence of a peak in power output located somewhere between 300–700 keV for the so-called low state. Such a spectrum resembles those in blazars or microquasars, and might suggest at work by the models accordingly. However, an analysis using a microquasar model consisting on a magnetized conical jet filled with relativistic electrons, shows that it is hard to comply with all observational constrains. This fact and the non-detection at hard X-rays introduce an a-posteriori question mark upon the physical reality of this source, what is discussed here.     

An improved spectral turning-bands algorithm for simulating stationary vector Gaussian random fields

We propose a spectral turning-bands approach for the simulation of second-order stationary vector Gaussian random fields. The approach improves existing spectral methods through coupling with importance sampling techniques. A notable insight is that one can simulate any vector random field whose direct and cross-covariance functions are continuous and absolutely integrable, provided that one knows the analytical expression of their spectral densities, without the need for these spectral densities to have a bounded support. The simulation algorithm is computationally faster than circulant-embedding techniques, lends itself to parallel computing and has a low memory storage requirement. Numerical examples with varied spatial correlation structures are presented to demonstrate the accuracy and versatility of the proposal.     

Estimates of groundwater discharge to a coastal wetland using multiple techniques: Taylor Slough, Everglades National Park, USA

Quantifying water exchange between a coastal wetland and the underlying groundwater is important for closing water, energy and chemical budgets. The coastal wetlands of the Florida Everglades (USA) are at the forefront of a large hydrologic restoration project, and understanding of groundwater/surface-water interactions is needed to comprehend the effects of the project. Four independent techniques were used to identify water exchange at varying spatial and temporal scales in Taylor Slough, Everglades National Park. The techniques included a water-budget study and measurements of hydraulic head gradients, geochemical tracers, and temperature. During the 18-month study, the four methods converged as to the timing of groundwater discharge, typically between June and September, contemporaneous with the wet season and increasing surface-water levels. These results were unexpected, as groundwater discharge was predicted to be greatest when surface-water levels were low, typically during the dry season. Either a time lag of 1?C5?months in the response of groundwater discharge to low surface-water levels or precipitation-induced groundwater discharge may explain the results. Groundwater discharge was a significant contributor (27?%) to the surface water in Taylor Slough with greater rates of discharge observed towards the coastline in response to seawater intrusion.