Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Curve matching approaches to waveform classification: a case study using ICESat data

Light Detection and Ranging (LiDAR) waveforms are being increasingly used in many forest and urban applications, especially for ground feature classification. However, most studies relied on either discretizing waveforms to multiple returns or extracting shape metrics from waveforms. The direct use of the full waveform, which contains the most comprehensive and accurate information has been scarcely explored. We proposed to utilize the complete waveform to test its ability to differentiate between objects having distinct vertical structures using curve matching approaches. Two groups of curve matching approaches were developed by extending methods originally designed for pixel-based hyperspectral image classification and object-based high spatial image classification. The first group is based on measuring the curve similarity between an unknown waveform and a reference waveform, including curve root sum squared differential area (CRSSDA), curve angle mapper (CAM), and Kullback–Leibler (KL) divergence. The second group assesses the curve similarity between an unknown and reference cumulative distribution functions (CDFs) of their waveforms, including cumulative curve root sum squared differential area (CCRSSDA), cumulative curve angle mapper (CCAM), and Kolmogorov–Smirnov (KS) distance. When employed to classify open space, trees, and buildings using ICESat waveform data, KL provided the highest average classification accuracy (87%), closely followed by CCRSSDA and CCAM, and they all significantly outperformed KS, CRSSDA, and CAM based on 15 randomized sample sets.     

Seasonal variation of carbon fluxes in a sparse savanna in semi arid Sudan

Background

Large spatial, seasonal and annual variability of major drivers of the carbon cycle (precipitation, temperature, fire regime and nutrient availability) are common in the Sahel region. This causes large variability in net ecosystem exchange and in vegetation productivity, the subsistence basis for a major part of the rural population in Sahel. This study compares the 2005 dry and wet season fluxes of CO₂ for a grass land/sparse savanna site in semi arid Sudan and relates these fluxes to water availability and incoming photosynthetic photon flux density (PPFD). Data from this site could complement the current sparse observation network in Africa, a continent where climatic change could significantly impact the future and which constitute a weak link in our understanding of the global carbon cycle.

Results

The dry season (represented by Julian day 35–46, February 2005) was characterized by low soil moisture availability, low evapotranspiration and a high vapor pressure deficit. The mean daily NEE (net ecosystem exchange, Eq. 1) was -14.7 mmol d^-1 for the 12 day period (negative numbers denote sinks, i.e. flux from the atmosphere to the biosphere). The water use efficiency (WUE) was 1.6 mmol CO₂ mol H₂O^-1 and the light use efficiency (LUE) was 0.95 mmol CO₂ mol PPFD^-1. Photosynthesis is a weak, but linear function of PPFD. The wet season (represented by Julian day 266–273, September 2005) was, compared to the dry season, characterized by slightly higher soil moisture availability, higher evapotranspiration and a slightly lower vapor pressure deficit. The mean daily NEE was -152 mmol d^-1 for the 8 day period. The WUE was lower, 0.97 mmol CO₂ mol H₂O^-1 and the LUE was higher, 7.2 μmol CO₂ mmol PPFD^-1 during the wet season compared to the dry season. During the wet season photosynthesis increases with PPFD to about 1600 μmol m^-2s^-1 and then levels off.

Conclusion

Based on data collected during two short periods, the studied ecosystem was a sink of carbon both during the dry and wet season 2005. The small sink during the dry season is surprising and similar dry season sinks have not to our knowledge been reported from other similar savanna ecosystems and could have potential management implications for agroforestry. A strong response of NEE versus small changes in plant available soil water content was found. Collection and analysis of flux data for several consecutive years including variations in precipitation, available soil moisture and labile soil carbon are needed for understanding the year to year variation of the carbon budget of this grass land/sparse savanna site in semi arid Sudan.     

Petroleum system analysis of the Khatatba Formation in the Shoushan Basin,north Western Desert,Egypt

The Middle Jurassic Khatatba Formation is an attractive petroleum exploration target in the Shoushan Basin, north Western Desert, Egypt. However, the Khatatba petroleum system with its essential elements and processes has not been assigned yet. This study throws the lights on the complete Khatatba petroleum system in the Shoushan Basin which has been evaluated and collectively named the Khatatba-Khatatba (!) petroleum system. To evaluate the remaining hydrocarbon potential of the Khatatba system, its essential elements were studied, in order to determine the timing of hydrocarbon generation, migration and accumulation. Systematic analysis of the petroleum system of the Khatatba Formation has identified that coaly shales and organic-rich shales are the most important source rocks. These sediments are characterised by high total organic matter content and have good to excellent hydrocarbon generative potential. Kerogen is predominantly types II–III with type III kerogen. The Khatatba source rocks are mature and, at the present time, are within the peak of the oil window with vitrinite reflectance values in the range of 0.81 to 1.08 % Ro. The remaining hydrocarbon potential is anticipated to exist mainly in stratigraphic traps in the Khatatba sandstones which are characterised by fine to coarse grain size, moderate to well sorted. It has good quality reservoir with relatively high porosity and permeability values ranging from 1 to 17 % and 0.05–1,000 mD, respectively. Modelling results indicated that hydrocarbon generation from the Khatatba source rocks began in the Late Cretaceous time and peak of hydrocarbon generation occurred during the end Tertiary time (Neogene). Hydrocarbon primarily migrated from the source rock via fractured pathways created by abnormally high pore pressures resulting from hydrocarbon generation. Hydrocarbon secondarily migrated from active Khatatba source rocks to traps side via vertical migration pathways through faults resulting from Tertiary tectonics during period from end Oligocene to Middle Miocene times.     

