Tiansheng Hong, Yan Su, Elena Pettinelli, Roberto Orosei, Sebastian Emanuel Lauro, Zhibin Li, Chunyu Ding, Elisabetta Mattei, Barbara Cosciotti, Chunlai Li
The radar surface echo can be separated into coherent and incoherent components by statistical approaches, and the coherent component can be described by a backscattering model related to the RMS height. According to backscattering models for fractal surfaces, the coherent power in decibels decreases with RMS height on a scale independent of the wavelength at a rate depending on the Hurst exponent and the roughness scale. We extract the coherent power in four research areas by fitting the amplitude distribution of the Martian surface echoes recorded by the SHARAD radar, and compare the coherent power with the RMS height derived from pulse width of the MOLA laser altimeter. Scatter plots of squared MOLA-derived RMS height-coherent power are drawn to estimate the rates of coherent power fall-off by linear fitting, and the fitting power fall-off rates are compared to the Hurst exponents derived from digital terrain models in those areas. The fitting rates decrease with the Hurst exponent, similar to the theoretical rates. However, the fitting rates decrease with the Hurst exponent more sharply than the theoretical prediction. We explain the mismatch with a linear assumption between different roughness parameters, which helps to estimate the Hurst exponent, and a significant discrepancy between the wavelength and the roughness scale might influence the estimation results due to the scaling dependence of the Hurst exponent. This paper offers an opportunity to learn about the Hurst exponent at a tens-of-meter scale.
{"title":"Investigation on the Coherent Component of SHARAD Surface Echo and Surface Roughness: Scaling Behavior and Influence of Hurst Exponent","authors":"Tiansheng Hong, Yan Su, Elena Pettinelli, Roberto Orosei, Sebastian Emanuel Lauro, Zhibin Li, Chunyu Ding, Elisabetta Mattei, Barbara Cosciotti, Chunlai Li","doi":"10.1029/2024EA004081","DOIUrl":"https://doi.org/10.1029/2024EA004081","url":null,"abstract":"<p>The radar surface echo can be separated into coherent and incoherent components by statistical approaches, and the coherent component can be described by a backscattering model related to the RMS height. According to backscattering models for fractal surfaces, the coherent power in decibels decreases with RMS height on a scale independent of the wavelength at a rate depending on the Hurst exponent and the roughness scale. We extract the coherent power in four research areas by fitting the amplitude distribution of the Martian surface echoes recorded by the SHARAD radar, and compare the coherent power with the RMS height derived from pulse width of the MOLA laser altimeter. Scatter plots of squared MOLA-derived RMS height-coherent power are drawn to estimate the rates of coherent power fall-off by linear fitting, and the fitting power fall-off rates are compared to the Hurst exponents derived from digital terrain models in those areas. The fitting rates decrease with the Hurst exponent, similar to the theoretical rates. However, the fitting rates decrease with the Hurst exponent more sharply than the theoretical prediction. We explain the mismatch with a linear assumption between different roughness parameters, which helps to estimate the Hurst exponent, and a significant discrepancy between the wavelength and the roughness scale might influence the estimation results due to the scaling dependence of the Hurst exponent. This paper offers an opportunity to learn about the Hurst exponent at a tens-of-meter scale.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 2","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2024EA004081","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146139147","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
R. Tilling, A. Egido, J. Harbeck, N. Kurtz, A. Petty, A. Ridout, A. Shepherd, J. Wimert, D. Yi
Differences in satellite sampling affect their ability to resolve small-scale features over Arctic sea ice. For CryoSat-2 (CS2) and ICESat-2 (IS2) these differences are driven by geometric (footprint resolution) and radiometric (radar or laser) sampling. Here we compare growth season (October-April) surface type densities (the detected densities of lead, floe, and ambiguous targets) from CS2 and IS2 products, Arctic-wide over a common mission period. We develop these products using standard and fully-focused CS2 sea ice processing, IS2 ATL07 sea ice height data, and IS2 ATL10 sea ice freeboard data. Our analysis shows agreement in the spatial distributions of lead and floe densities between products, but significant variations in magnitude. Average floe densities from CS2 standard and fully-focused processing are 40% and 41% respectively, but 91% for all IS2 products. The average lead density from CS2 standard processing is 45%, and below 10% for all other products. The factors causing ambiguous classifications and misclassifications differ between satellites; while CS2 is more susceptible to off-nadir ranging to leads, IS2 retrievals are complicated by variable apparent lead brightness at nadir, and the presence of ridged ice. We also investigate the impact of sampling on sea ice floe length estimates, which average 1.8–2.9 km. Finally, we assess the performance of CS2 and IS2 surface type classification along near-coincident CRYO2ICE orbits. Based on our results we encourage CS2 and IS2 data users to consider how satellite sampling impacts true geophysical retrievals, and to utilize both missions simultaneously to benefit from their complementary strengths.
{"title":"Lead and Floe Detection From CryoSat-2 Radar and ICESat-2 Laser Altimetry","authors":"R. Tilling, A. Egido, J. Harbeck, N. Kurtz, A. Petty, A. Ridout, A. Shepherd, J. Wimert, D. Yi","doi":"10.1029/2025EA004399","DOIUrl":"https://doi.org/10.1029/2025EA004399","url":null,"abstract":"<p>Differences in satellite sampling affect their ability to resolve small-scale features over Arctic sea ice. For CryoSat-2 (CS2) and ICESat-2 (IS2) these differences are driven by geometric (footprint resolution) and radiometric (radar or laser) sampling. Here we compare growth season (October-April) surface type densities (the detected densities of lead, floe, and ambiguous targets) from CS2 and IS2 products, Arctic-wide over a common mission period. We develop these products using standard and fully-focused CS2 sea ice processing, IS2 ATL07 sea ice height data, and IS2 ATL10 sea ice freeboard data. Our analysis shows agreement in the spatial distributions of lead and floe densities between products, but significant variations in magnitude. Average floe densities from CS2 standard and fully-focused processing are 40% and 41% respectively, but 91% for all IS2 products. The average lead density from CS2 standard processing is 45%, and below 10% for all other products. The factors causing ambiguous classifications and misclassifications differ between satellites; while CS2 is more susceptible to off-nadir ranging to leads, IS2 retrievals are complicated by variable apparent lead brightness at nadir, and the presence of ridged ice. We also investigate the impact of sampling on sea ice floe length estimates, which average 1.8–2.9 km. Finally, we assess the performance of CS2 and IS2 surface type classification along near-coincident CRYO2ICE orbits. Based on our results we encourage CS2 and IS2 data users to consider how satellite sampling impacts true geophysical retrievals, and to utilize both missions simultaneously to benefit from their complementary strengths.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 2","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004399","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146140239","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The longest known Aletai meteorite belt presents a unique phenomenon in meteoroid dynamics. To investigate its formation mechanism, this study introduces a bilobate-shaped meteoroid model, emphasizing aerodynamic interactions and structure evolution. The sintered bond model is applied to simulate the tensile, compressive, and shear strengths of the bilobate-shaped meteoroid. Its disintegration is analyzed under the combined effects of aerodynamic forces and self-rotation. After disintegration, the transverse velocity of the sub-spherical fragments is applied to track their dispersal trajectories and calculate the resulting strewn field of meteorites. The influence of aerodynamical shock wave and mass ablation is considered throughout the descent process. Numerical simulations are conducted with varying initial entry conditions, particularly focusing on the initial rotation of the bilobate-shaped meteoroid. The study focuses on the mechanism of the skipping trajectory and the associated strewn field during the meteoroid's dynamical evolution. The results highlight the critical role of bilobate-shaped meteoroids in generating skipping trajectories and provide new insights into the formation of Aletai-like ultra-long meteorite belt.
{"title":"Disintegration and Skipping Dynamics of Bilobate-Shaped Meteoroids for Generating Ultra-Long Strewn Fields","authors":"Haoyu Li, Ziwen Li, Qingbo Gan, Xiangyuan Zeng","doi":"10.1029/2025EA004383","DOIUrl":"https://doi.org/10.1029/2025EA004383","url":null,"abstract":"<p>The longest known Aletai meteorite belt presents a unique phenomenon in meteoroid dynamics. To investigate its formation mechanism, this study introduces a bilobate-shaped meteoroid model, emphasizing aerodynamic interactions and structure evolution. The sintered bond model is applied to simulate the tensile, compressive, and shear strengths of the bilobate-shaped meteoroid. Its disintegration is analyzed under the combined effects of aerodynamic forces and self-rotation. After disintegration, the transverse velocity of the sub-spherical fragments is applied to track their dispersal trajectories and calculate the resulting strewn field of meteorites. The influence of aerodynamical shock wave and mass ablation is considered throughout the descent process. Numerical simulations are conducted with varying initial entry conditions, particularly focusing on the initial rotation of the bilobate-shaped meteoroid. The study focuses on the mechanism of the skipping trajectory and the associated strewn field during the meteoroid's dynamical evolution. The results highlight the critical role of bilobate-shaped meteoroids in generating skipping trajectories and provide new insights into the formation of Aletai-like ultra-long meteorite belt.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 2","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004383","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146140096","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
With the increasing demands for global development, deep space exploration missions are taking place more frequently, rendering the traditional 24-hr Polar Motion (PM) predictions inadequate for supporting the frequent operations of spacecraft. Therefore, this study utilized the mainstream EOP C04 and EOP C04 spliced with the final.daily series to evaluate the performance of interpolating PM to 6-hr and then forecasting, under both ideal and real-world environments. Interpolation is performed using the Stable Ultra-high-Precision-Radial-Basis-Function (SURBF) and the Seasonal-Trend-Residual decomposition with adaptive interpolation (STR). Prediction is carried out using the Least Squares and Autoregressive (LS + AR). Our findings indicate that the two input formats, namely interpolating prior to column assignment in the idealized case and direct column assignment before interpolation in the real-world environments, greatly affect prediction results. The former takes full advantage of the disclosed subsequent information, achieving a first-day MAE of less than 0.1mas. The latter adheres to the actual conditions, resulting in a first-day MAE above 0.2mas and performing worse than traditional forecasting. Additionally, we evaluated the forecasting performance using 6-hr International GNSS Service Ultra-rapid (IGU) inputs and found that it not only surpasses traditional methods but also outperforms USNO forecasts over the 1–10 days. Therefore, we do not recommend interpolating 24-hr to 6-hr, as interpolation is essentially a guess and does not provide meaningful information. This study only recommends using IGU-6h derived from real observational data to obtain 6-hr PM predictions. These findings can assist practical deep space exploration projects in making informed decisions when selecting 6-hr PM predictions.
{"title":"Evaluate the Prediction Level of Polar Motion Interpolation for 6 hr in Ideal and Real Environments","authors":"Wei Miao, Xueqing Xu, Yonghong Zhou","doi":"10.1029/2025EA004866","DOIUrl":"https://doi.org/10.1029/2025EA004866","url":null,"abstract":"<p>With the increasing demands for global development, deep space exploration missions are taking place more frequently, rendering the traditional 24-hr Polar Motion (PM) predictions inadequate for supporting the frequent operations of spacecraft. Therefore, this study utilized the mainstream EOP C04 and EOP C04 spliced with the final.daily series to evaluate the performance of interpolating PM to 6-hr and then forecasting, under both ideal and real-world environments. Interpolation is performed using the Stable Ultra-high-Precision-Radial-Basis-Function (SURBF) and the Seasonal-Trend-Residual decomposition with adaptive interpolation (STR). Prediction is carried out using the Least Squares and Autoregressive (LS + AR). Our findings indicate that the two input formats, namely interpolating prior to column assignment in the idealized case and direct column assignment before interpolation in the real-world environments, greatly affect prediction results. The former takes full advantage of the disclosed subsequent information, achieving a first-day MAE of less than 0.1mas. The latter adheres to the actual conditions, resulting in a first-day MAE above 0.2mas and performing worse than traditional forecasting. Additionally, we evaluated the forecasting performance using 6-hr International GNSS Service Ultra-rapid (IGU) inputs and found that it not only surpasses traditional methods but also outperforms USNO forecasts over the 1–10 days. Therefore, we do not recommend interpolating 24-hr to 6-hr, as interpolation is essentially a guess and does not provide meaningful information. This study only recommends using IGU-6h derived from real observational data to obtain 6-hr PM predictions. These findings can assist practical deep space exploration projects in making informed decisions when selecting 6-hr PM predictions.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004866","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146057951","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
M. Hörlöck, B. Heber, S. Jensen, P. Kühl, H. Sierks
Context. The Electron Proton Helium INstrument (EPHIN) aboard SOlar and Heliospheric Observatory is a particle telescope that measures energetic protons and helium above 4 MeV/nuc and electrons above 150 keV. While a calibration of EPHIN has been performed before launch, it has only partially been assessed reducing the accuracy of the detector and dead layer thicknesses of each detector. However, these parameters are crucial for the correct representation of the instrument in simulation runs which are necessary for detailed analysis of the instruments measurements in space. Aims. An accurate representation of EPHINs detector geometry has been derived from the calibration results. These improvements compared to previous detector models will allow for lower systematic uncertainties as well as new data products for the EPHIN data at high energies. Methods. EPHIN was calibrated using a 4He beam aimed at a gold target. The various particle populations produced in the target were then filtered using a rigidity filter. The measured energy loss distributions in each individual detector have been compared to GEANT4 simulations. The thicknesses and dead layer thicknesses of each detector have been varied in the simulation setup to achieve best agreements between simulation and calibration. Results. Good agreements between simulation and calibration have been reached. Most energy loss distribution fall within a 3 % margin, compared to differences of more than 10 % in earlier models.
{"title":"Development and Validation of the GEANT4 Model of the Electron Proton Helium Instrument on Board the SOHO Spacecraft Utilizing the Accelerator Measurements at the Hahn-Meitner Institute Berlin","authors":"M. Hörlöck, B. Heber, S. Jensen, P. Kühl, H. Sierks","doi":"10.1029/2025EA004713","DOIUrl":"https://doi.org/10.1029/2025EA004713","url":null,"abstract":"<p>Context. The Electron Proton Helium INstrument (EPHIN) aboard SOlar and Heliospheric Observatory is a particle telescope that measures energetic protons and helium above 4 MeV/nuc and electrons above 150 keV. While a calibration of EPHIN has been performed before launch, it has only partially been assessed reducing the accuracy of the detector and dead layer thicknesses of each detector. However, these parameters are crucial for the correct representation of the instrument in simulation runs which are necessary for detailed analysis of the instruments measurements in space. Aims. An accurate representation of EPHINs detector geometry has been derived from the calibration results. These improvements compared to previous detector models will allow for lower systematic uncertainties as well as new data products for the EPHIN data at high energies. Methods. EPHIN was calibrated using a 4He beam aimed at a gold target. The various particle populations produced in the target were then filtered using a rigidity filter. The measured energy loss distributions in each individual detector have been compared to GEANT4 simulations. The thicknesses and dead layer thicknesses of each detector have been varied in the simulation setup to achieve best agreements between simulation and calibration. Results. Good agreements between simulation and calibration have been reached. Most energy loss distribution fall within a 3 % margin, compared to differences of more than 10 % in earlier models.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004713","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146058062","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
J. V. Clark, R. S. Jakubek, P. Conrad, E. Cardarelli, D. Buckner, A. Corpolongo, W. J. Abbey, K. Hand, M. Fries, P. Casbeer, S. Siljeström, F. M. McCubbin, S. Sharma, C. Lee, R. Bhartia, A. Steele
Oxychlorines (i.e., perchlorates (ClO4−) and chlorates (ClO3−)) have been detected by several landed missions on Mars at various locations. These missions have provided crucial information about the geographic distribution and abundances of oxychlorines on Mars but have not definitively identified the cation and anion type of in situ, solid oxychlorines. By speciating and precisely locating oxychlorines in Martian rocks, we may be able to better interpret the aqueous history of the rocks, understand the chlorine cycle on Mars, understand the chlorine isotope systematics on Mars, identify the potential for liquid brines on the surface, and advance in situ resource utilization activities for future robotic or landed missions. The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument on the Mars 2020 Perseverance rover may have the ability to identify oxychlorine species (i.e., cation and anion) in solid surface targets due to their characteristic Raman and fluorescence bands. Additionally, the location and distribution of oxychlorines within rocks can be determined using the Autofocus and Context Imager (ACI) or Wide Angle Topographic Sensor for Operations and eNgineering (WATSON) camera, a subsystem of SHERLOC that provides high-resolution, close-up images of targets analyzed by the SHERLOC Deep-Ultraviolet Raman spectrometer. The overarching goal of this work is to test SHERLOC's ability to identify and differentiate oxychlorine species in synthetic pure and natural mixed samples using a laboratory analog to the SHERLOC Raman and fluorescence spectrometer, identify instrumental limitations, and to further constrain potential detections made within Jezero crater, Mars.
{"title":"Deep-Ultraviolet (DUV) Raman and Fluorescence of Perchlorates and Chlorates: Implications for the Detection of Oxychlorines by the SHERLOC Instrument in Jezero Crater, Mars","authors":"J. V. Clark, R. S. Jakubek, P. Conrad, E. Cardarelli, D. Buckner, A. Corpolongo, W. J. Abbey, K. Hand, M. Fries, P. Casbeer, S. Siljeström, F. M. McCubbin, S. Sharma, C. Lee, R. Bhartia, A. Steele","doi":"10.1029/2025EA004512","DOIUrl":"https://doi.org/10.1029/2025EA004512","url":null,"abstract":"<p>Oxychlorines (i.e., perchlorates (ClO<sub>4</sub><sup>−</sup>) and chlorates (ClO<sub>3</sub><sup>−</sup>)) have been detected by several landed missions on Mars at various locations. These missions have provided crucial information about the geographic distribution and abundances of oxychlorines on Mars but have not definitively identified the cation and anion type of in situ, solid oxychlorines. By speciating and precisely locating oxychlorines in Martian rocks, we may be able to better interpret the aqueous history of the rocks, understand the chlorine cycle on Mars, understand the chlorine isotope systematics on Mars, identify the potential for liquid brines on the surface, and advance in situ resource utilization activities for future robotic or landed missions. The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument on the Mars 2020 <i>Perseverance</i> rover may have the ability to identify oxychlorine species (i.e., cation and anion) in solid surface targets due to their characteristic Raman and fluorescence bands. Additionally, the location and distribution of oxychlorines within rocks can be determined using the Autofocus and Context Imager (ACI) or Wide Angle Topographic Sensor for Operations and eNgineering (WATSON) camera, a subsystem of SHERLOC that provides high-resolution, close-up images of targets analyzed by the SHERLOC Deep-Ultraviolet Raman spectrometer. The overarching goal of this work is to test SHERLOC's ability to identify and differentiate oxychlorine species in synthetic pure and natural mixed samples using a laboratory analog to the SHERLOC Raman and fluorescence spectrometer, identify instrumental limitations, and to further constrain potential detections made within Jezero crater, Mars.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004512","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146057928","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
S. Scher, F. Ladstädter, M. Schwärz, J. Innerkofler, G. Kirchengast
Radio occultation is a well-established remote sensing method that provides reliable estimates of atmospheric profiles of diverse variables, including temperature and pressure. However, as with all indirect methods, radio occultation has some inherent systematic and random error effects, which lead to observational uncertainties. While propagation of uncertainties along the processing chain for individual radio occultation profiles was described in recent studies, this uncertainty information has not yet been carried forward to climatological fields. We close this gap and present an uncertainty propagation procedure that provides uncertainty estimates for aggregated means for climate applications. Estimated random uncertainties, basic and apparent systematic uncertainties and sampling uncertainties (due to the discrete sampling by profiles) are propagated through the aggregation process, resulting in uncertainty estimates for gridded fields. We demonstrate the new procedure for two test months and representative variables, inspecting monthly mean profiles for refractivity, dry temperature and physical temperature measurements. Results show that estimated random uncertainties and residual sampling uncertainties (after sampling bias correction) have similar magnitudes, both decreasing with increasing spatial aggregation sizes and corresponding increasing number of aggregated observations. At small aggregation they are the main contributors to uncertainty in refractivity, and important contributors to uncertainty of temperature. Systematic uncertainty, whose magnitude is independent of the number of profiles, is for refractivity the main source of uncertainty for larger aggregation sizes, and for pressure and dry temperature at all commonly used aggregation sizes. All uncertainty components exhibit pronounced spatial variation over the globe, with polar regions showing the greatest uncertainty.
{"title":"Uncertainty Propagation From Radio Occultation Profiles to Aggregated Atmospheric Gridded Fields","authors":"S. Scher, F. Ladstädter, M. Schwärz, J. Innerkofler, G. Kirchengast","doi":"10.1029/2025EA004389","DOIUrl":"https://doi.org/10.1029/2025EA004389","url":null,"abstract":"<p>Radio occultation is a well-established remote sensing method that provides reliable estimates of atmospheric profiles of diverse variables, including temperature and pressure. However, as with all indirect methods, radio occultation has some inherent systematic and random error effects, which lead to observational uncertainties. While propagation of uncertainties along the processing chain for individual radio occultation profiles was described in recent studies, this uncertainty information has not yet been carried forward to climatological fields. We close this gap and present an uncertainty propagation procedure that provides uncertainty estimates for aggregated means for climate applications. Estimated random uncertainties, basic and apparent systematic uncertainties and sampling uncertainties (due to the discrete sampling by profiles) are propagated through the aggregation process, resulting in uncertainty estimates for gridded fields. We demonstrate the new procedure for two test months and representative variables, inspecting monthly mean profiles for refractivity, dry temperature and physical temperature measurements. Results show that estimated random uncertainties and residual sampling uncertainties (after sampling bias correction) have similar magnitudes, both decreasing with increasing spatial aggregation sizes and corresponding increasing number of aggregated observations. At small aggregation they are the main contributors to uncertainty in refractivity, and important contributors to uncertainty of temperature. Systematic uncertainty, whose magnitude is independent of the number of profiles, is for refractivity the main source of uncertainty for larger aggregation sizes, and for pressure and dry temperature at all commonly used aggregation sizes. All uncertainty components exhibit pronounced spatial variation over the globe, with polar regions showing the greatest uncertainty.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004389","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146002270","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Elisa Bozzolan, Elisa Matteligh, Andrea Brenna, Martina Cecchetto, Nicola Surian, Patrice Carbonneau, Simone Bizzi
The active channel of alluvial rivers delineates areas of geomorphic activity over a defined time window. While increasing satellite data availability enables monthly active channel delineations, multi-year analyses often rely on temporal aggregates (e.g., annual medians) to reduce computational costs and intra-annual variability. The potential of monthly information to improve active channels delineation and geomorphic interpretation remains largely unexplored. In this work, we delineated active channels for the Po River (Italy) by aggregating monthly Sentinel-2 classifications of river water and sediment bars into annual frequency maps at 10 m resolution. Annual aggregation mitigated monthly sediment underestimation (12%) but also amplified model overestimation biases (15%). Monthly classification persistence (e.g., classified as active channel for more than N months/year) was then used to reduce these errors and produce active channel areas that closely match those manually delineated from 30 cm orthophotos. The spatiotemporal variability of monthly classifications also show that the active channel area of dynamic reaches can vary ∼50% over the year. These changes revealed areas most prone to water-stage fluctuations, sediment transport, as well as zones seasonally or progressively colonized by vegetation—patterns hidden in single orthophotos or annual medians. Less dynamic reaches, by contrast, showed minimal differences between annual and monthly-based delineation methods. These findings emphasize the importance of adapting temporal aggregation to the river type and process analysed, with sub-annual resolutions better capturing, in dynamic rivers, seasonal and progressive active channel reconfigurations, along with their interaction with sediment and vegetation dynamics.
{"title":"Enhancing Active Channel Delineation in Alluvial Rivers Using Monthly Aggregation of Sentinel-2 Imagery","authors":"Elisa Bozzolan, Elisa Matteligh, Andrea Brenna, Martina Cecchetto, Nicola Surian, Patrice Carbonneau, Simone Bizzi","doi":"10.1029/2025EA004642","DOIUrl":"https://doi.org/10.1029/2025EA004642","url":null,"abstract":"<p>The active channel of alluvial rivers delineates areas of geomorphic activity over a defined time window. While increasing satellite data availability enables monthly active channel delineations, multi-year analyses often rely on temporal aggregates (e.g., annual medians) to reduce computational costs and intra-annual variability. The potential of monthly information to improve active channels delineation and geomorphic interpretation remains largely unexplored. In this work, we delineated active channels for the Po River (Italy) by aggregating monthly Sentinel-2 classifications of river water and sediment bars into annual frequency maps at 10 m resolution. Annual aggregation mitigated monthly sediment underestimation (12%) but also amplified model overestimation biases (15%). Monthly classification persistence (e.g., classified as active channel for more than N months/year) was then used to reduce these errors and produce active channel areas that closely match those manually delineated from 30 cm orthophotos. The spatiotemporal variability of monthly classifications also show that the active channel area of dynamic reaches can vary ∼50% over the year. These changes revealed areas most prone to water-stage fluctuations, sediment transport, as well as zones seasonally or progressively colonized by vegetation—patterns hidden in single orthophotos or annual medians. Less dynamic reaches, by contrast, showed minimal differences between annual and monthly-based delineation methods. These findings emphasize the importance of adapting temporal aggregation to the river type and process analysed, with sub-annual resolutions better capturing, in dynamic rivers, seasonal and progressive active channel reconfigurations, along with their interaction with sediment and vegetation dynamics.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004642","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146002065","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
P. Bountzis, E. Lippiello, S. Baccari, G. Petrillo
In the widely adopted description of seismic occurrence, earthquakes are categorized as either background or triggered events. In this work, we present a fully automated, non-parametric algorithm for distinguishing between these two categories, a process known as seismic declustering, based on the widely used nearest-neighbor (NN) metric. We introduce a new measure, the susceptibility index, which identifies an optimal threshold to discriminate between background and triggered events within the NN metric. Through statistical testing on simulated epidemic type aftershock sequence catalogs, we demonstrate that our method yields classification metrics exceeding 90%, outperforming state-of-the art algorithms. Notably, we show that a single threshold is sufficient for reliable discrimination within a given data set. The identification of this threshold requires memory capacity and computational time that scale linearly and quadratically with the data set size, respectively, making the method particurarly suited for large earthquake catalogs. We also apply our method to the relocated Southern California catalog and the GeoNet catalog of New Zealand (NZ). Our method effectively adapts across the different tectonic settings, capturing the variability of background seismicity rates between the shallow crustal events of Southern California and the tectonically diverse seismicity of NZ.
{"title":"Automatic Earthquake Declustering Using the Nearest-Neighbor Distance","authors":"P. Bountzis, E. Lippiello, S. Baccari, G. Petrillo","doi":"10.1029/2025EA004539","DOIUrl":"https://doi.org/10.1029/2025EA004539","url":null,"abstract":"<p>In the widely adopted description of seismic occurrence, earthquakes are categorized as either background or triggered events. In this work, we present a fully automated, non-parametric algorithm for distinguishing between these two categories, a process known as seismic declustering, based on the widely used nearest-neighbor (NN) metric. We introduce a new measure, the susceptibility index, which identifies an optimal threshold to discriminate between background and triggered events within the NN metric. Through statistical testing on simulated epidemic type aftershock sequence catalogs, we demonstrate that our method yields classification metrics exceeding 90%, outperforming state-of-the art algorithms. Notably, we show that a single threshold is sufficient for reliable discrimination within a given data set. The identification of this threshold requires memory capacity and computational time that scale linearly and quadratically with the data set size, respectively, making the method particurarly suited for large earthquake catalogs. We also apply our method to the relocated Southern California catalog and the GeoNet catalog of New Zealand (NZ). Our method effectively adapts across the different tectonic settings, capturing the variability of background seismicity rates between the shallow crustal events of Southern California and the tectonically diverse seismicity of NZ.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004539","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145996669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The accurate co-registration of geospatial data is necessary to answer questions that cross-cut disciplines and are key to understanding fundamental questions about our Solar System. To address this need and provide an updated product for Mars that is tied to a common reference frame, we have photogrammetrically controlled Thermal Emission Imaging System (THEMIS) daytime and nighttime infrared IR images. Using this improved image position knowledge, we generated orthorectified daytime and nighttime IR mosaics of Mars at 100 m per pixel for the ±65° latitude region of Mars. The updated spacecraft position and pointing information for the images is also released as SPICE kernels. The co-registration between individual THEMIS images achieves sub-pixel precision, and the average accuracy with which we know the position of any feature within the THEMIS controlled products is approximately 200 m horizontally. A globally controlled image set, with quantified accuracy and precision, is necessary to facilitate exploration and discovery for all bodies in the Solar System. Controlling THEMIS data allows multi-instrument science to be performed with significantly higher confidence as precise co-registration, and the accuracy knowledge of that registration, is necessary for analyses designed to extract information from the subtle differences between multiple images. A global image mosaic of Mars where uncertainties in the absolute image position are well characterized serves a wide range of purposes, including landing site evaluations, providing an accurate base to which high-resolution images (e.g., CTX and HiRISE) can be tied, and enables the fusion of multiple data types within a single framework.
{"title":"The THEMIS Control Network of Mars","authors":"R. L. Fergason, L. Weller, M. T. Bland","doi":"10.1029/2025EA004758","DOIUrl":"https://doi.org/10.1029/2025EA004758","url":null,"abstract":"<p>The accurate co-registration of geospatial data is necessary to answer questions that cross-cut disciplines and are key to understanding fundamental questions about our Solar System. To address this need and provide an updated product for Mars that is tied to a common reference frame, we have photogrammetrically controlled Thermal Emission Imaging System (THEMIS) daytime and nighttime infrared IR images. Using this improved image position knowledge, we generated orthorectified daytime and nighttime IR mosaics of Mars at 100 m per pixel for the ±65° latitude region of Mars. The updated spacecraft position and pointing information for the images is also released as SPICE kernels. The co-registration between individual THEMIS images achieves sub-pixel precision, and the average accuracy with which we know the position of any feature within the THEMIS controlled products is approximately 200 m horizontally. A globally controlled image set, with quantified accuracy and precision, is necessary to facilitate exploration and discovery for all bodies in the Solar System. Controlling THEMIS data allows multi-instrument science to be performed with significantly higher confidence as precise co-registration, and the accuracy knowledge of that registration, is necessary for analyses designed to extract information from the subtle differences between multiple images. A global image mosaic of Mars where uncertainties in the absolute image position are well characterized serves a wide range of purposes, including landing site evaluations, providing an accurate base to which high-resolution images (e.g., CTX and HiRISE) can be tied, and enables the fusion of multiple data types within a single framework.</p>","PeriodicalId":54286,"journal":{"name":"Earth and Space Science","volume":"13 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025EA004758","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146007469","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}