HCMM energy budget data as a model input for assessing regions of high potential groundwater pollution

The author has identified the following significant results. In early April 1978, heavy spring runoff from snowmelt caused significant flooding along a portion of the Big Sioux River Basin in southeastern South Dakota. The flooded area was visible from surrounding areas on a May 15 HCMM IR test image. On May 15, the flood waters had receded but an area of anomalous residual high soil moisture remained. The high soil moisture area was not visible on a HCMM day visible test image of the same scene, or on LANDSAT imagery. To evaluate the effect of water table depth on surface temperatures, thermal scanner data collected on September 5 and 6, 1978 at approximate HCMM overpass times at an altitude of 3650 m were analyzed. Apparent surface temperatures measured by the scanner included emittance contributions from soil surface and the land cover. Results indicated that the shallow water tables produced a damping of the amplitude of the diurnal surface temperature wave

Pseudo-Stochastic Orbit Modeling Techniques for Low-Earth Orbiters

The Earth's non-spherical mass distribution and atmospheric drag cause the strongest perturbations on very low-Earth orbiting satellites (LEOs). Models of gravitational and non-gravitational accelerations are utilized in dynamic precise orbit determination (POD) with GPS data, but it is also possible to derive LEO positions based on GPS precise point positioning without dynamical information. We use the reduced-dynamic technique for LEO POD, which combines the geometric strength of the GPS observations with the force models, and investigate the performance of different pseudo-stochastic orbit parametrizations, such as instantaneous velocity changes (pulses), piecewise constant accelerations, and continuous piecewise linear accelerations. The estimation of such empirical orbit parameters in a standard least-squares adjustment process of GPS observations, together with other relevant parameters, strives for the highest precision in the computation of LEO trajectories. We used the procedures for the CHAMP satellite and found that the orbits may be validated by means of independent SLR measurements at the level of 3.2cm RMS. Validations with independent accelerometer data revealed correlations at the level of 95% in the along-track direction. As expected, the empirical parameters compensate to a certain extent for deficiencies in the dynamic models. We analyzed the capability of pseudo-stochastic parameters for deriving information about the mismodeled part of the force field and found evidence that the resulting orbits may be used to recover force field parameters, if the number of pseudo-stochastic parameters is large enough. Results based on simulations showed a significantly better performance of acceleration-based orbits for gravity field recovery than for pulse-based orbits, with a quality comparable to a direct estimation if unconstrained accelerations are set up every 30

Propagation of atmospheric model errors to gravity potential harmonics—impact on GRACE de-aliasing

High-frequency, time-varying mass redistributions in the ocean and atmosphere have an impact on GRACE gravity field solutions due to the space-time sampling characteristics of signal and orbit. Consequently, aliasing of these signals into the GRACE observations is present and needs to be taken into account during data analysis by applying atmospheric and oceanic model data (de-aliasing). As the accuracy predicted prior to launch could not yet be achieved in the analysis of real GRACE data, the de-aliasing process and related geophysical model uncertainties are regarded as a potential error source in GRACE gravity field determination. Therefore, this study aims to improve the de-aliasing process in order to obtain a more accurate GRACE gravity field time-series. As these time-series provide estimates for the integrated mass transport in the Earth system, like the global water cycle and solid Earth geophysical processes, any increase in accuracy will lead to improvements in the geophysical interpretation of the results. So in conclusion, improving the de-aliasing is of relevance for a better understanding of geophysical processes. By no longer regarding the atmosphere and ocean model output as error-free, deeper insight into the impact of such uncertainties on the de-aliasing and on the resulting GRACE gravity field models can be obtained. For this purpose, in a first step, a full error propagation of the atmospheric and oceanic model parameters up to the de-aliasing gravity field coefficients is performed and the GRACE K-Band-Satellite-to-Satellite Tracking (KBR-SST) residuals, as an intermediate gravity field result, are analysed. The paper reviews the standard GRACE de-aliasing process and presents the mathematical model applied for the error propagation. Specifically, the effect of uncertainties in the atmospheric input parameters (temperature, surface pressure, specific humidity, geopotential) on the gravity field potential coefficients used for de-aliasing is shown in several scenarios. Finally, the impact of de-aliasing products (with and without error propagation) on a GRACE gravity field solution is investigated on the level of observation residuals. From the results obtained in this study it can be concluded that with respect to the current GRACE error budget, atmospheric model uncertainties do not play a prominent role in the error budget of current GRACE gravity field solutions. Nevertheless, in order to fully exploit the GRACE measurements towards the baseline accuracy, an optimized de-aliasing is needed. In this case, GRACE gravity field solutions are sensitive to uncertainties in atmospheric and oceanic models. Thus, the associated geophysical model errors shall be taken into account in the de-aliasing proces

