Submarine groundwater discharge of nutrients to the ocean along a coastal lagoon barrier, Southern Brazil

The Patos–Mirim Lagoon system along the southern coast of Brazil is linked to the coastal ocean by a narrow mouth and by groundwater transport through a Holocene barrier. Although other groundwater systems are apparently active in this region, the hydraulic head of the lagoon, the largest in South America, drives groundwater transport to the coast. Water levels in wells placed in the barrier respond to changing water level in the lagoon. The wells also provide a measure of the nutrient concentrations of groundwater flowing toward the ocean. Additionally, temporary well points were used to obtain nutrient samples in groundwater on the beach face of the barrier. These samples revealed a subterranean freshwater–seawater mixing zone over a ca. 240 km shoreline. Previously published results of radium isotopic analyses of groundwater and of surface water from cross-shelf transects were used to estimate a water flux of submarine groundwater discharge (SGD) to nearshore surface waters of 8.5 × 10⁷ m³/day. Using this SGD and the nutrient concentrations in different compartments, nutrient fluxes between groundwater and surface water were estimated. Fluxes were computed using both average and median reservoir (i.e. groundwater and surface water) nutrient concentrations. The SGD total dissolved inorganic nitrogen, phosphate and silicate fluxes (2.42, 0.52, 5.92 × 10⁶ mol day^− 1, respectively) may represent as much as 55% (total N) to 10% (Si) of the nutrient fluxes to the adjacent shelf environment. Assuming nitrogen limitation, SGD may be capable of supporting a production rate of ca. 3000 g C m² year^− 1in the nearshore surf zone in this region.     

Simulation of Synthetic Exploration and Geometallurgical Database of Porphyry Copper Deposits for Educational Purposes

Natural Resources Research - The access to real geometallurgical data is very limited in practice, making it difficult for practitioners, researchers and students to test methods, models and...     

A locally conservative Eulerian–Lagrangian numerical method and its application to nonlinear transport in porous media

Eulerian-Lagrangian and Modified Method of Characteristics (MMOC) procedures provide computationally efficient techniques for approximating the solutions of transport-dominated diffusive systems. The original MMOC fails to preserve certain integral identities satisfied by the solution of the differential system; the recently introduced variant, called the MMOCAA, preserves the global form of the identity associated with conservation of mass in petroleum reservoir simulations, but it does not preserve a localized form of this identity. Here, we introduce an Eulerian-Lagrangian method related to these MMOC procedures that guarantees conservation of mass locally for the problem of two-phase, immiscible, incompressible flow in porous media. The computational efficiencies of the older procedures are maintained. Both the original MMOC and the MMOCAA procedures for this problem are derived from a nondivergence form of the saturation equation; the new method is based on the divergence form of the equation. A reasonably extensive set of computational experiments are presented to validate the new method and to show that it produces a more detailed picture of the local behavior in waterflooding a fractally heterogeneous medium. A brief discussion of the application of the new method to miscible flow in porous media is included. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of Explicit Urban-Canopy Representation on Local Circulations Above a Tropical Mega-City

The Advanced Regional Prediction System (ARPS) is coupled with the tropical town energy budget (tTEB) scheme to analyze the effects of the urban canopy circulation over the metropolitan area of São Paulo and its interactions with the sea breeze and mountain-valley circulation in the eastern state of São Paulo, Brazil. Two experiments are carried out for the typical sea-breeze event occurring on 22 August 2014 under weak synoptic forcing and clear-sky conditions: (a) a control run with the default semi-desert surface parametrization and; (b) a tTEB run for the urban canopy of São Paulo. A realistic land-use database over the south-eastern domain of Brazil is used in the downscaling simulation to a horizontal grid resolution of 3 km. Our results indicate that ARPS effectively simulates features of the nighttime and early morning land-breeze circulation, which is affected by the surrounding hills and the nocturnal heat island of São Paulo. By early afternoon, the south-eastern sea-breeze circulation moves inland perpendicular to the upslope of the Serra do Mar scarp, which generates a line of moisture convergence and updrafts further inland. Later, the convergence line reaches São Paulo and interacts with the circulation arising from the urban heat island (UHI), which increases the moisture convergence and strength of updrafts. The surface energy balance indicates that the UHI is caused by large sensible heat storage within the urban canopy during the day, which is later released in the afternoon and at night. The simulations are verified with available radiosonde and surface weather station data, land-surface-temperature estimates from the moderate resolution imaging spectroradiometer, as well as the National Center for Atmospheric Research reanalysis databases. The three-dimensional geometry of the urban canyons within the tTEB scheme consistently improves the thermodynamically-induced circulation over São Paulo.     

In situ Re–Os isotopic analysis of platinum-group minerals from the Mayarí-Cristal ophiolitic massif (Mayarí-Baracoa Ophiolitic Belt,eastern Cuba): implications for the origin of Os-isotope heterogeneities in podiform chromitites

