Requirements and Boundary Conditions for Near Surface Seismic Inversion

Author(s):  
N. Römer-Stange ◽  
V. Spieß
2006 ◽  
Vol 19 (15) ◽  
pp. 3751-3767 ◽  
Author(s):  
Véronique Bugnion ◽  
Chris Hill ◽  
Peter H. Stone

Abstract Multicentury sensitivities in a realistic geometry global ocean general circulation model are analyzed using an adjoint technique. This paper takes advantage of the adjoint model’s ability to generate maps of the sensitivity of a diagnostic (i.e., the meridional overturning’s strength) to all model parameters. This property of adjoints is used to review several theories, which have been elaborated to explain the strength of the North Atlantic’s meridional overturning. This paper demonstrates the profound impact of boundary conditions in permitting or suppressing mechanisms within a realistic model of the contemporary ocean circulation. For example, the so-called Drake Passage Effect in which wind stress in the Southern Ocean acts as the main driver of the overturning’s strength, is shown to be an artifact of boundary conditions that restore the ocean’s surface temperature and salinity toward prescribed climatologies. Advective transports from the Indian and Pacific basins play an important role in setting the strength of the overturning circulation under “mixed” boundary conditions, in which a flux of freshwater is specified at the ocean’s surface. The most “realistic” regime couples an atmospheric energy and moisture balance model to the ocean. In this configuration, inspection of the global maps of sensitivity to wind stress and diapycnal mixing suggests a significant role for near-surface Ekman processes in the Tropics. Buoyancy also plays an important role in setting the overturning’s strength, through direct thermal forcing near the sites of convection, or through the advection of salinity anomalies in the Atlantic basin.


2017 ◽  
Vol 30 (15) ◽  
pp. 6017-6036 ◽  
Author(s):  
Sachie Kanada ◽  
Tetsuya Takemi ◽  
Masaya Kato ◽  
Shota Yamasaki ◽  
Hironori Fudeyasu ◽  
...  

Intense tropical cyclones (TCs) sometimes cause huge disasters, so it is imperative to explore the impacts of climate change on such TCs. Therefore, the authors conducted numerical simulations of the most destructive historical TC in Japanese history, Typhoon Vera (1959), in the current climate and a global warming climate. The authors used four nonhydrostatic models with a horizontal resolution of 5 km: the cloud-resolving storm simulator, the fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model, the Japan Meteorological Agency (JMA) operational nonhydrostatic mesoscale model, and the Weather Research and Forecasting Model. Initial and boundary conditions for the control simulation were provided by the Japanese 55-year Reanalysis dataset. Changes between the periods of 1979–2003 and 2075–99 were estimated from climate runs of a 20-km-mesh atmospheric general circulation model, and these changes were added to the initial and boundary conditions of the control simulation to produce the future climate conditions. Although the representation of inner-core structures varies largely between the models, all models project an increase in the maximum intensity of future typhoons. It is found that structural changes only appeared around the storm center with sudden changes in precipitation and near-surface wind speeds as the radius of maximum wind speed (RMW) contracted. In the future climate, the water vapor mixing ratio in the lower troposphere increased by 3–4 g kg−1. The increased water vapor allowed the eyewall updrafts to form continuously inside the RMW and contributed to rapid condensation in the taller and more intense updrafts.


Author(s):  
Brittany N. Carson-Marquis ◽  
Jianglong Zhang ◽  
Peng Xian ◽  
Jeffrey S. Reid ◽  
Jared Marquis

AbstractWhen unaccounted for in numerical weather prediction (NWP) models, heavy aerosol events can cause significant unrealized biases in forecasted meteorological parameters such as surface temperature. To improve near-surface forecasting accuracies during heavy aerosol loadings, we demonstrate the feasibility of incorporating aerosol fields from a global chemical transport model as initial and boundary conditions into a higher resolution NWP model with aerosol-meteorological coupling. This concept is tested for a major biomass burning smoke event over the Northern Great Plains region of the United States that occurred during summer of 2015. Aerosol analyses from the global Navy Aerosol Analysis and Prediction System (NAAPS) are used as initial and boundary conditions for Weather Research and Forecasting with Chemistry (WRF-Chem) simulations. Through incorporating more realistic aerosol direct effects into the WRF-Chem simulations, errors in WRF-Chem simulated surface downward shortwave radiative fluxes and near-surface temperature are reduced compared with surface-based observations. This study confirms the ability to decrease biases induced by the aerosol direct effect for regional NWP forecasts during high-impact aerosol episodes through the incorporation of analyses and forecasts from a global aerosol transport model.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. T93-T108 ◽  
Author(s):  
Stig Hestholm

