On Possible Types of Magnetospheres of Hot Jupiters

As a rule, the orbits of “hot Jupiter” exoplanets are located close to the Alfven point of the stellar wind of the host star. Many hot Jupiters could be in the sub-Alfven zone, where the magnetic pressure of the stellar wind exceeds the dynamical pressure. Therefore, the magnetic field in the wind should play an extremely important role in the process of stellar wind flowing around the atmosphere of a hot Jupiter. This must be taken into account when constructing theoretical models and interpreting observational data. Analyses show that many typical hot Jupiters should have shockless induced magnetospheres, which have no analogs in the solar system. Such magnetospheres are characterized first and foremost by the fact that there is no bow shock, and the magnetic barrier (ionopause) is formed by induced currents in upper layers of the ionosphere. This conclusion is confirmed here using three-dimensional numerical simulations of the flow of the stellar wind from the host star around the hot Jupiter HD 209458b, taking into account both the intrinsic magnetic field of the planet and the magnetic field in the wind.

     

Types of gaseous envelopes of “hot Jupiter” exoplanets

As a rule, the orbital velocities of “hot Jupiters,” i.e., exoplanets with masses comparable to the mass of Jupiter and orbital semi-major axes less than 0.1 AU, are supersonic relative to the stellar wind, resulting in the formation of a bow shock. Gas-dynamical modeling shows that the gaseous envelopes around hot Jupiters can belong to two classes, depending on the position of the collision point. if the collision point is inside the Roche lobe of the planet, the envelopes have the almost spherical shapes of classical atmospheres, slightly distorted by the influence of the star and interactions with the stellar-wind gas; if the collision point is located outside the Roche lobe, outflows from the vicinity of the Lagrangian points L₁ and L₂ arise, and the envelope becomes substantially asymmetrical. The latter class of objects can also be divided into two types. If the dynamical pressure of the stellar-wind gas is high enough to stop the most powerful outflow from the vicinity of the inner Lagrangian point L₁, a closed quasi-spherical envelope with a complex shape forms in the system. If the wind is unable to stop the outflow from L₁, an open aspherical envelope forms. The possible existence of atmospheres of these three types is confirmed by 3D numerical modeling. Using the typical hot Jupiter HD 209458b as an example, it is shown that all three types of atmospheres could exist within the range of estimated parameters of this planet. Since different types of envelopes have different observational manifestations, determining the type of envelope in HD 209458b could apply additional constrains on the parameters of this exoplanet.     

Possible New Envelope Types of Hot Jupiters

Astronomy Reports - Abstract—For the case of the sub-Alfvén regime of the stellar wind flyby around the planet, the flow structure in the vicinity of a hot Jupiter have been studied by...     

Pulsations in the atmospheres of hot Jupiters possessing magnetic fields

The discovery of the possible existence of huge quasi-stationary envelopes around a number of hot Jupiters (i.e., with sizes appreciably exceeding their Roche lobes) and the need to correctly take into account their properties when interpreting observational data require a careful analysis of the main physical processes influencing their atmospheres. One important factor is the possibility that the planet has a magnetic field. It was shown earlier that the presence of even a modest dipolar magnetic field of a hot Jupiter (with a magnetic moment approximately 1/10 the magnetic moment of Jupiter) influences the properties of the planetary atmosphere, in particular, leading to expansion of the range of parameters for which a giant, quasi-closed envelope can form around the planet. It was also established that the presence of a planetary magnetic field reduced the mass-loss rate from the envelope, since matter flowing out from the inner Lagrange point moves perpendicular to the field lines. Three-dimensional magnetohydrodynamical (MHD) modeling on time scales appreciably exceeding the time for the formation of the envelope show that pulsations arise in the atmospheres of hot Jupiters possessing dipolar magnetic fields, with characteristic periods ~0.27P_orb. This behavior is easy to understand physically, since even in the case of a spherical atmosphere, the continuous expansion of the ionized atmsphere of a hot Jupiter can lead to the accumulation of matter in regions bounded by closed field lines, and to the periodic rupture of the atmosphere beyond the magnetic field. In the case considered, when the system contains a giant envelope fed by a stream of matter from the inner Lagrange point, the presence of such pulsations gives rise to appreciable variations in the gas-dynamical structure of the flow. In particular, pulsations of the atmosphere lead to tearing off of part of the flow and sharp fluctuations in the size of the envelope, leading to variations in the envelope’s observational properties.     

Mass-loss rates of “hot-Jupiter” exoplanets with various types of gaseous envelopes

