Characterization of crater morphometry on the Moon and Mercury from altimetry observations

Thesis (S.M.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, 2012.Cataloged from PDF version of thesis.Includes bibliographical references (p. 51-55).Recently acquired altimetry data from laser altimeters are used to assess the morphometry of impact craters. Data acquired by the Mercury Laser Altimeter on the MESSENGER spacecraft are used to measure the depths and diameters of 537 craters at the high northern latitudes on Mercury, including 42 polar-deposit-hosting craters (PDCs) which host material that is bright to earth-based radar observations. A comparative analysis suggests that the radar-bright material forms a thin (< 20 m) layer emplaced preferentially in comparatively young craters, contradicting an earlier morphometric study that indicated that PDCs contained a thick layer of water ice and dust. Topographic datasets from the lunar surface, collected by the Lunar Orbiter Laser Altimeter onboard the LRO spacecraft, are also used to evaluate the morphometry of 1,356 lunar craters. We study the morphologic change between the simple and complex crater regime, a manifestation of the transition between gravity-dominated and strength-dominated impact-forming processes, on the Moon and Mercury. The Moons transition diameter is near 16 km, in line with previous studies, while Mercurys is near 8 km, 2 km smaller than previously determined. The onset of gravity-dominated mechanisms scales inversely with gravity, which explains why Mercurys transition diameter is approximately half of the Moons.by Matthieu Jean Talpe.S.M

Developing a Complex Independent Component Analysis (CICA) technique to extract non-stationary patterns from geophysical time series

In recent decades, decomposition techniques have enabled increasingly more applications for dimension reduction, as well as extraction of additional information from geophysical time series. Traditionally, the principal component analysis (PCA)/empirical orthogonal function (EOF) method and more recently the independent component analysis (ICA) have been applied to extract, statistical orthogonal (uncorrelated), and independent modes that represent the maximum variance of time series, respectively. PCA and ICA can be classified as stationary signal decomposition techniques since they are based on decomposing the autocovariance matrix and diagonalizing higher (than two) order statistical tensors from centered time series, respectively. However, the stationarity assumption in these techniques is not justified for many geophysical and climate variables even after removing cyclic components, e.g., the commonly removed dominant seasonal cycles. In this paper, we present a novel decomposition method, the complex independent component analysis (CICA), which can be applied to extract non-stationary (changing in space and time) patterns from geophysical time series. Here, CICA is derived as an extension of real-valued ICA, where (a) we first define a new complex dataset that contains the observed time series in its real part, and their Hilbert transformed series as its imaginary part, (b) an ICA algorithm based on diagonalization of fourth-order cumulants is then applied to decompose the new complex dataset in (a), and finally, (c) the dominant independent complex modes are extracted and used to represent the dominant space and time amplitudes and associated phase propagation patterns. The performance of CICA is examined by analyzing synthetic data constructed from multiple physically meaningful modes in a simulation framework, with known truth. Next, global terrestrial water storage (TWS) data from the Gravity Recovery And Climate Experiment (GRACE) gravimetry mission (2003–2016), and satellite radiometric sea surface temperature (SST) data (1982–2016) over the Atlantic and Pacific Oceans are used with the aim of demonstrating signal separations of the North Atlantic Oscillation (NAO) from the Atlantic Multi-decadal Oscillation (AMO), and the El Niño Southern Oscillation (ENSO) from the Pacific Decadal Oscillation (PDO). CICA results indicate that ENSO-related patterns can be extracted from the Gravity Recovery And Climate Experiment Terrestrial Water Storage (GRACE TWS) with an accuracy of 0.5–1 cm in terms of equivalent water height (EWH). The magnitude of errors in extracting NAO or AMO from SST data using the complex EOF (CEOF) approach reaches up to ~50% of the signal itself, while it is reduced to ~16% when applying CICA. Larger errors with magnitudes of ~100% and ~30% of the signal itself are found while separating ENSO from PDO using CEOF and CICA, respectively. We thus conclude that the CICA is more effective than CEOF in separating non-stationary patterns

Initial Orbit Determination Results from the University of Luxembourg using Spire GNSS Tracking Data

CubeSats constellations using commercial off-the-shelf components have been studied for different applications, such as GNSS Radio Occultation (GNSS-RO). Furthermore, precise orbit determination of Low Earth Orbit (LEO) CubeSats based on multiple GNSS constellations would open new opportunities for scientific applications such as Earth’s gravity field measurements. In GNSS kinematic orbit determination, which is the common method used for small sats, the derived orbits are affected by noise, data gaps, outliers, measurement errors as well as poor geometry of the observations. Our work seeks to mitigate these issues and we present two areas of research: 1) GNSS network processing of GPS and Galileo constellations and 2) kinematic orbit determination of a set of Spire CubeSats that host a GNSS-RO payload. An initial architecture of kinematic orbit processing for the Spire GNSS-RO CubeSats constellation is obtained and the details on validations and limitations are discussed in more details. In addition, we showcase the agreement between the GNSS orbit products produced at the University of Luxembourg (UL) with those of the Center for Orbit Determination in Europe (CODE). Finally, the Spire kinematic orbits based on the raw observation approach are derived and compared to the L1B Spire orbit products

