Bulk characteristics of heat transfer in the unstable,baroclinic atmospheric boundary layer

Field data for the unstable, baroclinic, atmospheric boundary layer over land and over the sea are considered in the context of a general similarity theory of vertical heat transfer. The dependence of δθ/θ_* upon logarithmic functions of h _c z _T and stability (through the similarity function C) is clearly demonstrated in the data. The combined data support the conventional formulation for the heat transfer coefficient δθ/θ_* when,

1. the surface scaling length is z _T (« z ₀), the height at which the surface temperature over land is obtained by extrapolation of the temperature profile
2. the height scale is taken as the depth of convective mixing h _c
3. the temperature profile equivalent of the von Karman constant is taken as 0.41
4. areal average, rather than single point, values of δθ are employed in strongly baroclinic conditions. No significant effect of baroclinity or the height scale ratio as proposed in the general theory is found. Variations in C about a linear regression relation against stability are most probably due to uncertainties in the areal surface temperature and to experimental errors in general temperature measurements.

     

Atmospheric turbulence within and above a Douglas-fir stand. Part I: Statistical properties of the velocity field

This is the first of two papers reporting the results of a study of the turbulence regimes and exchange processes within and above an extensive Douglas-fir stand. The experiment was conducted on Vancouver Island during a two-week rainless period in July and August 1990. The experimental site was located on a 5^o slope. The stand, which was planted in 1962, and thinned and pruned uniformly in 1988, had a (projected) leaf area index of 5.4 and a heighth=16.7 m. Two eddy correlation units were operated in the daytime to measure the fluctuations in the three velocity components, air temperature and water vapour density, with one mounted permanently at a height of 23.0m (z/h=1.38) and the other at various heights in the stand with two to three 8-hour periods of measurement at each level. Humidity and radiation regimes both above and beneath the overstory and profiles of wind speed and air temperature were also measured. The most important findings are:

1. A marked secondary maximum in the wind speed profile occurred in the middle of the trunk space (aroundz/h=0.12). The turbulence intensities for the longitudinal and lateral velocity components increased with decreasing height, but the intensity for the vertical velocity component had a maximum atz/h=0.60 (middle of the canopy layer). Magnitudes of the higher order moments (skewness and kurtosis) for the three velocity components were higher in the canopy layer than in the trunk space and above the stand.
2. There was a 20% reduction in Reynolds stress fromz/h=1.00 to 1.38. Negative Reynolds stress or upward momentum flux perisistently occurred atz/h=0.12 and 0.42 (base of the canopy), and was correlated with negative wind speed gradients at the two heights. The longitudinal pressure gradient due to the land-sea/upslope-downslope circulations was believed to be the main factor responsible for the negative Reynolds stress.
3. Momentum transfer was highly intermittent. Sweep and ejection events dominated the transfer atz/h=0.60, 1.00 and 1.38, with sweeps playing the more important role of the two atz/h=0.60 and 1.00 and the less important role atz/h=1.38. But interaction events were of greater magnitude than sweep and ejection events atz/h=0.12 and 0.42.

     

Experiments on scalar dispersion within a model plant canopy part II: An elevated plane source

This paper describes a wind-tunnel experiment on the dispersion of trace heat from an effectively planar source within a model plant canopy, the source height being h _s = 0.80 h _c , where h _c is the canopy height. A sensor assembly consisting of three coplanar hot wires and one cold wire was used to make simultaneous measurements of the temperature and the streamwise and vertical velocity components. It was found that:

