Influence of subdaily polar motion model on nutation offsets estimated by very long baseline interferometry

This paper studies the connection between the subdaily model for polar motion used in the processing of very long baseline interferometry (VLBI) observations and the estimated nutation offsets. By convention accepted by the International Earth Rotation Service, the subdaily model for polar motion recommended for routine processing of geodetic observations does not contain any daily retrograde terms due to their one-to-one correlation with the nutation. Nevertheless, for a 24-h VLBI solution a part of the signal contained in the polar motion given by the used subdaily model is numerically mistaken for a retrograde daily sidereal signal. This fictitious retrograde daily signal contributes to the estimated nutation, leading to systematic differences between the nutation offsets from VLBI solutions computed with different subdaily polar motion models. We demonstrate this effect using solutions for all suitable 24-h VLBI sessions over a time span of 11 years (2000–2011). By changing the amplitudes of one tidal term in the underlying subdaily model for polar motion and comparing the estimated parameters to the solutions computed with the unchanged subdaily model, the paper shows and explains theoretically the effects produced by the individual subdaily terms on the VLBI nutation estimates.     

Altimetry-derived tide model for improved tide and water level forecasting along the European continental shelf

Ocean Dynamics - With the continued rise in global mean sea level, operational predictions of tidal height and total water levels have become crucial for accurate estimations and understanding of...     

A new concept for estimating the influence of vegetation on throughfall kinetic energy using aerial laser scanning

Soil loss caused by erosion has enormous economic and social impacts. Splash effects of rainfall are an important driver of erosion processes; however, effects of vegetation on splash erosion are still not fully understood. Splash erosion processes under vegetation are investigated by means of throughfall kinetic energy (TKE). Previous studies on TKE utilized a heterogeneous set of plant and canopy parameters to assess vegetation's influence on erosion by rain splash but remained on individual plant- or plot-levels. In the present study we developed a method for the area-wide estimation of the influence of vegetation on TKE using remote sensing methods. In a literature review we identified key vegetation variables influencing splash erosion and developed a conceptual model to describe the interaction of vegetation and raindrops. Our model considers both amplifying and protecting effect of vegetation layers according to their height above the ground and aggregates them into a new indicator: the Vegetation Splash Factor (VSF). It is based on the proportional contribution of drips per layer, which can be calculated via the vegetation cover profile from airborne LiDAR datasets. In a case study, we calculated the VSF using a LiDAR dataset for La Campana National Park in central Chile. The studied catchment comprises a heterogeneous mosaic of vegetation layer combinations and types and is hence well suited to test the approach. We calculated a VSF map showing the relation between vegetation structure and its expected influence on TKE. Mean VSF was 1.42, indicating amplifying overall effect of vegetation on TKE that was present in 81% of the area. Values below 1 indicating a protective effect were calculated for 19% of the area. For future work, we recommend refining the weighting factor by calibration to local conditions using field-reference data and comparing the VSF with TKE field measurements. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Optimized formulas for the gravitational field of a tesseroid

Various tasks in geodesy, geophysics, and related geosciences require precise information on the impact of mass distributions on gravity field-related quantities, such as the gravitational potential and its partial derivatives. Using forward modeling based on Newton’s integral, mass distributions are generally decomposed into regular elementary bodies. In classical approaches, prisms or point mass approximations are mostly utilized. Considering the effect of the sphericity of the Earth, alternative mass modeling methods based on tesseroid bodies (spherical prisms) should be taken into account, particularly in regional and global applications. Expressions for the gravitational field of a point mass are relatively simple when formulated in Cartesian coordinates. In the case of integrating over a tesseroid volume bounded by geocentric spherical coordinates, it will be shown that it is also beneficial to represent the integral kernel in terms of Cartesian coordinates. This considerably simplifies the determination of the tesseroid’s potential derivatives in comparison with previously published methodologies that make use of integral kernels expressed in spherical coordinates. Based on this idea, optimized formulas for the gravitational potential of a homogeneous tesseroid and its derivatives up to second-order are elaborated in this paper. These new formulas do not suffer from the polar singularity of the spherical coordinate system and can, therefore, be evaluated for any position on the globe. Since integrals over tesseroid volumes cannot be solved analytically, the numerical evaluation is achieved by means of expanding the integral kernel in a Taylor series with fourth-order error in the spatial coordinates of the integration point. As the structure of the Cartesian integral kernel is substantially simplified, Taylor coefficients can be represented in a compact and computationally attractive form. Thus, the use of the optimized tesseroid formulas particularly benefits from a significant decrease in computation time by about 45 % compared to previously used algorithms. In order to show the computational efficiency and to validate the mathematical derivations, the new tesseroid formulas are applied to two realistic numerical experiments and are compared to previously published tesseroid methods and the conventional prism approach.     

Bright water: hydrosols, water conservation and climate change

Because air?Cwater and water?Cair interfaces are equally refractive, cloud droplets and microbubbles dispersed in bodies of water reflect sunlight in much the same way. The lifetime of sunlight-reflecting microbubbles, and hence the scale on which they may be applied, depends on Stokes Law and the influence of ambient or added surfactants. Small bubbles backscatter light more efficiently than large ones, opening the possibility of using highly dilute micron-radius hydrosols to substantially brighten surface waters. Such microbubbles can noticeably increase water surface reflectivity, even at volume fractions of parts per million and such loadings can be created at an energy cost as low as J m^???2 to initiate and mW m^???2 to sustain. Increasing water albedo in this way can reduce solar energy absorption by as much as 100 W m^???2, potentially reducing equilibrium temperatures of standing water bodies by several Kelvins. While aerosols injected into the stratosphere tend to alter climate globally, hydrosols can be used to modulate surface albedo, locally and reversibly, without risk of degrading the ozone layer or altering the color of the sky. The low energy cost of microbubbles suggests a new approach to solar radiation management in water conservation and geoengineering: Don??t dim the Sun; Brighten the water.     