According to the compuations results obtained by Bisikalo et al. (2013) for the gas-dynamical effect of stellar winds on exoplanet atmospheres, three types of gaseous envelopes can form around hot Jupiters: closed, quasi-closed, and open. The type of envelope that forms depends on the position of the frontal collision point (where the dynamical pressure of the wind is equal to the pressure of the surrounding atmosphere) relative to the Roche-lobe boundaries. Closed envelopes are formed around planets whose atmospheres lie completely within their Roche lobes. If the frontal collision point is located outside the Roche lobe, the atmospheric material begins to flow out through the Lagrangian points L₁ and L₂, which can result in the formation of quasi-closed (if the dynamical pressure of the stellar wind stops the outflow through L₁) or open gaseous envelopes. The example of the typical hot Jupiter HD 209458b is considered for four sets of atmospheric parameters, to determine the mass-loss rates for the different types of envelopes arising with these parameters. The mass-loss rates based on the modeling results were estimated to be ? ≤ 10⁹ g/s for a closed atmosphere, ? ? 3 × 10⁹ g/s for a quasi-closed atmosphere, and ? ? 3 × 10¹⁰ g/s for an open atmosphere. The matter in the closed and quasi-closed atmospheres flows out mainly through L₂, and the matter in open envelopes primarily through L₁.     

Effect of Magnetic Field Diffusion on the Structure of Extended Envelopes of Hot Jupiters

Astronomy Reports - The paper presents a modification of the numerical MHD model of extended envelopes of hot Jupiters, in which the magnetic viscosity is considered. In the new version of the...     

Reduction of mass loss by the hot Jupiter WASP-12b due to its magnetic field

The influence of the dipolar magnetic field of a “hot Jupiter” with the parameters of the object WASP-12b on the mass-loss rate from its atmosphere is investigated. The results of three-dimensional gas-dynamical and magnetohydrodynamical computations show that the presence of a magnetic moment with a strength of ~0.1 the magnetic moment of Jupiter leads to appreciable variations of the matter flow structure. For example, in the case of the exoplanet WASP-12b with its specified set of atmospheric parameters, the stream from the vicinity of the Lagrange point L₁ is not stopped by the dynamical pressure of the stellar wind, and the envelope remains open. Including the effect of the magnetic field leads to a variation in this picture—the atmosphere becomes quasi-closed, with a characteristic size of order 14 planetary radii, which, in turn, substantially decreases the mass-loss rate by the exoplanet atmosphere (by~70%). This reduction of the mass-loss rate due to the influence of the magnetic fieldmakes it possible for exoplanets to form closed and quasi-closed envelopes in the presence of more strongly overflowing Roche lobes than is possible without a magnetic field.     

Effect of the Orbit Eccentricity on the Flow in the Envelope of a Hot Jupiter

Astronomy Reports - A three-dimensional numerical gas–dynamic model is presented. The aim of the model is to study the structure of the flow in the envelopes of hot Jupiters with elliptical...     

Possible scenarios for the formation of Ap/Bp stars

Possible paths for the formation of Ap/Bp stars—massive main-sequence stars with strong magnetic fields—are analyzed based on modern theories for the evolution of single and binary stars. Assuming that the strong magnetic fields of these stars are the main reason for their comparatively slow axial rotation and the observed anomalies in the chemical compositions of their atmospheres, possible origins for these high magnetic fields are considered. Analysis of several possible scenarios for the formation of these stars leads to the conclusion that their surface magnetic fields are probably generated in the convective envelopes of the precursor stars. These precursors may be young, single stars with masses 1.5–3 M _⊙ that formed at the peripheries of forming star clusters and ended their accretion at the Hayashi boundary, or alternatively close binaries whose components have convective envelopes, whose merger leads to the formation of an Ap/Bp star. Arguments are presented supporting the view that the merger of close binaries is the main channel for the formation of Ap/Bp stars, and a detailed analysis of this scenario is presented. The initial major axes of the merging binary systems must be in the range 6–12 R _⊙, and the masses of their components in the range 0.7–1.5 M _⊙. When the merging components possess developed convective envelopes and fairly strong initial magnetic fields, these can generate powerful magnetic fields “inherited” by the products of the merger—Ap/Bp stars. The reason the components of the close binaries merge is a loss of angular momentum via the magnetic stellar winds of the components.     

The influence of coronal mass ejections on the gas dynamics of the atmosphere of a “hot Jupiter” exoplanet

The results of three-dimensional numerical simulations of the gas dynamics of the atmosphere of a “hot Jupiter” exoplanet during the passage of a coronal mass ejection (CME) from the central star are presented. These computations assumed the parameters for the stellar wind and the CME to be typical of the solar values. The characteristic variations of the flow pattern are considered for quasi-closed and closed (but appreciably distorted by the gravitational influence of the star) gaseous envelopes of the exoplanet. It is shown that a typical CME is sufficient to tear off the outer part of an asymmetric envelope that is located beyond the Roche lobe and carry it away from the exoplanet. This leads to a substantial increase in the mass-loss rate from the exoplanet envelope during the passage of CMEs. The mass-loss rate grows by about a factor of 11 for a closed envelope, and by about a factor of 14 for a quasi-closed envelope. Possible evolutionary consequences of the loss of part of the atmosphere during the passage of CMEs are discussed.     

