Turkey Lakes Watershed,Ontario, Canada: 40 years of interdisciplinary whole-ecosystem research

The Turkey Lakes Watershed (TLW) study is a federal, interdepartmental study established in 1979 to investigate the effects of acid rain on terrestrial and aquatic ecosystems. The 10.5 km² watershed, located in the Eastern Temperate Mixed Forest on the Canadian Shield, has been the site of multidisciplinary studies on biogeochemical and ecological processes conducted across plot to catchment scales. The whole-ecosystem investigative approach was adopted from the outset and has allowed research to evolve from its original (and continuing) acidification focus to include investigations on the effects of climate change, forest harvesting and other forest ecosystem perturbations. The extensive scientific and support infrastructure allows for collection of a comprehensive data record essential for understanding long-term environmental trends. Data include atmospheric deposition, meteorology, stream hydrology and chemistry, soil, pore and ground water properties, understory and overstory vegetation, lake and outflow physical and chemical properties, and aquatic macroinvertebrate and fish community composition and abundance. These data have contributed to over 400 published research papers and graduate theses. The watershed has also figured prominently in many continent-wide comparisons advancing fundamental watershed theory. The knowledge gained at TLW has influenced pollutant emission and natural resource management policies provincially, nationally and internationally.     

New satellite derived sea ice motion tracks Arctic contamination

Sea ice has been reported to contain contaminants from atmospheric and nearshore sediment resuspension processes. In this study successive passive microwave images from the 85.5 GHz channels on the Special Sensor Microwave Imager (SSM/I) were merged with drifting buoy trajectories from the International Arctic Buoy Program to compute Arctic sea ice motion in the Russian Arctic between 1988 and 1994. Smooth daily motion fields were averaged to prepare monthly maps making it possible to compute the 7-year mean and mean seasonal ice motions as well as principal components of directional variability of sea ice motion for the entire Arctic and surrounding basins. These mean motion vectors are used to simulate the advection of contaminants deposited on or contained within the sea ice and subsequently transported into the Arctic Ocean in order to predict both their mean trajectories and dispersal over time. The 3-year displacement of contaminants from a number of Russian sites and one American site display various behaviours from substantial displacement and dispersal to almost no movement. This computational procedure could be applied to realtime SSM/I and ice buoy data to provide detailed, all-weather, vector motion maps of ice circulation to predict the path and dispersal of any new substance introduced to the sea ice and transported into the Arctic or Antarctic ocean surface.     

Submesoscale features and their interaction with fronts and internal tides in a high-resolution coupled atmosphere-ocean-wave model of the Bay of Bengal

Large freshwater fluxes into the Bay of Bengal by rainfall and river discharges result in strong salinity fronts in the bay. In this study, a high-resolution coupled atmosphere-ocean-wave model with comprehensive physics is used to model the weather, ocean circulation, and wave field in the Bay of Bengal. Our objective is to explore the submesoscale activity that occurs in a realistic coupled model that resolves mesoscales and allows part of the submesoscale field. Horizontal resolution in the atmosphere varies from 2 to 6 km and is 13 km for surface waves, while the ocean model is submesoscale permitting with resolutions as high as 1.5 km and a vertical resolution of 0.5 m in the upper 10 m. In this paper, three different cases of oceanic submesoscale features are discussed. In the first case, heavy rainfall and intense downdrafts produced by atmospheric convection are found to force submesoscale currents, temperature, and salinity anomalies in the oceanic mixed layer and impact the mesoscale flow. In a second case, strong solitary-like waves are generated by semidiurnal tides in the Andaman Sea and interact with mesoscale flows and fronts and affect submesoscale features generated along fronts. A third source of submesoscale variability is found further north in the Bay of Bengal where river outflows help maintain strong salinity gradients throughout the year. For that case, a comparison with satellite observations of sea surface height anomalies, sea surface temperature, and chlorophyll shows that the model captures the observed mesoscale eddy features of the flow field, but in addition, submesoscale upwelling and downwelling patterns associated with ageostrophic secondary circulations along density fronts are also captured by the model.     

Advancing tsunami risk assessment by improving spatio-temporal population exposure and evacuation modeling

