首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The current outbreak of mountain pine beetle (MPB) that started in the late 1990s in British Columbia, Canada, is the largest ever recorded in the north American native habitat of the beetle. The killing of trees is expected to change the vertical distribution of net radiation ( $Q^*$ Q ? ) and the partitioning of latent ( $Q_\mathrm{E}$ Q E ) and sensible ( $Q_\mathrm{H}$ Q H ) heat fluxes in the different layers of an attacked forest canopy. During an intensive observation period in the summer of 2010, eddy-covariance flux and radiation measurements were made at seven heights from ground level up to 1.34 times the canopy height in an MPB-attacked open-canopy forest stand $(\hbox {leaf area index} = 0.55~\mathrm{{m}}^{2}\ \mathrm{{m}}^{-2})$ ( leaf area index = 0.55 m 2 m - 2 ) in the interior of British Columbia, Canada. The lodgepole pine dominated stand with a rich secondary structure (trees and understorey not killed by the beetle) was first attacked by the MPB in 2003 and received no management. In this study, the vertical distribution of the energy balance components and their sources and sinks were analyzed and energy balance closure (EBC) was determined for various levels within the canopy. The low stand density resulted in approximately 60 % of the shortwave irradiance and 50 % of the daily total $Q^*$ Q ? reaching the ground. Flux divergence calculations indicated relatively strong sources of latent heat at the ground and where the secondary structure was located. Only very weak sources of latent heat were found in the upper part of the canopy, which was mainly occupied by dead lodgepole pine trees. $Q_\mathrm{H}$ Q H was the dominant term throughout the canopy, and the Bowen ratio ( $Q_\mathrm{H}/Q_\mathrm{E}$ Q H / Q E ) increased with height in the canopy. Soil heat flux ( $Q_\mathrm{G}$ Q G ) accounted for approximately 4 % of $Q^*$ Q ? . Sensible heat storage in the air ( $\Delta Q_\mathrm{S,H}$ Δ Q S , H ) was the largest of the energy balance storage components in the upper canopy during daytime, while in the lower canopy sensible heat storage in the boles ( $\Delta Q_\mathrm{S,B}$ Δ Q S , B ) and biochemical energy storage ( $\Delta Q_\mathrm{S,C}$ Δ Q S , C ) were the largest terms. $\Delta Q_\mathrm{S,H}$ Δ Q S , H was almost constant from the bottom to above the canopy. $\Delta Q_\mathrm{S,C}$ Δ Q S , C , $\Delta Q_\mathrm{S,B}$ Δ Q S , B and latent heat storage in the air ( $\Delta Q_\mathrm{S,E}$ Δ Q S , E ) varied more than $\Delta Q_\mathrm{S,H}$ Δ Q S , H throughout the canopy. During daytime, energy balance closure was high in and above the upper canopy, and in the lowest canopy level. However, where the secondary structure was most abundant, ${\textit{EBC}} \le 66\,\%$ EBC ≤ 66 % . During nighttime, the storage terms together with $Q_\mathrm{G}$ Q G made up the largest part of the energy balance, while $Q_\mathrm{H}$ Q H and $Q_\mathrm{E}$ Q E were relatively small. These radiation and energy balance measurements in an insect-attacked forest highlight the role of secondary structure in the recovery of attacked stands.  相似文献   