Spatial distribution and contamination assessment of heavy metals in surface soils of Hassi Messaoud, Algeria

An increase in heavy metal pollution in the soils of Hassi Messaoud (Algeria) due to intense industrialization and urbanization has become a serious environmental problem. There are three large industrial complexes that have been established in the region of Hassi Messaoud for petroleum extraction field and refinery. The region hosts several industrial facilities which are the main sources for hazardous wastes. Surface soil samples from 58 sampling sites (systematically sampled; 1 × 1 km regular grid), including different functional areas in Hassi Messaoud, were collected and analyzed. The results showed that the average concentrations of Cu, Ni, Mn, Pb and Zn in soil of Hassi Messaoud were up to 13.17, 35.78, 121.21, 130.97 and 61.08 mg/kg, respectively. Ni concentrations were comparable to background values, while Cu, Mn, Pb and Zn concentrations were higher than their corresponding background values. Among the functional areas, the industrial regions displayed the highest metal concentrations, while the lowest concentrations occurred in rural soil. Principal component analysis coupled with cluster analysis showed that: (1) Pb and Zn had anthropogenic sources; and (2) Ni, Cu and Mn were associated with parent materials. Contaminations in soils were classified as geoaccumulation index and enrichment factor. Pollution index values of Cu, Ni, Mn, Pb and Zn varied in the range of 0.04–5.41, 0.46–2.49, 0.01–5.73, 0.62–152.9 and 0.09–53.01, with mean values of 1.32, 1.08, 1.26, 5.64 and 3.1, respectively. The integrated pollution index (IPI) of all the analyzed samples varied from 0.42 to 31.59, with a mean of 2.48, and more than 5.45 % of samples are extremely contaminated; 18.18 % are heavily contaminated; 60 % are moderately contaminated; and others are low contaminated. The spatial distribution of IPI showed that desert and rural areas displayed relatively lower heavy metal contamination in comparison with other areas.     

Compressive and Tensile Behavior of Polymer Treated Sulfate Contaminated CL Soil

In this study, the compressive and tensile behavior of polymer treated sulfate contaminated CL soil was investigated. Based on the information in the literature, a field soil was contaminated with up to 4 % (40,000 ppm) of calcium sulfate in this study. In addition to characterizing the behavior of sulfate contaminated CL soil, the effect of treating the soil with a polymer solution was investigated and the performance was compared to 6 % lime treated soil. In treating the soil, acrylamide polymer solution (15 g of polymer dissolved in 85 g of water) content was varied up to 15 % (by dry soil weight). Addition of 4 % calcium sulfate to the soil decreased the compressive and tensile strengths of the compacted soils by 22 and 33 % respectively with the formation of calcium silicate sulfate [ternesite Ca₅(SiO₄)₂SO₄)], magnesium silicate sulfate (Mg₅(SiO₄)₂SO₄) and calcium-magnesium silicate (merwinite Ca₃Mg(SiO₄)₂). With the polymer treatment the strength properties of sulfate contaminated CL soil was substantially improved. Polymer treated sulfate soils had higher compressive and tensile strengths and enhanced compressive stress–strain relationships compared to the lime treated soils. Also polymer treated soils gained strength more rapidly than lime treated soil. With 10 % of polymer solution treatment, the maximum unconfined compressive and splitting tensile strengths for 4 % of calcium sulfate soil were 625 kPa (91 psi) and 131 kPa (19 psi) respectively in 1 day of curing. Similar improvement in the compressive modulus was observed with polymer treated sulfate contaminated CL soil. The variation of the compacted compressive strength and tensile strength with calcium sulfate concentrations for the treated soils were quantified and the parameters were related to calcium sulfate content in the soil and polymer content. Compressive stress–strain relationships of the sulfate soil, with and without lime and polymer treatment, have been quantified using two nonlinear constitutive models. The constitutive model parameters were sensitive to the calcium sulfate content and the type of treatment.     

Analysis of physical factors of climate and bioclimate and their effects on almonds production to increase the economy in Hebron area of Palestine

For the 1993–2009 period, we analyzed the relationship between almond yield and three climatic variables (mean annual temperature, soil water reserve, and precipitation), and four bioclimatic variables (annual ombrothermic index, water deficit, simple continentality index, and compensated thermicity index), for one major Hebron crop (soft and hard almonds). Moreover, we obtained data almond production from the Palestinian Central Bureau of Statistics, while the climate data from the Palestinian meteorological station during the study period from 1993 to 2009, and analysis is it by using bioclimatic classification of the Earth of Salvador Rivas-Martinez to study the relationship between the almond yield and climate and bioclimate factors (variables). The climatic and bioclimate variables of greatest importance to almond were used to develop regressions analysis relating yield to climatic conditions. Hebron was positively affected by annual ombrothermic index, simple continentality index, precipitation, water soil reserve, and mean annual temperature, but negatively affected by water deficit, with a large proportion of the variance explained by axis F1 (72.48%), F2 (22.38%), and axes F1and F2 (94. 86%). However, in order to produce a high amount of almonds and quality, it can be grown in the regions of the mesomediterranean region, with the value of annual ombrothemic index more than 3, compensated thermicity index between 220/220 to 350/350, simple continentality index between 14 and 20, and in areas where the average annual temperature is between 15 and 20 °C.