Ice mass change in Greenland and Antarctica between 1993 and 2013 from satellite gravity measurements

We construct long-term time series of Greenland and Antarctic ice sheet mass change from satellite gravity measurements. A statistical reconstruction approach is developed based on a Principal Component Analysis to combine high-resolution spatial modes from the Gravity Recovery and Climate Experiment (GRACE) mission with the gravity information from conventional satellite track-ing data. Uncertainties of this reconstruction are rigorously assessed; they include temporal limitations for short GRACE measurements, spatial limitations for the low-resolution conventional tracking data measurements, and limitations of the estimated statistical relationships between low and high degree potential coe�cients re ected in the PCA modes. Trends of mass variations in Greenland and Antarctica are assessed against a number of previous studies. The resulting time series for Greenland show a higher rate of mass loss than other methods before 2000, while the Antarctic ice sheet appears heavily in uenced by interannual variations

An iterative ICA-based reconstruction method to produce consistent time-variable total water storage fields using GRACE and swarm satellite data

Observing global terrestrial water storage changes (TWSCs) from (inter-)seasonal to (multi-)decade time-scales is very important to understand the Earth as a system under natural and anthropogenic climate change. The primary goal of the Gravity Recovery And Climate Experiment (GRACE) satellite mission (2002–2017) and its follow-on mission (GRACE-FO, 2018–onward) is to provide time-variable gravity fields, which can be converted to TWSCs with ∼300 km spatial resolution; however, the one year data gap between GRACE and GRACE-FO represents a critical discontinuity, which cannot be replaced by alternative data or model with the same quality. To fill this gap, we applied time-variable gravity fields (2013–onward) from the Swarm Earth explorer mission with low spatial resolution of ∼1500 km. A novel iterative reconstruction approach was formulated based on the independent component analysis (ICA) that combines the GRACE and Swarm fields. The reconstructed TWSC fields of 2003–2018 were compared with a commonly applied reconstruction technique and GRACE-FO TWSC fields, whose results indicate a considerable noise reduction and long-term consistency improvement of the iterative ICA reconstruction technique. They were applied to evaluate trends and seasonal mass changes (of 2003–2018) within the world’s 33 largest river basin

Characterization of the Morphometry of Impact Craters Hosting Polar Deposits in Mercury's North Polar Region