2.
The scalar flux–gradient relationships of temperature ( $\phi _{T}$ ? T ) and specific humidity ( $\phi _{q}$ ? q ) under unstable conditions are investigated using eddy-covariance measurements of air–sea turbulent fluxes and vertical profiles of temperature and specific humidity collected from a marine meteorological platform. The gradients of temperature and specific humidity are obtained from measurements at five heights above the sea surface using the log-square fitting method and the simpler first-order approximation method. The two methods yield similar results. The proposed flux–gradient relationships $\phi _{T}$ ? T and $\phi _{q}$ ? q covers a wide range of instability: the stability parameter $\zeta $ ζ ranges from $-$ ? 0.1 to $-$ ? 50. The functional form of the proposed flux–gradient relationships is an interpolation between the Businger–Dyer relation and the free convection relation, which includes the “ $-$ ? 1/2” and “ $-$ ? 1/3” scaling laws at two different stability regimes. The widely used COARE 3.0 algorithm, which is an interpolation between the integrals of the Businger–Dyer and the free convection relations, is also evaluated and compared. The analysis and comparisons show that both schemes generate reasonable values of $\phi _{q}$ ? q in the whole unstable regime. The COARE 3.0 algorithm, however, overestimates $\phi _{T}$ ? T values under very unstable conditions. The errors in the flux–gradient relationships induced by the random errors in the turbulence measurements are assessed. When the random errors are taken into account, the observations agree with predictions of various schemes fairly well, implying that the dominant transport mechanism is adequately captured by the Monin–Obukhov similarity theory. The study also shows that $\phi _{q}$ ? q is significantly ${>}\phi _{T}$ > ? T under unstable conditions and that the ratio $\phi _{q}/\phi _{T}$ ? q / ? T increases with $-\zeta $ ? ζ . The ratio of $\phi _{q}$ ? q to $\phi _{T}$ ? T and the ratio of turbulent transport efficiencies of heat and water vapour ( $R_{wT}/R_{wq}$ R wT / R wq ) suggest that heat is transported more efficiently than water vapour under unstable conditions.  相似文献   

3.
Investigation of the Stable Atmospheric Boundary Layer at Halley Antarctica   总被引：2，自引：2，他引：0  
Boundary-layer measurements from the Brunt Ice Shelf, Antarctica are analyzed to determine flux–profile relationships. Dimensionless quantities are derived in the standard approach from estimates of wind shear, potential temperature gradient, Richardson number, eddy diffusivities for momentum and heat, Prandtl number, mixing length and turbulent kinetic energy. Nieuwstadt local scaling theory for the stable atmospheric boundary-layer appears to work well departing only slightly from expressions found in mid-latitudes. An $E$ E – $l_{\mathrm{m}}$ l m single-column model of the stable boundary layer is implemented based on local scaling arguments. Simulations based on the first GEWEX Atmospheric Boundary-Layer Study case study are validated against ensemble-averaged profiles for various stability classes. A stability-dependent function of the dimensionless turbulent kinetic energy allows a better fit to the ensemble profiles.  相似文献   

4.
In a valley sheltered from strong synoptic effects, the dynamics of the valley atmosphere at night is dominated by katabatic winds. In a stably stratified atmosphere, these winds undergo temporal oscillations, whose frequency is given by $N \sin {\alpha }$ N sin α for an infinitely long slope of constant slope angle $\alpha $ α , $N$ N being the buoyancy frequency. Such an unsteady flow in a stably stratified atmosphere may also generate internal gravity waves (IGWs). The numerical study by Chemel et al. (Meteorol Atmos Phys 203:187–194, 2009) showed that, in the stable atmosphere of a deep valley, the oscillatory motions associated with the IGWs generated by katabatic winds are distinct from those of the katabatic winds. The IGW frequency was found to be independent of $\alpha $ α and about $0.8N$ 0.8 N . Their study did not consider the effects of the background stratification and valley geometry on these results. The present work extends this study by investigating those effects for a wide range of stratifications and slope angles, through numerical simulations for a deep valley. The two oscillatory systems are reproduced in the simulations. The frequency of the oscillations of the katabatic winds is found to be equal to $N$ N times the sine of the maximum slope angle. Remarkably, the IGW frequency is found to also vary as $C_\mathrm{w}N$ C w N , with $C_\mathrm{w}$ C w in the range $0.7$ 0.7 – $0.95$ 0.95 . These values for $C_\mathrm{w}$ C w are similar to those reported for IGWs radiated by any turbulent field with no dominant frequency component. Results suggest that the IGW wavelength is controlled by the valley depth.  相似文献   