Tsunamis are among the most destructive and lethal of coastal hazards. These are time-specific events, and despite directly affecting a narrow strip of coastline, a single occurrence can have devastating effects and cause massive loss of life, especially in urbanized coastal areas. In this work, in order to consider the time dependence of population exposure to tsunami threat, the variation of spatio-temporal population distribution in the daily cycle is mapped and analyzed in the Lisbon Metropolitan Area. High-resolution daytime and nighttime population distribution maps are developed using ‘intelligent dasymetric mapping,’ that is, applying areal interpolation to combine best-available census data and statistics with land use and land cover data. Workplace information and mobility statistics are considered for mapping daytime distribution. In combination with a tsunami hazard map, information on infrastructure, land use and terrain slope, the modeled population distribution is used to assess people’s evacuation speed, applying a geospatial evacuation modeling approach to the city of Lisbon. The detailed dynamic population exposure assessment allows producing both daytime and nighttime evacuation time maps, which provide valuable input for evacuation planning and management. Results show that a significant amount of population is at risk, and its numbers increase dramatically from nighttime to daytime, especially in the zones of high tsunami flooding susceptibility. Also, full evacuation can be problematic in the daytime period, even if initiated immediately after a major tsunami-triggering earthquake. The presented approach greatly improves tsunami risk assessment and can benefit all phases of the disaster management process.     

Deep versus shallow origin of gravity anomalies, topography and volcanism on Earth, Venus and Mars

The relation between gravity anomalies, topography and volcanism can yield important insights about the internal dynamics of planets. From the power spectra of gravity and topography on Earth, Venus and Mars we infer that gravity anomalies have likely predominantly sources below the lithosphere up to about spherical harmonic degree l=30 for Earth, 40 for Venus and 5 for Mars. To interpret the low-degree part of the gravity spectrum in terms of possible sublithospheric density anomalies we derive radial mantle viscosity profiles consistent with mineral physics. For these viscosity profiles we then compute gravity and topography kernels, which indicate how much gravity anomaly and how much topography is caused by a density anomaly at a given depth. With these kernels, we firstly compute an expected gravity-topography ratio. Good agreement with the observed ratio indicates that for Venus, in contrast to Earth and Mars, long-wavelength topography is largely dynamically supported from the sublithospheric mantle. Secondly, we combine an empirical power spectrum of density anomalies inferred from seismic tomography in Earth’s mantle with gravity kernels to model the gravity power spectrum. We find a good match between modeled and observed gravity power spectrum for all three planets, except for 2?l?4 on Venus. Density anomalies in the Venusian mantle for these low degrees thus appear to be very small. We combine gravity kernels and the gravity field to derive radially averaged density anomaly models for the Martian and Venusian mantles. Gravity kernels for l?5 are very small on Venus below ≈800 km depth. Thus our inferences on Venusian mantle density are basically restricted to the upper 800 km. On Mars, gravity anomalies for 2?l?5 may originate from density anomalies anywhere within its mantle. For Mars as for Earth, inferred density anomalies are dominated by l=2 structure, but we cannot infer whether there are features in the lowermost mantle of Mars that correspond to Earth’s Large Low Shear Velocity Provinces (LLSVPs). We find that volcanism on Mars tends to occur primarily in regions above inferred low mantle density, but our model cannot distinguish whether or not there is a Martian analog for the finding that Earth’s Large Igneous Provinces mainly originate above the margins of LLSVPs.     

Detection of Earth-mass and super-Earth Trojan planets using transit timing variation method

We have carried out an extensive study of the possibility of the detection of Earth-mass and super-Earth Trojan planets using transit timing variation method with the Kepler space telescope. We have considered a system consisting of a transiting Jovian-type planet in a short period orbit, and determined the induced variations in its transit timing due to an Earth-mass/super-Earth Trojan planet. We mapped a large section of the phase space around the 1:1 mean-motion resonance and identified regions corresponding to several other mean-motion resonances where the orbit of the planet would be stable. We calculated transit timing variations (TTVs) for different values of the mass and orbital elements of the transiting and perturbing bodies as well as the mass of central star, and identified orbital configurations of these objects (ranges of their orbital elements and masses) for which the resulted TTVs would be within the range of the variations of the transit timing of Kepler’s planetary candidates. Results of our study indicate that in general, the amplitudes of the TTVs fall within the detectable range of timing precision obtained from the Kepler’s long-cadence data, and depending on the parameters of the system, their magnitudes may become as large as a few hours. The probability of detection is higher for super-Earth Trojans with slightly eccentric orbits around short-period Jovian-type planets with masses slightly smaller than Jupiter. We present the details of our study and discuss the implications of its results.     

Geochemistry of fine-grained sediments of the upper Cretaceous to Paleogene Gosau Group （Austria, Slovakia） Implications for paleoenvironmental and provenance studies

