The WIM method for refraction statics

Geophysics ◽  
1996 ◽  
Vol 61 (6) ◽  
pp. 1859-1870 ◽  
Author(s):  
Luigi Zanzi
Keyword(s):  
Real Data ◽  
Linear Operators ◽  
Surface Model ◽  
Near Surface ◽  
Model Driven ◽  
Iterative Model ◽  
Interpolation Problems ◽  
Final Optimization ◽  
Irregular Spacing ◽  
Reciprocal Method

We know that first breaks may be efficiently inverted in the wavenumber‐offset domain to determine a robust estimate of the near‐surface model. This domain is indicated particularly for regularly covered surveys so that data can be Fourier transformed without interpolation problems. However, it usually happens that data are collected with irregular spacing of the source points. In this case the initial interpolation may bias the solution. An efficient method to prevent bias consists of an iterative model‐driven interpolation. The database is decomposed and recomposed iteratively according to a linear approach to fill in the data space with model‐consistent first arrivals. The method is proved to converge to the unbiased solution and will be called the wavenumber iterative modeling (WIM) method. Experiments, both on synthetic and real data, show that convergence occurs after few iterations so that the procedure is still quite efficient. Moreover, the WIM method can be extended easily to automatically detect and remove the mispicks and suppress coherent and incoherent noise. Simple linear operators are also available for the conversion of the solution to the generalized reciprocal method (GRM) if the interpreter intends to guide the final optimization of the model.

scholarly journals Near-surface attenuation estimation using wave-propagation modeling

Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. U27-U37 ◽  
Author(s):  
Nizare El Yadari ◽  
Fabian Ernst ◽  
Wim Mulder
Keyword(s):  
Wave Attenuation ◽  
Waveform Inversion ◽  
Real Data ◽  
P Wave ◽  
Surface Model ◽  
Near Surface ◽  
Propagation Modeling ◽  
Static Corrections ◽  
Starting Model ◽  
Wave Propagation Modeling

The effect of the near surface on seismic land data can be so severe that static corrections are insufficient. Full-waveform inversion followed by redatuming may be an alternative, but inversion will work only if the starting model is sufficiently close to the true model. As a first step toward determining a viscoelastic near-surface model, we assume that existing methods can provide a horizontally layered velocity and density model. Because near-surface attenuation is strongest, we propose a method to estimate the P-wave attenuation based on viscoacoustic finite-difference modeling. We compare energy decay along traveltime curves of reflection and refraction events in the modeled and observed seismic data for a range of attenuation parameters. The best match provides an estimate of the attenuation. First, we estimate only the attenuation of the top layer and study the sensitivity to depth and velocity perturbations. Then, we consider multiple layers. We apply the method to synthetic and real data and investigate the effect of source wavelet and topography. The method is robust against depth and velocity perturbations smaller than 10%. The results are sensitive to the source wavelet. Incorporating the surface topography in the computed traveltimes reduces the uncertainty of the attenuation estimates, especially for deeper layers.

Regularization and datuming of seismic data by weighted, damped least squares

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. U67-U76 ◽  
Author(s):  
Robert J. Ferguson
Keyword(s):  
Least Squares ◽  
Seismic Data ◽  
Synthetic Data ◽  
Real Data ◽  
Structural Data ◽  
Irregular Sampling ◽  
Data Set ◽  
Near Surface ◽  
Damped Least Squares ◽  
Wavefield Extrapolation

The possibility of improving regularization/datuming of seismic data is investigated by treating wavefield extrapolation as an inversion problem. Weighted, damped least squares is then used to produce the regularized/datumed wavefield. Regularization/datuming is extremely costly because of computing the Hessian, so an efficient approximation is introduced. Approximation is achieved by computing a limited number of diagonals in the operators involved. Real and synthetic data examples demonstrate the utility of this approach. For synthetic data, regularization/datuming is demonstrated for large extrapolation distances using a highly irregular recording array. Without approximation, regularization/datuming returns a regularized wavefield with reduced operator artifacts when compared to a nonregularizing method such as generalized phase shift plus interpolation (PSPI). Approximate regularization/datuming returns a regularized wavefield for approximately two orders of magnitude less in cost; but it is dip limited, though in a controllable way, compared to the full method. The Foothills structural data set, a freely available data set from the Rocky Mountains of Canada, demonstrates application to real data. The data have highly irregular sampling along the shot coordinate, and they suffer from significant near-surface effects. Approximate regularization/datuming returns common receiver data that are superior in appearance compared to conventional datuming.