5.
In October 2012 Hurricane Sandy devastated New York City and its vicinity caused mainly by the storm surge, which is the water height above normal astronomical tide level. The meteorological conditions were as follows: minimum central pressure, 962 hPa, highest sustained wind speed 27.1 m s $^{-1}$ ? 1 and maximum gust 37.8 m s $^{-1}$ ? 1 . The peak storm surge was at 3.9 m and the peak storm tide at 4.4 m (which is referenced above mean lower low water). The wind-stress tide relation shows that $S=K\,V^{2}$ S = K V 2 , where $S$ S is the storm surge, $V$ V is the wind speed and $K$ K is the coefficient. It is found that with $S$ S in units of m, and $V$ V in  m s $^{-1}$ ? 1 , $K = 0.0051$ K = 0.0051 with $R^{2}= 0.91$ R 2 = 0.91 ( $R$ R is the correlation coefficient) indicating that 91 % of the total variation of the storm surge can be explained by variations in the wind stress, which is proportional to $V^{2}$ V 2 . Similar results were obtained during Hurricane Irene in 2011, which also affected the New York area. Therefore, this simple wind stress-tide relation should be useful in coastal engineering, urban planning, and emergency management.  相似文献   

6.
Methane ( ${\mathrm {CH}}_{4}$ ) fluxes observed with the eddy-covariance technique using an open-path ${\mathrm {CH}}_{4}$ analyzer and a closed-path ${\mathrm {CH}}_{4}$ analyzer in a rice paddy field were evaluated with an emphasis on the flux correction methodology. A comparison of the fluxes obtained by the analyzers revealed that both the open-path and closed-path techniques were reliable, provided that appropriate corrections were applied. For the open-path approach, the influence of fluctuations in air density and the line shape variation in laser absorption spectroscopy (hereafter, spectroscopic effect) was significant, and the relative importance of these corrections would increase when observing small ${\mathrm {CH}}_{4}$ fluxes. A new procedure proposed by Li-Cor Inc. enabled us to accurately adjust for these effects. The high-frequency loss of the open-path ${\mathrm {CH}}_{4}$ analyzer was relatively large (11 % of the uncorrected covariance) at an observation height of 2.5 m above the canopy owing to its longer physical path length, and this correction should be carefully applied before correcting for the influence of fluctuations in air density and the spectroscopic effect. Uncorrected ${\mathrm {CH}}_{4}$ fluxes observed with the closed-path analyzer were substantially underestimated (37 %) due to high-frequency loss because an undersized pump was used in the observation. Both the bandpass and transfer function approaches successfully corrected this flux loss. Careful determination of the bandpass frequency range or the transfer function and the cospectral model is required for the accurate calculation of ${\mathrm {CH}}_{4}$ fluxes with the closed-path technique.  相似文献   

7.
The characteristics of the temporal and height variations of the temperature structure parameter $C_\mathrm{T}^{2}$ in strongly convective situations derived from the sodar echo-signal intensity measurements were analyzed for the first 100 m. It was corroborated that the probability density function (pdf) of the logarithm of $C_\mathrm{T}^{2}$ in the lower convective boundary layer is markedly non-Gaussian, whereas turbulence theory predicts it to be normal. It was also corroborated that the sum of two weighted Gaussians, which characterize the statistics of $C_\mathrm{T}^{2}$ within convective plumes and in their environment and the probability of plume occurrence, well approximates the observed pdfs. It was shown that the height behaviour of the arithmetic mean of $ C_\mathrm{T}^{2}$ (both total and within plumes) follows well a power law $C_\mathrm{T}^{2} (z) \sim z^{-q}$ with the exponent $q$ close to the theoretically predicted value of 4/3. But for the geometrical means of $C_\mathrm{T}^{2}$ (both total and within the plumes), $q$ is close to 1. The difference between arithmetically and geometrically averaged $C_\mathrm{T}^{2}$ profiles was analyzed. The vertical profiles of the standard deviation, skewness and kurtosis of $\hbox {ln}C_\mathrm{T}^{2}$ pdfs were analyzed to show their steady behaviour with height. The standard deviations of the logarithm of $C_\mathrm{T}^{2}$ within the plumes and between them are similar and are 1.5 times less than the total standard deviation. The estimate of the variability index $F_\mathrm{T}$ and its height behaviour were obtained, which can be useful to validate some theoretical and modelling predictions. The vertical profiles of the skewness and kurtosis show the negative asymmetry of pdfs and their flatness, respectively. The spectra of variations in $\hbox {ln}C_\mathrm{T}^{2}$ are shown to be satisfactorily fitted by the power law $f^{-\gamma } $ in the frequency range 0.02 and 0.2 Hz, with the average exponent $\approx $ 1.27  $\pm $  0.22.  相似文献   