Chromitite pods in the Mayarí-Cristal ophiolitic massif (eastern Cuba) were formed in the Late Cretaceous when island arc tholeiites and MORB-like back-arc basin basalts reacted with residual mantle peridotites and generated chromite-rich bodies enclosed in dunite envelopes. Platinum-group minerals (PGM) in the podiform chromitites exhibit important Os-isotope heterogeneities at the kilometric, hand sample and thin section scales. ¹⁸⁷Os/¹⁸⁸Os calculated at the time of chromitite crystallization (~90 Ma) ranges between 0.1185 and 0.1295 (γOs = −7.1 to +1.6, relative to enstatite chondrite), and all but one PGM have subchondritic ¹⁸⁷Os/¹⁸⁸Os. Grains in a single hand sample have initial ¹⁸⁷Os/¹⁸⁸Os that spans from 0.1185 to 0.1274, and in one thin section it varies between 0.1185 and 0.1232 in two PGM included in chromite which are only several millimeters apart. As the Os budget of a single micrometric grain derives from a mantle region that was at least several m³ in size, the variable Os isotopic composition of PGM in the Mayarí-Cristal chromitites probably reflects the heterogeneity of their mantle sources on the 10–100 m scale. Our results show that this heterogeneity was not erased by pooling and mingling of individual melt batches during chromitite crystallization but was transferred to the ore deposits on mineral scale. The distribution of the Os model ages calculated for PGM shows four main peaks, at ~100, 500, 750 and 1,000 Ma. These variable Os model ages reflect the presence of different depleted domains in the oceanic (Pacific-related) upper mantle of the Greater Antilles paleo-subduction zone. The concordance between the age of crystallization of the Mayarí-Cristal chromitites and the most recent peak of the Os model age distribution in PGM supports that Os in several grains was derived from fertile domains of the upper mantle, whose bulk Os isotopic composition is best approximated by that of enstatite chondrites; on the other hand, most PGM are crystallized by melts that tapped highly refractory mantle sources.     

Periodic orbits for a class of galactic potentials

In this work, applying general results from averaging theory, we find periodic orbits for a class of Hamiltonian systems H whose potential models the motion of elliptic galaxies.     

Stochastic computational modelling of highly heterogeneous poroelastic media with long‐range correlations

The compaction of highly heterogeneous poroelastic reservoirs with the geology characterized by long‐range correlations displaying fractal character is investigated within the framework of the stochastic computational modelling. The influence of reservoir heterogeneity upon the magnitude of the stresses induced in the porous matrix during fluid withdrawal and rock consolidation is analysed by performing ensemble averages over realizations of a log‐normally distributed stationary random hydraulic conductivity field. Considering the statistical distribution of this parameter characterized by a coefficient of variation governing the magnitude of heterogeneity and a correlation function which decays with a power‐law scaling behaviour we show that the combination of these two effects result in an increase in the magnitude of effective stresses of the rock during reservoir depletion. Further, within the framework of a perturbation analysis we show that the randomness in the hydraulic conductivity gives rise to non‐linear corrections in the upscaled poroelastic equations. These corrections are illustrated by a self‐consistent recursive hierarchy of solutions of the stochastic poroelastic equations parametrized by a scale parameter representing the fluctuating log‐conductivity standard deviation. A classical example of land subsidence caused by fluid extraction of a weak reservoir is numerically simulated by performing Monte Carlo simulations in conjunction with finite elements discretizations of the poroelastic equations associated with an ensemble of geologies. Numerical results illustrate the effects of the spatial variability and fractal character of the permeability distribution upon the evolution of the Mohr–Coulomb function of the rock. Copyright © 2004 John Wiley & Sons, Ltd.     

Granada oscillation code (GraCo)

Granada oscillation code (GraCo) is a software constructed to compute adiabatic and non-adiabatic oscillation eigenfunctions and eigenvalues. The adiabatic version gives the standard numerical resolution, and also the Richardson extrapolation, different sets of eigenfunctions, different outer mechanical boundary conditions or different integration variables. The non-adiabatic version can include the atmosphere-pulsation interaction. The code has been used for intensive studies of δ Scuti, γ Doradus, β Ceph., SdO and, SdB stars. The non adiabatic observables “phase-lag” (the phase between the effective temperature variations and the radial displacement) and ${\delta T_{\mathrm{eff}}\over T_{\mathrm{eff}}}$ (relative surface temperature variation) can help to the modal identification. These quantities together with the energy balance (“growth rate”) provide useful additional information to the adiabatic resolution (eigenfrequencies and eigenfunctions).     

Evolution of hard gluons in a medium with q^(t)∼tn and the thermalization problem

We study the energy distribution of hard gluons traversing a dense quark-gluon plasma by comparing various transverse momentum broadening rates , using a probabilistic perturbative approach. These results were applied to address the thermalization problem in heavy ion collisions. Within the weak coupling model, thermalization follows a “bottom-up” process: early-formed high-energy partons emit low-energy gluons, leading to their equilibrium formation, creating a thermal bath that facilitates equilibrium in the high-energy sector. Under this scenario, we model the time dependencies of $$ as a power-law $$, and assess the impact of $$ on the distribution of hard gluons passing through the medium.