Integrability of motion around galactic razor-thin disks

We consider the three-dimensional bounded motion of a test particle around razor-thin disk configurations, by focusing on the adiabatic invariance of the vertical action associated with disk-crossing orbits. We find that it leads to an approximate third integral of motion predicting envelopes of the form \(Z(R)\propto [\varSigma (R)]^{-1/3}\), where R is the radial galactocentric coordinate, Z is the z-amplitude (vertical amplitude) of the orbit and \(\varSigma \) represents the surface mass density of the thin disk. This third integral, which was previously formulated for the case of flattened 3D configurations, is tested for a variety of trajectories in different thin-disk models.     

Inversion of two-dimensional tidal open boundary conditions of M<Subscript>2</Subscript> constituent in the Bohai and Yellow Seas

Two-dimensional tidal open boundary conditions of the M2 constituent in the Bohai and Yellow Seas(BYS) have been estimated by assimilating T/P altimeter data.During inversion,independent point(IP) strategy was used,in which several IPs on the open boundary is assumed,values at these IPs can be optimized with an adjoint method,and those at other grid points are determined by linearly interpolating the values at IPs.The reasonability and feasibility of the model are tested by ideal twin experiments.In the practical experiment(PE) after assimilation,the cost function may reach 1% or less of its initial value.Mean absolute errors in amplitude and phase can be less than 5 cm and 5°,respectively,and the obtained co-chart can show the character of the M2 constituent in the BYS.The results of the PE indicate that using only two IPs on the open boundary can yield better simulated results.     

Carbon costs and benefits of France’s biomass energy production targets

Background

Concern about climate change has motivated France to reduce its reliance on fossil fuel by setting targets for increased biomass-based renewable energy production. This study quantifies the carbon costs and benefits for the French forestry sector in meeting these targets. A forest growth and harvest simulator was developed for French forests using recent forest inventory data, and the wood-use chain was reconstructed from national wood product statistics. We then projected wood production, bioenergy production, and carbon balance for three realistic intensification scenarios and a business-as-usual scenario. These intensification scenarios targeted either overstocked, harvest-delayed or currently actively managed stands.

Results

All three intensification strategies produced 11.6–12.4 million tonnes of oil equivalent per year of wood-based energy by 2026, which corresponds to the target assigned to French wood-energy to meet the EU 2020 renewable energy target. Sustaining this level past 2026 will be challenging, let alone further increasing it. Although energy production targets can be reached, the management intensification required will degrade the near-term carbon balance of the forestry sector, compared to continuing present-day management. Even for the best-performing intensification strategy, i.e., reducing the harvest diameter of actively managed stands, the carbon benefits would only become apparent after 2040. The carbon balance of a strategy putting abandoned forests back into production would only break even by 2055; the carbon balance from increasing thinning in managed but untended stands would not break even within the studied time periods, i.e. 2015–2045 and 2046–2100. Owing to the temporal dynamics in the components of the carbon balance, i.e., the biomass stock in the forest, the carbon stock in wood products, and substitution benefits, the merit order of the examined strategies varies over time.

Conclusions

No single solution was found to improve the carbon balance of the forestry sector by 2040 in a way that also met energy targets. We therefore searched for the intensification scenario that produces energy at the lowest carbon cost. Reducing rotation time of actively managed stands is slightly more efficient than targeting harvest-delayed stands, but in both cases, each unit of energy produced has a carbon cost that only turns into a benefit between 2060 and 2080.

     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.     

A new method for specular and diffuse pseudorange multipath error extraction using wavelet analysis

