Coupled convection and tidal dissipation in Europa’s ice shell

We performed 2D numerical simulations of oscillatory tidal flexing to study the interrelationship between tidal dissipation (calculated using the Maxwell model) and a heterogeneous temperature structure in Europa’s ice shell. Our 2D simulations show that, if the temperature is spatially uniform, the tidal dissipation rate peaks when the Maxwell time is close to the tidal period, consistent with previous studies. The tidal dissipation rate in a convective plume encased in a different background temperature depends on both the plume and background temperature. At a fixed background temperature, the dissipation increases strongly with plume temperature at low temperatures, peaks, and then decreases with temperature near the melting point when a melting-temperature viscosity of 10¹³ Pa s is used; however, the peak occurs at significantly higher temperature in this heterogeneous case than in a homogeneous medium for equivalent rheology. For constant plume temperature, the dissipation rate in a plume decreases as the surrounding temperature increases; plumes that are warmer than their surroundings can exhibit enhanced heating not only relative to their surroundings but relative to the Maxwell-model prediction for a homogeneous medium at the plume temperature. These results have important implications for thermal feedbacks in Europa’s ice shell.To self-consistently determine how convection interacts with tidal heating that is correctly calculated from the time-evolving heterogeneous temperature field, we coupled viscoelastic simulations of oscillatory tidal flexing (using Tekton) to long-term simulations of the convective evolution (using ConMan). Our simulations show that the tidal dissipation rate resulting from heterogeneous temperature can have a strong impact on thermal convection in Europa’s ice shell. Temperatures within upwelling plumes are greatly enhanced and can reach the melting temperature under plausible tidal-flexing amplitude for Europa. A pre-existing fracture zone (at least 6 km deep) promotes the concentration of tidal dissipation (up to ∼20 times more than that in the surroundings), leading to lithospheric thinning. This supports the idea that spatially variable tidal dissipation could lead locally to high temperatures, partial melting, and play an important role in the formation of ridges, chaos, or other features.     

Efficacy of a multi-metric fish index as an analysis tool for the transitional fish component of the Water Framework Directive

The WFD has introduced an international commitment to assess the ecological status of transitional waters (TWs), within which fish communities are a key biological monitoring component. The Transitional Fish Classification Index (TFCI) outlined in this paper uses 10 ecological measures to analyse fish populations caught from various ecological niches using a variety of gear types within the Thames estuary. These reach and method-specific communities are then compared to a reference population created from a 'healthy' population from TWs of a similar type. The results indicate a progressive downstream increase the quality of fish communities, consistent with previous work; variation between methods can be accounted for by gear selectivity. Overall, the TFCI is an effective communication tool for converting ecological information into an easily understood format for managers, policy makers and the general public.     

Comparison of P- and S-waves velocities estimated from Biot-Gassmann and Kuster-Toksöz models with results obtained from acoustic wavetrains interpretation

Kuster-Toksöz and Biot-Gassmann models for estimating velocities of longitudinal and shear waves on the basis of well-logging data were analysed. P-wave and S-wave velocity models are crucial for interpretation of seismic data. Discussed models enable determination with quite good accuracy, in some cases higher than the acoustic full wavetrains interpretation. Because velocity strongly depends on lithology and saturation of pore space, the selection of parameters of rock matrix, hydrocarbons and formation waters has a strong effect on the quality of velocities estimation.     

HST observations of azimuthal asymmetry in Saturn's rings

