Derivatives of the Gravity Potential with Respect to Rectangular Coordinates

Starting from a complex operator of derivation, we give expressions for derivatives of arbitrary order of the gravity potential with respect to rectangular coordinates. These expressions have a form similar to the original potential expanded in spherical harmonics and are free of singularity at the poles. Computing sets of numerical coefficients once for all, we can compute the derivatives with a very limited work: the same functions are used to compute all derivatives by means of a unique parametrized formula. This is very comfortable for further algebraic manipulations. Numerical tests prove the accuracy and the efficiency of the algorithm derived from our formula to compute the gravity acceleration vector and the gravity gradient tensor.     

Heat blanketing envelopes and thermal radiation of strongly magnetized neutron stars

Strong (B?10⁹ G) and superstrong (B?10¹⁴ G) magnetic fields profoundly affect many thermodynamic and kinetic characteristics of dense plasmas in neutron star envelopes. In particular, they produce strongly anisotropic thermal conductivity in the neutron star crust and modify the equation of state and radiative opacities in the atmosphere, which are major ingredients of the cooling theory and spectral atmosphere models. As a result, both the radiation spectrum and the thermal luminosity of a neutron star can be affected by the magnetic field. We briefly review these effects and demonstrate the influence of magnetic field strength on the thermal structure of an isolated neutron star, putting emphasis on the differences brought about by the superstrong fields and high temperatures of magnetars. For the latter objects, it is important to take proper account of a combined effect of the magnetic field on thermal conduction and neutrino emission at densities ρ?10¹⁰ g?cm^?3. We show that the neutrino emission puts a B-dependent upper limit on the effective surface temperature of a cooling neutron star.     

A rare reflection effect eclipsing sdB+dM binary: 2M 1533+3759

In this paper, we report a rare reflection effect eclipsing sdB+dM binary, 2M?1533+3759. It is the seventh eclipsing sdB+dM binary that has been discovered to date. This system has an orbital period of 0.16177042 day and a velocity semi-amplitude of 71.1 km?s^?1. Using a grid of zero-metallicity NLTE model atmospheres, we derived T _eff=29250 K, log?g=5.58 and [He/H]=?2.37 from spectra taken near the reflection effection minimum. Lightcurve modeling resulted in a system mass ratio of 0.301 and an orbital inclination angle of 86.6°. The derived primary mass for 2M?1533+3759, 0.376±0.055 M _⊙, is significantly lower than the canonical mass (0.48 M _⊙) found for most previously investigated sdB stars. This implies an initial progenitor mass >1.8 M _⊙, at least a main sequence A star and perhaps even one massive enough to undergo non-degenerate helium ignition.     

Evolving metasomatic agent in the Northern Andean subduction zone, deduced from magma composition of the long-lived Pichincha volcanic complex (Ecuador)

Geochemical studies of long-lived volcanic complexes are crucial for the understanding of the nature and composition of the subduction component of arc magmatism. The Pichincha Volcanic Complex (Northern Andean Volcanic Zone) consists of: (1) an old, highly eroded edifice, the Rucu Pichincha, whose lavas are mostly andesites, erupted from 1,100 to 150 ka; and (2) a younger, essentially dacitic, Guagua Pichincha composite edifice, with three main construction phases (Basal Guagua Pichincha, Toaza, and Cristal) which developed over the last 60 ka. This structural evolution was accompanied by a progressive increase of most incompatible trace element abundances and ratios, as well as by a sharp decrease of fluid-mobile to fluid-immobile element ratios. Geochemical data indicate that fractional crystallization of an amphibole-rich cumulate may account for the evolution from the Guagua Pichincha andesites to dacites. However, in order to explain the transition between the Rucu Pichincha andesites and Guagua Pichincha dacites, the mineralogical and geochemical data indicate the predominance of magma mixing processes between a mafic, trace-element depleted, mantle-derived end-member, and a siliceous, trace-element enriched, adakitic end-member. The systematic variation of trace element abundances and ratios in primitive samples leads us to propose that the Rucu Pichincha magmas came from a hydrous-fluid metasomatized mantle wedge, whereas Guagua Pichincha magmas are related to partial melting of a siliceous-melt metasomatized mantle. This temporal evolution implies a change from dehydration to partial melting of the slab, which may be associated with an increase in the geothermal gradient along the slab due to the presence of the subducted Carnegie Ridge at the subduction system. This work emphasizes the importance of studying arc-magma systems over long periods of time (of at least 1 million of years), in order to evaluate the potential variations of the slab contribution into the mantle source of the arc magmatism.     