Computational resources have increased in capacity over time — mostly by speed, partly by memory. Consequently, people have continuously explored the possibilities of performing wave modeling and inversion of increasing physical complexity. Achieving a detailed as possible image of the earth’s subsurface improves the success of hydrocarbon exploration, and it is important for other applications, such as archeology, mining, and engineering. I have developed an accurate computational method for elastic wave modeling up to tilted orthorhombic symmetry of anisotropy. The model may be covered by an arbitrary topographic function along the free surface. Through snapshots and seismograms of the wavefield, I confirm known effects from applying the code to plane, free surfaces (horizontal or tilted) as well as more complex topographies. The method is based on adapting a curved grid to a free-surface topography at hand, and transforming the wave equations and the topography free-surface boundary conditions from this grid to a rectangular grid, where finite-difference (FD) calculations can be performed. Free-surface topography boundary conditions for the particle velocities originate from locally setting the normal stress components to zero at the curved grid free surface. Vanishing normal traction is achieved by additionally imposing mirror conditions on stresses across the free surface. This leads me to achieve a more accurate modeling of free-surface waves (Rayleigh — Rg-waves in particular), using either FDs or any other numerical discretization method. Statics correction, muting, and destructive processing, which all consider free-surface effects as noise, can hence be avoided in inversion/imaging because surface effects can be more accurately simulated. By including near-surface effects in the full wavefield, we ultimately obtain superior inversion for interior earth materials, also for deeper physical medium properties.


2013 ◽  
Vol 134 ◽  
pp. 116-130 ◽  
Author(s):  
D. Akritidis ◽  
P. Zanis ◽  
E. Katragkou ◽  
M.G. Schultz ◽  
I. Tegoulias ◽  
...  

2018 ◽  
Author(s):  
Thomas Schneider von Deimling ◽  
Thomas Kleinen ◽  
Gustaf Hugelius ◽  
Christian Knoblauch ◽  
Christian Beer ◽  
...  

Abstract. We have developed a new module to calculate soil organic carbon (SOC) accumulation in perennially frozen ground in the land surface model JSBACH. Running this offline version of MPI-ESM we have modelled permafrost carbon accumulation and release from the Last Glacial Maximum (LGM) to the Pre-industrial (PI). Our simulated near-surface PI permafrost extent of 16.9 Mio km2 is close to observational evidence. Glacial boundary conditions, especially ice sheet coverage, result in profoundly different spatial patterns of glacial permafrost extent. Deglacial warming leads to large-scale changes in soil temperatures, manifested in permafrost disappearance in southerly regions, and permafrost aggregation in formerly glaciated grid cells. In contrast to the large spatial shift in simulated permafrost occurrence, we infer an only moderate increase of total LGM permafrost area (18.3 Mio km2) – together with pronounced changes in the depth of seasonal thaw. Reconstructions suggest a larger spread of glacial permafrost towards more southerly regions, but with a highly uncertain extent of non-continuous permafrost. Compared to a control simulation without describing the transport of SOC into perennially frozen ground, the implementation of our newly developed module for simulating permafrost SOC accumulation leads to a doubling of simulated LGM permafrost SOC storage (amounting to a total of ~ 150 PgC). Despite LGM temperatures favouring a larger permafrost extent, simulated cold glacial temperatures – together with low precipitation and low CO2 levels – limit vegetation productivity and therefore prevent a larger glacial SOC build-up in our model. Changes in physical and biogeochemical boundary conditions during deglacial warming lead to an increase in mineral SOC storage towards the Holocene (168 PgC at PI), which is below observational estimates (575 PgC in continuous and discontinuous permafrost). Additional model experiments clarified the sensitivity of simulated SOC storage to model parameters, affecting long-term soil carbon respiration rates and simulated active layer depths. Rather than a steady increase in carbon release from the LGM to PI as a consequence of deglacial permafrost degradation, our results suggest alternating phases of soil carbon accumulation and loss as an effect of dynamic changes in permafrost extent, active layer depths, soil litter input, and heterotrophic respiration.


2021 ◽  
Author(s):  
Kassem Asfour ◽  
Roland Martin ◽  
Ludovic Bodet ◽  
Didier El Baz ◽  
Bastien Plazolles ◽  
...  