scholarly journals Climate pattern-scaling set for an ensemble of 22 GCMs – adding uncertainty to the IMOGEN version 2.0 impact system

2018 ◽  
Vol 11 (2) ◽  
pp. 541-560 ◽  
Author(s):  
Przemyslaw Zelazowski ◽  
Chris Huntingford ◽  
Lina M. Mercado ◽  
Nathalie Schaller
Keyword(s):  
Climate Change ◽  
Relative Humidity ◽  
Land Surface ◽  
Land Surface Model ◽  
Individual Performance ◽  
Balance Model ◽  
Surface Model ◽  
Near Surface ◽  
Global Circulation Models ◽  
Version 2.0

Abstract. Global circulation models (GCMs) are the best tool to understand climate change, as they attempt to represent all the important Earth system processes, including anthropogenic perturbation through fossil fuel burning. However, GCMs are computationally very expensive, which limits the number of simulations that can be made. Pattern scaling is an emulation technique that takes advantage of the fact that local and seasonal changes in surface climate are often approximately linear in the rate of warming over land and across the globe. This allows interpolation away from a limited number of available GCM simulations, to assess alternative future emissions scenarios. In this paper, we present a climate pattern-scaling set consisting of spatial climate change patterns along with parameters for an energy-balance model that calculates the amount of global warming. The set, available for download, is derived from 22 GCMs of the WCRP CMIP3 database, setting the basis for similar eventual pattern development for the CMIP5 and forthcoming CMIP6 ensemble. Critically, it extends the use of the IMOGEN (Integrated Model Of Global Effects of climatic aNomalies) framework to enable scanning across full uncertainty in GCMs for impact studies. Across models, the presented climate patterns represent consistent global mean trends, with a maximum of 4 (out of 22) GCMs exhibiting the opposite sign to the global trend per variable (relative humidity). The described new climate regimes are generally warmer, wetter (but with less snowfall), cloudier and windier, and have decreased relative humidity. Overall, when averaging individual performance across all variables, and without considering co-variance, the patterns explain one-third of regional change in decadal averages (mean percentage variance explained, PVE, 34.25±5.21), but the signal in some models exhibits much more linearity (e.g. MIROC3.2(hires): 41.53) than in others (GISS_ER: 22.67). The two most often considered variables, near-surface temperature and precipitation, have a PVE of 85.44±4.37 and 14.98±4.61, respectively. We also provide an example assessment of a terrestrial impact (changes in mean runoff) and compare projections by the IMOGEN system, which has one land surface model, against direct GCM outputs, which all have alternative representations of land functioning. The latter is noted as an additional source of uncertainty. Finally, current and potential future applications of the IMOGEN version 2.0 modelling system in the areas of ecosystem modelling and climate change impact assessment are presented and discussed.

Development of stochastically perturbed parameterization scheme for the Noah Land Surface Model with the optimized random forcing parameters using the micro-genetic algorithm

2021 ◽  
Author(s):  
Sujeong Lim ◽  
Claudio Cassardo ◽  
Seon Ki Park
Keyword(s):  
Genetic Algorithm ◽  
Land Surface ◽  
Weather Forecasting ◽  
Land Surface Model ◽  
Surface Model ◽  
Ensemble Spread ◽  
Near Surface ◽  
Random Forcing ◽  
Micro Genetic Algorithm ◽  
Stochastically Perturbed