1. The thermal layer consisted of two parts with different length scales, an inner sublayer (scaling with h _s and h _c ) which quickly reached streamwise equilibrium downstream of the leading edge of the source, and an outer sublayer which was self-preserving with a length scale proportional to the depth of the thermal layer.
2. Below 2h _c , the vertical eddy diffusivity for heat from the plane source (K _HP ) was substantially less than the far-field limit of the corresponding diffusivity for heat from a lateral line source at the same height as the plane source. This shows that dispersion from plane or other distributed sources in canopies is influenced, near the canopy, by turbulence ‘memory’ and must be considered as a superposition of both near-field and far-field processes. Hence, one-dimensional models for scalar transport from distributed sources in canopies are wrong in principle, irrespective of the order of closure.
3. In the budgets for temperature variance, and for the vertical and streamwise components of the turbulent heat flux, turbulent transport was a major loss between h _s and h _c and a principal gain mechanism below h _s , as also observed in the budgets for turbulent energy and shear stress.
4. Quadrant analysis of the vertical heat flux showed that sweeps and ejections contributed about equal amounts to the heat flux between h _s and h _c , though among the more intense events, sweeps were dominant. Below h _s , almost all the heat was transported by sweeps.

     

A comparison of ABL heights inferred routinely from lidar and radiosondes at noontime

The height of the atmospheric boundary layer (ABL) obtained with lidar and radiosondes is compared for a data set of 43 noon (12.00 GMT) cases in 1984. The data were selected to represent the synoptic circulation types appropriately. Lidar vertical profiles at 1064 nm were used to obtain three estimates for the ABL height (h _lid), based on the first gradient in the back-scatter profile, namely, at the beginning, middle and top of the gradient. The boundary-layer height obtained with the radiosondes (h _s) was determined with the dry-parcel-intersection method in unstable conditions. As a first guess for near-neutral and stable conditions, the height of the first significant level in the potential temperature profile was taken. Overall, the boundary-layer thickness estimates agree surprisingly well (regression lineh _lidb=h_s:cc.=0.93 and the standard error=121 m). However, in 10% of the cases, the lidar estimate was significantly lower (difference>400 m) than the routinely inferredh _s. These outliers are discussed separately. For stable conditions, an estimate of ABL height (h _N) is also made based on the friction velocity and the Brunt-Väisälä frequency. The agreement betweenh _Nandh _lidbis good. Discrepancies between the two methods are caused by:

1. rapid growth of the boundary layer arround the measurement time;
2. the presence of a deep entrainment layer leading to a large zone in which quantities are not well mixed;
3. a large systematic error of 100–200 m in the estimate of boundary-layer height obtained from the radiosonde due to the way that profiles are recorded, as a series of significant points.

     

Atmospheric boundary layer research at Cabauw

At Cabauw, The Netherlands, a 213 m high mast specifically built for meteorological research has been operational since 1973. Its site, construction, instrumentation and observation programs are reviewed. Regarding analysis of the boundary layer at Cabauw, the following subjects are discussed:

- terrain roughness;

- Monin-Obukhov theory in practice;

- the structure of stable boundary layers;

- observed evolution of fog layers;

- inversion rise and early morning entrainment;

- use of the geostrophic wind as a predictor for wind profiles;

- height variation of wind climate statistics;

- air pollution applications: long range transport and short range dispersion;

- dependence of sound wave propagation on boundary-layer structure;

- testing of weather and climate models.

     

Wind profiles at the Boulder Tower

Analysis of wind profiles at the Boulder Tower (BAO) leads to these conclusions:

1. The variation of roughness with wind direction found earlier is confirmed. Roughness lengths measured on the tower are larger than those measured close to the surface.
2. The profiles and measurements of Reynolds stress are consistent with a von-Karman constant of 0.35.
3. The form φ_m=(1?15z/L)^-1/3 fits best in the range -0.6 < z/L < 0. In the range 0 < z/L < 0.5, θ _m ～ 1 + 4.7z/L provides a good fit to the observations. For z/L < 0.1, φ _m also depends on h, the thickness of the PBL. For z/L < -0.6, Φ _m approaches the constant 0.5, in contrast to all previous suggestions. For larger stabilities, the upper level is usually not in the surface layer, and wind ratios become independent of z/L.
4. With snow cover, the effective roughness diminishes to about 1 cm, even for directions for which the roughness length without snow is large.
5. Estimation of winds at 100 or 150 m from information near the surface is best for similarity theory provided that the ratio of height to Monin-Obukhov L is less than 0.1. For larger z/L, simple power laws seem more appropriate.

     

Experiments on scalar dispersion within a model plant canopy part I: The turbulence structure