From 378 Hubble Space Telescope WFPC2 images obtained between 1996-2004, we have measured the detailed nature of azimuthal brightness variations in Saturn's rings. The extensive geometric coverage, high spatial resolution (), and photometric precision of the UBVRI images have enabled us to determine the dependence of the asymmetry amplitude and longitude of minimum brightness on orbital radius, ring elevation, wavelength, solar phase angle, and solar longitude. We explore a suite of dynamical models of self-gravity wakes for two particle size distributions: a single size and a power law distribution spanning a decade in particle radius. From these N-body simulations, we calculate the resultant wake-driven brightness asymmetry for any given illumination and viewing geometry. The models reproduce many of the observed properties of the asymmetry, including the shape and location of the brightness minimum and the trends with ring elevation and solar longitude. They also account for the “tilt effect” in the A and B rings: the change in mean ring brightness with effective ring opening angle, |Beff|. The predicted asymmetry depends sensitively on dynamical ring particle properties such as the coefficient of restitution and internal mass density, and relatively weakly on photometric parameters such as albedo and scattering phase function. The asymmetry is strongest in the A ring, reaching a maximum amplitude A∼25% near a=128,000 km. Here, the observations are well-matched by an internal particle density near 450 kg m⁻³ and a narrow particle size distribution. The B ring shows significant asymmetry (∼5%) in regions of relatively low optical depth (τ∼0.7). In the middle and outer B ring, where τ?1, the asymmetry is much weaker (∼1%), and in the C ring, A<0.5%. The asymmetry diminishes near opposition and at shorter wavelengths, where the albedo of the ring particles is lower and multiple-scattering effects are diminished. The asymmetry amplitude varies strongly with ring elevation angle, reaching a peak near |Beff|=10° in the A ring and at |Beff|=15-20° in the B ring. These trends provide an estimate of the thickness of the self-gravity wakes responsible for the asymmetry. Local radial variations in the amplitude of the asymmetry within both the A and B rings are probably caused by regional differences in the particle size distribution.     

Solar shield: forecasting and mitigating space weather effects on high-voltage power transmission systems

In this paper, central elements of the Solar Shield project, launched to design and establish an experimental system capable of forecasting the space weather effects on high-voltage power transmission system, are described. It will be shown how Sun–Earth system data and models hosted at the Community Coordinated Modeling Center (CCMC) are used to generate two-level magnetohydrodynamics-based forecasts providing 1–2 day and 30–60 min lead-times. The Electric Power Research Institute (EPRI) represents the end-user, the power transmission industry, in the project. EPRI integrates the forecast products to an online display tool providing information about space weather conditions to the member power utilities. EPRI also evaluates the economic impacts of severe storms on power transmission systems. The economic analysis will quantify the economic value of the generated forecasting system. The first version of the two-level forecasting system is currently running in real-time at CCMC. An initial analysis of the system’s capabilities has been completed, and further analysis is being carried out to optimize the performance of the system. Although the initial results are encouraging, definite conclusions about system’s performance can be given only after more extensive analysis, and implementation of an automatic evaluation process using forecasted and observed geomagnetically induced currents from different nodes of the North American power transmission system. The final output of the Solar Shield will be a recommendation for an optimal forecasting system that may be transitioned into space weather operations.     

Investigating the igneous petrogenesis of Martian volcanic rocks using augite quantitative textural analysis of the Yamato nakhlites

To better understand volcanism on planetary bodies other than the Earth, the quantification of physical processes is needed. Here, the petrogenesis of the achondrite Martian Yamato (Y) nakhlites (Y 000593, Y 000749, and Y 000802) is reinvestigated via quantitative analysis of augite (high-Ca clinopyroxene) phenocrysts: crystal size distribution (CSD), spatial distribution patterns (SDP), and electron backscatter diffraction (EBSD). Results from CSD and EBSD quantitative data sets show augite to have continuous uninterrupted growth resulting in calculated minimum magma chamber residence times of either 88–117 ± 6 yr or 9–12 yr. All samples exhibit low-intensity S-LS type crystallographic preferred orientation. Directional strain is observed across all samples with intracrystalline misorientation patterns indicative of (100)[001]:(001)[100] (Y 000593 and Y 000802) and {110}<001>or {110}¹/₂<110> (Y 000749) slip systems. SDP results indicate phenocryst-bearing crystal-clustered rock signatures. Combined findings from this work show that the Yamato nakhlites formed on Mars as individual low-viscosity lava flows or sills. This study shows that through combining these different quantitative techniques over multiple samples, one can more effectively compare and interpret resulting data to gain a more robust, geologically contextualized petrogenetic understanding of the rock suite being studied. The techniques used in this study should be equally applicable to igneous achondrites from other parent bodies.     

Realizable accessibility: evaluating the reliability of public transit accessibility using high-resolution real-time data

Journal of Geographical Systems - The widespread availability of high spatial and temporal resolution public transit data is improving the measurement and analysis of public transit-based...     

Diatoms indicate that calcium decline,not acidification,explains recent cladoceran assemblage changes in south-central Ontario softwater lakes