Intercombinational line profiles in the UV spectra of T Tauri stars and analysis of the accretion zone

We have analyzed for the first time profiles of the SiIII 1892 Å and CIII 1909 Å intercombinational lines in HST spectra of the stars RY Tau and RU Lup. The widths of these optically thin lines exceeded 400 km/s, ruling out formation in the stellar chromosphere. Since the intensity of the Si line exceeds that of the C line, it is unlikely that a large fraction of the observed line flux is formed in a stellar wind. The observed profiles can be reproduced in the framework of an accretion shock model if the velocity field in the accretion zone is appreciably nonaxisymmetric. In this case, the line profiles should display periodic variations, which can be used to determine the accretion zone geometry and the topology of the magnetic field near the stellar surface; corresponding formulas are presented. In addition, periodic variations of the 0.3–0.7 keV X-ray flux should be observed.     

Groundwater flow modelling of the regional aquifer of the Pampa del Tamarugal,northern Chile

The Pampa del Tamarugal Aquifer (PTA) is an important source of groundwater in northern Chile. In this study, a groundwater flow model of this aquifer is developed and calibrated for the period 1983–2004. The model reproduces the observed flow-field and the water balance components reasonably well. Five scenarios are defined to evaluate the response to different pumping situations. These scenarios show that groundwater heads will continue to decrease with the present pumping discharge rates. To account for variations in the model results due to uncertainties in average recharge rates, randomly generated recharge realizations with different levels of uncertainty are simulated. Evaporation flow rates and groundwater flowing out of the modelled area seem unaffected by the recharge uncertainty, whereas the storage terms can vary considerably. For the most intensive pumping scenario under the generated random recharge rates, it is unlikely that the cumulative discharged volume from the aquifer, at the end of the simulation period, will be larger than 12% of the estimated groundwater reserve. Fluctuations in simulated groundwater heads due to uncertainties in the average recharge values are more noticeable in certain areas. These fluctuations could explain unusual behaviour in the observed groundwater heads in these areas.

---

Resumen El Acuífero de la Pampa del Tamarugal (PTA) es una fuente importante de agua subterránea en el norte de Chile. En este estudio se desarrolla y calibra un modelo de flujo de agua subterránea para el periodo 1983–2004. El modelo reproduce razonablemente bien el campo de flujo observado y los componentes del balance hídrico. Se definen cinco escenarios para evaluar la respuesta a diferentes situaciones de bombeo. Estos escenarios muestran que con las tasas de descarga de bombeo actuales las presiones de agua subterránea continuarán en descenso. Para explicar las variaciones en los resultados del modelo debidas a incertidumbres en tasas de recarga promedio se han simulado realizaciones de recarga generadas aleatoriamente con diferentes niveles de incertidumbre. Las tasa de flujo por evaporación y el agua subterránea que fluye fuera del área modelizada parecen no haber sido afectadas por la incertidumbre en recarga mientras que los términos de almacenamiento pueden variar considerablemente. Para el escenario de bombeo más intenso, bajo las tasas de recarga generadas aleatoriamente, no es probable que el volumen descargado acumulado del acuífero al final del periodo de simulación sea mayor del 12% de la reserva de agua subterránea estimada. Las fluctuaciones en las presiones simuladas de agua subterránea debido a las incertidumbres en los valores de recarga promedio son más notorios en ciertas áreas. Las fluctuaciones podrían explicar el comportamiento inusual en las presiones de agua subterránea observadas en estas áreas.

---

Résumé L’aquifère de la Pampa del Tamarugal (PTA) représente une importante source d’eau souterraine dans le Nord du Chili. Dans cette étude, un modèle d’écoulement souterrain de cet aquifère est réalisé et calibré pour la période 1983–2004. Celui-ci reproduit le champs d’écoulement observé et les composantes du bilan d’eau raisonnablement bien. Cinq scénarios sont définis pour évaluer la réponse du système à des conditions de pompage différentes. D’après ces scénarios, les niveaux piézomètriques vont continuer à baisser sous l’effet des taux de pompage actuels. Pour tenir compte de la variabilité des résultats du modèle due aux incertitudes sur les taux de recharge moyens, des épisodes de recharge, générés aléatoirement et présentant différents niveaux d’incertitude sont simulés. Les taux d’évaporation et l’eau souterraine s’écoulant hors de la zone modélisée semblent non affectés par l’incertitude sur la recharge, les termes d’emmagasinement en revanche peuvent varier considérablement. Dans le cas du scénario avec le pompage le plus important et des taux de recharge aléatoires, il est peu probable que le volume cumulé pompé hors de l’aquifère à la fin de la simulation atteigne une valeur supérieure à 12% de la réserve d’eau souterraine estimée. Les fluctuations des niveaux piézomètriques simulés liées aux incertitudes sur les recharges moyennes sont plus perceptibles dans certaines zones. Ces fluctuations pourraient expliquer un comportement inhabituel des niveaux piézomètriques observés dans ces zones.

     