This is the first of a series of three papers describing experiments on the dispersion of trace heat from elevated line and plane sources within a model plant canopy in a wind tunnel. Here we consider the wind field and turbulence structure. The model canopy consisted of bluff elements 60 mm high and 10 mm wide in a diamond array with frontal area index 0.23; streamwise and vertical velocity components were measured with a special three-hot-wire anemometer designed for optimum performance in flows of high turbulence intensity. We found that:

1. The momentum flux due to spatial correlations between time-averaged streamwise and vertical velocity components (the dispersive flux) was negligible, at heights near and above the top of the canopy.
2. In the turbulent energy budget, turbulent transport was a major loss (of about one-third of local production) near the top of the canopy, and was the principal gain mechanism lower down. Wake production was greater than shear production throughout the canopy. Pressure transport just above the canopy, inferred by difference, appeared to be a gain in approximate balance with the turbulent transport loss.
3. In the shear stress budget, wake production was negligible. The role of turbulent transport was equivalent to that in the turbulent energy budget, though smaller.
4. Velocity spectra above and within the canopy showed the dominance of large eddies occupying much of the boundary layer and moving downstream with a height-independent convection velocity. Within the canopy, much of the vertical but relatively little of the streamwise variance occurred at frequencies characteristic of wake turbulence.
5. Quadrant analysis of the shear stress showed only a slight excess of sweeps over ejections near the top of the canopy, in contrast with previous studies. This is a result of improved measurement techniques; it suggests some reappraisal of inferences previously drawn from quadrant analysis.

     

Lake Kinneret (Sea of Galilee) and its exceptional wind system

From 1973–1976, research was performed around the Sea of Galilee, aimed at examining the wind regime in the area and whether the area develops a land-sea breeze despite its particular topographical location.

The main conclusions were:
1. During the summer mornings a lake breeze develops, blowing towards the shores of the lake. It ceases at the peak of its development when a westerly wind, originating in the development of a breeze along the Israeli Mediterranean coast, plunges towards the lake.
2. Late at night, a wind flow develops from the land towards the lake, which combines with the katabatic winds that blow along the steep slopes surrounding the Kinneret.
3. The stations at the upper level, at a height of 400–500 m above the Kinneret, are not affected by the lake breeze during the day or by the land breeze at night.
4. In winter, the Kinneret lake breeze is almost as developed as in summer, because the westerly winds, originating in the Mediterranean sea breeze which hardly develops in this season, do not plunge into the Kinneret.

     

Measurements of internal waves and turbulence in two-dimensional stratified shear flows

A turbulent stratified shear flow is generated in a towing tank by towing a grid or a circular cylinder through a tank of stratified salt water. The internal waves and turbulence generated in these flows are visualized with shadowgraphs and measured with quartz-coated hot-film probes (up to four probes for velocity fluctuations) and single-electrode conductivity probes (up to four probes for salinity fluctuations) which are towed at the same speed as the obstacle. The velocity and salinity signals are recorded on magnetic tapes. A portion of these signals is processed directly-on-line with a digital computer. From these shadowgraphs and probe measurements, we observe that

1. Far downstream of the obstacle where the turbulence has already subsided, the stratified fluid always has a layered structure. This layered structure persists for a long time, and is a result of the convection of turbulently mixed layers by the mean flow. These results indicate that in the regions of a stably stratified atmosphere and ocean where the turbulence has subsided, one could often find layered structure.
2. There are spectral peaks and valleys in the measured velocity and salinity autospectra when the stratifications are sufficiently strong. Under certain conditions, these spectral peaks tend to lift up the spectral curves to show substantialf ^?5/3 subranges, although the turbulence Reynolds numbers are too low for the flows to have recognizable inertial subranges. This anomalousf ^?5/3 subrange demonstrates the pitfalls of using spectral measurements in thef ^?5/3 subrange to predict the turbulent energy dissipation rate through the Kolmogorov hypothesis.
3. A diagnostic method is developed for distinguishing internal waves from turbulence, utilizing their phase characteristics. The phase characteristics can be conveniently examined from the cospectra and quadrature spectra measurements of: (a), two vertically separated velocity probes; (b), two vertically separated density probes; and (c), a velocity probe and a density probe. This method is demonstrated to be useful in the laboratory and can be applied directly to atmospheric and oceanic measurements to distinguish internal waves from turbulence.
4. From the coherency measurements, it is found that the entire turbulent stratified wake is actually whipping up and down at a frequency corresponding to the Brunt-Väisälä frequency. This indicates that similar stratified shear flows in the atmosphere and in the ocean, such as the jet streams in the atmosphere and the Cromwell current in the ocean, may oscillate vertically, which in turn can induce horizontal oscillation and meandering.

     

