Scattering and reflection of SH waves around a slope on an elastic wedged space

The scattering and reflection of SH waves by a slope on an elastic wedged space is investigated. A series solution is obtained by using the wave function expansion method. The slope on a wedged space is divided into two subregions by an artificial, auxiliary circular arc. The wave fields with unknown complex coefficients within each sub-region are derived. Applying Graf addition theorem, the scattered waves in the sub-regions are expressed in a global coordinate system. Fourier transform is adopted to derive a consistent form of standing waves in the inner region using the orthogonality of the cosine functions. The boundary-valued problem is solved by stress and displacement continuity along the artificial, auxiliary arc to obtain the unknown complex coefficients. Parametric studies are next performed to investigate how the topography from the slope on the wedged space will affect the scattering and diffraction, and hence the amplification and de-amplification of the SH waves. Numerical results show that the surface motions on the slope of the wedged space is influenced greatly by the topography. Amplification of the surface motions near the slope vertex is significant. The corresponding phases along the wedged space surfaces are consistent with the direction that the SH waves are propagating.     

The role of hydrogeological setting in two Canadian peatlands investigated through 2D steady-state groundwater flow modelling

This study investigated how hydrogeological setting influences aquifer–peatland connections in slope and basin peatlands. Steady-state groundwater flow was simulated using Modflow on 2D transects for an esker slope peatland and for a basin peatland in southern Quebec (Canada). Simulations investigated how hydraulic heads and groundwater flow exported toward runoff from the peatland can be influenced by recharge, hydraulic properties, and heterogeneity. The slope peatland model was strongly dominated by horizontal flow from the esker. This suggests that slope peatlands are dependent on the hydrogeological conditions of the adjacent aquifer reservoir, but are resilient to hydrological changes. The basin peatland produced groundwater outflow to the surface aquifer. Lateral and vertical peat heterogeneity due to peat decomposition or compaction were identified as having a significant influence on fluxes. These results suggest that basin peatlands are more dependent on recharge conditions, and could be more susceptible to land use and climate changes.     

Performance of the LandSoil expert-based model to map erosion and sedimentation: application to a cultivated catchment in central Belgium

Intensive agricultural practices on sensitive soils induce high erosion rates in central Belgium. Expert-rules models quantify runoff and erosion at catchment scale, avoiding over-parameterization, and can include some direct or indirect connectivity features. The aim of this article is to test the ability of an expert-based model, LandSoil, to quantify runoff and to locate erosion and sedimentation areas in a small cultivated loamy catchment in Belgium during the years 2014, 2015 and 2016. Spatialized data are important for assessing model outputs and the erosive response. Measurements of runoff and observation of spatial erosion/deposition patterns, especially around major connectivity points, permitted an assessment of the reliability of the model results. Runoff modelling gave contrasting results (good linear adjustment at the outlet of the 83 ha sub-catchment (point 1): r² of 0.96, Nash–Sutcliffe criterion of 0.95; less good at the outlet of the 3.9 ha sub-catchment (point 2): r² of 0.28, Nash–Sutcliffe criterion of –0.47). For point 2 the poor results are explained by the very few runoff events observed, a scaling effect and the small area with a single land use. Graduated rulers demonstrate that the model is able to provide a coherent pattern of erosion/deposition. The study highlights great sensitivity to the effect of land use, land allocation, landscape design and slope gradients. Grass strips induce deposition of eroded particles when slopes are gentle (< 2%). Woodland strips decrease connectivity by being in the stream but deposit thinner sediment layers. Field boundaries have a role in the transport, but not really the quantity, of sediments. This model validation in the Belgian loess context allows us to use LandSoil in other similar environments in order to estimate the effects of landscape management scenarios. © 2020 John Wiley & Sons, Ltd.     

Geochemical characterization of soil organic matter and variability of a postglacial detrital organic supply (chaillexon lake,france)

Geochemical analysis of sedimentary organic matter in recent lacustrine sediments appears to be a useful tool in providing information concerning past environmental conditions. However, such analysis is often made without knowing the geochemical characteristics of the organic matter derived from the watershed and, more explicitly, its soils. The present work deals with (i) a geochemical investigation (Rock-Eval pyrolysis) of soil organic matter sampled in a lake watershed, and (ii) the study of the sedimentary organic matter trapped in the lake deposits. The research was conducted on Chaillexon Lake which was created by a rock collapse that dammed the palaeovalley of the Doubs River about 12 000 years ago. Since this event, the sediment trap provides a continuous palaeoclimatic record for the Postglacial period. Results obtained lead to two main conclusions. First, the variability of Rock-Eval pyrolysis values observed in soils modifies the common interpretation given to these parameters in the characterization of sedimentary organic matter. Indeed, variations in these parameters point not only to varying proportions of terrestrial and lacustrine organic matter in a lacustrine infilling but also to variations of the terrestrial supply linked with the evolution of vegetal cover in the catchment. The second conclusion is that the story of the Chaillexon lacustrine system is marked by a rather sudden soil and forest development at the Preboreal–Boreal transition (9000 BP ). Copyright © 1998 John Wiley & Sons, Ltd.     

The effects of resolving the Canadian Arctic Archipelago in a finite element sea ice model