Groundwater risk intensity mapping in semi-arid regions using optical remote sensing data as an additional tool

Due to semi-arid to arid climatic conditions, Jordan has limited groundwater resources. As a result of agricultural activities and untreated wastewater, most of the groundwater in the karstic areas is microbiologically contaminated. Groundwater vulnerability, hazards, and risk intensity were mapped (scale 1:50,000) at a test site near the Jordan Rift Valley. The mapping included the use of optical remote sensing to complement conventional data in areas with poor data. LANDSAT ETM+ data, and colour and panchromatic aerial photographs at different scales were incorporated using visual image interpretation and digital image processing. The applicability of the different remote sensing data sources is discussed and recommendations for their usage are given. Information derived from digital images offers new opportunities for vulnerability and hazard assessment, particularly when related to land use, vegetation cover, urbanisation and infrastructure. The resulting maps indicate clearly the vulnerable areas and the “hot spots” of potential contamination in the test site and form an important basis for integrated groundwater management studies and the long-term planning of protective measures. The application and transferability of the European vulnerability approach (COST Action 620) to the test site in Jordan proved to be good, in general, although modifications were necessary to suit local conditions.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorised users.

---

Résumé La Jordanie ne dispose que de ressources limitées en eaux souterraines, étant donné le climat semi-aride à aride prévalant dans ce pays. La contamination microbiologique de la plupart des eaux souterraines des régions karstiques est due aux activités agricoles et l’absence de traitement des eaux usées. La vulnérabilité des eaux souterraines, et le danger et l’intensité du risque de leur contamination ont été cartographiés (à l’échelle de 1:50 000) pour un site d’étude situé près de la vallée du Rift du Jourdain. La cartographie inclut l’utilisation de données optiques de télédétection pour complémenter les sources de données conventionnelles, dans les zones où la couverture est déficiente. Des données de LANDSAT ETM+ et des photographies aériennes en couleur et panchromatiques à différentes échelles ont été incorporées en utilisant des méthodes d’interprétation visuelle et informatique de l’image. La qualité des différentes sources de données de télédétection est discutée et évaluée, et des recommandations pour leur utilisation sont suggérées. L’informations obtenue des images digitales offre de nouvelles opportunités pour l’évaluation des vulnérabilités et des dangers, surtout reliés à l’usage du sol, la couverture végétale, l’urbanisation et les diverses infrastructures. Les cartes ainsi obtenues indiquent clairement les zones vulnérables et les points critiques potentiels pour la contamination du site d’étude et forment une base importante pour l’étude de la gestion intégrée des eaux souterraines et pour la planification à long terme des mesures de protection. L’applicabilité de l’approche européenne de cartographie de la vulnérabilité et des risques (COST action 620) au site d’étude jordanien est généralement approprié, même si certaines modifications ont été nécessaires pour s’ajuster aux particularités régionales.

---

Resumen Los recursos hídricos subterráneos en Jordania son limitados debido principalmente a su clima semi árido. En la mayoría de las áreas kársticas, las aguas subterráneas están contaminadas microbiológicamente a consecuencia de la actividad agrícola y de las aguas residuales. Por ello, se han realizado mapas de peligro, vulnerabilidad e intensidad del riesgo (a escala 1/50.000) de un área piloto próxima al Rift en Jordania. Los mapas han sido elaborados con el apoyo de técnicas de teledetección en aquellas zonas donde la información era escasa. Se han utilizado imágenes del LANDSAT ETM+ a y fotografía aérea, de color y pancromática, a diferentes escalas para el procesamiento de las imágenes digitales y la fotointerpretación. En este trabajo se discute y se dan recomendaciones sobre la aplicabilidad de las diferentes imágenes de satélite en los análisis de teledetección. Los datos que se pueden derivar de las imágenes digitales son un avance para la evaluación de la vulnerabilidad y los peligros de contaminación, particularmente aquellos referidos a usos del suelos, cubierta vegetal e infraestructuras urbanas. Los mapas resultantes muestran claramente la vulnerabilidad y los peligros potenciales de contaminación en la zona de estudio, lo que genera una importante base para los estudios de gestión integrada de las aguas subterráneas y de planificación de las medidas protectoras. La aplicabilidad de la aproximación europea a los estudios de vulnerabilidad en karst (COST Action 620) en la zona de estudio (Jordania) es buena, aunque fueron necesarias algunas leves modificaciones de la metodología para ajustarla a las condiciones locales.

     