<p>The ensemble data assimilation system is beneficial to represent the initial uncertainties and flow-dependent background error covariance (BEC). In particular, the inevitable model uncertainties can be expressed by ensemble spread, that is the standard deviation of ensemble BEC. However, the ensemble spread generally suffers from under-estimated problems. To alleviate this problem, recent studies employed stochastic perturbation schemes to increases the ensemble spreads by adding the random forcing in the model tendencies (i.e., physical or dynamical tendencies) or parameterization schemes (i.e., PBL, convective scheme, etc.). In this study, we focus on the near-surface uncertainties which are affected by the interactions between the land and atmosphere process. The land surface model (LSM) provides various fluxes as the lower boundary condition to the atmosphere, influencing the accuracy of hourly-to-seasonal scale weather forecasting, but the surface uncertainties were not much addressed yet. In this study, we developed the stochastically perturbed parameterization (SPP) scheme for the Noah LSM. The Weather Research and Forecasting (WRF) ensemble system is used for regional weather forecasting over East Asia, especially over the Korean Peninsula. As a testbed experiment with the newly-developed Noah LSM-SPP system, we first perturbed the soil temperature — a crucial variable for the near-surface forecasts by affecting sensible heat fluxes, land surface skin temperature and surface air temperature, and hence lower-tropospheric temperature. Here, the random forcing used in perturbation is made by the tuning parameters for amplitude, length scale, and time scales: they are commonly determined empirically by trial and error. In order to find optimal tuning parameter values, we applied a global optimization algorithm — the micro-genetic algorithm (micro-GA) — to achieve the smallest root-mean-squared errors. Our results indicate that optimization of the random forcing parameters contributes to an increase in the ensemble spread and a decrease in the ensemble mean errors in the near-surface and lower-troposphere uncertainties. Further experiments will be conducted by including soil moisture in the testbed.</p>

FULL-WAVEFORM INVERSION WITH SOURCE AND RECEIVER COUPLING EFFECTS CORRECTION

Geophysics ◽  
2021 ◽  
pp. 1-37
Author(s):  
Xinhai Hu ◽  
Wei Guoqi ◽  
Jianyong Song ◽  
Zhifang Yang ◽  
Minghui Lu ◽  
...  
Keyword(s):  
Waveform Inversion ◽  
Nonlinear Least Squares ◽  
Real Data ◽  
Model Parameters ◽  
Full Waveform Inversion ◽  
Linear Inversion ◽  
Near Surface ◽  
Model Parameter ◽  
Full Waveform ◽  
Coupling Factors

Coupling factors of sources and receivers vary dramatically due to the strong heterogeneity of near surface, which are as important as the model parameters for the inversion success. We propose a full waveform inversion (FWI) scheme that corrects for variable coupling factors while updating the model parameter. A linear inversion is embedded into the scheme to estimate the source and receiver factors and compute the amplitude weights according to the acquisition geometry. After the weights are introduced in the objective function, the inversion falls into the category of separable nonlinear least-squares problems. Hence, we could use the variable projection technique widely used in source estimation problem to invert the model parameter without the knowledge of source and receiver factors. The efficacy of the inversion scheme is demonstrated with two synthetic examples and one real data test.

scholarly journals High-resolution fully-coupled atmospheric–hydrological modeling: a cross-compartment regional water and energy cycle evaluation

2020 ◽  
Author(s):  
Benjamin Fersch ◽  
Alfonso Senatore ◽  
Bianca Adler ◽  
Joël Arnault ◽  
Matthias Mauder ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Land Surface ◽  
Hydrological Modeling ◽  
Weather Prediction ◽  
Surface Model ◽  
Near Surface ◽  
Ground Heat Flux ◽  
The Impact ◽  
Fully Coupled ◽  
Regional Water

