Precise orbit determination for the FORMOSAT-3/COSMIC satellite mission using GPS

The joint Taiwan–US mission FORMOSAT-3/ COSMIC (COSMIC) was launched on April 17, 2006. Each of the six satellites is equipped with two POD antennas. The orbits of the six satellites are determined from GPS data using zero-difference carrier-phase measurements by the reduced dynamic and kinematic methods. The effects of satellite center of mass (COM) variation, satellite attitude, GPS antenna phase center variation (PCV), and cable delay difference on the COSMIC orbit determination are studied. Nominal attitudes estimated from satellite state vectors deliver a better orbit accuracy when compared to observed attitude. Numerical tests show that the COSMIC COM must be precisely calibrated in order not to corrupt orbit determination. Based on the analyses of the 5 and 6-h orbit overlaps of two 30-h arcs, orbit accuracies from the reduced dynamic and kinematic solutions are nearly identical and are at the 2–3 cm level. The mean RMS difference between the orbits from this paper and those from UCAR (near real-time) and WHU (post-processed) is about 10 cm, which is largely due to different uses of GPS ephemerides, high-rate GPS clocks and force models. The kinematic orbits of COSMIC are expected to be used for recovery of temporal variations in the gravity field.     

Factors controlling seasonal groundwater and solute flux from snow‐dominated basins

Critical zone influences on hydrologic partitioning, subsurface flow paths and reactions along these flow paths dictate the timing and magnitude of groundwater and solute flux to streams. To isolate first‐order controls on seasonal streamflow generation within highly heterogeneous, snow‐dominated basins of the Colorado River, we employ a multivariate statistical approach of end‐member mixing analysis using a suite of daily chemical and isotopic observations. Mixing models are developed across 11 nested basins (0.4 to 85 km²) spanning a gradient of climatological, physical, and geological characteristics. Hydrograph separation using rain, snow, and groundwater as end‐members indicates that seasonal contributions of groundwater to streams is significant. Mean annual groundwater flux ranges from 12% to 33% whereas maximum groundwater contributions of 17% to 50% occur during baseflow. The direct relationship between snow water equivalent and groundwater flux to streams is scale dependent with a trend toward self‐similarity when basins exceed 5.5 km². We find groundwater recharge increases in basins of high relief and within the upper subalpine where maximum snow accumulation is coincident with reduced conifer cover and lower canopy densities. The mixing model developed for the furthest downstream site did not transfer to upstream basins. The resulting error in predicted stream concentrations points toward weathering reactions as a function of source rock and seasonal shifts in flow path. Additionally, the potential for microbial sulfate reduction in floodplain sediments along a low‐gradient, meandering portion of the river is sufficient to modify hillslope contributions and alter mixing ratios in the analysis. Soil flushing in response to snowmelt is not included as an end‐member but is identified as an important mechanism for release of solutes from these mountainous watersheds. End‐member mixing analysis used in combination with high‐frequency observations reveals important aspects of catchment hydrodynamics across scale.     

Proper orthogonal decomposition reduced model for mass transport in heterogenous media

Numerical models with fine discretization normally demand large computational time and space, which lead to computational burden for state estimations or model parameter inversion calculation. This article presented a reduced implicit finite difference scheme that based on proper orthogonal decomposition (POD) for two-dimensional transient mass transport in heterogeneous media. The reduction of the original full model was achieved by projecting the high-dimension full model to a low-dimension space created by POD bases, and the bases are derived from the snapshots generated from the model solutions of the forward simulations. The POD bases were extracted from the ensemble of snapshots by singular value decomposition. The dimension of the Jacobian matrix was then reduced after Galerkin projection. Thus, the reduced model can accurately reproduce and predict the original model’s transport process with significantly decreased computational time. This scheme is practicable with easy implementation of the partial differential equations. The POD method is illustrated and validated through synthetic cases with various heterogeneous permeability field scenarios. The accuracy and efficiency of the reduced model are determined by the optimal selection of the snapshots and POD bases.     

Interplanetary dust detected by the Cassini CDA Chemical Analyser