Bulk rock geochemistry of 169 fine-grained sediment samples of the upper Cretaceous to Paleogene Gosau Group（Northern Calcareous Alps,Austria and Slovakia） from borehole and outcrop localities was performed to separate non-marine and marine deposits.Geochemical characteristics of different Gosau depositional systems,basins and sediment provenance using major-,trace-,and rare earth elements were also investigated.Geochemical proxies such as boron concentrations were tested for seeking the possibilities of paleosalinity indicators.Due to the fact that several pelagic sections are represented by extremely low boron contents.B/Al^* ratios are recognized as more robust and differentiate reliably between marine（mean:160±34） and non-marine（mean:133±33） samples.Using statistical factor analysis,hemipelagic to pelagic samples from the Gieβhbl Syncline and Slovakian equivalents can be differentiated from marginal-marine to non-marine samples from the Grnbach and Glinzendorf Syncline related to terrigenous（SiCh.Al₂O₃,K₂O,Th,Rb,Zr and others） and pelagic indicative elements （CaO,Sr,TOT/C and B/Al^*）.A clear indication for ophiolitic provenance is traced by high amounts of chromium and nickel.Only non-marine successions of the Glinzendorf Syncline show higher Cr and Ni concentrations（up to 250 and 400 ppm,respectively） and enriched Cr/V and Y/Ni ratios trending to an ultramafic source.     

Comparison of rates of ureolysis between Sporosarcina pasteurii and an indigenous groundwater community under conditions required to precipitate large volumes of calcite

Ureolysis-driven calcite precipitation has potential to seal porosity and fracture networks in rocks thus preventing groundwater flow and contaminant transport. In this study urea hydrolysis and calcite precipitation rates for the model bacterium Sporosarcina pasteurii were compared with those of indigenous groundwater communities under conditions required to precipitate large volumes of calcite (up to 50 g L⁻¹). We conducted microcosm experiments in oxic artificial and anoxic natural groundwaters (collected from the Permo-Triassic sandstone aquifer at Birmingham, UK) that were inoculated with aerobically grown S. pasteurii. The rate constants for urea hydrolysis, k_urea, ranged between 0.06 and 3.29 d⁻¹ and were only affected by inoculum density. Higher Ca²⁺ concentration (50-500 mM Ca²⁺) as well as differences in fO₂ did not inhibit the ureolytic activity of S. pasteurii and did not significantly impact k_urea. These results demonstrate that S. pasteurii has potential to improve calcite precipitation in both oxic and anoxic groundwaters, especially if indigenous communities lack ureolytic activity. Urea hydrolysis by indigenous groundwater communities was investigated in anoxic, natural groundwaters amended with urea and CaCl₂. A notable increase in ureolysis rates was measured only when these communities were stimulated with dilute nutrients (with best results from blackstrap molasses). Furthermore, there was a considerable lag time (12-20 days) before ureolysis and calcite precipitation began. Calculated ureolysis rate constants, k_urea, ranged between 0.03 and 0.05 d⁻¹ and were similar to k_urea values produced by S. pasteurii at low inoculum densities. Overall, this comparative study revealed that the growth of ureolytic microorganisms present within groundwaters can easily be stimulated to enhance rates of urea hydrolysis in the subsurface, and thus can be used to induce calcite precipitation in these environments. The time required for urea hydrolysis to begin is almost instantaneous if an inoculum of S. pasteurii is included, while it may take several weeks for ureolytic groundwater communities to grow and become ureolytically active.     

Gas migration through Opouawe Bank at the Hikurangi margin offshore New Zealand

This study presents 2D seismic reflection data, seismic velocity analysis, as well as geochemical and isotopic porewater compositions from Opouawe Bank on New Zealand’s Hikurangi subduction margin, providing evidence for essentially pure methane gas seepage. The combination of geochemical information and seismic reflection images is an effective way to investigate the nature of gas migration beneath the seafloor, and to distinguish between water advection and gas ascent. The maximum source depth of the methane that migrates to the seep sites on Opouawe Bank is 1,500–2,100 m below seafloor, generated by low-temperature degradation of organic matter via microbial CO₂ reduction. Seismic velocity analysis enabled identifying a zone of gas accumulation underneath the base of gas hydrate stability (BGHS) below the bank. Besides structurally controlled gas migration along conduits, gas migration also takes place along dipping strata across the BGHS. Gas migration on Opouawe Bank is influenced by anticlinal focusing and by several focusing levels within the gas hydrate stability zone.