Scattering of cyclotron radiation in a moving plasma

We consider the scattering of cyclotron radiation in a plasma moving along a homogeneous magnetic field. The equation of radiation transfer in a co-moving frame is derived and two limiting cases are pointed out. In the first case of a “small” velocity gradient, the total Doppler frequency shift due to variations in the plasma velocity over the flow is much smaller than the width of the line. The second, opposite, case of a “large” velocity gradient is analogous to the Sobolev approximation in the theory of moving stellar envelopes. The solution of the transfer equation for a wind-type flow illuminated by radiation of a given intensity is obtained in the latter case, when the influence of the plasma motion on cyclotron scattering is most important. It is shown that cyclotron scattering in a moving plasma differs from the known (and qualitatively similar) problems of resonance scattering in moving stellar envelopes and cyclotron scattering in a motionless plasma permeated by an inhomogeneous magnetic field. In particular, a symmetric absorption band with residual intensity proportional to the velocity gradient appears in the spectrum of the outgoing radiation, while in these two other problems, the depth of the corresponding spectral features cannot exceed half the continuum level. Detailed qualitative analysis reveals that this difference is due to the particular form of the frequency redistribution for cyclotron scattering.     

Investigating pesticide transport in the León-Chinandega aquifer, Nicaragua

Nicaragua’s León-Chinandega aquifer has seen extensive contamination by persistent organochlorine pesticides applied over decades of intensive agricultural activity. Models of flow and transport of a 330 km² sub-region of the aquifer were developed to test conceptual models of contaminant transport, to constrain the value of certain key transport parameters, and to investigate contamination-related concerns raised in past studies. To support these models, a variety of hydrogeologic and geochemical data were collected. It was concluded that the organochlorine pesticides seen in groundwater originate in soils, and are transported to the water table through widespread preferential flow, through shortcutting around wells, or through wind-blown particles delivered to poorly protected hand-dug wells. The distribution coefficient (K _d) of these pesticides is estimated to be between 0.1 and 100 ml/g and the concentration of pesticides being delivered to the water table is estimated to be between 10² and 10⁵ ng/L. It was found that the distribution and concentration of pesticides in the aquifer would be affected by an increase in groundwater abstraction within the region.

---

Résumé L’aquifère du León-Chinandega au Nicaragua a été soumis à des contaminations étendues par des pesticides organochlorés persistants qui ont été répandus pendant des dizaines d’années d’agriculture intensive. Des modèles d’écoulement et de transport d’une sous-région de 330 km² de l’aquifère ont été construits pour tester des modèles conceptuels du transport de contaminants, pour contraindre la valeur de certains paramètres de transport clefs, et pour étudier les problèmes liés au transport de contaminations qui avaient été rencontrés lors d’études précédentes. Différents types de données hydrogéologiques et géochimiques ont été collectés pour réaliser ces modèles. La conclusion montre que les pesticides organochlorés rencontrés dans les eaux souterraines proviennent des sols et qu’ils sont transportés jusqu’à l’aquifère via de grandes zones d’écoulement préférentiel, via des raccourcis hydrauliques autour des puits, ou encore via des particules emportées par le vent et tombées dans des puits creusés à la main et mal protégés. Le coefficient de distribution (K _d) de ces pesticides est estimé entre 0.1 et 100 ml/g et la concentration de pesticides arrivant dans la nappe est estimée entre 10² et 10⁵ ng/L. L’étude a montré qu’une augmentation des prélèvements d’eau souterraine dans la région aurait un impact sur la distribution et la concentration des pesticides dans l’aquifère.

---

Resumen El acuifero León-Chinandega (Nicaragua) presenta una contaminación extensiva por la aplicación de plaguicidas organoclorados en la agricultura intensiva durante décadas. Se han desarrollado modelos de flujo y transporte de una subregión del acuífero de 330 km² para comprobar los modelos conceptuales del transporte de contaminantes, estimar los valores de ciertos parámetros clave de transporte e investigar los temas relacionados con la contaminación en estudios anteriores. Para elaborar estos modelos, se recogió una variedad de datos hidrogeológicos y geoquímicos. Se ha concluido que los plaguicidas organoclorados presentes en el agua subterránea se originan en el suelo y son transportados hacia el nivel freático a través de flujos preferentes generalizados, a través de cortes alrededor de los pozos o a través de partículas transportadas por el viento depositadas en pozos excavados a mano con poca protección. El coeficiente de distribución (K _d) de estos plaguicidas se ha estimado entre 0.1 y 100 ml/g y la concentración de plaguicidas que sale hacia el nivel freático se estima que puede estar entre 10² y 10⁵ ng/L. Se ha concluido que el incremento de la extracción de aguas subterráneas dentro de la región afectaría la distribución y concentración de plaguicidas en el acuífero.

     