During its cruise phase, prior to encountering Jupiter, the Cosmic Dust Analyser (CDA) onboard the Cassini spacecraft returned time of flight mass spectra (TOF MS) of two interplanetary dust particles. Both particles were found to be iron-rich, with possible traces of hydrogen, carbon, nickel, chromium, manganese, titanium, vanadium and minor silicates. Carbon, hydrogen, oxygen and potassium are also present as possible contaminants of the impact target of CDA. Silicates and magnesium do not feature predominantly in the spectra; this is surprising considering the expected dominance of silicate-rich minerals in interplanetary dust particles. The particle masses are and . The corresponding radii ranges for the particles, assuming densities from 7874-2500 kg m⁻³ are 0.7-4 μm and 2.6-6.8 μm, respectively. With the same density assumptions the β values (ratio of radiation pressure to gravitational force) are estimated as 0.027-0.21 and 0.016-0.06 respectively, allowing possible orbits to be calculated. The resulting orbits are bound and prograde with semi-major axes, eccentricities and inclinations in the region of 0.3-1.26 AU, 0.4-1.0 and 0-60° for the first particle and 0.8-2.5 AU, 0.2-0.9 and 0-30° for the second. The more probable orbits within these ranges indicate that the first particle is in an Aten-like orbit, whilst the second particle is in an Apollo-like orbit, despite both grains having very similar, predominantly metallic compositions. Other possible orbital solutions for both particles encompass orbits which more closely resemble those of Jupiter-family comets.     

Lapse Rate Adjustments of Gridded Surface Temperature Normals in an Area of Complex Terrain: Atmospheric Reanalysis versus Statistical Up-Sampling

The applicability of elevation-regression based interpolation methods for long-term temperature normals, for example the Parameter-elevation Regressions on Independent Slopes Model (PRISM), becomes increasingly limited in data sparse, complex terrain such as that found in mountainous British Columbia (BC), Canada. Recent methods to improve both the resolution and accuracy of interpolation models have focused on the development of “up-sampling” algorithms based on local lapse rate adjustments to the original interpolated surfaces. Lapse rates can be derived from statistical models (e.g., elevation-based polynomial regression equations) or dynamical models (e.g., vertical temperature profiles from numerical weather prediction (NWP) models). This study compares a widely used statistical up-sampling algorithm, ClimateBC, with two NWP reanalysis products, the National Centers for Environmental Prediction/National Corporation for Atmospheric Research, Reanalysis 1 (NCEP1) and the more modern European Centre for Medium-range Weather Forecasts (ECMWF) Reanalysis Interim (ERA-Interim). Thirty-year climate normals for maximum and minimum temperatures were calculated using statistical up-sampling and NWP lapse rate adjustments to existing PRISM-based climate normals at a subset of stations in BC. Specifically, up-sampling model evaluation was performed using 1951–80 climate normals from an independent set of 54 surface stations (1 m to 2347 m) which were not included in the PRISM interpolation or assimilated into the NWP reanalysis products. All models performed similarly for minimum temperature, which showed only a slight improvement over PRISM. For maximum temperature, ClimateBC, NCEP1 and ERA-Interim all performed significantly better than PRISM, in particular during spring and summer. The ERA-Interim reanalysis outperformed NCEP1 in almost all months. The results suggest that lapse rate adjustment algorithms based on reanalysis products will have greater potential as progress continues on developing NWP components.