The recording of floods and earthquakes in Lake Chichój,Guatemala during the twentieth century

Laguna Chichój (Lake Chichój) is the only deep permanent lake in the central highlands of Guatemala. The lake is located in the boundary zone between the North American and Caribbean plates. The lake has been struck by devastating earthquakes and tropical cyclones in historical times. We investigated the imprint of twentieth century extreme events on the sedimentary record of this tropical lake using a bathymetric survey of the lake, coring the lake floor, and providing a chronology of sediment accumulation. The lake occupies a series of circular depressions likely formed by the rapid dissolution of a buried body of gypsum. ²¹⁰Pb and ¹³⁷Cs inventories and varve counting indicate high rates of sedimentation (1–2 cm year^?1). The annually layered sediment is interrupted by turbidites of two types: a darker-colored turbidite, enriched in lake-derived biogenic constituents, and interpreted as a seismite, and a lighter-colored type, enriched in catchment-derived constituents, interpreted as a flood layer. Comparison of our ¹³⁷Cs-determined layer ages with a catalog of twentieth century earthquakes shows that an earthquake on the Motagua fault in 1976 generated a conspicuous darker-colored turbidite and slumped deposits in separate parts of the lake. The entire earthquake inventory further reveals that mass movements in the lake are triggered at Modified Mercalli Intensities higher than V. Tropical cyclonic depressions known to have affected the lake area had limited effect on the lake, including Hurricane Mitch in 1998. One storm however produced a significantly thicker flood layer in the 1940s. This storm is reportedly the only event to have generated widespread slope failures in the lake catchment. It is thus inferred that abundant landsliding provided large amounts of concentrated sediment to the lake, through hyperpycnal flows.     

Climatological study of ionospheric irregularities over the European mid-latitude sector with GPS

High-frequency variability of the ionosphere, or irregularities, constitutes the main threat for real-time precise positioning techniques based on Global Navigation Satellite Systems (GNSS) measurements. Indeed, during periods of enhanced ionospheric variability, GNSS users in the field—who cannot verify the integrity of their measurements—will experience positioning errors that can reach several decimeters, while the nominal accuracy of the technique is cm-level. In the frame of this paper, a climatological analysis of irregularities over the European mid-latitude region is presented. Based on a 10 years GPS dataset over Belgium, the work analyzes the occurrence rate (as a function of the solar cycle, season and local time) as well as the amplitude of ionospheric irregularities observed at a single GPS station. The study covers irregularities either due to space weather events (solar origin) or of terrestrial origin. If space weather irregularities are responsible for the largest effects in terms of ionospheric error, their occurrence rate highly depends on solar activity. Indeed, the occurrence rate of ionospheric irregularities is about 9 % during solar maximum, whereas it drops to about 0 % during medium or low solar activity periods. Medium-scale ionospheric disturbances (MSTIDs) occurring during daytime in autumn/winter are the most recurrent pattern of the time series, with yearly proportions slightly varying with the solar cycle and an amplitude of about 10 % of the TEC background. Another recurrent irregularity type, though less frequent than MSTIDs, is the noise-like variability in TEC observed during summer nighttime, under quiet geomagnetic conditions. These summer nighttime irregularities exhibit amplitudes ranging between 8 and 15 % of the TEC background.     