8.
The air–sea transfer velocity of $\mathrm{CO}_{2}\, (k_{\mathrm{CO}_{2}})$ was investigated in a shallow estuary in March to July 2012, using eddy-covariance measurements of $\mathrm{CO}_{2}$ fluxes and measured air–sea $\mathrm{CO}_{2}$ partial-pressure differences. A data evaluation method that eliminates data by nine rejection criteria in order to heighten parametrization certainty is proposed. We tested the data evaluation method by comparing two datasets: one derived using quality criteria related solely to the eddy-covariance method, and the other derived using quality criteria based on both eddy-covariance and cospectral peak methods. The best parametrization of transfer velocity normalized to a Schmidt number of 600 $(k_{600})$ was determined to be: $k_{600} = 0.3\,{U_{10}}^{2.5}$ where $U_{10}$ is the wind speed in m $\mathrm{s}^{-1}$ at 10 m; $k_{600}$ is based on $\mathrm{CO}_{2}$ fluxes calculated by the eddy-covariance method and including the cospectral peak method criteria. At low wind speeds, the transfer velocity in the shallow water estuary was lower than in other coastal waters, possibly a symptom of low tidal amplitude leading to low intensity water turbulence. High transfer velocities were recorded above wind speeds of 5 m $\mathrm{s}^{-1}$ , believed to be caused by early-breaking waves and the large fetch (6.5 km) of the estuary. These findings indicate that turbulence in both air and water influences the transfer velocity.  相似文献   

9.
The effects of abrupt streamwise transitions of the aerodynamic roughness length ( $z_\mathrm{o}$ z o ) on the stable atmospheric boundary layer are evaluated using a series of large-eddy simulations based on the first Global Energy and Water Cycle Experiment Atmospheric Boundary Layer intercomparison study (GABLS1). Four $z_\mathrm{o}$ z o values spanning three orders of magnitude are used to create all possible binary distributions with each arranged into patches of characteristic length scales equal to roughly one-half, one, and two times the equivalent homogeneous boundary-layer height. The impact of the heterogeneity on mean profiles of wind speed and temperature, on surface fluxes of heat and momentum, and on internal boundary-layer dynamics are considered. It is found that $z_\mathrm{o}$ z o transitions do not significantly alter the functional relationship between the average surface fluxes and the mean profiles of wind speed and potential temperature. Although this suggests that bulk similarity theory is applicable for modelling the stable boundary layer over $z_\mathrm{o}$ z o heterogeneity, effective surface parameters must still be specified. Existing models that solve for effective roughness lengths of momentum and heat are evaluated and compared to values derived from the simulation data. The existing models are unable to accurately reproduce both the values of the effective aerodynamic roughness lengths and their trends as functions of patch length scale and stability. A new model for the effective aerodynamic roughness length is developed to exploit the benefits of the other models tested. It accurately accounts for the effects of the heterogeneity and stratification on the blending height and effective aerodynamic roughness length. The new model provides improved average surface fluxes when used with bulk similarity.  相似文献   

10.
We examine the influence of a modern multi-megawatt wind turbine on wind and turbulence profiles three rotor diameters ( $D$ D ) downwind of the turbine. Light detection and ranging (lidar) wind-profile observations were collected during summer 2011 in an operating wind farm in central Iowa at 20-m vertical intervals from 40 to 220 m above the surface. After a calibration period during which two lidars were operated next to each other, one lidar was located approximately $2D$ 2 D directly south of a wind turbine; the other lidar was moved approximately $3D$ 3 D north of the same wind turbine. Data from the two lidars during southerly flow conditions enabled the simultaneous capture of inflow and wake conditions. The inflow wind and turbulence profiles exhibit strong variability with atmospheric stability: daytime profiles are well-mixed with little shear and strong turbulence, while nighttime profiles exhibit minimal turbulence and considerable shear across the rotor disk region and above. Consistent with the observations available from other studies and with wind-tunnel and large-eddy simulation studies, measurable reductions in wake wind-speeds occur at heights spanning the wind turbine rotor (43–117 m), and turbulent quantities increase in the wake. In generalizing these results as a function of inflow wind speed, we find the wind-speed deficit in the wake is largest at hub height or just above, and the maximum deficit occurs when wind speeds are below the rated speed for the turbine. Similarly, the maximum enhancement of turbulence kinetic energy and turbulence intensity occurs at hub height, although observations at the top of the rotor disk do not allow assessment of turbulence in that region. The wind shear below turbine hub height (quantified here with the power-law coefficient) is found to be a useful parameter to identify whether a downwind lidar observes turbine wake or free-flow conditions. These field observations provide data for validating turbine-wake models and wind-tunnel observations, and for guiding assessments of the impacts of wakes on surface turbulent fluxes or surface temperatures downwind of turbines.  相似文献   