R ésumé ?[Traduit par la rédaction] L'application des techniques d'interpolation par régression en fonction de l'altitude pour les normales de température à long terme, comme le Parameter-elevation Regressions on Independent Slopes Model (PRISM), devient très difficile dans les régions accidentées pour lesquelles on dispose de données insuffisantes, par exemple les secteurs montagneux de la Colombie-Britannique (C.-B.) au Canada. Les toutes dernières méthodes destinées à augmenter le degré de résolution des modèles d'interpolation et leur précision reposent sur la conception d'algorithmes d’échantillonnage vertical fondés sur l'ajustement des surfaces interpolées originales au moyen du gradient vertical local. Nous pouvons établir les gradients verticaux à partir de modèles statistiques (p. ex., des équations de régression polynomiales en fonction de l'altitude) ou de modèles dynamiques (p. ex., des profils verticaux de température à partir de modèles de prévision numérique du temps (PNT)). Dans la présente étude, nous comparons un algorithme d’échantillonnage vertical statistique communément utilisé, le programme ClimateBC, à deux produits de réanalyse de PNT, celle des National Centres for Environmental Prediction/National Corporation for Atmospheric Research Reanalysis 1 (NCEP1), et la réanalyse provisoire (ERA-Interim) du Centre européen pour les prévisions météorologiques à moyen terme (ECMWF). Les normales climatiques de trente ans pour les températures maximums et minimums ont été calculées en appliquant la méthode d’échantillonnage vertical statistique et l'ajustement du gradient obtenu par PNT aux normales climatiques établies à partir du PRISM pour un sous-ensemble de stations en Colombie-Britannique. Plus particulièrement, nous avons procédé à l’évaluation du modèle d’échantillonnage vertical en nous servant des normales climatiques (1951–1980), pour un ensemble de 54 stations d'observation en surface indépendantes (1?m à 2347?m), exclues du modèle d'interpolation PRISM et des produits de réanalyse de PNT. Pour tous les modèles, nous avons obtenu des résultats comparables pour la température minimum, soit une légère amélioration seulement par rapport au PRISM. Pour la température maximum, nous avons obtenu avec ClimateBC, NCEP1 et ERA-Interim, des résultats nettement plus probants qu'avec PRISM, notamment au printemps et en été. Les réanalyses ERA-Interim ont donné de meilleurs résultats que NCEP1 pour pratiquement tous les mois. D'après ces résultats, le potentiel des algorithmes d'ajustements des gradients verticaux de température, établis à partir de produits de réanalyse se renforcera à mesure que les composantes de PNT se développeront.     

Placing dots in dot maps

Dot mapping is a cartographic representation method to visualise discrete absolute values and their spatial distribution. To achieve this, dots equal in size and represented value are used. According to the dot value, a certain number of dots are used to depict a data value. These dots usually form dot clusters. The data value needs to be rounded to a multiple of the dot value. It is possible to roughly determine the visualised data value by counting the dots and multiplying this number with the dot value. As there are many parameters – dot size, dot value, map scale – to consider when designing a dot map, the manual way is very complex and time consuming. This paper presents a method to automatically create a dot representation of a dot map from given statistical data that needs no cartographic expertise. The dot representation may be combined with other elements, such as a topographic background, to form a complete map. So the algorithm can easily be integrated into the map design process. The paper refines the basic approach of automated dot mapping published earlier. The dot placement and arrangement have been improved compared to the basic method.     

Long-term transit timing monitoring and homogenous study of WASP-32

We report new photometric observations of the transiting exoplanetary system WASP-32 made by using CCD cameras at Yunnan Observatories and Ho Koon Nature Education cum Astronomical Centre, China from 2010 to 2012. Following our usual procedure, the observed data are corrected for systematic errors according to the coarse decorrelation and SYSREM algorithms so as to enhance the signal of the transit events. Combined with radial velocity data presented in the literature, our newly observed data and earlier photometric data in the literature are simultaneously analyzed to derive the physical parameters describing the system by employing the Markov chain Monte Carlo technique. The derived parameters are consistent with the result published in the original paper about WASP-32b, but the uncertainties of the new parameters are smaller than those in the original paper. Moreover, our modeling result supports a circular orbit for WASP-32b. Through the analysis of all available mid-transit times, we have refined the orbital period of WASP-32b; no evident transit timing variation is found in these transit events.     

Micron‐scale hypervelocity impact craters: Dependence of crater ellipticity and rim morphology on impact trajectory,projectile size,velocity, and shape

The interstellar collector on NASA's Stardust mission captured many particles from sources other than the interstellar dust stream. Impact trajectory may provide a means of discriminating between these different sources, and thus identifying/eliminating candidate interstellar particles. The collector's aerogel preserved a clear record of particle impact trajectory from the inclination and direction of the resultant tracks. However, the collector also contained aluminum foils and, although impact crater studies to date suggest only the most inclined impacts (>45° from normal) produce crater morphologies that indicate trajectory (i.e., distinctly elliptical), these studies have been restricted to much larger (mm and above) scales than are relevant for Stardust (μm). It is unknown how oblique impact crater morphology varies as a function of length scale, and therefore how well Stardust craters preserve details of impactor trajectory. Here, we present data from a series of impact experiments, together with complementary hydrocode modeling, that examine how crater morphology changes with impact angles for different‐sized projectiles. We find that, for our smallest spherical projectiles (2 μm diameter), the ellipticity and rim morphology provide evidence of their inclined trajectory from as little as 15° from normal incidence. This is most likely a result of strain rate hardening in the target metal. Further experiments and models find that variation in velocity and impactor shape complicate these trends, but that rim morphology remains useful in determining impact direction (where the angle of impact is >20° from normal) and may help identify candidate interstellar particle craters on the Stardust collector.