Experiments on scalar dispersion within a model plant canopy,part III: An elevated line source

An experiment is reported in which heat was released as a passive tracer from an elevated lateral line source within a model plant canopy, with h _s = 0.85 h _c (h _s and h _c being the source and canopy heights, respectively). A sensor assembly consisting of three coplanar hot wires and one cold wire was used to measure profiles of mean temperature % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikamaana% aabaGaeqiUdehaaiaacMcaaaa!390C!\[(\overline \theta )\], temperature variance (Σ_θ ²), vertical and streamwise turbulent heat fluxes, and third moments of wind and temperature fluctuations. Conclusions were:

1. Despite the very heterogeneous flow within the canopy, the observed dispersive heat flux (due to spatial correlation between time-averaged temperature and vertical velocity) was small. However, there is evidence from the plume centroid (which was lower than h _s at the source) of systematic recirculating motions within the canopy.
2. The ratio % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS% baaSqaaiabeI7aXjaab2gacaqGHbGaaeiEaaqabaGccaGGVaWaa0aa% aeaacqaH4oqCaaWaaSbaaSqaaiaab2gacaqGHbGaaeiEaaqabaaaaa!41DF!\[\sigma _{\theta {\text{max}}} /\overline \theta _{{\text{max}}} \] (of maximum values on vertical profiles) decreased from 1 near the source to an asymptotic value of 0.4 far downstream, in good agreement with previous experimental and theoretical work for concentration fluctuations in the surface layer well above the canopy.
3. The eddy diffusivity for heat from the line source (K _HL ) increased, downstream of the source, to a nearly constant ‘far-field’ vertical profile. Within the canopy, the far-field K _HL was an order of magnitude larger than K _HP , the equivalent diffusivity for a plane source; well above the canopy, the two were equal. The time scale defined by (far-field K _HL )/(vertical velocity variance) was independent of height within the canopy.
4. Budgets for temperature variance, vertical heat flux and streamwise heat flux are remarkably similar to the equivalent budgets for an elevated line source in the surface layer well above the canopy, except in the lower part of the canopy in the far field, where vertical transport is much more important than in the surface layer.
5. A random flight simulation of the mean height and depth of the temperature plume was generally in good agreement with experiment. However, details of the temperature and streamwise turbulent heat flux profiles were not correct, suggesting that the model formulation needs to be improved.

     

Vertical cross-spectra of horizontal velocity components at the Boulder observatory

Cross-spectra between horizontal wind components at different levels of the Boulder Atmospheric Observatory (BAO) tower lead to the following conclusions:

1. Davenport's hypothesis is satisfied that coherence decays exponentially with the ratio of vertical separation to horizontal wave length, at least to very small values of coherence.
2. The decay coefficients increase with z/L for z/L < 0.5. For larger stabilities, irregular fluctuations with periods of order 10–20 min have considerable vertical coherence. Results at BAO are quite consistent with those elsewhere.
3. Eddy slopes in vertical planes increase with wind shear up to a point where the slope (horizontal delay over vertical separation) is just above 2. Beyond that point, the systematic increase of slopes with shear ceases. Since wind shear decreases upward, slopes tend to decrease upward. Slopes for lateral components are significantly larger than those for u-components.

     

Inversions,and fog,stratus and cumulus formation in warm air over cooler water

Two aspects of convection over oceans are discussed and the following conclusions are derived from theoretical considerations.

