JPEG2000 Image Compression on Solar EUV Images

For future solar missions as well as ground-based telescopes, efficient ways to return and process data have become increasingly important. Solar Orbiter, which is the next ESA/NASA mission to explore the Sun and the heliosphere, is a deep-space mission, which implies a limited telemetry rate that makes efficient onboard data compression a necessity to achieve the mission science goals. Missions like the Solar Dynamics Observatory (SDO) and future ground-based telescopes such as the Daniel K. Inouye Solar Telescope, on the other hand, face the challenge of making petabyte-sized solar data archives accessible to the solar community. New image compression standards address these challenges by implementing efficient and flexible compression algorithms that can be tailored to user requirements. We analyse solar images from the Atmospheric Imaging Assembly (AIA) instrument onboard SDO to study the effect of lossy JPEG2000 (from the Joint Photographic Experts Group 2000) image compression at different bitrates. To assess the quality of compressed images, we use the mean structural similarity (MSSIM) index as well as the widely used peak signal-to-noise ratio (PSNR) as metrics and compare the two in the context of solar EUV images. In addition, we perform tests to validate the scientific use of the lossily compressed images by analysing examples of an on-disc and off-limb coronal-loop oscillation time-series observed by AIA/SDO.     

The sedimentology of an ephemeral fluvial–aeolian succession

Ephemeral fluvial systems are commonly associated with arid to semi-arid climates. Although their complex sedimentology and depositional settings have been described in much detail, depositional models depicting detailed lateral and vertical relationships, and interactions with coeval depositional environments, are lacking compared to well-recognized meandering and braided fluvial systems. This study critically evaluates the applicability of current models for ephemeral fluvial systems to an ancient arid fluvial example of the Lower Jurassic Kayenta Formation of the Colorado Plateau, USA. The study employs detailed sedimentary logging, palaeocurrent analysis and photogrammetric panels across the regional extent of the Kayenta. A generic model that accounts for the detailed sedimentology of a sandy arid ephemeral fluvial system (drawing upon both ancient and geomorphological studies) is developed, along with analysis of the spatial and temporal interactions with the aeolian setting. Results show that the ephemeral system is dominated by laterally and vertically amalgamated, poorly channelized to sheet-like elements, with abundant upper flow regime flat beds and high sediment load structures formed between periods of lower flow regime conditions. Through interaction with a coeval aeolian system, most of the fluvial deposits are dominated by sand-grade sediment, unlike many modern ephemeral fluvial systems that contain a high proportion of conglomeratic and/or finer grained mudstone and siltstone deposits. During dominantly fluvial deposition, high width to thickness ratios are observed for channelized and sheet-like elements. However, with increasing aridity, the aeolian environment becomes dominant and fluvial deposition is restricted to interdune corridors, resulting in lower width to thickness ratio channels dominated by flash-flood and debris-flow facies. The data presented here, coupled with modern examples of ephemeral systems and flood regimes, suggest that ephemeral flow produces and preserves distinctive sedimentological traits that can not only be recognized in outcrops, but also within core.     

Numerical implementation and sensitivity analysis of a wave energy converter in a time-dependent mild-slope equation model

Several Wave Energy Converters (abbreviated as WECs) have intensively been studied and developed during the last decade and currently small farms of WECs are getting installed. WECs in a farm are partly absorbing, partly redistributing the incident wave power. Consequently, the power absorption of each individual WEC in a farm is affected by its neighbouring WECs. The knowledge of the wave climate around the WEC is needed to predict its performance in the farm. In this paper a technique is developed to implement a single and multiple WECs based on the overtopping principle in a time-dependent mild-slope equation model. So far, the mild-slope equations have been widely used to study wave transformations around coastal and offshore structures, such as breakwaters, piles of windmills and offshore platforms. First the limitations of the WEC implementation are discussed through a sensitivity analysis. Next the developed approach is applied to study the wave height reduction behind a single WEC and a farm. The wake behind an isolated WEC is investigated for uni- and multidirectional waves; it is observed that an increase of the directional spread leads to a faster wave redistribution behind the WEC. Further the wake in the lee of multiple WECs is calculated for two different farm lay-outs, i.e. an aligned grid and a staggered grid, by adapting the performance of each WEC to its incident wave power. The evolved technique is a fast tool to find the optimal lay-out of WECs in a farm and to study the possible influence on surrounding activities in the sea.     