Effect of the processing methodology on satellite altimetry-based global mean sea level rise over the Jason-1 operating period

Determining how the global mean sea level (GMSL) evolves with time is of primary importance to understand one of the main consequences of global warming and its potential impact on populations living near coasts or in low-lying islands. Five groups are routinely providing satellite altimetry-based estimates of the GMSL over the altimetry era (since late 1992). Because each group developed its own approach to compute the GMSL time series, this leads to some differences in the GMSL interannual variability and linear trend. While over the whole high-precision altimetry time span (1993–2012), good agreement is noticed for the computed GMSL linear trend (of $3.1\pm 0.4$  mm/year), on shorter time spans (e.g., ${<}10~\hbox {years}$ ), trend differences are significantly larger than the 0.4 mm/year uncertainty. Here we investigate the sources of the trend differences, focusing on the averaging methods used to generate the GMSL. For that purpose, we consider outputs from two different groups: the Colorado University (CU) and Archiving, Validation and Interpretation of Satellite Oceanographic Data (AVISO) because associated processing of each group is largely representative of all other groups. For this investigation, we use the high-resolution MERCATOR ocean circulation model with data assimilation (version Glorys2-v1) and compute synthetic sea surface height (SSH) data by interpolating the model grids at the time and location of “true” along-track satellite altimetry measurements, focusing on the Jason-1 operating period (i.e., 2002–2009). These synthetic SSH data are then treated as “real” altimetry measurements, allowing us to test the different averaging methods used by the two processing groups for computing the GMSL: (1) averaging along-track altimetry data (as done by CU) or (2) gridding the along-track data into $2^{\circ }\times 2^{\circ }$ meshes and then geographical averaging of the gridded data (as done by AVISO). We also investigate the effect of considering or not SSH data at shallow depths $({<}120~\hbox {m})$ as well as the editing procedure. We find that the main difference comes from the averaging method with significant differences depending on latitude. In the tropics, the $2^{\circ }\times 2^{\circ }$ gridding method used by AVISO overestimates by 11 % the GMSL trend. At high latitudes (above $60^{\circ }\hbox {N}/\hbox {S}$ ), both methods underestimate the GMSL trend. Our calculation shows that the CU method (along-track averaging) and AVISO gridding process underestimate the trend in high latitudes of the northern hemisphere by 0.9 and 1.2 mm/year, respectively. While we were able to attribute the AVISO trend overestimation in the tropics to grid cells with too few data, the cause of underestimation at high latitudes remains unclear and needs further investigation.     

Calcium,Na, K and Mg Concentrations in Seawater by Inductively Coupled Plasma‐Atomic Emission Spectrometry: Applications to IAPSO Seawater Reference Material,Hydrothermal Fluids and Synthetic Seawater Solutions

We report a measurement procedure to determine simultaneously the major cation concentrations (Na, Ca, K and Mg) of seawater‐derived solutions by inductively coupled plasma‐atomic emission spectrometry. The best results were obtained when the IAPSO (‘standard’) seawater reference material was diluted by thirty times with Milli‐Q^® water. We obtained an average reference value r_K (the ratio of the mass fraction of potassium to that of chlorine, i.e., (g kg^?1)/(g kg^?1)) for IAPSO seawater of 0.0205 ± 0.0006 (2.9% RSD), not significantly different from 0.0206 ± 0.0005 (2.4% RSD) for seawater composition reported in the literature. The measured Na, Ca and Mg concentrations correspond to r_Na, r_Ca and r_Mg values of 0.5406 ± 0.0026 (0.5% RSD), 0.02192 ± 0.00048 (2.2% RSD) and 0.06830 ± 0.00047 (0.7% RSD), respectively, in line with previous values measured by wet‐chemistry and atomic absorption spectrophotometry or wet‐chemical titration. Our measurement procedure was used successfully on synthetic seawater solutions and high‐temperature hydrothermal fluids.