11.
Scintillometer measurements of the turbulence inner-scale length $l_\mathrm{o }$ l o and refractive index structure function $C_n^2$ C n 2 allow for the retrieval of large-scale area-averaged turbulent fluxes in the atmospheric surface layer. This retrieval involves the solution of the non-linear set of equations defined by the Monin–Obukhov similarity hypothesis. A new method that uses an analytic solution to the set of equations is presented, which leads to a stable and efficient numerical method of computation that has the potential of eliminating computational error. Mathematical expressions are derived that map out the sensitivity of the turbulent flux measurements to uncertainties in source measurements such as $l_\mathrm{o }$ l o . These sensitivity functions differ from results in the previous literature; the reasons for the differences are explored.  相似文献   

12.
In studies of lake–atmosphere interactions, the fluxes of momentum, water vapour and sensible heat are often parametrized as being proportional to the differences in wind, humidity and air temperature between the water surface and a reference height above the surface. Here, the proportionality via transfer coefficients in these relationships was investigated with the eddy-covariance method at three sites within an eddy-covariance mesonet across Lake Taihu, China. The results indicate that the transfer coefficients decreased with increasing wind speed for weak winds and approached constant values for strong winds. The presence of submerged macrophytes reduced the momentum transfer (drag) coefficient significantly. At the two sites free of submerged macrophytes, the 10-m drag coefficients under neutral stability were 1.8 $(\pm \,0.4) \times \,10^{-3}$ ( ± 0.4 ) × 10 ? 3 and $1.7\,(\pm \,0.3) \times \,10^{-3 }$ 1.7 ( ± 0.3 ) × 10 ? 3 at the wind speed of $9\,\text{ m } \text{ s }^{-1}$ 9 m s ? 1 , which are 38 and 34 % greater than the prediction by the Garratt model for the marine environment.  相似文献   

13.
The bulk Richardson number ( $Ri_{Bh}$ ; defined over the entire stable boundary layer) is commonly utilized in observational and modelling studies for the estimation of the boundary-layer height. Traditionally, $Ri_{Bh}$ is assumed to be a quasi-universal constant. Recently, based on large-eddy simulation and wind-tunnel data, a stability-dependent relationship has been proposed for $Ri_{Bh}$ . In this study, we analyze extensive observational data from several field campaigns and provide further support for this newly proposed relationship.  相似文献   

14.
Airflow over two-dimensional hills was investigated in a wind tunnel using particle image velocimetry. We focus on the flow separation behaviour. A trapezoidal hill shape was used in most of the experimental runs, but the critical slope angle for flow separation was approximately the same as that established for smooth hill shapes. The re-attachment point of the separated flow became farther from the hill as the slope angle $\theta $ increased, reaching a saturation of about seven times the hill height for $\theta \gtrsim 60^\circ $ . Increasing the upwind surface roughness length was found to suppress flow separation. This tendency is analogous to the previous experimental results for turbulent boundary layers on flat plates. The boundary-layer thickness varied by the presence or absence of Counihan-type spires and a castellated fence at the test-section entrance had negligible effect on the flow separation.  相似文献   

