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.     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

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.

     

Using a regional numerical weather prediction model for GNSS positioning over Brazil

The global navigation satellite system (GNSS) can provide centimeter positioning accuracy at low costs. However, in order to obtain the desired high accuracy, it is necessary to use high-quality atmospheric models. We focus on the troposphere, which is an important topic of research in Brazil where the tropospheric characteristics are unique, both spatially and temporally. There are dry regions, which lie mainly in the central part of the country. However, the most interesting area for the investigation of tropospheric models is the wet region which is located in the Amazon forest. This region substantially affects the variability of humidity over other regions of Brazil. It provides a large quantity of water vapor through the humidity convergence zone, especially for the southeast region. The interconnection and large fluxes of water vapor can generate serious deficiencies in tropospheric modeling. The CPTEC/INPE (Center for Weather Forecasting and Climate Studies/Brazilian Institute for Space Research) has been providing since July 2012 a numerical weather prediction (NWP) model for South America, known as Eta. It has yield excellent results in weather prediction but has not been used in GNSS positioning. This NWP model was evaluated in precise point positioning (PPP) and network-based positioning. Concerning PPP, the best positioning results were obtained for the station SAGA, located in Amazon region. Using the NWP model, the 3D RMS are less than 10 cm for all 24 h of data, whereas the values reach approximately 60 cm for the Hopfield model. For network-based positioning, the best results were obtained mainly when the tropospheric characteristics are critical, in which case an improvement of up to 7.2 % was obtained in 3D RMS using NWP models.     

User Contribution Patterns and Completeness Evaluation of Mapillary,a Crowdsourced Street Level Photo Service

Mapillary is a Web 2.0 application which allows users to contribute crowdsourced street level photographs from all over the world. In the first part of the analysis this article reviews Mapillary data growth for continents and countries as well as the contribution behavior of individual mappers, such as the number of days of active mapping. In the second part of the analysis the study assesses Mapillary data completeness relative to a reference road network dataset at the country level. In addition, a more detailed completeness analysis is conducted for selected urban and rural areas in the US and part of northern Europe for which the completeness of Mapillary data will also be compared with that of Google Street View. Results show that Street View provides generally a better coverage on almost all road categories with some exceptions for pedestrian and cycle paths in selected cities. However, Mapillary data can be conveniently collected from any mobile device that is equipped with a photo camera. This gives Mapillary the potential to reach better coverage along off‐road segments than Google Street View.     

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...     

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.