Calibration of a frontal ablation parameterization applied to Greenland's peripheral calving glaciers

<p>Greenland's Peripheral Glaciers (PGs) are glaciers that are weakly or not connected to the Ice Sheet. Many are tidewater, losing mass via frontal ablation. Without comprehensive regional observations or enough individual estimates of frontal ablation, constraining model parameters remains a challenging task in this region. We present three independent ways to calibrate the calving parameterization implemented in the Open Global Glacier Model (OGGM) and asses the impact of accounting for frontal ablation on the estimate of ice stored in PGs. We estimate an average regional frontal ablation flux for PGs of 7.94±4.15 Gtyr<sup>-1</sup> after calibrating the model with two different satellite velocity products, and of 0.75±0.55 Gt yr<sup>-1</sup> if the model is constrained using frontal ablation fluxes derived from independent modelled Surface Mass Balance (SMB) averaged over an equilibrium reference period (1961-1990). This second method is based on the assumption that most PGs during that time have an equilibrium between mass gain via SMB and mass loss via frontal ablation. This assumption can serve as a basis to assess the order of magnitude of dynamic mass loss of glaciers when compared to the SMB imbalance. By comparing the model output after applying both calibration methods, we find that the model is not able to predict individual tidewater glacier dynamics if it relies only on SMB estimates and the assumption of a closed budget to constrain the model. The differences between the results from both calibration methods serve as an indication of how strong the dynamic imbalance might have been for PGs during that reference period.</p>

Brief communication: Heterogenous thinning and subglacial lake activity on Thwaites Glacier, West Antarctica

A system of subglacial lakes drained on Thwaites Glacier from 2012–2014. To improve coverage for subsequent drainage events, we extended the elevation and ice-velocity time series on Thwaites Glacier through austral winter 2019. These new observations document a second drainage cycle in 2017/18 and identified two new lake systems located in the western tributaries of Thwaites and Haynes glaciers. In situ and satellite velocity observations show temporary < 3 % speed fluctuations associated with lake drainages. In agreement with previous studies, these observations suggest that active subglacial hydrology has little influence on thinning and retreat of Thwaites Glacier on decadal to centennial timescales.

Wave Orbital Velocity Effects on Radar Doppler Altimeter for Sea Monitoring

The orbital velocity of sea wave particles affects the value of sea surface parameters as measured by radar Doppler altimeters (also known as delay Doppler altimeter (DDA)). In DDA systems, the along-track resolution is attained by algorithms that take into account the Doppler shift induced by the component along the Earth/antenna direction of the satellite velocity, VS. Since the vertical component of the wave particle orbital velocity also induces an additional Doppler effect (in the following R-effect), an error arises on the positioning of the target on the sea surface. A numerical investigation shows that when the wavelength of sea waves is of the same order of magnitude of the altimeter resolution, the shape of the waveform might be significantly influenced by the R-effect. The phenomenon can be particularly important for the monitoring of long swells, such as those that often take place in the oceans.

Brief Communication: Heterogenous thinning and subglacial lake activity on Thwaites Glacier, West Antarctica

A system of subglacial lakes drained on Thwaites Glacier from 2012–2014. To improve coverage for subsequent drainage events, we extended the elevation and ice velocity time series on Thwaites Glacier through austral winter 2019. These new observations document a second drainage cycle and identified two new lake systems located in the western tributaries of Thwaites and Haynes Glaciers. In situ and satellite velocity observations show temporary

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

Basilisk: Bayesian hierarchical inference of the galaxy–halo connection using satellite kinematics – I. Method and validation

ABSTRACT We present a Bayesian hierarchical inference formalism (Basilisk) to constrain the galaxy–halo connection using satellite kinematics. Unlike traditional methods, Basilisk does not resort to stacking the kinematics of satellite galaxies in bins of central luminosity, and does not make use of summary statistics, such as satellite velocity dispersion. Rather, Basilisk leaves the data in its raw form and computes the corresponding likelihood. In addition, Basilisk can be applied to flux-limited, rather than volume-limited samples, greatly enhancing the quantity and dynamic range of the data. And finally, Basilisk is the only available method that simultaneously solves for halo mass and orbital anisotropy of the satellite galaxies, while properly accounting for scatter in the galaxy–halo connection. Basilisk uses the conditional luminosity function to model halo occupation statistics, and assumes that satellite galaxies are a relaxed tracer population of the host halo’s potential with kinematics that obey the spherical Jeans equation. We test and validate Basilisk using mocks of varying complexity, and demonstrate that it yields unbiased constraints on the galaxy–halo connection and at a precision that rivals galaxy–galaxy lensing. In particular, Basilisk accurately recovers the full PDF of the relation between halo mass and central galaxy luminosity, and simultaneously constrains the orbital anisotropy of the satellite galaxies. Basilisk ’s inference is not affected by potential velocity bias of the central galaxies, or by slight errors in the inferred, radial profile of satellite galaxies that arise as a consequence of interlopers and sample impurity.

Differential astrometric framework for the Jupiter relativistic experiment with Gaia

We employ differential astrometric methods to establish a small field reference frame stable at the microarcsecond (μas) level on short time-scales using high-cadence simulated observations taken by Gaia in 2017 February of a bright star close to the limb of Jupiter, as part of the relativistic experiment on Jupiter’s quadrupole. We achieve subμas-level precision along scan through a suitable transformation of the field angles into a small-field tangent plane and a least-squares fit over several overlapping frames for estimating the plate and geometric calibration parameters with tens of reference stars that lie within ∼0.5 deg from the target star, assuming perfect knowledge of stellar proper motions and parallaxes. Furthermore, we study the effects of unmodelled astrometric parameters on the residuals and find that proper motions have a stronger effect than unmodelled parallaxes, e.g. unmodelled Gaia DR2 proper motions introduce extra residuals of ∼23 μas (AL) and 69 μas (AC) versus the ∼5 μas (AL) and 17 μas (AC) due to unmodelled parallaxes. On the other hand, assuming catalogue errors in the proper motions and parallaxes such as those from Gaia DR2 has a minimal impact on the stability introducing subμas and μas level residuals in the along and across scanning direction, respectively. Finally, the effect of a coarse knowledge in the satellite velocity components (with time-dependent errors of 10 μas s−1) is capable of enlarging the size of the residuals to roughly 0.2 mas.

Performance Assessment of BDS and GPS/BDS Velocity Estimation with Stand-alone Receiver

Accurate velocity estimates are critical in highly dynamic positioning, airborne gravimetry, and geophysics applications. This paper focuses on the evaluation of the performance of velocity estimation using the BeiDou navigation satellite system (BDS) alone and integrated Global Positioning System (GPS)/BDS. Firstly, we analyse and compare the position-derivation method and analytical method which are used to calculate BDS satellite velocity from broadcast ephemeris. Results show that the accuracy of the estimated velocity by position-derivation method can be within 1 mm/s and better than that of the analytical method. Secondly, velocity estimation tests were carried out both in static and kinematic modes. The results show that: 1) the accuracy of BDS velocity estimation is in the same order of magnitude to that of GPS; 2) Compared with a single navigation system, the stability and accuracy of velocity estimation can be remarkably improved by integrated GPS/BDS, especially under conditions of poor observation; 3) Compared with Helmert variance component estimation, it is more appropriate and efficient to assign the weights of different types of observations using equivalent weight ratio. Finally, the ionospheric influence on velocity estimation with single-frequency observations can reach several mm/s; this influence can be significantly mitigated by using ionosphere-free combination observations.