15.
Vertical mixing of the nocturnal stable boundary layer (SBL) over a complex land surface is investigated for a range of stabilities, using a decoupling index ( $0 < D_{rb} < 1$ ) based on the 2–50 m bulk gradient of the ubiquitous natural trace gas radon-222. The relationship between $D_{rb}$ and the bulk Richardson number ( $R_{ib}$ ) exhibits three broad regions: (1) a well-mixed region ( $D_{rb} \approx 0.05$ ) in weakly stable conditions ( $R_{ib} < 0.03$ ); (2) a steeply increasing region ( $0.05 < D_{rb} < 0.9$ ) for “transitional” stabilities ( $0.03 < R_{ib} < 1$ ); and (3) a decoupled region ( $D_{rb} \approx 0.9$ –1.0) in very stable conditions ( $R_{ib} > 1$ ). $D_{rb}$ exhibits a large variability within individual $R_{ib}$ bins, however, due to a range of competing processes influencing bulk mixing under different conditions. To explore these processes in $R_{ib}$ – $D_{rb}$ space, we perform a bivariate analysis of the bulk thermodynamic gradients, various indicators of external influences, and key turbulence quantities at 10 and 50 m. Strong and consistent patterns are found, and five distinct regions in $R_{ib}$ – $D_{rb}$ space are identified and associated with archetypal stable boundary-layer regimes. Results demonstrate that the introduction of a scalar decoupling index yields valuable information about turbulent mixing in the SBL that cannot be gained directly from a single bulk thermodynamic stability parameter. A significant part of the high variability observed in turbulence statistics during very stable conditions is attributable to changes in the degree of decoupling of the SBL from the residual layer above. When examined in $R_{ib}$ – $D_{rb}$ space, it is seen that very different turbulence regimes can occur for the same value of $R_{ib}$ , depending on the particular combination of values for the bulk temperature gradient and wind shear, together with external factors. Extremely low turbulent variances and fluxes are found at 50 m height when $R_{ib} > 1$ and $D_{rb} \approx 1$ (fully decoupled). These “quiescent” cases tend to occur when geostrophic forcing is very weak and subsidence is present, but are not associated with the largest bulk temperature gradients. Humidity and net radiation data indicate the presence of low cloud, patchy fog or dew, any of which may aid decoupling in these cases by preventing temperature gradients from increasing sufficiently to favour gravity wave activity. The largest temperature gradients in our dataset are actually associated with smaller values of the decoupling index ( $D_{rb} < 0.7$ ), indicating the presence of mixing. Strong evidence is seen from enhanced turbulence levels, fluxes and submeso activity at 50 m, as well as high temperature variances and heat flux intermittencies at 10 m, suggesting this region of the $R_{ib}$ – $D_{rb}$ distribution can be identified as a top-down mixing regime. This may indicate an important role for gravity waves and other wave-like phenomena in providing the energy required for sporadic mixing at this complex terrain site.  相似文献   

16.
Large-eddy simulations of the atmospheric boundary layer (ABL) under a wide range of stabilities are conducted to educe very-large-scale motions and then to study their dynamics and how they are influenced by buoyancy. Preliminary flow visualizations suggest that smaller-scale motions that resemble hairpins are embedded in much larger scale streamwise meandering rolls. Using simulations that represent more than 150 h of physical time, many snapshots in the \(xy\) -, \(yz\) - and \(xz\) -planes are then collected to perform snapshot proper orthogonal decomposition and further investigate the large structures. These analyses confirm that large streamwise rolls that share several features with the very-large-scale motions observed in laboratory studies arise as the dominant modes under most stabilities, but the effect of the surface kinematic buoyancy flux on the energy content of these dominant modes is very significant. The first two modes in the \(yz\) -plane in the neutral case contain up to 3 % of the total turbulent kinetic energy; they also have a vertical tilt angle in the \(yz\) -plane of about 0 to 30 \(^\circ \) due to the turning effect associated with the Coriolis force. Unstable cases also feature streamwise rolls, but in the convective ABL they are strengthened by rising plumes in between them, with two to four rolls spanning the whole domain in the first few modes; the Coriolis effect is much weaker in the unstable ABL. These rolls are no longer the dominant modes under stable conditions where the first mode is observed to contain sheet-like motions with high turbulent kinetic energy. Using these proper orthogonal decomposition modes, we are also able to extract the vertical velocity fields corresponding to individual modes and then to correlate them with the horizontal velocity or temperature fields to obtain the momentum and heat flux carried by individual modes. Structurally, the fluxes are explained by the topology of their corresponding modes. However, the fraction of the fluxes produced by the modes is invariably smaller than the fraction of energy they contain, particularly under stable conditions where the first modes are found to perform weak counter-gradient fluxes.  相似文献   