1. The air layer over the sea will usually convect even when the water surface is ten degrees or more colder than the initial air temperature.
2. An inversion at stratus cloud tops is created by the stratus, and is not a necessary preexisting condition. Such inversions persist after subsidence evaporates the cloud.
3. Radiation heat exchange does not play an essential role in stratus formation or maintenance, and can either heat or cool the cloud.
4. Dry air convection does not erode inversions at the top of the convecting layer. Examples of soundings are discussed.
5. Fogs are most likely to form at sea where the water is coolest, and need no radiation effects to initiate cooling, or a boost from patches of warmer water, to begin convection.
6. Both stratus cloud growth, and the evaporation of clouds by cloud top entrainment, readjust the vertical structure of the air to leave a constant wet-bulb potential temperature with height.

These conclusions are supported by, firstly, a convective model which has been developed and which shows that vapor-driven convection over the ocean will proceed with zero or negative heat fluxes, at rates which saturate the lowest layer of the atmosphere in a few hours to altitudes of many tens of meters. Secondly, the availability of condensed moisture at the top of the surface layer cools the warmer entrained overlying dry air parcels so that when they descend they are no warmer than the sea surface temperature, and this induces downward moving plumes. This occurs if the wet-bulb potential temperature of the overlying air is less than the sea surface temperature, even if it is ten degrees C, or more, warmer in actual temperature.     

Numerical study on the bulk heat transfer coefficient for a variety of vegetation types and densities

The relationship between the geometrical structure of a canopy layer and the bulk transfer coefficient was investigated using a numerical canopy model. The following results were obtained:

1. The bulk transfer coefficients for momentum and heat, C _M and C _H , change with non-dimensional canopy density C _* each has a maximum.
2. The value of C _M is always larger than the value of C _H for a canopy with c _m > c _h , c _m and c _h being the drag coefficient and the heat transfer coefficient of an individual canopy element, respectively.
3. The value of C _* at which C _H has its maximum value is larger than the value of C _* at which C _M has its maximum. Therefore, the reciprocal of the sublayer Stanton number b _h ^?1 ranges between 50 and 65 for C _* around 0.1 while it ranges between 0 and 30 for C _* < 10^?2 and C _* > 2 (when c _m = 0.5).
4. The value of B _H ^?1 in the present study is consistent with most available observations, except for canopies of medium density (when C _* is around 0.1) for which no observational value has been obtained.

     

A wind tunnel study of turbulent flow around single and multiple windbreaks, part I: Velocity fields

This paper describes wind-tunnel experiments on the flow around single and multiple porous windbreaks (height H), sheltering a model plant canopy (height H/3). The mean wind is normal to the windbreaks, which span the width of the wind tunnel. The incident turbulent flow simulates the adiabatic atmospheric surface layer. Five configurations are examined: single breaks of three solidities (low, medium, high; solidity = 1 - porosity), and medium-solidity multiple breaks of streamwise spacing 12H and 6H. The experimental emphases are on the interactions of the windbreak flow with the underlying plant canopy; the effects of solidity; the differences in shelter between single and multiple windbreaks; and the scaling properties of the flow. Principal results are: (1) the "quiet zones" behind each windbreak are smaller in multiple than single arrays, because of the higher turbulence level in the very rough-wall internal boundary layer which develops over the multiple arrays. Nevertheless, the overall shelter effectiveness is higher for multiple arrays than single windbreaks because of the "nonlocal shelter" induced by the array as a whole. (2) The flow approaching the windbreak decelerates above the canopy but accelerates within the canopy, particularly when the windbreak solidity is high. (3) A strong mixing layer forms just downwind of the top of each windbreak, showing some of the turbulence and scaling properties of the classical mixing layer formed between uniform, coflowing streams. (4) No dramatic increase in turbulence levels in the canopy is evident at the point where the deepening mixing layer contacts the canopy (around x/H = 3) but the characteristic inflection in the canopy wind profile is eliminated at this point.     

Turbulence Anisotropy in the Near-Surface Atmosphere and the Evaluation of Multiple Outer Length Scales