Earth-based radar images dating back two decades show that the floors of some polar craters on Mercury host radar-bright deposits that have been proposed to consist of frozen volatiles. Several hypotheses have been put forth to explain their source, including volcanic outgassing, chemical sputtering, and deposition of exogenous water ice. Calculations show that volatiles are thermally stable in permanently shadowed areas. An earlier study of the depths of north polar craters determined with photoclinometric techniques applied to Mariner 10 images yielded the conclusion that the mean ratio of crater depth d to rim-crest diameter D for craters hosting polar deposits is two-thirds that of the mean ratio for a comparable population of neighboring craters lacking such deposits. This result could be explained by (though doesn't require) the presence of a thick layer of volatiles within the polar deposit-hosting craters. Here we use altimetric profiles and topographic maps obtained by the Mercury Laser Altimeter (MLA) to revisit this analysis. MLA is an instrument on the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft, which has been orbiting Mercury since March 2011. MLA transmits a 1064-nm laser pulse at 8 Hz during MESSENGER's trajectory over Mercury s surface. The MLA illuminates surface areas averaging between 15 m and 100 m in diameter, spaced approx 400 m apart along the spacecraft ground track. The radial precision of individual measurements is <1 m, and the current accuracy with respect to Mercury s center of mass is better than 20 m. As of mid-December 2011, MLA coverage had reached to 15 S and has yielded a comprehensive map of the topography of Mercury s northern hemisphere. The MLA data are used here to quantify the shapes of craters in the north polar region and to avoid the shadowing bias of photoclinometric techniques

Mass balance of the Greenland Ice Sheet from 1992 to 2018

In recent decades, the Greenland Ice Sheet has been a major contributor to global sea-level rise1,2, and it is expected to be so in the future3. Although increases in glacier flow4–6 and surface melting7–9 have been driven by oceanic10–12 and atmospheric13,14 warming, the degree and trajectory of today’s imbalance remain uncertain. Here we compare and combine 26 individual satellite measurements of changes in the ice sheet’s volume, flow and gravitational potential to produce a reconciled estimate of its mass balance. Although the ice sheet was close to a state of balance in the 1990s, annual losses have risen since then, peaking at 335 ± 62 billion tonnes per year in 2011. In all, Greenland lost 3,800 ± 339 billion tonnes of ice between 1992 and 2018, causing the mean sea level to rise by 10.6 ± 0.9 millimetres. Using three regional climate models, we show that reduced surface mass balance has driven 1,971 ± 555 billion tonnes (52%) of the ice loss owing to increased meltwater runoff. The remaining 1,827 ± 538 billion tonnes (48%) of ice loss was due to increased glacier discharge, which rose from 41 ± 37 billion tonnes per year in the 1990s to 87 ± 25 billion tonnes per year since then. Between 2013 and 2017, the total rate of ice loss slowed to 217 ± 32 billion tonnes per year, on average, as atmospheric circulation favoured cooler conditions15 and as ocean temperatures fell at the terminus of Jakobshavn Isbræ16. Cumulative ice losses from Greenland as a whole have been close to the IPCC’s predicted rates for their high-end climate warming scenario17, which forecast an additional 50 to 120 millimetres of global sea-level rise by 2100 when compared to their central estimate

Extending the Record of Time-Variable Gravity by Combining Grace and Conventional Tracking Data

The Gravity Recovery and Climate Experiment (GRACE) mission has provided invaluable insight on the behavior of mass redistribution on the surface of the Earth since its launch in March 2002. GRACE measures the range between its two satellites to obtain a new, high-resolution, global gravity field of the Earth every month. This dissertation presents the first record of high-resolution time-variable gravity that extends to the 1990s. This gravity reconstruction is obtained by combining GRACE with gravity solutions generated from conventional tracking of satellites that start in November 1992. These tracking observations originate from two ground-based systems, Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellites (DORIS). The statistical framework for the reconstructions of the two fields can be described in three major steps. First, high-resolution spatial patterns are obtained from a Principal Component Analysis (PCA) decomposition of GRACE gravity coefficients. In the second step, the top four patterns are used as the model relating the SLR/DORIS coefficients to their temporal modulation, which are solved for in a Least Squares Adjustment. Lastly, the full patterns and their temporal modulation are linearly combined to form gravity fields that span the SLR/DORIS time frame and therefore add a decade of perspective before GRACE. An error budget is also designed to quantify the influence of assumptions made in the reconstruction. The primary scientific goal was to add perspective on the behavior of polar ice sheet melt in the 1990s. The second goal was to extend the record of terrestrial water storage. One last application examined the possibility of using this technique as a way to bridge the upcoming multi-month gap between GRACE and its successor mission, GRACE Follow-On

Investigation of regional variation in Lunar crater morphometry from (Lunar Orbiter Laser Altimeter) LOLA observations

Thesis: S.B., Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences, 2011.Cataloged from PDF version of thesis.Includes bibliographical references (pages 41-42).The advent of global Digital Elevation Models of the lunar surface, obtained from the Lunar Orbiter Laser Altimeter (LOLA), has allowed for a quantitative assessment of crater morphometry. 351 simple and complex craters in the Mare Serenitatis, far side highlands, near side highlands, and South Pole-Aitken basin are decomposed into 50 elevation profiles, from which key geometric crater properties are extracted. The geometric properties and their respective standard variation, such as height-to-diameter ratios, and average elevation profile are compared on a global level to investigate regional differences in terrain rheology and study the transition between the simple and complex crater regime. Furthermore, the relationship between known degradation mechanisms and crater morphometry is discussed, as well as the current state of quantitative methods to assess crater degradation. The resulting regional differences observed in crater morphometry are explained in the context of lunar geologic history. Finally, the addition of other crater geometric properties in future quantitative assessments will broaden the study of crater morphometry, and improvements to current methods are necessary to conclusively define degradation states in terms of quantitative factors.by Matthieu Jean Talpe.S.B

An Improved Ionospheric Correction Model for Grazing Angle GNSS-R Altimetry

peer reviewedIn this work, Spire Global, Inc. carrier phase of GNSS dual-frequency reflected signals under grazing angle (GG-R) are used for sea surface altimetry over the Java Sea, Indonesia. Processing of 714 GG-R coherent profiles from January 2020 to October 2021 is carried out to retrieve the relative surface height from WGS84 applying two methods to remove the ionospheric delay: an ionosphere-free combination and smoothing of the frequency-specific correction based on dual-frequency carrier phase observations. A substantial improvement of ∼ 8 cm is observed in the root mean square error (RMSE) between the GG-R retrievals and the reference surface model using the frequency-specific ionospheric correction method.13. Climate actio