Though narrow straits may have a strong influence on the large-scale sea ice mass balance, they are often crudely represented in coarse resolution sea ice models. Unstructured meshes, with their natural ability to fit boundaries and locally increase the mesh resolution, propose an alternative framework to capture the complex oceanic areas formed by coasts and islands. In this paper, we develop a finite element sea ice model to investigate the sensitivity of the Arctic sea ice cover features to the resolution of the narrow straits constituting the Canadian Arctic Archipelago. The model is a two-level dynamic-thermodynamic sea ice model, including a viscous-plastic rheology. It is run over 1979–2005, forced by daily NCEP/NCAR reanalysis data. Confronting qualitatively numerical experiments with observations shows a good agreement with satellite and buoys measurements. Due to its simple representation of the oceanic interactions, the model overestimates the sea ice extent during winter in the southernmost parts of the Arctic, while the Baffin Bay and Kara Sea remain ice-covered during summer. In order to isolate the benefits from resolving the Canadian Arctic Archipelago, a numerical experiment is performed where we artificially close the archipelago. Focusing on the large-scale sea ice thickness pattern, no significant change is found in our model, except in the close surroundings of the archipelago. However, the local and short-term influences of the ice exchanges are nonnegligible. In particular, we show that the ice volume associated to the Canadian Arctic Archipelago represents 10% of the Northern Hemisphere sea ice volume and that the annual mean ice export towards Baffin Bay amounts to 125 km³ yr⁻¹, which may play an important role on the convective overturning in the Labrador Sea.     

Preliminary 10Be chronology for the last deglaciation of the western margin of the Greenland Ice Sheet

The now acknowledged thinning of the Greenland Ice Sheet raises concerns about its potential contribution to future sea level rise. In order to appreciate the full extent of its contribution to sea level rise, reconstruction of the ice sheet's most recent last deglaciation could provide key information on the timing and the height of the ice sheet at a time of rapid climate readjustment. We measured ¹⁰Be concentrations in 12 samples collected along longitudinal and altitudinal transects from Sisimiut to within 10 km of the Isunguata Sermia Glacier ice margin on the western coast of Greenland. Along the longitudinal transect, we collected three perched boulders and two bedrocks. In addition, we sampled seven perched boulders along a vertical transect in a valley within 10 km of the Isunguata Sermia Glacier ice margin. Our pilot dataset constrains the height of the ice sheet during the Last Glacial Maximum (LGM) between 500 m and 840 m (including the 120 m relative sea level depression at the time of the LGM, 21 ka BP). From the transect we estimate the thinning of the ice sheet at the end of the deglaciation between 12.3 ± 1.5 ¹⁰Be ka (n = 2) and 8.3 ± 1.2 ¹⁰Be ka (n = 3) to be ～6 cm a^?1 over this time period. Direct dating of the retreat of the western margin of the Greenland Ice Sheet has the potential to better constrain the retreat rate of the ice margin, the thickness of the former ice sheet as well as its response to climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

An improved geomagnetic data selection algorithm for global geomagnetic field modelling

A spherical harmonic degrees 60, global internal field model is described (called BGS/G/L/0706). This model includes a degree 15 core and piecewise-linear secular variation model and is derived from quiet-time Ørsted and Champ satellite data sampled between 2001.0 and 2005.0. For the satellite data selection, a wide range of geomagnetic index and other data selection filters have been used to best isolate suitably quiet magnetospheric and ionospheric conditions. Only a relatively simple, degree one spherical harmonic, external field model is then required. It is found that a new 'Vector Magnetic Disturbance' index ( VMD ), the existing longitude sector A indices, the auroral zone index IE , and the polar cap index PC are better than Kp and Dst at rejecting rapidly varying external field signals at low, middle, auroral and polar latitudes. The model quality is further enhanced by filling spatial and temporal gaps in the quiet data selection with a second selection containing slightly more disturbed data. It is shown that VMD provides a better parametrization than Dst of the large-scale, rapidly changing, external field. The lithospheric field model between degrees 16 and 50 is robust and displays good coherence with other recently published models for this epoch. BGS/G/L/0706 also shows crustal anomalies consistent with other studies, although agreement is poorer in the southern polar cap. Intermodel coherency reduces above about degree 40, most likely due to incompletely filtered signals from polar ionospheric currents and auroral field aligned currents. The absence of the PC index for the southern hemisphere for 2003 onwards is a particular concern.     

平面P波在饱和半空间中洞室周围的散射(I):解析解

利用波函数展开法给出了入射平面P波在饱和半空间中圆形洞室周围散射问题的一个解析解。半空间假定为无粘性流体饱和介质,满足Biot理论。采用一种基于实验数据的孔隙率和模量之间的线性关系来确定Biot模型中的介质参数。解答考虑了透水边界和非透水边界两种情况。对边界条件进行了数值检验,结果表明,随着级数截断项数的增大,边界残量衰减很快。解答为进一步研究入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响奠定了基础。     

Sensitivity equations for hyperbolic conservation law-based flow models

The present paper focuses on the governing equations for the sensitivity of the variables to the parameters in flow models that can be described by one-dimensional scalar, hyperbolic conservation laws. The sensitivity is shown to obey a hyperbolic, scalar conservation law. The sensitivity is a conserved scalar except in the case of discontinuous flow solutions, where an extra, point source term must be added to the equations in order to enforce conservation. The propagation speed of the sensitivity waves being identical to that of the conserved variable in the original conservation law, the system of conservation laws formed by the original hyperbolic equation and the equation satisfied by the sensitivity is linearly degenerate. A consequence on the solution of the Riemann problem is that rarefaction waves for the variable of the original equation result in vacuum regions for the sensitivity. The numerical solution of the hyperbolic conservation law for the sensitivity by finite volume methods requires the implementation of a specific shock detection procedure. A set of necessary conditions is defined for the discretisation of the source term in the sensitivity equation. An application to the one-dimensional kinematic wave equation shows that the proposed numerical technique allows analytical solutions to be reproduced correctly. The computational examples show that first-order numerical schemes do not yield satisfactory numerical solutions in the neighbourhood of moving shocks and that higher-order schemes, such as the MUSCL scheme, should be used for sharp transients.