<p>The land surface and the atmospheric boundary layer are closely intertwined with respect to the exchange of water, trace gases and energy. Nonlinear feedback and scale dependent mechanisms are obvious by observations and theories. Modeling instead is often narrowed to single compartments of the terrestrial system or bound to traditional viewpoints of definite scientific disciplines. Coupled terrestrial hydrometeorological modeling systems attempt to overcome these limitations to achieve a better integration of the processes relevant for regional climate studies and local area weather prediction. We examine the ability of the hydrologically enhanced version of the Weather Research and Forecasting Model (WRF-Hydro) to reproduce the regional water cycle by means of a two-way coupled approach and assess the impact of hydrological coupling with respect to a traditional regional atmospheric model setting. It includes the observation-based calibration of the hydrological model component (offline WRF-Hydro) and a comparison of the classic WRF and the fully coupled WRF-Hydro models both with identical calibrated parameter settings for the land surface model (Noah-MP). The simulations are evaluated based on extensive observations at the pre-Alpine Terrestrial Environmental Observatory (TERENO Pre-Alpine) for the Ammer (600 km²) and Rott (55 km²) river catchments in southern Germany, covering a five month period (Jun–Oct 2016).</p><p>The sensitivity of 7 land surface parameters is tested using the <em>Latin-Hypercube One-factor-At-a-Time</em> (LH-OAT) method and 6 sensitive parameters are subsequently optimized for 6 different subcatchments, using the Model-Independent <em>Parameter Estimation and Uncertainty Analysis software</em> (PEST).</p><p>The calibration of the offline WRF-Hydro leads to Nash-Sutcliffe efficiencies between 0.56 and 0.64 and volumetric efficiencies between 0.46 and 0.81 for the six subcatchments. The comparison of classic WRF and fully coupled WRF-Hydro shows only tiny alterations for radiation and precipitation but considerable changes for moisture- and energy fluxes. By comparison with TERENO Pre-Alpine observations, the fully coupled model slightly outperforms the classic WRF with respect to evapotranspiration, sensible and ground heat flux, near surface mixing ratio, temperature, and boundary layer profiles of air temperature. The subcatchment-based water budgets show uniformly directed variations for evapotranspiration, infiltration excess and percolation whereas soil moisture and precipitation change randomly.</p>

scholarly journals Sensitivity of Near-Surface Temperature Forecasts to Soil Properties over a Sparsely Vegetated Dryland Region

2014 ◽  
Vol 53 (8) ◽  
pp. 1976-1995 ◽  
Author(s):  
Jeffrey D. Massey ◽  
W. James Steenburgh ◽  
Sebastian W. Hoch ◽  
Jason C. Knievel
Keyword(s):  
Thermal Conductivity ◽  
Soil Moisture ◽  
Sandy Loam ◽  
Early Morning ◽  
Surface Model ◽  
Silt Loam ◽  
Soil Thermal Conductivity ◽  
Near Surface ◽  
Ground Heat Flux ◽  
Loam Soil

AbstractWeather Research and Forecasting Model forecasts over the Great Salt Lake Desert erroneously underpredict nocturnal cooling over the sparsely vegetated silt loam soil area of Dugway Proving Ground in northern Utah, with a mean positive bias error in temperature at 2 m AGL of 3.4°C in the early morning [1200 UTC (0500 LST)]. Positive early-morning bias errors also exist in nearby sandy loam soil areas. These biases are related to the improper initialization of soil moisture and parameterization of soil thermal conductivity in silt loam and sandy loam soils. Forecasts of 2-m temperature can be improved by initializing with observed soil moisture and by replacing Johansen's 1975 parameterization of soil thermal conductivity in the Noah land surface model with that proposed by McCumber and Pielke in 1981 for silt loam and sandy loam soils. Case studies illustrate that this change can dramatically reduce nighttime warm biases in 2-m temperature over silt loam and sandy loam soils, with the greatest improvement during periods of low soil moisture. Predicted ground heat flux, soil thermal conductivity, near-surface radiative fluxes, and low-level thermal profiles also more closely match observations. Similar results are anticipated in other dryland regions with analogous soil types, sparse vegetation, and low soil moisture.

scholarly journals The impact of near-surface soil moisture assimilation at subseasonal, seasonal, and inter-annual timescales

