Refined modeling and movement characteristics analyses of irregularly shaped particles

Irregularly shaped (IRS) particles widely exist in many engineering and industrial fields. The macro physical and mechanical properties of the particle system are governed by the interaction between the particles in the system. The interaction between IRS particles is more complicated because of their complex geometric shape with extremely irregular and co‐existed concave and convex surfaces. These particles may interlock each other, making the sliding and friction of IRS particles more complex than that of particles with regular shape. In order to study the interaction of IRS particles more efficiently, a refined method of constructing discrete element model based on computed tomography scanning of IRS particles is proposed. Three parameters were introduced to control the accuracy and the number of packing spheres. Subsequently, the inertia tensor of the IRS particle model was optimized. Finally, laboratory and numerical open bottom cylinder tests were carried out to verify the refined modeling method. The influence of particle shape, particle position, and mesoscopic friction coefficient on the interaction of particles was also simulated. It is noteworthy that with the increase of mesoscopic friction coefficient, the fluidity of IRS particle assembly decreases, and intermittent limit equilibrium state may appear. Copyright © 2014 John Wiley & Sons, Ltd.     

Mercury and methylmercury biogeochemistry in a thawing permafrost wetland complex,Northwest Territories,Canada

In arctic and sub‐arctic environments, mercury (Hg), more specifically toxic methylmercury (MeHg), is of growing concern to local communities because of its accumulation in fish. In these regions, there is particular interest in the potential mobilization of atmospherically deposited Hg sequestered in permafrost that is thawing at unprecedented rates. Permafrost thaw and the resulting ground surface subsidence transforms forested peat plateaus into treeless and permafrost‐free thermokarst wetlands where inorganic Hg released from the thawed permafrost and draining from the surrounding peat plateaus may be transformed to MeHg. This study begins to characterize the spatial distribution of MeHg in a peat plateau–thermokarst wetland complex, a feature that prevails throughout the wetland‐dominated southern margin of thawing discontinuous permafrost in Canada's Northwest Territories. We measured pore water total Hg, MeHg, dissolved organic matter characteristics and general water chemistry parameters to evaluate the role of permafrost thaw on the pattern of water chemistry. A gradient in vegetation composition, water chemistry and dissolved organic matter characteristics followed a toposequence from the ombrotrophic bogs near the crest of the complex to poor fens at its downslope margins. We found that pore waters in poor fens contained elevated levels of MeHg, and the water draining from these features had dissolved MeHg concentrations 4.5 to 14.5 times higher than the water draining from the bogs. It was determined through analysis of historical aerial images that the poor fens in the toposequence had formed relatively recently (early 1970s) as a result of permafrost thaw. Differences between the fens and bogs are likely to be a result of their differences in groundwater function, and this suggests that permafrost thaw in this landscape can result in hotspots for Hg methylation that are hydrologically connected to downstream ecosystems. Copyright © 2016 John Wiley & Sons, Ltd.     

A multi‐proxy record of changing environments from ca. 30 000 to 9000 cal. a BP: Onepoto maar palaeolake,Auckland, New Zealand

We present a high‐resolution record of lacustrine sedimentation spanning ca. 30 000 to 9000 cal. a BP from Onepoto maar, northern North Island, New Zealand. The multi‐proxy record of environmental change is constrained by tephrochronology and accelerator mass spectrometric ¹⁴C ages and provides evidence for episodes of rapid environmental change during the Last Glacial Coldest Period (LGCP) and Last Glacial–Interglacial Transition (LGIT) from northern New Zealand. The multi‐proxy palaeoenvironmental record from Onepoto indicates that the LGCP was cold, dry and windy in the Auckland region, with vegetation dominated by herb and grass in a beech forest mosaic between ca. 28 500 and 18 000 cal. a BP. The LGCP was accompanied by more frequent fires and influx of clastic sediment indicating increased erosion during the LGCP, with a mid‐LGCP interstadial identified between ca. 25 000 and 23 000 cal. a BP. Rapid climate amelioration at ca. 18 000 cal. a BP was accompanied by increased terrestrial biomass exemplified by the expansion of lowland podocarp forest, especially Dacrydium cupressinum. Increasing biomass production is reversed briefly by LGIT perturbations which are apparent in many of the proxies that span ca. 14 000–10 500 cal. a BP, suggesting generally increased wetness and higher in situ aquatic plant productivity with reduced terrestrial organic matter and terrigenous detrital influx. Furthermore, conditions at that time were probably warmer and frosts rare based on the increasing importance of Ascarina. The subsequent early Holocene is characterised by podocarp conifer forest and moist mild conditions. Postglacial sea‐level rise breached the crater rim and deposited 36 m of estuarine mud after ca. 9000 cal. a BP. Copyright © 2011 John Wiley & Sons, Ltd.     

Fracture mechanics analysis of multiple cracks in anisotropic media

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.     

基于垂直剖面函数的天顶对流层延迟插值算法

通过分析由ERA-Interim气象再分析资料积分方法得到的天顶对流层总延迟随高程变化的规律,提出一种基于垂直剖面函数的天顶对流层延迟(ZTD)插值算法。该算法以ZTD的垂直分布规律为基础,通过垂直剖面函数实现ZTD在高程方向上的精准投影延拓,可以避免因高差较大造成的空间内插结构畸形。采用IGS站提供的高精度对流层产品进行实验验证表明,该算法相对于传统算法能够有效提高ZTD改正值的精度,尤其在高差超过1 km的情况下,相对于反距离加权法精度提升了96%,相对于空间回归法精度提升了79%。     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

Vegetation history of the penultimate glacial period (Marine isotope stage 6) at Ioannina,north‐west Greece

The vegetational history of the penultimate glacial period, Marine Isotope Stage (MIS) 6 (c. 185–135 ka), has remained relatively unexplored. Here we present a new record from the Ioannina basin, north‐west Greece, which constitutes the highest‐resolution terrestrial pollen record for this interval produced to date. It shows that the vegetation history of MIS 6 in this region can be divided into two parts: an early period (185–155 ka) with pronounced oscillations in tree population extent, and a later period (155–135 ka) with much smaller tree populations and subdued oscillations. This pattern is analogous to the MIS 3/MIS 2 division during the last glacial in the same sequence, although the early part of MIS 6 had larger Pinus populations and fewer temperate trees relative to the equivalent interval in MIS 3. This implies cooler and wetter conditions, which is somewhat counterintuitive given the high summer insolation during MIS 6e, but is in line with other palaeoclimatic evidence from the Mediterranean. Comparison with North Atlantic records suggests that despite the absence of pronounced iceberg discharges during MIS 6, North Atlantic millennial‐scale variability had a significant downstream impact on tree populations in north‐west Greece. Copyright © 2011 John Wiley & Sons, Ltd.