C‐band SAR backscatter characteristics of Arctic sea and land ice during winter

Abstract

Synthetic Aperture Radar (SAR) data has become an important tool for studies of polar regions, due to high spatial resolution even during the polar night and under cloudy skies. We have studied the temporal variation of sea and land ice backscatter of twenty‐four SAR images from the European Remote Sensing satellite (ERS‐1) covering an area in Lady Ann Strait and Jones Sound, Nunavut, from January to March 1992. The presence of fast ice in Jones Sound and glaciers and ice caps on the surrounding islands provides an ideal setting for temporal backscatter studies of ice surfaces. Sample regions for eight different ice types were selected and the temporal backscatter variation was studied. The observed backscatter values for each ice type characterize the radar signatures of the ice surfaces. This time series of twenty‐four SAR images over a 3‐month period provides new insights into the degree of temporal variability of each surface. Ice caps exhibit the highest backscatter value of ‐3.9 dB with high temporal variability. Valley glacier ice backscatter values decrease with decreasing altitude, and are temporally the most stable, with standard deviations of 0.08–0.10 dB over the 90‐day period.

First‐year ice and lead ice show a negative trend in backscatter values in time and a positive correlation of up to 0.59 with air temperature over the 90‐day period. For first‐year ice and lead ice, episodes of large temperature fluctuations (±12°C) are associated with rapid changes in backscatter values (±2 dB). We attribute the backscatter increase to a temperature‐induced increase in brine volume at the base of the snow pack. Multi‐year ice, conglomerate ice and shore ice are relatively stable over the 3‐month period, with a backscatter variation of only a few dBs. An observed lag time of up to three days between backscatter increase/decrease and air temperature can be attributed to the insulation effect of the snow cover over sea ice. The net range of the backscatter values observed on the most temporally stable surface, valley glacier ice, of about 0.30 dB indicates that the ERS‐1 SAR instrument exceeds the 1 dB calibration accuracy specified for the Alaska SAR Facility processor for the three winter months.     

Predictability of large interannual Arctic sea-ice anomalies

In projections of twenty-first century climate, Arctic sea ice declines and at the same time exhibits strong interannual anomalies. Here, we investigate the potential to predict these strong sea-ice anomalies under a perfect-model assumption, using the Max-Planck-Institute Earth System Model in the same setup as in the Coupled Model Intercomparison Project Phase 5 (CMIP5). We study two cases of strong negative sea-ice anomalies: a 5-year-long anomaly for present-day conditions, and a 10-year-long anomaly for conditions projected for the middle of the twenty-first century. We treat these anomalies in the CMIP5 projections as the truth, and use exactly the same model configuration for predictions of this synthetic truth. We start ensemble predictions at different times during the anomalies, considering lagged-perfect and sea-ice-assimilated initial conditions. We find that the onset and amplitude of the interannual anomalies are not predictable. However, the further deepening of the anomaly can be predicted for typically 1 year lead time if predictions start after the onset but before the maximal amplitude of the anomaly. The magnitude of an extremely low summer sea-ice minimum is hard to predict: the skill of the prediction ensemble is not better than a damped-persistence forecast for lead times of more than a few months, and is not better than a climatology forecast for lead times of two or more years. Predictions of the present-day anomaly are more skillful than predictions of the mid-century anomaly. Predictions using sea-ice-assimilated initial conditions are competitive with those using lagged-perfect initial conditions for lead times of a year or less, but yield degraded skill for longer lead times. The results presented here suggest that there is limited prospect of predicting the large interannual sea-ice anomalies expected to occur throughout the twenty-first century.     

Spatial and temporal changes in dry-snow line altitude on the Antarctic Peninsula

Drastic changes were detected in glacial systems of the Antarctic Peninsula in the last decades. The observed phenomena comprise the disintegration of ice shelves, acceleration and thinning of glaciers, and retreat of glacier fronts. However, due to the lack of consistent systematic observations in particular of the higher parts of the glacial systems, it is difficult to predict further responses of the Antarctic Peninsula glaciers to climatic change. The present paper analyses spatial and temporal variations of changes in the dry-snow line altitude on the Antarctic Peninsula as extracted from a time series (1992–2005) of ERS-1/2 SAR and Envisat ASAR data. Upward changes in dry-snow line altitude were observed in general, and are attributed to extreme high-temperature events impacting the central plateaus of the Antarctic Peninsula and the increasing duration of warming periods. A mean decrease in dry-snow line altitude was detected on the west side of the peninsula and is identified as a response to recorded increase in precipitation and accumulation. These results validate the capability of SAR data for deriving superficial parameters of glaciers to be used as indicators of climatic changes in high-latitude regions where operational restrictions limit conventional meteorological observations.     

Grid-size dependence of Cauchy boundary conditions used to simulate stream–aquifer interactions

This work examines the simulation of stream–aquifer interactions as grids are refined vertically and horizontally and suggests that traditional methods for calculating conductance can produce inappropriate values when the grid size is changed. Instead, different grid resolutions require different estimated values. Grid refinement strategies considered include global refinement of the entire model and local refinement of part of the stream. Three methods of calculating the conductance of the Cauchy boundary conditions are investigated. Single- and multi-layer models with narrow and wide streams produced stream leakages that differ by as much as 122% as the grid is refined. Similar results occur for globally and locally refined grids, but the latter required as little as one-quarter the computer execution time and memory and thus are useful for addressing some scale issues of stream–aquifer interactions. Results suggest that existing grid-size criteria for simulating stream–aquifer interactions are useful for one-layer models, but inadequate for three-dimensional models. The grid dependence of the conductance terms suggests that values for refined models using, for example, finite difference or finite-element methods, cannot be determined from previous coarse-grid models or field measurements. Our examples demonstrate the need for a method of obtaining conductances that can be translated to different grid resolutions and provide definitive test cases for investigating alternative conductance formulations.     