Application of Large‐Scale Inversion Algorithms to Hydraulic Tomography in an Alluvial Aquifer

Large‐scale inversion methods have been recently developed and permitted now to considerably reduce the computation time and memory needed for inversions of models with a large amount of parameters and data. In this work, we have applied a deterministic geostatistical inversion algorithm to a hydraulic tomography investigation conducted in an experimental field site situated within an alluvial aquifer in Southern France. This application aims to achieve a 2‐D large‐scale modeling of the spatial transmissivity distribution of the site. The inversion algorithm uses a quasi‐Newton iterative process based on a Bayesian approach. We compared the results obtained by using three different methodologies for sensitivity analysis: an adjoint‐state method, a finite‐difference method, and a principal component geostatistical approach (PCGA). The PCGA is a large‐scale adapted method which was developed for inversions with a large number of parameters by using an approximation of the covariance matrix, and by avoiding the calculation of the full Jacobian sensitivity matrix. We reconstructed high‐resolution transmissivity fields (composed of up to 25,600 cells) which generated good correlations between the measured and computed hydraulic heads. In particular, we show that, by combining the PCGA inversion method and the hydraulic tomography method, we are able to substantially reduce the computation time of the inversions, while still producing high‐quality inversion results as those obtained from the other sensitivity analysis methodologies.     

A new mathematical model of asymmetric hydraulic fracture growth

Hydraulic fractures generated by fluid injection in rock formations are often mapped by seismic monitoring. In many cases, the microseismicity is asymmetric relative to the injection well, which has been interpreted by stress gradient along the direction of the hydraulic fracture. We present a mathematical model of asymmetric hydrofracture growth based on relations between the solid‐phase stress and the fracture hydraulics. For single fracture and single injection point, the model has three parameters, hydraulic conductivities of the fracture wings, and normalised stress gradient and predicts the positions of the fracture tips as functions of time. The model is applied to a set of microseismic event locations that occurred during and after an injection process. Two different methods are suggested that make it possible to delineate the fracture tips from the set of microseismic events. This makes it possible to determine the model parameters and to check the agreement between the model prediction and the measured data. The comparison of the measured and modelled growth of fracture wings supports both the assumption of the non‐zero stress gradient and the existence of the post‐injection unilateral growth.     

Validation of official erosion modelling based on high‐resolution radar rain data by aerial photo erosion classification

The universal soil loss equation (USLE) is the most frequently applied erosion prediction model and it is also implemented as an official decision‐making instrument for agricultural regulations. The USLE itself has been already validated using different approaches. Additional errors, however, arise from input data and interpolation procedures that become necessary for field‐specific predictions on a national scale for administrative purposes. In this study, predicted event soil loss using the official prediction system in Bavaria (Germany) was validated by comparison with aerial photo erosion classifications of 8100 fields. Values for the USLE factors were mainly taken from the official Bavarian high‐resolution (5 × 5 m²) erosion cadastre. As series of erosion events were examined, the cover and management factor was replaced by the soil loss ratio. The event erosivity factor was calculated from high‐resolution (1 × 1 km², 5 min), rain gauge‐adjusted radar rain data (RADOLAN). Aerial photo erosion interpretation worked sufficiently well and average erosion predictions and visual classifications correlated closely. This was also true for data broken down to individual factors and different crops. There was no reason to assume a general invalidity of the USLE and the official parametrization procedures. Event predictions mainly suffered from errors in the assumed crop stage period and tillage practices, which do not reflect interannual and farm‐specific variation. In addition, the resolution of radar data (1 km²) did not seem to be sufficient to predict short‐term erosion on individual fields given the strong spatial gradients within individual rains. The quality of the input data clearly determined prediction quality. Differences between USLE predictions and observations are most likely caused by parametrization weaknesses but not by a failure of the model itself. Copyright © 2017 John Wiley & Sons, Ltd.     

Sulfur Hexafluoride and Potassium Bromide as Groundwater Tracers for Managed Aquifer Recharge