In recent decades, softwater lakes across Canada have experienced a wide array of anthropogenic influences, with acidification and climate warming of particular concern. Here, we compare modern and pre-industrial sedimentary diatom assemblages from 36 softwater lakes located on the Canadian Shield in south-central Ontario to determine whether lake acidification or reduced calcium availability was the main stressor responsible for recent declines in Ca-sensitive cladoceran taxa. Regional surveys of south-central Ontario water chemistry have identified the pH recovery of many formerly acidified lakes, and our fossil diatom-inferred pH analyses indicate that modern lakewater pH in the 36 study lakes is similar to (or higher than) pre-industrial levels, with diatom assemblages from both time periods dominated by taxa with similar pH preferences. In addition, modern diatom assemblages compared to pre-industrial assemblages contained higher relative abundances of planktonic diatom taxa (e.g. Asterionella formosa and the Discostella stelligera complex) and lower relative abundances of heavily silicified diatoms (e.g. Aulacoseira taxa) and benthic fragilarioid taxa. These taxonomic shifts are consistent with warming-induced changes in lake properties including a longer ice-free period, decreased wind speed and/or increased thermal stability. We conclude that recent changes observed within the cladoceran assemblages of these lakes are not a response to acidification, but are likely a consequence of Ca declines. In addition, our data suggest that regional climate warming is now responsible for the diatom changes observed in this region.     

The formation of Ganymede's grooved terrain: Numerical modeling of extensional necking instabilities

Ganymede's grooved terrain likely formed during an epoch of global expansion, when unstable extension of the lithosphere resulted in the development of periodic necking instabilities. Linear, infinitesimal-strain models of extensional necking support this model of groove formation, finding that the fastest growing modes of an instability have wavelengths and growth rates consistent with Ganymede's grooves. However, several questions remain unanswered, including how nonlinearities affect instability growth at large strains, and what role instabilities play in tectonically resurfacing preexisting terrain. To address these questions we numerically model the extension of an icy lithosphere to examine the growth of periodic necking instabilities over a broad range of strain rates and temperature gradients. We explored thermal gradients up to 45 K km⁻¹ and found that, at infinitesimal strain, maximum growth rates occur at high temperature gradients (45 K km⁻¹) and moderate strain rates (10⁻¹³ s⁻¹). Dominant wavelengths range from 1.8 to 16.4 km (post extension). Our infinitesimal growth rates are qualitatively consistent with, but an order of magnitude lower than, previous linearized calculations. When strain exceeds ∼10% growth rates decrease, limiting the total amount of amplification that can result from unstable extension. This fall-off in growth occurs at lower groove amplitudes for high-temperature-gradient, thin-lithosphere simulations than for low-temperature-gradient, thick-lithosphere simulations. At large strains, this shifts the ideal conditions for producing large amplitude grooves from high temperature gradients to more moderate temperature gradients (15 K km⁻¹). We find that the formation of periodic necking instabilities can modify preexisting terrain, replacing semi-random topography up to 100 m in amplitude with periodic ridges and troughs, assisting the tectonic resurfacing process. Despite this success, the small topographic amplification produced by our model presents a formidable challenge to the necking instability mechanism for groove formation. Success of the necking instability mechanism may require rheological weakening or strain localization by faulting, effects not included in our analysis.     

Meridiani Planum hematite deposit and the search for evidence of life on Mars—iron mineralization of microorganisms in rock varnish

The extensive hematite deposit in Meridiani Planum was selected as the landing site for the Mars Exploration Rover Opportunity because the site may have been favorable to the preservation of evidence of possible prebiotic or biotic processes. One of the proposed mechanisms for formation of this deposit involves surface weathering and coatings, exemplified on Earth by rock varnish. Microbial life, including microcolonial fungi and bacteria, is documented in rock varnish matrices from the southwestern United States and Australia. Limited evidence of this life is preserved as cells and cell molds mineralized by iron oxides and hydroxides, as well as by manganese oxides. Such mineralization of microbial cells has previously been demonstrated experimentally and documented in banded iron formations, hot spring deposits, and ferricrete soils. These types of deposits are examples of the four “water-rock interaction” scenarios proposed for formation of the hematite deposit on Mars. The instrument suite on Opportunity has the capability to distinguish among these proposed formation scenarios and, possibly, to detect traces that are suggestive of preserved martian microbiota. However, the confirmation of microfossils or preserved biosignatures will likely require the return of samples to terrestrial laboratories.