One Million Years tephra record at IODP Sites U1436 and U1437: Insights into explosive volcanism from the Japan and Izu arcs

The 1 Myr tephra records of IODP (International Ocean Discovery Program) Holes U1436A and U1437B in the Izu‐Bonin fore‐ and reararc were investigated in order to assess provenance and eruptive volumes, respectively. In total, 304 tephra samples were examined and 260 primary tephra layers were identified. Tephra provenance was determined by means of major and trace element compositions of glass shards and distinguished between Japan and Izu‐Bonin arc origin of the tephra layers. A total of 33 marine tephra compositions were correlated to the Japan arc and 227 to the Izu arc. Twenty marine tephra layers were correlated between the two drilling sites. Additionally, we defined eleven correlations of marine tephra deposits to major widespread Japanese eruptions; from the 1.05 Ma Shishimuta‐Pink Tephra to the 30 ka Aira‐Tn Tephra, both from Kyushu Island. These eruptions provide independent time markers within the sediment record and six correlations were used to date tephra layers from Japan in Hole U1436A to establish an alternative age model for this hole. Furthermore, the minimum distal tephra volumes of all detected events were calculated, which enabled the comparison of the tephra volumes that derived from the Japan and the Izu‐Bonin arcs. For some of the major Japanese eruptions these are the first volume estimations that also include distal deposits. All of the Japanese tephras derived from events with eruption magnitude M_v ≥ 5.6 and three of the investigated eruptions reach magnitudes M_v ≥ 7. Volcanic events of the Izu‐Bonin arc have mostly eruption magnitudes M_v ≤ 5.     

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     

Karst spring responses examined by process-based modeling

Ground water in karst terrains is highly vulnerable to contamination due to the rapid transport of contaminants through the highly conductive conduit system. For contamination risk assessment purposes, information about hydraulic and geometric characteristics of the conduits and their hydraulic interaction with the fissured porous rock is an important prerequisite. The relationship between aquifer characteristics and short-term responses to recharge events of both spring discharge and physicochemical parameters of the discharged water was examined using a process-based flow and transport model. In the respective software, a pipe-network model, representing fast conduit flow, is coupled to MODFLOW, which simulates flow in the fissured porous rock. This hybrid flow model was extended to include modules simulating heat and reactive solute transport in conduits. The application of this modeling tool demonstrates that variations of physicochemical parameters, such as solute concentration and water temperature, depend to a large extent on the intensity and duration of recharge events and provide information about the structure and geometry of the conduit system as well as about the interaction between conduits and fissured porous rock. Moreover, the responses of solute concentration and temperature of spring discharge appear to reflect different processes, thus complementing each other in the aquifer characterization.     

Walking through volcanic mud: the 2,100-year-old Acahualinca footprints (Nicaragua)

We present the stratigraphy, lithology, volcanology, and age of the Acahualinca section in Managua, including a famous footprint layer exposed in two museum pits. The ca. 4-m-high walls of the main northern pit (Pit I) expose excellent cross sections of Late Holocene volcaniclastic deposits in northern Managua. We have subdivided the section into six lithostratigraphic units, some of which we correlate to Late Holocene eruptions. Unit I (1.2 m thick), chiefly of hydroclastic origin, begins with the footprint layer. The bulk is dominated by mostly massive basaltic-andesitic tephra layers, interpreted to represent separate pulses of a basically phreatomagmatic eruptive episode. We correlate these deposits based on compositional and stratigraphic evidence to the Masaya Triple Layer erupted at Masaya volcano ca. 2,120 ± 120 a B.P.. The eruption occurred during the dry season. A major erosional channel unconformity up to 1 m deep in the western half of Pit I separates Units II and I. Unit II begins with basal dacitic pumice lapilli up to 10 cm thick overlain by a massive to bedded fine-grained dacitic tuff including a layer of accretionary lapilli and pockets of well-rounded pumice lapilli. Angular nonvesicular glass shards are interpreted to represent hydroclastic fragmentation. The dacitic tephra is correlated unequivocally with the ca. 1.9-ka-Plinian dacitic Chiltepe eruption. Unit III, a lithified basaltic-andesitic deposit up to 50 cm thick and extremely rich in branch molds and excellent leaf impressions, is correlated with the Masaya Tuff erupted ca. 1.8 ka ago. Unit IV, a reworked massive basaltic-andesitic deposit, rich in brown tuff clasts and well bedded and cross bedded in the northwestern corner of Pit I, cuts erosionally down as far as Unit I. A poorly defined, pale brown mass flow deposit up to 1 m thick (Unit V) is overlain by 1–1.5 m of dominantly reworked, chiefly basaltic tephra topped by soil (Unit VI). A major erosional channel carved chiefly between deposition of Units II and I may have existed as a shallow drainage channel even prior to deposition of the footprint layer. The swath of the footprints is oriented NNW, roughly parallel to, and just east of, the axis of the channel. The interpretation of the footprint layer as the initial product of a powerful eruption at Masaya volcano followed without erosional breaks by additional layers of the same eruptive phase is strong evidence that the group of 15 or 16 people tried to escape from an eruption.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.