Circumstellar nebulae of WR and OB stars with strong winds: Simultaneous excitation by stellar radiation and shocks

We use a two-phase model for the structure of the circumstellar nebulae of hot stars to analyze the radiative cooling of a dense, compact cloud behind the shock produced by the compression of the cloud by hot gas from the stellar wind, taking into account ionization and heating by radiation from the central star. We can distinguish three stages of the evolution of the cloud during its compression. In the first stage, relevant for the entire cloud before compression and the gas ahead of the shock front, the state of the gas is determined purely by ionization by the stellar radiation. The next stage is characterized by the simultaneous action of two gas excitation mechanisms—photoionization by the stellar radiation and shock heating. In this stage, the gas intensively radiates thermal energy received at the shock front. After radiative cooling, in the final stage, ionization and heating of the gas are again determined mainly by the star. To compute the spectrum of the cloud radiation, we solved for the propagation of a plane-parallel, homogeneous flux through the shock front in the radiation field of the hot star. The computations show that a combination of two excitation mechanisms considerably enriches the theoretical spectrum. The relative intensities of emission lines of a single cloud may resemble either those for an HII region or of a supernova remnant.     

Hydrogeology and hydrochemistry of a shallow alluvial aquifer,western Saudi Arabia

A hydrogeological and hydrochemical study was conducted on a shallow alluvial aquifer, Wadi Wajj, in western Saudi Arabia to assess the influence of protection measures on groundwater quality. The hydrochemistry was assessed up-gradient and down-gradient from potential contamination sources in the main city in dry and wet seasons prior to and after the installation of major drainage and wastewater facilities. Wadi Wajj is an unconfined aquifer where water is stored and transmitted through fractured and weathered bedrock and the overlying alluvial sediments. Natural recharge to the aquifer is about 5% of rainfall-runoff. Hydrochemistry of the aquifer shows temporal and seasonal changes as influenced by protection measures and rainfall runoff. Both groundwater and runoff showed similar chemical signature, which is mostly of chloride-sulfate-bicarbonate and sodium-calcium type. Groundwater downstream of the city, though of poorer quality than upstream, showed significant improvement after the installation of a concrete runoff tunnel and a wastewater treatment plant. Concentrations of many of the groundwater quality indicators (e.g., TDS, coliform bacteria, and nitrate) exceed US Environmental Protection Agency drinking-water standards. Heavy metal content is, however, within allowable limits by local and international standards. The chemical analyses also suggest the strong influence of stream runoff and sewage water on the groundwater quality.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé Une étude hydrogéologique et géochimique a été menée sur l’aquifère phréatique alluviale Wadi Wajj dans l’Ouest de l’Arabie Saoudite afin d’évaluer l’influence de mesures de protection sur la qualité de l’eau souterraine. L’hydrogéochimie de l’eau a été étudiée en amont et en aval de sources potentielles de contamination dans la ville principale, pendant les saisons sèches et humides, avant et après l’installation de réseaux majeurs de drainage et d’eaux usées. L’aquifère Wadi Wajj est libre, l’eau est stockée et s’écoule dans les sédiments alluviaux et dans le socle fracturé et altéré sous-jacent. La recharge naturelle de l’aquifère représente 5% des eaux de pluie et de ruissellement. L’hydrogéochimie de l’eau de l’aquifère montre que les changements saisonniers et temporaires sont influencés par les mesures de protection et par le ruissellement des eaux pluviales. L’eau souterraine et l’eau de ruissellement ont présenté le même faciès chimique, de type bicarbonatée-sulfatée-chlorée et calco-sodique. En aval de la ville, l’eau souterraine, bien que de moins bonne qualité qu’en amont, a présenté une nette amélioration après l’installation d’un système de récupération et d’une station de traitement des eaux usées. Les concentrations de plusieurs paramètres indicateurs de la qualité de l’eau (tels que la charge totale dissoute, les coliformes, et les nitrates) dépassent les normes de potabilité de l’eau de consommation de l’agence américaine de la protection de l’environnement. Les teneurs en métaux lourds en revanche n’excèdent pas les normes locales et internationales. Les analyses chimiques indiquent aussi l’influence importante de l’écoulement par ruissellement et des eaux usées sur la qualité de l’eau souterraine.

---