17.
Analytical Lagrangian equations capable of predicting concentration profiles from known source distributions offer the opportunity to calculate source/sink distributions through inverted forms of these equations. Inverse analytical Lagrangian equations provide a practical means of estimating source profiles using concentration and turbulence measurements. Uncertainty concerning estimates of the essentially immeasurable Lagrangian length scale ( ${\mathcal{L}}$ ), a key input, impedes the operational practicality of this method. The present study evaluates ${\mathcal{L}}$ within a corn canopy by using field measurements to constrain an analytical Lagrangian equation. Measurements of net CO2 flux, soil-to-atmosphere CO2 flux, and in-canopy profiles of CO2 concentration provided the information required to solve for ${\mathcal{L}}$ in a global optimization algorithm for 30-min time intervals. For days when the canopy was a strong CO2 sink, the optimization frequently located ${\mathcal{L}}$ profiles that follow a convex shape. A constrained optimization then fit the profile shape to a smooth sigmoidal equation. Inputting the optimized ${\mathcal{L}}$ profiles in the forward and inverse Lagrangian equations leads to strong correlations between measured and calculated concentrations and fluxes. Coefficients of the sigmoidal equation were specific to each 30-min period and did not scale with any measured variable. Plausible looking ${\mathcal{L}}$ profiles were associated with negative bulk Richardson number values. Once the canopy senesced, a simple eddy diffusivity profile sufficed to relate concentrations and sources in the analytical Lagrangian equations.  相似文献   

18.
Using large-eddy simulation (LES), the effects of mesoscale local surface heterogeneity on the temporal evolution of low-level flows in the convective boundary layer driven by two-dimensional surface heat-flux variations are investigated at a height of about 100 m over flat terrain. The surface variations are prescribed with sinusoids of wavelength 32 km and varying amplitudes of 0, 50, 100, and 200 W m $^{-2}$ . The Weather Research and Forecasting numerical model is used as a mesoscale-domain LES model that has a grid spacing fine enough to explicitly resolve energy-containing turbulent eddies and a model domain large enough to include mesoscale circulations. Mesoscale circulations induced by the two-dimensional surface heterogeneity may undergo a flow transition and an associated spectral energy cascade, which has been found previously but only with one-dimensional surface heat-flux variations. Over a strongly heterogeneous surface prescribed with a two-dimensional sinusoid of amplitude 200 W m $^{-2}$ , the domain-averaged variance of the horizontal wind component initially grows rapidly, then undergoes a flow transition and subsequently rapidly decays. With a background wind, the induced mesoscale circulations are inhibited in the streamwise direction. However in the spanwise direction, somewhat stronger mesoscale circulations are induced, compared with those with no background wind. The background wind attenuates the significant reduction of the low-level temperature gradient by the fully-developed mesoscale horizontal flow. Spectral decomposition reveals that this rapid transition also exists in the mesoscale horizontal flows induced by the intermediate surface heterogeneity prescribed with a sinusoid of amplitude 100 W m $^{-2}$ . However the transition is masked by continuously growing turbulence.  相似文献   

19.
Measurements of vertical fluxes and concentration differences above a spring wheat crop (height $h=0.9$ – $0.95$  m, row spacing 0.25 m, displacement height $d=0.5$ – $0.6$  m) were analyzed to determine the Schmidt numbers for water vapour ( $S^\mathrm{v}$ ) and carbon dioxide ( $S^\mathrm{c}$ ) based on concentration differences between intakes 2.55 and 3.54 m above the ground. During nearly-neutral stratification $S^\mathrm{v}(0) = 0.68 \pm 0.1$ while $S^\mathrm{c} = 0.78 \pm 0.2$ , implying that the roughness sublayer extended above $2.5 h$ .  相似文献   

20.
Sensible heat fluxes ( \(Q_{H}\) ) are determined using scintillometry and eddy covariance over a suburban area. Two large-aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2–4 \(\text {km}^{2}\) , whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5–10 \(\text {km}^{2}\) . These large-scale heat fluxes are compared with local-scale eddy-covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly \(Q_{H}\) is related to the meteorological conditions. At shorter time scales the response of \(Q_{H}\) to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation ( \(Q^{*}\) ) coincide with greater differences between heat fluxes. For clear days \(Q_{H}\) lags \(Q^{*}\) , thus the ratio of \(Q_{H}\) to \(Q^{*}\) increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号