Multipath remains one of the major challenges in Global Navigation Satellite System (GNSS) positioning because it is considered the dominant source of ranging errors, which can be classified into specular and diffuse types. We present a new method using wavelets to extract the pseudorange multipath in the time domain and breaking it down into the two components. The main idea is an analysis-reconstruction approach based on application of both continuous wavelet transform (CWT) and discrete wavelet transform (DWT). The proposed procedure involves the use of L1 code-minus-carrier (CMC) observable where higher-frequency terms are isolated as residuals. CMC residuals are analyzed by applying the CWT, and we propose the scalogram as a technique for discerning time–frequency variations of the multipath signal. Unlike Fourier transform, the potential of the CWT scalogram for examining the non-stationary and multifrequency nature of the multipath is confirmed as it simultaneously allows fine detection and time localization of the most representative frequencies of the signal. This interpretation of the CWT scalogram is relevant when choosing the levels of reconstruction with DWT, allowing accurate time domain extraction of both the specular and diffuse multipath. The performance and robustness of the method and its boundary applicability are assessed. The experiment was carried out using a receiver of Campania GNSS Network. The results are given in which specular multipath error is achieved using DWT level 7 approximation component and diffuse multipath error is achieved using DWT level 6 denoised detail component.     

Impact of different sampling rates on precise point positioning performance using online processing service

In this study, the effect of different sampling rates (i.e. observation recording interval) on the Precise Point Positioning (PPP) solutions in terms of accurac...     

Mullite in Libyan Desert Glass: Evidence for high-temperature/low-pressure formation

Libyan Desert Glass (LDG) is a SiO₂-rich natural glass whose origin, formation mechanism, and target material are highly debated. We here report on the finding of a lens-shaped whitish inclusion within LDG. The object is dominantly composed of siliceous glass and separated from the surrounding LDG by numerous cristobalite grains. Within cristobalite, several regions rich in mullite often associated with ilmenite were detected. Mineral assemblage, chemical composition, and grain morphologies suggest that mullite was formed by thermal decomposition of kaolinitic clay at atmospheric pressure and T ≥ 1600 °C and also attested to high cooling rates under nonequilibrium conditions. Cristobalite contains concentric and irregular internal cracks and is intensely twinned, indicating that first crystallized β-cristobalite inverted to α-cristobalite during cooling of the SiO₂-rich melt. The accompanied volume reduction of 4% induced the high density of defects. The whitish inclusion also contains several partly molten rutile grains evidencing that at least locally the LDG melt was at T ≥ 1800 °C. Based on these observations, it is concluded that LDG was formed by high-temperature melting of kaolinitic clay-, rutile-, and ilmenite-bearing Cenozoic sandstone or sand very likely during an asteroid or comet impact onto Earth. While melting and ejection occurred at high pressures, the melt solidified quickly at atmospheric pressure.     

Observations of the sodium emission from Jupiter,region C

Measurements made when testing an interferometric system for imaging the sodium emission indicate that, between 3 and 7.5 arc min (9 and 23R_J) from Jupiter, the observed emissions can probably be attributed to terrestrial airglow.     

Surface gravitational field and topography changes induced by the Earth’s fluid core motions

Using a Love number formalism, the elastic deformations of the mantle and the mass redistribution gravitational potential within the Earth induced by the fluid pressure acting at the core–mantle boundary (CMB) are computed. This pressure field changes at a decadal time scale and may be estimated from observations of the surface magnetic field and its secular variation. First, using a spherical harmonic expansion, the poloidal and toroidal part of the fluid velocity field at the CMB for the last 40 years is computed, under the hypothesis of tangential geostrophy. Then the associated geostrophic pressure, whose order of magnitude is about 1000 Pa, is computed. The surface topography induced by this pressure field is computed using Love numbers, and is a few millimetres. The mass redistribution gravitational potential induced by these deformations and, in particular, the zonal components of the related surface gravitational potential perturbation (J₂, J₃ and J₄ coefficients), are calculated. Overall perturbations for the J₂ coefficient of about 10^–10, for J₃ of about 10^–11 and for J₄ are found of about 0.3×10^–11. Finally, these theoretical results are compared with recent observations of the decadal variation of J₂ from satellite laser ranging. Results concerning J₂ can be described as follows: first, they are one order of magnitude too small to explain the observed decadal variation of J₂ and, second, they show a significant linear trend over the last 40 years, whose rate of decrease amounts to 7% of the observed value.     

Computerised development control and approval system for City Hall of Kuala Lumpur

This paper will discuss the computerised development control and approval system being developed for the Planning and Development Control Department, City Hall of Kuala Lumpur, with stress on the GIS architecture developed within the system. The prospects and challenges towards implementation of the system are also discussed.