Resumen Un estudio hidrogeológico e hidroquímico fue hecho en un acuífero somero de Wadi Wajj, en Arabia Saudi oeste para evaluar la influecia de medidas de protección en la calidad del agua subterránea. La hidroquímica fue evaluada gradiente-arriba y gradiente-abajo de las fuentes potenciales de contaminación de la ciudad principal durante las estaciones seca y lluviosa, antes y después de la instalación de sistemas principales de drenaje y aguas servidas. Wadi Wajj es un acuífero no-confinado donde el agua es almacenada y transmitida a través de roca fracturada y meteorizada, y los sedimentos aluviales que le sobreyacen. La recarga natural del acuífero es de cerca del 5% de la precipitación-escorrentía. La hidroquímica del acuífero muestra cambios temporales y estacionales influenciados por las medidas de protección y la escorrentía de precipitación. Ambas, agua subterránea y escorrentía mostraron composición química similar, siendo mayoritariamente de tipos cloruro-sulfato-bicarbonato y sodio-calcio. El agua subterránea aguas arriba de la ciudad, aunque de calidad más pobre que aguas abajo, mostró significante mejoría después de la instalación de un tunel de concreto para escorrentía y una planta de tratamiento de aguas servidas. Las concentraciones de muchos de los indicadores de calidad de agua subterránea (e.g., STD, coliformes, y nitrato) exceden los estándares de la Agencia de Protección Ambiental USA para agua potable. El contenido de metales pesados está, sin embargo, dentro de los límites permisibles de los estándares locales e internacionales. Los análisis químicos también sugieren la fuerte influencia de la escorrentía y aguas residuales en la calidad del agua subterránea.

     

Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---

     

Multiple-method estimation of recharge rates at diverse locations in the North Carolina Coastal Plain, USA

Recharge rates determined at diverse study sites in a shallow, unconfined aquifer differed from one another depending on the analytical method used and on each method’s applicability and limitations. Total recharge was quantified with saturated-zone methods using water-table fluctuations at seven sites in North Carolina, USA and using groundwater-age dating at three of the seven sites; at two of the sites, potential recharge was quantified with an unsaturated-zone method using Darcy’s law; and at five of the sites, net recharge was quantified with a stream hydrograph-separation method using streamflow-recession curves. Historical mean net recharge was 25 to 69% of the historical total recharge rates. The large disparity is attributed to groundwater losses between recharge and discharge areas, primarily by evapotranspiration and seepage to underlying aquifers. The spatial distribution of historical mean annual total recharge did not vary between landscape units, as suggested in a previous study. Similarly, total recharge did not correlate significantly with mean annual rainfall, mean annual water table depth, or the surficial soil properties of percent clay and bulk density. Total recharge did correlate significantly with the surficial soil properties of percent sand and percent silt.

---

Résumé Les taux de recharge déterminés sur divers sites d’études dans un aquifère phréatique libre, varient les uns par rapport aux autres selon la méthode analytique utilisée et les limites et l’applicabilité de chaque méthode. La recharge totale a été quantifiée avec des méthodes en zone saturée utilisant les fluctuations des nappes sur sept sites de la Caroline du Nord, USA, et en utilisant les datations des eaux souterraines sur trois des sept sites ; sur deux des sites, la recharge potentielle a été quantifiée avec une méthode en zone non-saturée utilisant la loi de Darcy, et sur deux sites, la recharge nette a été quantifiée suivant une méthode de séparation des hydrographes de cours d’eau utilisant les courbes de récession du débit. La moyenne historique de la recharge nette est comprise entre 25 et 69% des taux historiques de recharge totale. L’importante disparité est attribuée aux pertes en eaux souterraines entre les zones de recharge et les zones d’exutoire, d’abord par évapotranspiration et ensuite par infiltration vers les aquifères sous-jacents. La distribution spatiale de la recharge historique moyenne annuelle et totale ne varie pas selon les unités paysagères, comme cela était suggéré dans une étude antérieure. De même, la recharge totale n’est pas corrélée significativement avec la moyenne annuelle des précipitations, la moyenne annuelle de la profondeur de la nappe, ou les propriétés de surface des sols que sont le pourcentage d’argile ou la densité apparente. La recharge totale est corrélée significativement avec les propriétés de surface du sol et les teneurs en sable et en silt.

---