Efficient satellite orbit modelling using pseudo-stochastic parameters

If the force field acting on an artificial Earth satellite is not known a priori with sufficient accuracy to represent its observations on their accuracy level, one may introduce so-called pseudo-stochastic parameters into an orbit determination process, e.g. instantaneous velocity changes at user-defined epochs or piecewise constant accelerations in user-defined adjacent time subintervals or piecewise linear and continuous accelerations in adjacent time subintervals. The procedures, based on standard least-squares, associated with such parameterizations are well established, but they become inefficient (slow) if the number of pseudo-stochastic parameters becomes large. We develop two efficient methods to solve the orbit determination problem in the presence of pseudo-stochastic parameters. The results of the methods are identical to those obtained with conventional least-squares algorithms. The first efficient algorithm also provides the full variance-covariance matrix; the second, even more efficient algorithm, only parts of i

Evaluation of HCMM data for assessing soil moisture and water table depth

Soil moisture in the 0-cm to 4-cm layer could be estimated with 1-mm soil temperatures throughout the growing season of a rainfed barley crop in eastern South Dakota. Empirical equations were developed to reduce the effect of canopy cover when radiometrically estimating the soil temperature. Corrective equations were applied to an aircraft simulation of HCMM data for a diversity of crop types and land cover conditions to estimate the soil moisture. The average difference between observed and measured soil moisture was 1.6% of field capacity. Shallow alluvial aquifers were located with HCMM predawn data. After correcting the data for vegetation differences, equations were developed for predicting water table depths within the aquifer. A finite difference code simulating soil moisture and soil temperature shows that soils with different moisture profiles differed in soil temperatures in a well defined functional manner. A significant surface thermal anomaly was found to be associated with shallow water tables

Validating ocean tide loading models using GPS

Abstract.: Ocean tides cause periodic deformations of the Earth's surface, also referred to as ocean tide loading (OTL). Tide-induced displacements of the Earth's crust relying on OTL models are usually taken into account in GPS (Global Positioning System) data analyses. On the other hand, it is also possible to validate OTL models using GPS analyses. The following simple approach is used to validate OTL models. Based on a particular model, instantaneous corrections of the site coordinates due to OTL are computed. Site-specific scale factors, f, for these corrections are estimated in a standard least-squares adjustment process of GPS observations together with other relevant parameters. A resulting value of f close to unity indicates a good agreement of the model with the actual site displacements. Such scale factors are computed for about 140 globally distributed IGS (International GPS Service) tracking sites. Three OTL models derived from the ocean tide models FES95.2.1, FES99, and GOT00.2 are analyzed. As expected, the most reliable factors are estimated for sites with a large loading effect. In general, the scaling factors have a value close to unity and no significant differences between the three ocean tide models could be observed. It is found that the validation approach is easy to apply. Without requiring much additional effort for a global and self-consistent GPS data analysis, it allows detection of general model misfits on the basis of a large number of globally distributed sites. For detailed validation studies on OTL models, the simultaneous estimation of amplitudes and phases for the main contributing partial tides within a GPS parameter adjustment process would provide more detailed answer

Osculating orbits in Schwarzschild spacetime, with an application to extreme mass-ratio inspirals

We present a method to integrate the equations of motion that govern bound, accelerated orbits in Schwarzschild spacetime. At each instant the true worldline is assumed to lie tangent to a reference geodesic, called an osculating orbit, such that the worldline evolves smoothly from one such geodesic to the next. Because a geodesic is uniquely identified by a set of constant orbital elements, the transition between osculating orbits corresponds to an evolution of the elements. In this paper we derive the evolution equations for a convenient set of orbital elements, assuming that the force acts only within the orbital plane; this is the only restriction that we impose on the formalism, and we do not assume that the force must be small. As an application of our method, we analyze the relative motion of two massive bodies, assuming that one body is much smaller than the other. Using the hybrid Schwarzschild/post-Newtonian equations of motion formulated by Kidder, Will, and Wiseman, we treat the unperturbed motion as geodesic in a Schwarzschild spacetime whose mass parameter is equal to the system's total mass. The force then consists of terms that depend on the system's reduced mass. We highlight the importance of conservative terms in this force, which cause significant long-term changes in the time-dependence and phase of the relative orbit. From our results we infer some general limitations of the radiative approximation to the gravitational self-force, which uses only the dissipative terms in the force.Comment: 18 pages, 6 figures, final version to be published in Physical Review