<div><span>In order to accurately study the properties of partially saturated unconsolidated media at the near surface scale, or be able to image deeper structures through them, accurate 2D and 3D wave propagation numerical modelling tools are required. The rheology/mechanical properties of such media are frequently extremely complex (nonlinear, anisotropic, ... ), even when considered at dry state and of homogeneous mixture. Experimental observations (both at the laboratory and field scales) show that the seismic wave-field in unconsolidated granular materials remains difficult to interpret within standard methodological frameworks. We present here a numerical study aiming at exploring possible alternative forward modelling approaches to better extract information from recorded signals. <br></span></div><div><span> </span></div><div><span>We first present a finite volume method (Asfour et al. 2021) in which exact Riemann solvers are introduced. Solutions are compared to high-order finite-differences (Seismic_CPML code) and spectral finite element (SPECFEM code) solutions. A first series of synthetic cases is shown to benchmark the code at the hundred meters scale with a 100-300Hz wavelet source content. Another synthetic and more reallistic case is then presented with a medium affected by a steep nonlinear velocity gradient in depth, typical of an unconsolidated granular medium (as previously considered  at laboratory scale). For this model, a 1500Hz dominant frequency point source wavelet is considered and fluid saturation is also tested by applying a fluid-solid coupling. First arrival times and PSV-wave dispersion obtained from the different codes are compared. <br></span></div><div><span> </span></div><div><span>In a second step, and considering the real data recorded at the laboratory, we apply a more realistic source wavelet (obtained through signal spectrum ratio) and we perform parallelized high-order finite difference simulations (UNISOLVER code) to compare 2D and 3D elastic as well as poroelastic solutions on finely discretized meshes. Computed and observed data are compared. The poroelastic rheology provides better amplitudes in the seismograms and better exhibits some PSV modes in the phase velocity dispersion observations. Sensitivity kernels are also shown for the different rheologies. The different results obtained are now paving the way to seismic inversion at the near surface scale and to image shallow fluid/water saturated layers.</span></div>


2021 ◽  
Author(s):  
Hannah Watts ◽  
Adam Booth ◽  
Benedict Reinardy ◽  
Siobhan Killingbeck ◽  
Peter Jansson ◽  
...  

<p>Glacier forelands contain valuable information on past glacier dynamics and associated climatic conditions, particularly at small mountain glaciers where responses to climate change are rapid. To maximize the potential of glacial landforms as palaeoclimate indicators, a thorough understanding of the controls on landform genesis and subsequent evolution is required. Traditionally, such landforms have been studied using glacial geological techniques such as sedimentary logging. While these provide valuable in situ information they have numerous limitations, namely poor availability and spatial extent of exposures. Near-surface geophysics provides an efficient and non-invasive means of studying subsurface conditions in numerous sedimentary settings, offering spatially extensive information on substrate material properties and architecture. However, the logistically challenging terrain, remote location and complex structure of proglacial environments has limited the development of geophysical techniques for studying the internal architecture of glacial landforms.</p><p>Here, we explore the application of three geophysical methods to investigate proglacial substrates: ground penetrating radar (GPR), seismic refraction and multi-channel analysis of surface waves (MASW). Three sites with contrasting sediment properties were surveyed at the foreland of Midtdalsbreen glacier in southern Norway; (a) a 100 m<sup>2</sup> area of glaciotectonised sandy sediments, (b) a ~2 m high lateral moraine ridge containing stratified silts, sands, and gravel and (c) a terminal moraine ridge with a peak crest height of ~5 m and an open blockwork of cobbles and boulders at its surface. At all sites, we deployed 25 MHz and 100 MHz GPR antennas and undertook seismic surveys with 50−75 m long geophone spreads and a sledge-hammer source to sample to target depths of around 10−15 m. Through comparing the results from sites (a) to (c), we assess the capabilities and limitations of each of the aforementioned techniques for proglacial substrate imaging and characterisation, we analyse how their performances vary across these settings and outline factors that contribute to a successful geophysical investigation. </p><p>The ease of analysis and achievable investigation depths of the geophysical data and the applicability of seismic interpretation methods varied considerably depending on the surface terrain and structural complexity of the site. Our results show how the combination of GPR and seismic data can assist with the internal characterisation of glacial moraines when a relatively simple subsurface structure is present. However, basic seismic inversions likely lack the sophistication to resolve seismic structure in all but the simplest of layered models. We offer suggestions on how to optimise field time in more complex settings, where more sophisticated seismic inversion algorithms (e.g. tomography) or 3-D GPR surveys could be better-suited.</p><p>Our experience should help advance the use of geophysics in proglacial studies. It should serve as a guide for future survey planning, and help avoid typical pitfalls such that field time can be optimised.  It is hoped that geophysical survey methods will play an increasing role in the understanding of proglacial sedimentary landforms and their associated palaeoenvironments.</p>


Sign in / Sign up

Export Citation Format

Share Document