Resumen Las tasas de recarga determinadas en diversos sitios de estudio en un acuífero somero no confinado difieren entre sí dependiendo del método analítico usado y de las limitaciones y aplicabilidad de cada método. Se cuantificó la recarga total con métodos de zona saturada usando fluctuaciones del nivel freático en siete sitios en Carolina del Norte, Estados Unidos, y usando datación de edades de agua subterránea en tres de los siete sitios; en dos de los sitios se cuantificó la recarga potencial con un método de zona no saturada usando la Ley de Darcy y en cinco de los sitios se cuantificó la recarga neta con el método de separación hidrográfica usando curvas de recesión de flujo en arroyos. La recarga neta media histórica varió entre 25 y 69% de las tasas de recarga total histórica. La diferencia tan grande se atribuye a pérdidas de agua subterránea entre áreas de descarga y recarga, principalmente por evapotranspiración y escurrimiento hacia acuíferos subyacentes. La distribución espacial de la recarga total anual media histórica no varió entre las unidades de paisaje, como ya se había sugerido en un estudio previo. De modo similar, la recarga total no guarda correlación significativa con la lluvia media anual, la profundidad media anual del nivel freático, o las propiedades superficiales del suelo, particularmente porcentaje de arcilla y densidad volumétrica. La recarga total si tuvo una correlación significativa con las propiedades superficiales del suelo tal como porcentaje de arena y porcentaje de limo.

     

Development of a GIS-based model for extrapolating mesoscale groundwater recharge estimates using integrated geospatial data sets

A geographic information systems (GIS)-based model was used to derive spatially explicit estimates of recharge using the elevation-dependent chloride mass-balance (EDCMB) approach in a 14-basin area in southwestern Nevada, USA. For each of the 14 basins, a non-linear regression equation was developed relating chloride enrichment to the elevation of spring watersheds. The ratio of chloride enrichment as a function of elevation was used in a GIS framework to transform continuous precipitation data to recharge. Spatial masks that represented two definitions of the lower limits of recharge—one definition based on alluvium/non-alluvium boundaries, the other based on both alluvium/non-alluvium boundaries and an elevation threshold—were applied to each basin. Resultant recharge estimates were then summed for each basin. When compared to summaries of previous recharge estimates for the study area, the EDCMB approach produced results that were within 14 and 3% of two studies, but were significantly greater (31%) than a third. GIS proved to be a very effective tool for combining large spatial data with widely different resolutions into an integrated data set, and also proved to be an efficient mechanism for implementing robust statistical models to estimate recharge.

---

Résumé Un modèle basé sur les systèmes d’informations géographiques (SIG) a été utilisé pour obtenir des estimations explicites de la recharge dans l’espace, en utilisant l’approche du bilan des chlorures en fonction de l’altitude (EDCMB en anglais) dans une zone de 14 bassins, dans le Sud-Ouest du Nevada aux Etats-Unis. Pour chacun des 14 bassins, une régression non-linéaire a été établie, reliant l’enrichissement des chlorures à l’altitude des bassins d’alimentation des sources. Le rapport d’enrichissement des chlorures en fonction de l’altitude a été utilisé dans une approche SIG pour transformer des données pluviométriques continues en terme de recharge. Des masques spatiaux représentant deux définitions des limites de recharge les plus basses ont été appliquées à chaque bassin ; l’une basée sur les limites alluvions/autres formations, l’autre sur ces mêmes limites ainsi que sur un seuil altitudinal. Les recharges estimées résultantes ont ensuite été additionnées pour chaque bassin. Comparée aux précédents bilans de recharge estimés dans la zone d’étude, l’approche EDCMB conduit à des résultats se situant entre 14 et 3% de ceux de deux autres études et largement supérieurs (31%) à ceux d’une troisième étude. Le SIG s’est révélé être un outil très efficace pour combiner dans une base de données intégrée, des données spatiales étendues avec des résolutions largement différentes, ainsi que pour implémenter des modèles statistiques conséquents permettant d’estimer la recharge.

---

Resumen Se utilizó un modelo basado en un Sistema de Información Geográfico (SIG) para derivar espacialmente estimados explícitos de recarga usando el enfoque de balance de masa de cloruro dependiente de elevación (EDCMB) en un área de 14 cuencas en el suroeste de Nevada, Estados Unidos. Para cada una de las 14 cuencas se desarrolló una ecuación de regresión no linear que relaciona el enriquecimiento de cloruro con la elevación de las cuencas con manantiales. La relación de enriquecimiento de cloruro en función de la elevación se usó en un marco SIG para transformar datos continuos de precipitación a recarga. Se aplicaron máscaras espaciales a cada cuenca las cuales representaron dos definiciones de los límites inferiores de recarga, una definición basada en límites de aluvión/sin aluvión, el otro basado en ambos límites aluvión/sin aluvión y una elevación normal de fondo. Los estimados resultantes de recarga fueron sumados para cada cuenca. Al compararse con estimados de recarga previos para el área de estudio se encontró que el enfoque EDCMB produjo resultados que estaban entre el 14 y 3% de dos estudios pero fueron significativamente mayores (31%) que un tercer estudio. El SIG ha demostrado ser una herramienta muy efectiva para combinar datos espaciales grandes con resoluciones muy distintas en un grupo integrado de datos y también ha mostrado ser un mecanismo eficiente para implementar modelos estadísticos robustos para estimar recarga.