Length scales determined by maximum turbulent kinetic energy (TKE), the integral scale, and two length scales based on Reynolds stress-tensor anisotropy are compared to the often stated outer length scales of boundary-layer depth and distance from the earth’s surface, $z$ . The scales are calculated using sonic anemometer data from two elevations, 5 and 50 m above the ground at the main tower site of the CASES-99 field campaign. In general, none of these scales agrees with the other, although the scale of maximum TKE is often similar to the boundary-layer depth during daytime hours, and the length scales derived from anisotropy characteristics are sometimes similar to $\kappa \!z, z$ , and $2z$ depending on scale definition and thermal stability. Except for the scale with the strictest isotropy threshold, the turbulence is anisotropic for each of the various candidates for the outer scale. Length scales for maximum buoyancy flux and temperature variance are evaluated and the turbulence characteristics at these scales are almost always found to be anisotropic.     

The internal boundary layer — A review

A review is given of relevant work on the internal boundary layer (IBL) associated with:

1. Small-scale flow in neutral conditions across an abrupt change in surface roughness,
2. Small-scale flow in non-neutral conditions across an abrupt change in surface roughness, temperature or heat/moisture flux,
3. Mesoscale flow, with emphasis on flow across the coastline for both convective and stably stratified conditions.

The major theme in all cases is on the downstream, modified profile form (wind and temperature), and on the growth relations for IBL depth.     

Wind characteristics at the Boulder Atmospheric Observatory

Wind speeds at the 300 m tower at the Boulder Atmospheric Observatory have been analyzed. This tower is located in slightly rolling farmland. The following conclusions have been reached:

1. For west winds, the terrain is sufficiently uniform for simple surface-layer theory to be adequate without modification even though the air has moved up a small slope to reach the tower. For south and southeast winds, ‘effective’ roughness lengths must be introduced, which are significantly larger than the ‘true’ roughness length.
2. Useful wind estimates up to 150 m can be made from winds at 10 m and stability information, provided the ‘effective’ roughness length is known.
3. The observations are consistent with a von Kármán constant of 0.35.

     

Aircraft measurements of oxides of nitrogen along the eastern rim of the Asian continent: Winter observations

Aircraft observations of oxides of nitrogen (NO _y ), measured with a ferrous sulfate converter, over the sea surrounding the Japanese islands (30–43° N, 131–141° E) were carried out in the winter of 1983 and 1984 at altitudes mostly between 3 and 8 km. NO _y defined here is the sum of NO, NO₂, and other unstable oxides of nitrogen that are converted to NO by ferrous sulfate. The main observations were:

1. Over the Pacific Ocean between the latitudes of 30–35° N, the observed NO _y mixing ratio between 3 and 8 km was a fairly constant 200 pptv. The NO mixing ratio increased with altitude from 15 pptv at 3 km to 35 pptv at 7 km.
2. Over the Sea of Japan, tropospheric NO _y mesured between 1 and 6 km started increasing with latitude North of 35° N and reached about 1000 pptv at 40° N.
3. NO _y was measured in an air mass transported from the stratosphere near a tropopause fold region. When the ozone mixing ratio was between 80 and 140 ppbv, the NO _y mixing ratio was about 200 pptv.

     

Limitations of an eddy-correlation technique for the determination of the carbon dioxide and sensible heat fluxes

A modified infrared CO₂ gas analyzer, a small thermocouple assembly, a heated-thermocouple anemometer for horizontal wind, and a propeller-type vertical wind sensor were used to measure the eddy fluxes of heat and CO₂ above a corn crop. Experimental results of these fluxes are discussed. The main sources of errors of the eddy fluxes using these instruments were estimated:

1. Sensors with a time constant of 0.5 s appear to be fast enough to detect most of the vertical CO₂ transfer as long as the sensors are located at least one meter above the crop surface.
2. The deviation from steady-state conditions for 10-min periods was found to have a significant effect on the eddy flux estimates.
3. Temperature fluctuations of the air sample passing through the CO₂ infrared gas analyzer were found to be non-negligible but could be easily corrected.
4. A 1° misalignment of the vertical anemometer affected these eddy fluxes by less than 10% under all circumstances studied.