A modified landform development model for the topography of drained thermokarst lake basins in fine‐grained sediments

Permafrost degradation associated with the expansion of thermokarst lakes is commonly interrupted by catastrophic drainage. Subsequently, in tundra areas, permafrost aggradation in drained basins leads to uneven topography characterized by raised centres and wet, depressed margins. The genesis of such topography has been investigated in Old Crow Flats (OCF), a glaciolacustrine plain in the continuous permafrost of northern Yukon. The thermokarst lakes of OCF have a mean depth of only 1.5 m because excess ice is dominantly found only in the uppermost 10 m of the ground. Surface conditions were measured in three drained thermokarst lake basins, including relief, snow conditions, ground temperatures, near‐surface ground ice, and sediment stratigraphy. Four nearby lakes provided information on wave base, shore recession patterns, and bathymetry before drainage: the bottoms of these lakes were not raised in the centre. An elevation difference of up to 2 m was recorded between drained basin margins and centres but was not associated with variations in ice‐wedge density or segregated ice content. Hence basin topography was not controlled by differences in volumetric ground‐ice content between margins and centres. We propose that transport of fine sediment away from eroding lake margins during lake development is the primary mechanism for the genesis of depressed margins and raised centres in drained basins of OCF. Over time, the transport results in the deposition of more and finer sediment in the central parts of lakes, where the lake bottom has subsided below wave base, than at the shallow margins, where resuspension by wave action occurs frequently. This difference in sediment volume is revealed in the topography after drainage, when permafrost aggrades in the lake‐bottom sediment and underlying talik. Copyright © 2016 John Wiley & Sons, Ltd.     

Geochemistry and quality of the groundwater from the karstic and coastal aquifer of Geropotamos River Basin at north-central Crete, Greece

In Geropotamos River Basin, located on the north-central part of Crete, Greece, two main factors were believed to be affecting the geochemistry of the groundwater with high salt contents: seawater intrusion and/or Miocene evaporates. To identify the origin of the high salinity in groundwater, a hydrogeochemical and isotopic study was performed. Water samples from 22 wells and 2 springs were analyzed for physico-chemical parameters, major ions analysis, as well as stable isotopes (??¹⁸O, ??D). From the present survey, in which detailed hydrogeochemical investigation was conducted, the uncertainty of the contamination sources was decreased in the northern part of Geropotamos Basin. The results complement the scenario in which seawater and the widespread human activities are the principal sources of groundwater contamination. Moreover, the results of the stable isotopes analyses (??¹⁸O and ??D) support the same hypothesis and make seawater intrusion the most probable cause for the highest salinity waters. It is indicated that saline intrusion is likely to occur along fractures in a fault zone through otherwise low-permeability phyllite?Cquartzite bedrock, which demonstrates the critical role of fracture pathways in salination problems of coastal aquifers.     

Ck-Reconstruction of Surfaces from Partial Data

In this paper, we study the problem of constructing a smooth approximant of a surface defined by the equation z = f(x ₁, x ₂), the data being a finite set of patches on this surface. This problem occurs, for example, after geophysical processing such as migration of time-maps or depth-maps. The usual algorithms to solve this problem are picking points on the patches to get Lagrange's data or trying to get local junctions on patches. But the first method does not use the continuous aspect of the data and the second one does not perform well to get a global regular approximant (C ¹ or more). As an approximant of f, a discrete smoothing spline belonging to a suitable piecewise polynomial space is proposed. The originality of the method consists in the fidelity criterion used to fit the data, which takes into account their particular aspect (surface's patches): the idea is to define a function that minimizes the volume located between the data patches and the function, and which is globally C ^k. We first demonstrate the new method on a theoretical aspect and numerical results on real data are given.     

Dissolution/precipitation kinetics of boehmite and gibbsite: Application of a pH-relaxation technique to study near-equilibrium rates

The dissolution and precipitation rates of boehmite, AlOOH, at 100.3 °C and limited precipitation kinetics of gibbsite, Al(OH)₃, at 50.0 °C were measured in neutral to basic solutions at 0.1 molal ionic strength (NaCl + NaOH + NaAl(OH)₄) near-equilibrium using a pH-jump technique with a hydrogen-electrode concentration cell. This approach allowed relatively rapid reactions to be studied from under- and over-saturation by continuous in situ pH monitoring after addition of basic or acidic titrant, respectively, to a pre-equilibrated, well-stirred suspension of the solid powder. The magnitude of each perturbation was kept small to maintain near-equilibrium conditions. For the case of boehmite, multiple pH-jumps at different starting pHs from over- and under-saturated solutions gave the same observed, first order rate constant consistent with the simple or elementary reaction: .

This relaxation technique allowed us to apply a steady-state approximation to the change in aluminum concentration within the overall principle of detailed balancing and gave a resulting mean rate constant, (2.2 ± 0.3) × 10⁻⁵ kg m⁻² s⁻¹, corresponding to a 1σ uncertainty of 15%, in good agreement with those obtained from the traditional approach of considering the rate of reaction as a function of saturation index. Using the more traditional treatment, all dissolution and precipitation data for boehmite at 100.3 °C were found to follow closely the simple rate expression:

R_net,boehmite=10^-5.485{m_OH^-}{1-exp(ΔG_r/RT)}, with R_net in units of mol m⁻² s⁻¹. This is consistent with Transition State Theory for a reversible elementary reaction that is first order in OH⁻ concentration involving a single critical activated complex. The relationship applies over the experimental ΔG_r range of 0.4–5.5 kJ mol⁻¹ for precipitation and −0.1 to −1.9 kJ mol⁻¹ for dissolution, and the pH_m ≡ −log(mH⁺) range of 6–9.6. The gibbsite precipitation data at 50 °C could also be treated adequately with the same model:R_net,gibbsite=10^-5.86{m_OH^-}{1-exp(ΔG_r/RT)}, over a more limited experimental range of ΔG_r (0.7–3.7 kJ mol⁻¹) and pH_m (8.2–9.7).     

The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection

The LMDZ4 general circulation model is the atmospheric component of the IPSL–CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley–Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke’s convection scheme, used in previous versions, the Emanuel’s scheme improves the representation of the Hadley–Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke’s parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.     

Orbital imprint on Holocene palaeohydrological variations in west‐central Europe as reflected by lake‐level changes at Cerin (Jura Mountains,eastern France)

This paper presents a lake‐level record for the Holocene at Lake Cerin (Jura Mountains, eastern France). It is based on a range of sedimentological techniques validated in previous studies, with a combination of systematic lithostratigraphic investigations of the infillings accumulated in the lacustrine basin, and sediment analyses of two selected cores. The chronology is based on 10 radiocarbon dates and pollen stratigraphy. On a millennial scale, the Cerin lake‐level record shows three distinct successive phases characterised by higher lake‐level conditions until ca. 9000 cal. a BP, followed by a maximal lowering at ca. 9000–8500 cal. a BP, and a progressive rise until the present. This rise was punctuated by centennial‐scale fluctuations, with major events around 4000, 2800 and after 1500 cal. a BP. Considered on a multimillennial scale, the general pattern of palaeohydrological changes reconstructed at Cerin reflects the impact of orbitally driven summer insolation. This is in agreement with other regional and extra‐regional palaeoclimatic records, although every record shows peculiarities in timing and shape depending on the proxy used for reconstruction. In this general context, centennial to multicentennial oscillations appear to have been second‐order events in comparison with the major influence of the orbital factor. Copyright © 2011 John Wiley & Sons, Ltd.     

An assessment of large-scale flood modelling based on LiDAR data

Large-scale flood modelling approaches designed for regional to continental scales usually rely on relatively simple assumptions to represent the potentially highly complex river bathymetry at the watershed scale based on digital elevation models (DEMs) with a resolution in the range of 25–30 m. Here, high-resolution (1 m) LiDAR DEMs are employed to present a novel large-scale methodology using a more realistic estimation of bathymetry based on hydrogeomorphological GIS tools to extract water surface slope. The large-scale 1D/2D flood model LISFLOOD-FP is applied to validate the simulated flood levels using detailed water level data in four different watersheds in Quebec (Canada), including continuous profiles over extensive distances measured with the HydroBall technology. A GIS-automated procedure allows to obtain the average width required to run LISFLOOD-FP. The GIS-automated procedure to estimate bathymetry from LiDAR water surface data uses a hydraulic inverse problem based on discharge at the time of acquisition of LiDAR data. A tiling approach, allowing several small independent hydraulic simulations to cover an entire watershed, greatly improves processing time to simulate large watersheds with a 10-m resampled LiDAR DEM. Results show significant improvements to large-scale flood modelling at the watershed scale with standard deviation in the range of 0.30 m and an average fit of around 90%. The main advantage of the proposed approach is to avoid the need to collect expensive bathymetry data to efficiently and accurately simulate flood levels over extensive areas.     

Data analysis pipeline for EChO end-to-end simulations

Atmospheric spectroscopy of extrasolar planets is an intricate business. Atmospheric signatures typically require a photometric precision of 1×10^?4 in flux over several hours. Such precision demands high instrument stability as well as an understanding of stellar variability and an optimal data reduction and removal of systematic noise. In the context of the EChO mission concept, we here discuss the data reduction and analysis pipeline developed for the EChO end-to-end simulator EChOSim. We present and discuss the step by step procedures required in order to obtain the final exoplanetary spectrum from the EChOSim ‘raw data’ using a simulated observation of the secondary eclipse of the hot-Neptune 55 Cnc e.     

Reply to “Reply to comments on defining biominerals and organominerals: Direct and indirect indicators of life [Perry et al., Sedimentary Geology, 201, 157–179]” by R.S. Perry and M.A. Sephton: [Sedimentary Geology 213 (2009) 156]

This is a reply to R.S. Perry and M.A. Sephton's “Reply to comments on defining biominerals and organominerals: direct and indirect indicators of life [Perry et al., Sedimentary Geology, 201, 157–179]” [Sedimentary Geology 213 (2009) 156].