Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.

Silicate mineralogy of basalts at DSDP Leg 54 Holes

Basalts drilled from the East Pacific Rise, OCP Ridge, and Siqueiros fracture zone during Leg 54 are texturally diverse. Dolerites are equigranular at Sites 422 and 428 and porphyritic, with phenocrysts of plagioclase (An69.73) and Ca-rich clinopyroxene (Ca42Mg48Fe10) at Site 427. The East Pacific Rise lavas and some of those from the OCP Ridge are fine-grained and porphyritic. The majority of the large crystals are clustered skeletal glomerocrysts of plagioclase An64-77), together with olivine (Fo80-87), Ca-rich clinopyroxene, or both. Euhedral phenocrysts of plagioclase, together with olivine, Carich clinopyroxene, and Cr-Al spinel in some cases, occur in most of the fine-grained lavas. These phenocrysts are small (maximum dimension <1 mm in all but one sample), sparse (combined modal amount <1% in all samples), and distinctive from the megacrysts which characterize many ocean-floor lavas. In two East Pacific Rise lavas, zoned plagioclase (An83 cores) is the sole phenocryst phase. In other porphyritic lavas from all the main East Pacific Rise and OCP Ridge units drilled during Leg 54, the plagioclase phenocrysts contain cores of bytownite (An79-87) surrounded by more-sodic feldspar (An67-77). Core/rim relationships vary from continuous normal zoning, through discontinuous zoning, to extensive resorption of the calcic cores in some samples. The compositions of the plagioclase calcic cores are systematically related to those of the glomerophyric plagioclase and olivine in the lavas containing them. Furthermore, only one compositional population of calcic cores occurs in each rock. The possible causes of these relationships are far from clear. Magma mixing, although superficially applicable, is inconsistent with important aspects of the phenocryst mineralogy of these particular lavas. A more satisfactory model to explain both phenocryst zoning and rapid glomerocryst growth immediately before extrusion may be constructed by postulating influx of water into the upwelling magmas within Layer 3 of the oceanic crust beneath the East Pacific Rise, and subsequent loss of part of this water during effervescence within feeder dykes between Layer 3 and the ocean floor. It is shown that this model is fully consistent with published data on water and carbon dioxide contents and ratios in the pillow-margin glasses, vesicles, and phenocryst inclusions of ocean-floor basalts. The evidence for the precipitation of plagioclase- dominated crystalline assemblages from these magmas in the upper part of Layer 3 is concordant with recent geophysically based modeling of the structure of the East Pacific Rise. Calcium-rich clinopyroxenes in dolerites from the OCP Ridge and Siqueiros fracture zone show radial, oscillatory, and sector-zoning. In Sample 428A-5-2 (Piece 5a), the compositional trends resulting from this zoning closely resemble those of the pyroxenes in some lunar lavas. The controls on crystallization of interstitial pigeonite - epitaxial upon augite - in this rock are discussed. Both sector-zoning of the augite and nucleation of pigeonite within microvolumes of magma with a low Ca(Mg + Fe) ratio appear to be important factors.

Downscaled land surface temperatures obtained by a particle smoothing approach, with links to NetCDF files

This dataset contains continuous time series of land surface temperature (LST) at spatial resolution of 300m around the 12 experimental sites of the PAGE21 project (grant agreement number 282700, funded by the EC seventh Framework Program theme FP7-ENV-2011). This dataset was produced from hourly LST time series at 25km scale, retrieved from SSM/I data (André et al., 2015, doi:10.1016/j.rse.2015.01.028) and downscaled to 300m using a dynamic model and a particle smoothing approach. This methodology is based on two main assumptions. First, LST spatial variability is mostly explained by land cover and soil hydric state. Second, LST is unique for a land cover class within the low resolution pixel. Given these hypotheses, this variable can be estimated using a land cover map and a physically based land surface model constrained with observations using a data assimilation process. This methodology described in Mechri et al. (2014, doi:10.1002/2013JD020354) was applied to the ORCHIDEE land surface model (Krinner et al., 2005, doi:10.1029/2003GB002199) to estimate prior values of each land cover class provided by the ESA CCI-Land Cover product (Bontemps et al., 2013) at 300m resolution . The assimilation process (particle smoother) consists in simulating ensemble of LST time series for each land cover class and for a large number of parameter sets. For each parameter set, the resulting temperatures are aggregated considering the grid fraction of each land cover and compared to the coarse observations. Miniminizing the distance between the aggregated model solutions and the observations allow us to select the simulated LST and the corresponding parameter sets which fit the observations most closely. The retained parameter sets are then duplicated and randomly perturbed before simulating the next time window. At the end, the most likely LST of each land cover class are estimated and used to reconstruct LST maps at 300m resolution using ESA CCI-Land Cover. The resulting temperature maps on which ice pixels were masked, are provided at daily time step during the nine-year analysis period (2000-2009).