Combination of multi-mission altimetry data along the Mekong River with spatio-temporal kriging

River water-level time series at fixed geographical locations, so-called virtual stations, have been computed from single altimeter crossings for many years. Their temporal resolution is limited by the repeat cycle of the individual altimetry missions. The combination of all altimetry measurements along a river enables computing a water-level time series with improved temporal and spatial resolutions. This study uses the geostatistical method of spatio-temporal ordinary kriging to link multi-mission altimetry data along the Mekong River. The required covariance models reflecting the water flow are estimated based on empirical covariance values between altimetry observations at various locations. In this study, two covariance models are developed and tested in the case of the Mekong River: a stationary and a non-stationary covariance model. The proposed approach predicts water-level time series at different locations along the Mekong River with a temporal resolution of 5 days. Validation is performed against in situ data from four gauging stations, yielding RMS differences between 0.82 and 1.29 m and squared correlation coefficients between 0.89 and 0.94. Both models produce comparable results when used for combining data from Envisat, Jason-1, and SARAL for the time period between 2002 and 2015. The quality of the predicted time series turns out to be robust against a possibly decreasing availability of altimetry mission data. This demonstrates that our method is able to close the data gap between the end of the Envisat and the launch of the SARAL mission with interpolated time series.     

Post-injection trapping of mobile CO<Subscript><Emphasis Type="Bold">2</Emphasis></Subscript> in deep aquifers: Assessing the importance of model and parameter uncertainties

Predicting the fate of the injected CO₂ is crucial for the safety of carbon storage operations in deep saline aquifers: especially the evolution of the position, the spreading and the quantity of the mobile CO₂ plume during and after the injection has to be understood to prevent any loss of containment. Fluid flow modelling is challenging not only given the uncertainties on subsurface formation intrinsic properties (parameter uncertainty) but also on the modelling choices/assumptions for representing and numerically implementing the processes occurring when CO₂ displaces the native brine (model uncertainty). Sensitivity analysis is needed to identify the group of factors which contributes the most to the uncertainties in the predictions. In this paper, we present an approach for assessing the importance of model and parameter uncertainties regarding post-injection trapping of mobile CO₂. This approach includes the representation of input parameters, the choice of relevant simulation outputs, the assessment of the mobile plume evolution with a flow simulator and the importance ranking for input parameters. A variance-based sensitivity analysis is proposed, associated with the ACOSSO-like meta-modelling technique to tackle the issues linked with the computational burden posed by the use of long-running simulations and with the different types of uncertainties to be accounted for (model and parameter). The approach is tested on a potential site for CO₂ storage in the Paris basin (France) representative of a project in preliminary stage of development. The approach provides physically sound outcomes despite the challenging context of the case study. In addition, these outcomes appear very helpful for prioritizing the future characterisation efforts and monitoring requirements, and for simplifying the modelling exercise.     

Bücher- und Zeitschriftenschau

Ohne Zusammenfassung     

The dust distribution within the inner coma of comet P/Halley 1982i: encounter by Giotto's impact detectors

Analysis of the data from Giotto's Dust Impact Detection System experiment (DIDSY) is presented. These data represent measurement of the size of dust grains incident on the Giotto dust shield along its trajectory through the coma of comet P/Halley on 1986 March 13/14. First detection occurred at some 287000 km distance from the nucleus on the inbound leg; the majority of the DIDSY subsystems remained operational after closest approach (604 km) yielding the last detection at about 202000 km from the nucleus. In order to improve the data coverage (and especially for the smallest grains, to approximately 10(-19) kg particle mass), data from the PIA instrument has been combined with DIDSY data. Flux profiles are presented for the various mass channels showing, to a first approximation, a 1/R2 flux dependence, where R is the distance of the detection point from the cometary nucleus, although significant differences are noted. Deviations from this dependence are observed, particularly close to the nucleus. From the flux profiles, mass and geometrical area distributions for the dust grains are derived for the trajectory through the coma. Groundbased CCD imaging of the dust continuum in the inner coma at the time of encounter is also used to derive the area of grains intercepted by Giotto. The results are consistent with the area functions derived by Giotto data and the low albedo of the grains deduced from infrared emission. For the close encounter period (-5 min to +5 min), the cumulative mass distribution function has been investigated, initially in 20 second periods; there is strong evidence from the data for a steepening of the index of the mass distribution for masses greater than 10(-13) kg during passage through dust jets which is not within the error limits of statistical uncertainty. The fluences for dust grains along the entire trajectory is calculated; it is found that extrapolation of the spectrum determined at intermediate masses (cumulative mass index alpha = 0.85) is not able to account for the spacecraft deceleration as observed by the Giotto Radio Science Experiment and by ESOC tracking operations. Data at large masses (>10(-8) kg) recently analysed from the DIDSY data set show clear evidence of a decrease in the mass distribution index at these masses within the coma, and it is shown that such a value of the mass index can provide sufficient mass for consistency with the observed deceleration. The total particulate mass output from the nucleus of comet P/Halley at the time of encounter would be dependent on the maximum mass emitted if this change in slope observed in the coma were also applicable to the emission from the nucleus; this matter is discussed in the text. The flux time profiles have been converted through a simple approach to modeling of the particle trajectories to yield an indication of nucleus surface activity. There is indication of an enhancement in flux at t approximately -29 s corresponding to crossing of the dawn terminator, but the flux detected prior to crossing of the dawn terminator is shown to be higher than predicted by simple modelling. Further enhancements corresponding to jet activity are detected around +190 s and +270 s.