Sulfur hexafluoride (SF₆) is an established tracer for use in managed aquifer recharge projects. SF₆ exsolves from groundwater when it encounters trapped air according to Henry's law. This results in its retardation relative to groundwater flow, which can help determine porous media saturation and flow dynamics. SF₆ and the conservative, nonpartitioning tracer, bromide (Br⁻ added as KBr), were introduced to recharge water infiltrated into stacked glacial aquifers in Thurston County, Washington, providing the opportunity to observe SF₆ partitioning. Br⁻, which is assumed to travel at the same velocity as the groundwater, precedes SF₆ at most monitoring wells (MWs). Average groundwater velocity in the unconfined aquifer in the study area ranges from 3.9 to 40 m/d, except in the southwestern corner where it is slower. SF₆ in the shallow aquifer exhibits an average retardation factor of 2.5 ± 3.8, suggesting an air-to-water ratio on the order of 10⁻³ to 10⁻² in the pore space. Notable differences in tracer arrival times at adjacent wells indicate very heterogeneous conductivity. One MW exhibits double peaks in concentrations of both tracers with different degrees of retardation for the first and second peaks. This suggests multiple flowpaths to the well with variable saturation. The confining layer between the upper two aquifers appears to allow intermittent connection between aquifers but serves as an aquitard in most areas. This study demonstrates the utility of SF₆ partitioning for evaluating hydrologic conditions at prospective recharge sites.     

The Holme Pierrepont sand and gravel and the timing of Middle and Late Devensian floodplain aggradation in the English Midlands

This paper provides the first detailed, multi-proxy environmental record for the character of Lateglacial conditions in the lowest Pleistocene terrace of the main valley floor of the River Trent at Holme Pierrepont near Nottingham. The analysis of pollen, plant, insect and mollusc remains preserved within organic channels near the base of the terrace, named the Holme Pierrepont Sand and Gravel by the British Geological Survey (historically known as the Floodplain Terrace), provided evidence of cool, though not fully arctic climatic conditions and a largely treeless landscape, roamed by large herbivores. Radiocarbon dating indicates that these sediments were deposited during the Loch Lomond Stadial (Younger Dryas GS-1). Comparison of these dates from Holme Pierrepont with those from morphostratigraphically similar deposits in the wider Trent catchment suggests that the Holme Pierrepont Sand and Gravel may have been laid down in two separate pulses of braidplain aggradation either side of the ‘Last Glacial Maximum’. It has been demonstrated from historical documentation and riverine archaeological evidence that the middle Trent has been particularly sensitive to changing flood frequency and magnitude associated with climatic oscillations during the late Holocene; this study demonstrates that such sensitivity appears to extend back into the late Pleistocene. The timing of fluvial aggradation recorded at Holme Pierrepont agrees broadly with that recorded from other sites across England and north-west Europe.     

Quantifying uncertainties in projections of extremes��a perturbed land surface parameter experiment

Uncertainties in the climate response to a doubling of atmospheric CO₂ concentrations are quantified in a perturbed land surface parameter experiment. The ensemble of 108 members is constructed by systematically perturbing five poorly constrained land surface parameters of global climate model individually and in all possible combinations. The land surface parameters induce small uncertainties at global scale, substantial uncertainties at regional and seasonal scale and very large uncertainties in the tails of the distribution, the climate extremes. Climate sensitivity varies across the ensemble mainly due to the perturbation of the snow albedo parameterization, which controls the snow albedo feedback strength. The uncertainty range in the global response is small relative to perturbed physics experiments focusing on atmospheric parameters. However, land surface parameters are revealed to control the response not only of the mean but also of the variability of temperature. Major uncertainties are identified in the response of climate extremes to a doubling of CO₂. During winter the response both of temperature mean and daily variability relates to fractional snow cover. Cold extremes over high latitudes warm disproportionately in ensemble members with strong snow albedo feedback and large snow cover reduction. Reduced snow cover leads to more winter warming and stronger variability decrease. As a result uncertainties in mean and variability response line up, with some members showing weak and others very strong warming of the cold tail of the distribution, depending on the snow albedo parametrization. The uncertainty across the ensemble regionally exceeds the CMIP3 multi-model range. Regarding summer hot extremes, the uncertainties are larger than for mean summer warming but smaller than in multi-model experiments. The summer precipitation response to a doubling of CO₂ is not robust over many regions. Land surface parameter perturbations and natural variability alter the sign of the response even over subtropical regions.