2015 ◽  
Vol 19 (12) ◽  
pp. 4831-4844 ◽  
Author(s):  
C. Draper ◽  
R. Reichle
Keyword(s):  
Soil Moisture ◽  
Data Assimilation ◽  
Surface Soil ◽  
Root Zone ◽  
Surface Model ◽  
Surface Soil Moisture ◽  
Near Surface ◽  
In Situ Observations ◽  
The Impact

Abstract. A 9 year record of Advanced Microwave Scanning Radiometer – Earth Observing System (AMSR-E) soil moisture retrievals are assimilated into the Catchment land surface model at four locations in the US. The assimilation is evaluated using the unbiased mean square error (ubMSE) relative to watershed-scale in situ observations, with the ubMSE separated into contributions from the subseasonal (SMshort), mean seasonal (SMseas), and inter-annual (SMlong) soil moisture dynamics. For near-surface soil moisture, the average ubMSE for Catchment without assimilation was (1.8 × 10−3 m3 m−3)2, of which 19 % was in SMlong, 26 % in SMseas, and 55 % in SMshort. The AMSR-E assimilation significantly reduced the total ubMSE at every site, with an average reduction of 33 %. Of this ubMSE reduction, 37 % occurred in SMlong, 24 % in SMseas, and 38 % in SMshort. For root-zone soil moisture, in situ observations were available at one site only, and the near-surface and root-zone results were very similar at this site. These results suggest that, in addition to the well-reported improvements in SMshort, assimilating a sufficiently long soil moisture data record can also improve the model representation of important long-term events, such as droughts. The improved agreement between the modeled and in situ SMseas is harder to interpret, given that mean seasonal cycle errors are systematic, and systematic errors are not typically targeted by (bias-blind) data assimilation. Finally, the use of 1-year subsets of the AMSR-E and Catchment soil moisture for estimating the observation-bias correction (rescaling) parameters is investigated. It is concluded that when only 1 year of data are available, the associated uncertainty in the rescaling parameters should not greatly reduce the average benefit gained from data assimilation, although locally and in extreme years there is a risk of increased errors.

scholarly journals Effects of Historical Urbanization in the Brussels Capital Region on Surface Air Temperature Time Series: A Model Study

2009 ◽  
Vol 48 (10) ◽  
pp. 2181-2196 ◽  
Author(s):  
R. Hamdi ◽  
A. Deckmyn ◽  
P. Termonia ◽  
G. R. Demarée ◽  
P. Baguis ◽  
...  
Keyword(s):  
Land Surface ◽  
Minimum Temperature ◽  
Computational Cost ◽  
Land Surface Model ◽  
Maximum Temperature ◽  
Surface Model ◽  
Impervious Surfaces ◽  
Urban Bias ◽  
Near Surface ◽  
Capital Region

Abstract The authors examine the local impact of change in impervious surfaces in the Brussels capital region (BCR), Belgium, on trends in maximum, minimum, and mean temperatures between 1960 and 1999. Specifically, data are combined from remote sensing imagery and a land surface model including state-of-the-art urban parameterization—the Town Energy Balance scheme. To (i) isolate effects of urban growth on near-surface temperature independent of atmospheric circulations and (ii) be able to run the model over a very long period without any computational cost restrictions, the land surface model is run in a stand-alone mode coupled to downscaled 40-yr ECMWF reanalysis data. BCR was considered a lumped urban volume and the rate of urbanization was assessed by estimating the percentage of impervious surfaces from Landsat images acquired for various years. Model simulations show that (i) the annual mean urban bias (AMUB) on minimum temperature is rising at a higher rate (almost 3 times more) than on maximum temperature, with a linear trend of 0.14° and 0.05°C (10 yr)−1, respectively, (ii) the 40-yr AMUB on mean temperature is estimated to be 0.62°C, (iii) 45% of the overall warming trend is attributed to intensifying urban heat island effects rather than to changes in local–regional climate, and (iv) during summertime, a stronger dependence between the increase of urban bias on minimum temperature and the change in percentage of impervious surfaces is found.

Sign in / Sign up

Export Citation Format

Share Document