Home About us Contact
|
Home About us Contact | ||
|
|
||
Velocity Models (velocity + models)
Selected AbstractsSedimentary and crustal structure from the Ellesmere Island and Greenland continental shelves onto the Lomonosov Ridge, Arctic OceanGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2010H. Ruth Jackson SUMMARY On the northern passive margin of Ellesmere Island and Greenland, two long wide-angle seismic reflection/refraction (WAR) profiles and a short vertical incident reflection profile were acquired. The WAR seismic source was explosives and the receivers were vertical geophones placed on the sea ice. A 440 km long North-South profile that crossed the shelf, a bathymetric trough and onto the Lomonosov Ridge was completed. In addition, a 110 km long profile along the trough was completed. P -wave velocity models were created by forward and inverse modelling. On the shelf modelling indicates a 12 km deep sedimentary basin consisting of three layers with velocities of 2.1,2.2, 3.1,3.2 and 4.3,5.2 km s,1. Between the 3.1,3.2 km s,1 and 4.3,5.2 km s,1 layers there is a velocity discontinuity that dips seaward, consistent with a regional unconformity. The 4.3,5.2 km s,1 layer is interpreted to be Palaeozoic to Mesozoic age strata, based on local and regional geological constraints. Beneath these layers, velocities of 5.4,5.9 km s,1 are correlated with metasedimentary rocks that outcrop along the coast. These four layers continue from the shelf onto the Lomonosov Ridge. On the Ridge, the bathymetric contours define a plateau 220 km across. The plateau is a basement high, confirmed by short reflection profiles and the velocities of 5.9,6.5 km s,1. Radial magnetic anomalies emanate from the plateau indicating the volcanic nature of this feature. A lower crustal velocity of 6.2,6.7 km s,1, within the range identified on the Lomonosov Ridge near the Pole and typical of rifted continental crust, is interpreted along the entire line. The Moho, based on the WAR data, has significant relief from 17 to 27 km that is confirmed by gravity modelling and consistent with the regional tectonics. In the trough, Moho shallows eastward from a maximum depth of 19,16 km. No indication of oceanic crust was found in the bathymetric trough. [source] Waveform modelling of teleseismic S, Sp, SsPmP, and shear-coupled PL waves for crust- and upper-mantle velocity structure beneath AfricaGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2007Abhijit Gangopadhyay SUMMARY We describe a waveform modelling technique and demonstrate its application to determine the crust- and upper-mantle velocity structure beneath Africa. Our technique uses a parallelized reflectivity method to compute synthetic seismograms and fits the observed waveforms by a global optimization technique based on a Very Fast Simulated Annealing (VFSA). We match the S, Sp, SsPmP and shear-coupled PL phases in seismograms of deep (200,800 km), moderate-to-large magnitude (5.5,7.0) earthquakes recorded teleseismically at permanent broad-band seismic stations in Africa. Using our technique we produce P - and S -wave velocity models of crust and upper mantle beneath Africa. Additionally, our use of the shear-coupled PL phase, wherever observed, improves the constraints for lower crust- and upper-mantle velocity structure beneath the corresponding seismic stations. Our technique retains the advantages of receiver function methods, uses a different part of the seismogram, is sensitive to both P - and S -wave velocities directly, and obtains helpful constraints in model parameters in the vicinity of the Moho. The resulting range of crustal thicknesses beneath Africa (21,46 km) indicates that the crust is thicker in south Africa, thinner in east Africa and intermediate in north and west Africa. Crustal P - (4.7,8 km s,1) and S -wave velocities (2.5,4.7 km s,1) obtained in this study show that in some parts of the models, these are slower in east Africa and faster in north, west and south Africa. Anomalous crustal low-velocity zones are also observed in the models for seismic stations in the cratonic regions of north, west and south Africa. Overall, the results of our study are consistent with earlier models and regional tectonics of Africa. [source] P - and S -velocity images of the lithosphere,asthenosphere system in the Central Andes from local-source tomographic inversionGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2006Ivan Koulakov SUMMARY About 50 000 P and S arrival times and 25 000 values of t* recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P, S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp, Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting. [source] Crosswell seismic waveguide phenomenology of reservoir sands & shales at offsets >600 m, Liaohe Oil Field, NE ChinaGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2005P. C. Leary SUMMARY Crosswell seismic data recorded at 620,650 m offsets in an oil-bearing sand/shale reservoir formation at the Liaohe Oil Field, northeast China, provide robust evidence for waveguide action by low-velocity reservoir layers. Crosswell-section velocity models derived from survey-well sonic logs and further constrained by observed waveguide seismic wavegroup amplitudes and phases yield plausible evidence for interwell reservoir,sand continuity and discontinuity. A pair of back-to-back Liaohe crosswell vector-seismic surveys were conducted using a source well between two sensor wells at 650 and 620 m offsets along a 200-m-thick reservoir formation dipping 7° down-to-east between depths of 2.5 and 3 km. A downhole orbital vibrator generated seismic correlation wavelets with frequency range 50,350 Hz and signal/noise ratio up to 5:1 over local downhole ambient noise. The sensor wells were instrumented with a mobile 12- to 16-level string of clamped vector-motion sensor modules at 5 m intervals. Using 5 m source depth increments, crosswell Surveys 1 and 2 cover source/sensor well intervals above and through the reservoir of, respectively, 600 m/600 m (13 000 vector traces in 9 common sensor fans) and 300 m/560 m (7000 vector traces in 7 common sensor fans). Survey 1 common sensor gathers show clear, consistent high-amplitude 20 ms waveletgroup lags behind the first-arrival traveltime envelope. Such arrivals are diagnostic of seismic low-velocity waveguides connecting the source and sensor wells. Observed Survey 1 retarded wavegroup depths tally with source and sensor depths in low-velocity layers identified in sonic well logs. Finite-difference acoustic model wavefields computed for waveguide acoustic layers constrained by well-log sonic velocity data match the observed waveguide traveltime and amplitude systematics. Model waveforms duplicate the observed m-scale and ms-scale sensitivity of waveguide spatio-temporal energy localization. Survey 2 crosswell data, in contrast, provide no comparable evidence for waveguide action despite a sensor-well sonic well log similar to that of Survey 1. Instead, acoustic wavefield modelling of Survey 2 data clearly favours an interpreted waveguide model with 10° downdip interrupted by a 75,100 m throw down-fault near the sensor well. The absence of clear waveguide arrivals is adequately explained by dispersal of waveguide energy at the fault discontinuity. Auxiliary well sonic velocity and lithologic logs confirm the model-implied 75,100 m of down-throw faulting near the sensor well. [source] Streamer tomography velocity models for the Gulf of Corinth and Gulf of Itea, GreeceGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2004Barry C. Zelt SUMMARY The Gulf of Corinth (GOC), Greece is a rapidly extending, active continental rift with a record of large, damaging earthquakes. An extensive multichannel seismic (MCS) survey of the GOC conducted in 2001 provided, in addition to the processed MCS images, the opportunity to constrain velocity structure using refracted arrivals recorded along a 6-km-long streamer. We use first-arrival traveltimes to derive tomographic P -wave velocity models for several profiles collected in the central portion of the GOC. Eight of the profiles are closely spaced, north,south lines crossing the GOC and extending into the Gulf of Itea (GOI); a ninth profile is an east,west-oriented tie line. The N,S profiles image the relatively simple velocity structure of the deep Corinth rift basin and more complicated structure of the northern margin of the currently active rift. Integration of the velocity models with migrated MCS sections shows that south of the GOI the basement, which comprises Mesozoic nappes, occurs at a velocity of 4.5 km s,1 in the velocity models, although the actual velocity at, or just below, the top of basement is probably closer to 5,5.5 km s,1. The maximum sediment thickness in the Corinth basin is 2.2 km. The basement shallows to the north into a fault-bounded terrace in the central region between the two gulfs. Sediment cover in this central region decreases in thickness from west to east. Beneath the GOI, low average velocities beneath the rift-onset reflector indicate the presence of pre-rift sediments. The pre-rift velocity structure in the GOI is complex, with significant lateral variation from west to east. The E,W line shows that high-velocity basement is shallow (,1 km depth) and flat to the west of the GOI but dips ,20° east down to ,1.5 km beneath the pre-rift sediments of the GOI. [source] Constraints on earthquake epicentres independent of seismic velocity modelsGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2004T. Nicholson SUMMARY We investigate the constraints that may be placed on earthquake epicentres without assuming a model for seismic wave speed variation within the Earth. This allows location improvements achieved using 1-D or 3-D models to be put into perspective. A simple, arrival order misfit criterion is proposed that may be used in standard location schemes. The arrival order misfit criterion does not use a seismic velocity model but simply assumes that the traveltime curve for a particular phase is monotonic with distance. Greater robustness is achieved by including a contribution from every possible pairing of stations and the effect of timing inconsistencies reduced by smoothing. An expression is found that relates the smoothing parameter to the number of observations. A typical event is studied in detail to demonstrate the properties of the misfit function. A pathological case is shown that illustrates that, like other location methods, the arrival order misfit is susceptible to poor station distribution. 25 ground truth and 5000 other teleseismically observed events are relocated and the arrival order locations compared to those found using a least-squares approach and a 1-D earth model. The arrival order misfit is found to be surprisingly accurate when more than 50 observations are used and may be useful in obtaining a model independent epicentre estimate in regions of poorly known velocity structure or the starting point for another location scheme. [source] Sequential integrated inversion of refraction and wide-angle reflection traveltimes and gravity data for two-dimensional velocity structuresGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2000Rosaria Tondi A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the ,maximum likelihood' regularization scheme. We construct a ,probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity,density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion. We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently. An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area. [source] Improved microseismic event location by inclusion of a priori dip particle motion: a case study from EkofiskGEOPHYSICAL PROSPECTING, Issue 5 2010G.A. Jones ABSTRACT Microseismic monitoring in petroleum settings provides insights into induced and naturally occurring stress changes. Such data are commonly acquired using an array of sensors in a borehole, providing measures of arrival times and polarizations. Events are located using 1D velocity models, P- and S-wave arrival times and the azimuths of P-wave particle motions. However in the case of all the sensors being deployed in a vertical or near-vertical borehole, such analysis leads to an inherent 180° ambiguity in the source location. Here we present a location procedure that removes this ambiguity by using the dip of the particle motion as an a priori information to constrain the initial source location. The new procedure is demonstrated with a dataset acquired during hydraulic fracture stimulation, where we know which side of the monitoring well the events are located. Using a 5 -step location procedure, we then reinvestigate a microseismic data set acquired in April 1997 at the Ekofisk oilfield in the North Sea. Traveltimes for 2683 candidate events are manually picked. A noise-weighted analytic-signal polarization analysis is used to estimate the dip and azimuth of P-wave particle motions. A modified t-test is used to statistically assess the reliability of event location. As a result, 1462 events are located but 627 are deemed to be statistically reliable. The application of a hierarchal cluster analysis highlights coherent structures that cluster around wells and inferred faults. Most events cluster at a depth of roughly 3km in the Ekofisk chalk formation but very little seismicity is observed from the underlying Tor chalk formation, which is separated from the Ekofisk formation by an impermeable layer. We see no evidence for seismicity in the overburden but such events may be too distant to detect. The resulting picture of microseismicity at Ekofisk is very different from those presented in previous studies. [source] Constrained tomography of realistic velocity models in microseismic monitoring using calibration shotsGEOPHYSICAL PROSPECTING, Issue 5 2010T. Bardainne ABSTRACT The knowledge of the velocity model in microseismic jobs is critical to achieving statistically reliable microseismic event locations. The design of microseismic networks and the limited sources for calibration do not allow for a full tomographic inversion. We propose optimizing a priori velocity models using a few active shots and a non-linear inversion, suitable to poorly constrained systems. The considered models can be described by several layers with different P- and S-wave velocities. The velocities may be constant or have 3D gradients; the layer interfaces may be simple dipping planes or more complex 3D surfaces. In this process the P- and S- wave arrival times and polarizations measured on the seismograms constitute the observed data set. They are used to estimate two misfit functions: i) one based on the measurement residuals and ii) one based on the inaccuracy of the source relocation. These two functions are minimized thanks to a simulated annealing scheme, which decreases the risk of converging to a local solution within the velocity model. The case study used to illustrate this methodology highlights the ability of this technique to constrain a velocity model with dipping layers. This was performed by jointly using sixteen perforation shots recorded during a multi-stage fracturing operation from a single string of 3C-receivers. This decreased the location inaccuracies and the residuals by a factor of six. In addition, the retrieved layer dip was consistent with the pseudo-horizontal trajectories of the wells and the background information provided by the customer. Finally, the theoretical position of each calibration shot was contained in the uncertainty domain of the relocation of each shot. In contrast, single-stage inversions provided different velocity models that were neither consistent between each other nor with the well trajectories. This example showed that it is essential to perform a multi-stage inversion to derive a better updated velocity model. [source] Migration velocity analysis for tilted transversely isotropic mediaGEOPHYSICAL PROSPECTING, Issue 1 2009Laxmidhar Behera ABSTRACT Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold-and-thrust belts) and in subsalt exploration. Here, we introduce a methodology for P-wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters , and , and linearly varying symmetry-direction velocity VP0 controlled by the vertical (kz) and lateral (kx) gradients. Since determination of tilt from P-wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity VP0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients kz and kx and the anisotropy parameters , and , in the layer-stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry-axis direction is fixed and VP0 is known, the parameters kz, kx, , and , can be resolved from reflection data. It should be emphasized that estimation of , in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas. [source] Velocity analysis based on data correlationGEOPHYSICAL PROSPECTING, Issue 6 2008T. Van Leeuwen ABSTRACT Several methods exist to automatically obtain a velocity model from seismic data via optimization. Migration velocity analysis relies on an imaging condition and seeks the velocity model that optimally focuses the migrated image. This approach has been proven to be very successful. However, most migration methods use simplified physics to make them computationally feasible and herein lies the restriction of migration velocity analysis. Waveform inversion methods use the full wave equation to model the observed data and more complicated physics can be incorporated. Unfortunately, due to the band-limited nature of the data, the resulting inverse problem is highly nonlinear. Simply fitting the data in a least-squares sense by using a gradient-based optimization method is sometimes problematic. In this paper, we propose a novel method that measures the amount of focusing in the data domain rather than the image domain. As a first test of the method, we include some examples for 1D velocity models and the convolutional model. [source] Near-surface models in Saudi ArabiaGEOPHYSICAL PROSPECTING, Issue 6 2007Ralph Bridle ABSTRACT A single-layer model of the near surface throughout the Kingdom of Saudi Arabia is available. While this simple model suffices for most areas and large subsurface structures, it fails in situations where the surface topography is complex, the base of weathering is below the datum, or where the time structural closure is less than the uncertainty in the static correction. In such cases, multiple-layered models that incorporate velocities derived from analysis of first arrivals picked from seismic shot records have proved to be successful in defining the lateral heterogeneity of the near surface. The additional velocity information obtained from this first-arrival analysis (direct as well as refracted arrivals) vastly improves the velocity,depth model of the near surface, regardless of the topography. Static corrections computed from these detailed near-surface velocity models have significantly enhanced subsurface image focusing, thereby reducing the uncertainty in the closure of target structures. Other non-seismic methods have been used either to confirm qualitatively or to enhance the layer models previously mentioned. Gravity data may be particularly useful in sandy areas to confirm general structure, while geostatistical modelling of vibrator base-plate attributes has yielded information that enhances the velocity field. In the global context, exploration targets of the oil and gas industry are seeking smaller and lower relief-time structures. Thus, near-surface models will need to enhance and integrate these methods, particularly in areas where the assumption of flat-lying near-surface layers cannot be met. [source] Artificial neural networks for parameter estimation in geophysicsGEOPHYSICAL PROSPECTING, Issue 1 2000Carlos Calderón-Macías Artificial neural systems have been used in a variety of problems in the fields of science and engineering. Here we describe a study of the applicability of neural networks to solving some geophysical inverse problems. In particular, we study the problem of obtaining formation resistivities and layer thicknesses from vertical electrical sounding (VES) data and that of obtaining 1D velocity models from seismic waveform data. We use a two-layer feedforward neural network (FNN) that is trained to predict earth models from measured data. Part of the interest in using FNNs for geophysical inversion is that they are adaptive systems that perform a non-linear mapping between two sets of data from a given domain. In both of our applications, we train FNNs using synthetic data as input to the networks and a layer parametrization of the models as the network output. The earth models used for network training are drawn from an ensemble of random models within some prespecified parameter limits. For network training we use the back-propagation algorithm and a hybrid back-propagation,simulated-annealing method for the VES and seismic inverse problems, respectively. Other fundamental issues for obtaining accurate model parameter estimates using trained FNNs are the size of the training data, the network configuration, the description of the data and the model parametrization. Our simulations indicate that FNNs, if adequately trained, produce reasonably accurate earth models when observed data are input to the FNNs. [source] Simulation of shockwave propagation with a thermal lattice Boltzmann modelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2003ShiDe Feng Abstract A two-dimensional 19-velocity (D2Q19) lattice Boltzmann model which satisfies the conservation laws governing the macroscopic and microscopic mass, momentum and energy with local equilibrium distribution order O(u4) rather than the usual O(u3) has been developed. This model is applied to simulate the reflection of shockwaves on the surface of a triangular obstacle. Good qualitative agreement between the numerical predictions and experimental measurements is obtained. As the model contains the higher-order terms in the local equilibrium distribution, it performs much better in terms of numerical accuracy and stability than the earlier 13-velocity models with the local equilibrium distribution accurate only up to the second order in the velocity u. Copyright © 2003 John Wiley & Sons, Ltd. [source] Recent modelling of sedimentation of suspended particles: a survey,IRRIGATION AND DRAINAGE, Issue 2 2001P. Boogerd modélisation de la sédimentation; dépostition granulaire; transport sédimentaire Abstract Recent literature on modelling of sedimentation was studied. Attention was paid to hydrodynamics, numerical simulation, settling velocity models, sediment and velocity distribution functions, and sediment transport equations. Many popular theories, e.g. those regarding stratification and preferential sweeping, are under discussion. The traditional view that large-scale, energy-containing fluid motions dominate the transport of particles is found to be under attack, as is the modification of the von Karman coefficient to account for the presence of sediment. It is unclear which model for hindered settling should be used under what circumstances, and the effect of particle distribution cannot yet be calculated. Even the most basic problems, such as settling of multiple and/or non-spherical particles in a quiescent liquid, still require research. In the field of sediment distribution functions new solutions are still not entirely satisfactory. Furthermore, the predictive value of transport rate models is still rather low, and several popular sediment transport functions consistently allow more degradation than aggradation. Copyright © 2001 John Wiley & Sons, Ltd. Ce document est une étude de la littérature récente de la modélisation de la sédimentation, en prêtant attention à la hydrodynamique, les simulations numériques, les modèles de la vitesse de déposer, les fonctions de la répartition du sédiment et de la vitesse, et les équations du transport sédimentaire. Beaucoup de théories modernes, celles concernant la stratification et l'entraînage préférentiel par exemple, sont en cours de discussion. L'idée traditionnelle que les mouvements fluides au champ extensif et contenant de l'energie dominent le transport granulaire est attaquée, ainsi que la modification du coefficient von Karman pour tenir compte de la présence du sédiment. Quel modèle à user pour la déposition gênée, et en quelle situation, n'est pas évident, et l'effet de la répartition granulométrique est incalculable. Même les problèmes les plus fondamentaux, comme la déposition d'e multiples granules ou des granules non-sphériques dans un liquide quiescent, ont besoin de recherche. Dans le domaine des fonctions de la répartition du sédiment, les nouvelles solutions ne sont pas encore entièrement satisfaisantes. En outre la valeur prédictionnaire des modèles de la vitesse de transport est encore assez basse, et plusieurs fonctions de transport populaires permettent constamment plus de dégradation que d'accroissement. Copyright © 2001 John Wiley & Sons, Ltd. [source] A coupled dispersion and exchange model for short-range dry deposition of atmospheric ammoniaTHE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 618 2006Benjamin Loubet Abstract The MODDAS-2D model (MOdel of Dispersion and Deposition of Ammonia over the Short-range in two dimensions) is presented. This stationary model couples a two-dimensional Lagrangian stochastic model for short-range dispersion, with a leaf-scale bi-directional exchange model for ammonia (NH3), which includes cuticular uptake and a stomatal compensation point. The coupling is obtained by splitting the upward and downward components of the flux, which can be generalized for any trace gas, and hence provides a way of simply incorporating bi-directional exchanges in existing deposition velocity models. The leaf boundary-layer resistance is parametrized to account for mixed convection in the canopy, and the model incorporates a stability correction for the Lagrangian time-scale for vertical velocity, which tends to increase the Lagrangian time-scale in very stable conditions compared with usual parametrizations. The model is validated against three datasets, where concentrations of atmospheric NH3 were measured at several distances from a line source. Two datasets are over grassland and one is over maize, giving a range of canopy structure. The model correctly simulates the concentration in one situation, but consistently overestimates it at further distances or underestimates it at small distances in the two other situations. It is argued that these discrepancies are mainly due to the lack of length of one of the line sources and non-aligned winds. Analysis shows that the surface exchange parameters and the turbulent mixing at the source level are the predominant factors controlling short-range deposition of NH3. Copyright © 2006 Royal Meteorological Society [source] | |||||||||||||