首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Neutral surface layer flow over low hills and varying surface roughness is considered with emphasis on closure schemes in relation to the prediction of turbulence quantities. The equations are linearised, Fourier transformed in the two horizontal directions and solved by means of a finite difference method in the vertical. Three closure schemes are. employed, namely mixing length, E- and e-- closure where E, and indicate that differential equations are used for turbulent kinetic energy, dissipation rate and shear stress. Model calculations are compared with experimental data for the step in roughness problem and for the Askervein hill. The mean flow results turn out to be relatively insensitive to the closure scheme. The shear stress and the dimensionless shear, however, are much better predicted with the E- equations than with mixing length closure. In the outer layer of the hill problem, advection of shear stress becomes important. An equation for is needed here.  相似文献   

2.
In a recent paper, the author introduced a new viscous boundary layer, called the mesolayer, in turbulent shear flow. Its importance stems from its location between the inner and outer regions which are controlled by the law of the wall and Reynolds number similarity, respectively. This intrusion prevents the classical overlap assumption which appears to be fundamental in the derivation of the classical logarithmic behavior. The mesolayer has a thickness proportional to Taylor's microscale . This, and the analogy between the energy equation for the spectrum function of isotropic turbulence and the momentum equation for shear flow, suggest the existence of a similar region in wavenumber space with wavenumber k ~ -1. This mesoregion separates the inner region k ~ k s(where k s-1 and is the Kolmogorov length) and the outer region k k e(where k e -1 and l is the energy-containing eddy size) and again invalidates the overlap assumption which appears to be fundamental in the derivation of the classical k -5/3-behavior of the inertial subrange.Incorporation of the mesoregion into the argument leads to a new theory with k -5/3-behavior in two regions (-1 k k s) and (k e k -1) although with two different coefficients of proportionality (Kolmogorov constants). This leads to a wandering of the spectrum curve about the classical k -5/3 line similar to a wandering in turbulent shear flow about the logarithmic curve. This is clearly indicated by the data for the variation of the Kolmogorov constant.Other data support the new theory. In particular, the location of the point k mwhere the curve of the nonlinear energy-transfer function goes through zero shows agreement with the theory, i.e., k m-1.  相似文献   

3.
From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - Ci constants i=u, v, w, (see Equation (5) - kme,kmT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - kmi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - ksb, kc characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - m nondimensional mean shear - H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - D nondimensional divergence of vertical transports of turbulent linetic energy  相似文献   

4.
Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T * ). Normal atmospheric monsoon circulations are defined by long-term average winds (u n,v n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u L,v L); 7–20 day filtered winds (u M,v M); and 2–6 day filtered winds (u S,v S). Three-month running mean anomaly kinetic energy (signified asK L , K M , andK S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK L , K M , andK S perturbations propagate slowly eastward with same phase speed (0.3 ms–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K L , K M , K S ) and interannualu (OLR) modes. Namely, (K L , K M , K S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K L , K M , K S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T * 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT * variations on interannual time scales. HighestT * occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT * variations are also controlled by the prevailing monsoonal zonal windsu n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.  相似文献   

5.
Analytical solutions for the Ekman layer   总被引：1，自引：0，他引：1  
The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U *z (0.4,U *=friction velocity,z=elevation) to the geostrophic asymptoteK GU *d forzd is solved through an expansion in fd/U *1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU *d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U *d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U * (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U * and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U *d[1–exp(–z/d)] are developed in appendices.  相似文献   

6.
Effect of finite sampling on atmospheric spectra   总被引：2，自引：0，他引：2  
The effect of a finite averaging time on variances is well known, but its effect on power spectra is less clearly understood. We present numerical solutions for the spectral distortion arising from sampling over a finite time interval T and show that the commonly used filter function (1 – sinc2f T), valid for variances, is a reasonable approximation for power spectra only when T 10 m , where f is the cyclic frequency, and m is the dominant time scale of the process. Our results exhibit an increasingly steeper low-frequency roll-off as T decreases relative to m , indicating that the measured spectrum is subject to a greater suppression of the lower frequencies (f > 1/T) than predicted by (1 – sinc2f T). This suppression is, in a sense, compensated by an overestimation of spectral estimates in the frequency range f 1/T.  相似文献   

7.
Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO2 concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO2, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO2. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= cw for CO2, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux c w. Thus, when the net flux of CO2 is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., c w = – c w. Corrections for the mean convective flux are particularly significant for CO2 because cw and c w are often of similar magnitude. The correct measurement of the net CO2 flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.  相似文献   

8.
An E- turbulence model is used to study air-sea interaction characteristics and turbulence structure using a coupled model for air-sea boundary layers. The E- turbulence model consists of equations for the turbulent kinetic energy, the energy-dissipation, and for the turbulent exchange coefficient expressed in terms of turbulent kinetic energy and energy-dissipation. The energy-dissipation equations for the air-sea interface are solved analytically to obtain boundary conditions for energy-dissipation at the interface. The air-sea interaction and turbulence characteristics of the E- model are compared with those of the mixing-length model and with available observations.The simulations demonstrate that the air-sea interaction parameters obtained by the E- model agree well with observations. The numerical studies also show that the E- turbulence model with appropriate constants can give good results in modeling coupled air-sea boundary-layer flows.  相似文献   

9.
Summary The effect of the Alpine orography on prototype cold fronts approaching from the west is investigated by three-dimensional numerical model simulations. The numerical experiments cover a range of parameter constellations which govern the prefrontal environment of the front. Especially, the appearance and intensity of prefrontal northern Alpine foehn varies from case to case.The behaviour of a cold front north of the Alps depends much on the prefrontal condition it encounters. It is found that prefrontal foehn can either accelerate or retard the approaching front.An important feature is the pressure depression along the northern Alpine rim that results from the southerly foehn flow. In cases where this depression compensates the eastward directed pressure gradient associated with the largescale flow, the front tends to accelerate and the foehn breaks down as soon as the front passes. In contrast, the foehn prevents the front from a rapid eastward propagation if it is connected with a strong southerly wind component.No-foehn experiments are performed for comparison, where either the mountains are removed, or the static stability is set to neutral. Also shown are effects of different crossfrontal temperature contrasts.List of Symbols c F propagation speed of a front - x, y horizontal grid spacing (cartesian system) - , horizontal grid spacing (geographic system) - t time step - z vertical grid spacing (cartesian system) - cross-frontal potential temperature difference - i potential temperature step at an inversion - E turbulent kinetic energy - f Coriolis parameter - FGP frontogenesis parameter (see section 2.2) - g gravity acceleration (g=9.81 m s–2) - vertical gradient of potential temperature - h terrain elevation (above MSL) - h i height of an inversion (h i =1000 m MSL) - H height of model lid (H=9000 m MSL) - K M exchange coefficient of momentum - K H exchange coefficient of heat and moisture - longitude - N Brunt-Väisäla-frequency - p pressure - Exner function (=T/) - latitude - q v specific humidity - R d gas constant of dry air (R d =287.06 J kg–1 K–1) - density of dry air - t time - T temperature - potential temperature - TFP thermal front parameter (see section 2.2) - u, v, w cartesian wind components - u g ,v g geostrophic wind components - horizontal wind vector - x, y, z cartesian coordinates Abbreviations GND (above) ground level - MSL (above) mean sea level - UTC universal time coordinated With 20 Figures  相似文献   

10.
Summary The integral aerosol optical depths (k ) at the hour of 08:20 Local Standard Time (LST), are compared with those calculated previously at 11:20 and 14:20 LST, for clear days during summer in Athens over the period 1962–1988. The mean values at 08:20 LST were consistently lower than the values at 11:20 and 14:20 LST. The influence of the vertical wind profile on the values ofk was also investigated. A comparison was made of the wind profiles at 02:00 and 14:00 LST, for days in which the 11:20 and 14:20 LST values ofk were 0.200 andk 0.350, respectively. The corresponding bulk wind shear s was also found for the period 1980–1988. The most significant results occurred with the first category of days. The resultant wind velocities from the surface to the 900 hPa level, in each hour were higher by 2–4 m·s–1 with respect to the corresponding values for the second category. At 02:00 LST the bulk wind shear showed a considerable difference (1.8) between the two categories of days in the surface to 700 hPa layer at 02:00 LST. Finally, the associated weather conditions that appear to initiate a period of low values ofk (k 0.200) at 11:20 and 14:20 LST were examined for the period 1980–1988. Fifteen such cases were identified and it was found that they all occurred after the passage of weak cold fronts.With 6 Figures  相似文献   

11.
The characteristics of a Lyman-alpha humidiometer have been carefully examined in an air-conditioned test chamber. The results confirm that when carefully used, this humidiometer is suitable for measurements of turbulent humidity fluctuations. Measurements with a Lyman-alpha humidiometer were carried out in the surface boundary layer over the ocean. The relation between turbulent intensity ( a = a ov2) and the friction humidity (a *) can be expressed as a = l.6a *. The spectrum of turbulent humidity for wind speeds larger than 3 m s –1 conforms to the similarity law in the surface boundary layer. The spectrum has two characteristic normalized frequencies, namely, a higher peak and a secondary peak (or a shoulder).  相似文献   

12.
Turbulence measurements performed in a stable boundary layer over the sloping ice surface of the Vatnajökull in Iceland are described. The boundary layer, in which katabatic forces are stronger than the large-scale forces, has a structure that closely resembles that of a stable boundary layer overlying a flat land surface, although there are some important differences. In order to compare the two situations the set-up of the instruments on an ice cap in Iceland was reproduced on a flat grass surface at Cabauw, the Netherlands. Wind speed and temperature gradients were calculated and combined with flux measurements made with a sonic anemometer in order to obtain the local stability functions m and h as a function of the local stability parameter z/L. Unlike the situation at Cabauw, where m was linear as a function of z/L, in the katabatically forced boundary layer, the dependence of m on stability was found to be non-linear and related to the height of the wind maximum. Thermal stratification and the depth of the stable boundary layer however seem to be rather similar under these two different forcing conditions.Furthermore, measurements on the ice were used to construct the energy balance. These showed good agreement between observed melt and components contributing to the energy balance: net radiation (supplying 55% of the energy), sensible heat flux (30%) and latent heat flux (15%).Local sources and sinks in the turbulent kinetic energy budget are summed and indicate a reasonable balance in near-neutral conditions but not in more stable situations. The standard deviation of the velocity fluctuations u, v, and w, can be scaled satisfactorily with the local friction velocity u* and the standard deviation of the temperature fluctuation with the local temperature scale *.  相似文献   

13.
Summary A radiative transfer model has been used to determine the large scale effective 6.6 GHz and 37 GHz optical depths of the vegetation cover. Knowledge of the vegetation optical depth is important for satellite-based large scale soil moisture monitoring using microwave radiometry. The study is based on actual observed large scale surface soil moisture data and observed dual polarization 6.6 and 37 GHz Nimbus/SMMR brightness temperatures over a 3-year period. The derived optical depths have been compared with microwave polarization differences and polarization ratios in both frequencies and with Normalized Difference Vegetation Index (NDVI) values from NOAA/AVHRR. A synergistic approach to derive surface soil emissivity from satellite observed brightness temperatures by inverse modelling is described. This approach improves the relationship between satellite derived surface emissivity and large scale top soil moisture fromR 2=0.45 (no correction for vegetation) toR 2=0.72 (after correction for vegetation). This study also confirms the relationship between the microwave-based MPDI and NDVI earlier described and explained in the literature.List of Symbols f frequency [Hz] - f i(p) fractional absorption at polarizationp - h surface roughness - h h cos2 - H horizontal polarization - n i complex index of refraction - p polarization (H orV) - R s microwave surface reflectivity - T B(p) brightness temperature at polarizationp - T * normalized brightness temperature - T polarization difference (T v-T H) - T s temperature of soil surface - T c temperature of canopy - T max daily maximum air temperature - T min daily minimum air temperature - V vertical polarization - soil moisture distribution factor; also used for the constant to partition the influence of bound and free water components to the dielectric constant of the mixture - empirical complex constant related to soil texture - microwave transmissivity of vegetation (=e –) - * effective transmissivity of vegetation (assuming =0) - microwave emissivity - s emissivity of smooth soil surface - rs emissivity of rough soil surface - vs emissivity of vegetated surface - soil moisture content (% vol.) - K dielectric constant [F·m–1] - K fw dielectric constant of free water [F·m–1] - K ss dielectric constant of soil solids [F·m–1] - K m dielectric constant of mixture [F·m–1] - K o permittivity of free space [8.854·10–12 F·m–1] - high frequency limit ofK wf [F·m–1] - wavelength [m] - incidence angle [degrees from nadir] - polarization ratio (T H/T V) - b soil bulk density [gr·cm–3] - s soil particle density [gr·cm–3] - R surface reflectivity in red portion of spectrum - NIR surface reflectivity in near infrared portion of spectrum - eff effective conductivity of soil extract [mS·cm–1] - vegetation optical depth - 6.6 vegetation optical depth at 6.6 GHz - 37 vegetation optical depth at 37 GHz - * effective vegetation optical depth (assuming =0) - single scattering albedo of vegetation With 12 Figures  相似文献   

14.
Previous results of non-dimensional wind and temperature profiles as functions of ( = z/L) show systematic deviations between different experiments. These discrepancies are generally believed not to reflect real differences but rather instrumental shortcomings. In particular, it is clear that flow distortion has not been adequately treated in most previous experiments. In the present paper, results are presented from a surface-layer field experiment where great care was taken to remove any effects from this kind of error and also to minimize other measuring errors. Data from about 90 30-min runs with turbulence measurements at three levels (3, 6, and 14 m) and simultaneous profile data have been analysed to yield information on flux-gradient relationships for wind and temperature.The flux measurements themselves show that the fluxes of momentum and sensible heat are constant within ± 7% on average for the entire 14 m layer in daytime conditions and when the stratification is slightly stable. For more stable conditions, the flux starts to decrease systematically somewhere in the layer 6 to 14 m. From a large body of data for near-neutral conditions (¦¦ 0.1), values are derived for von Kármán's constant: 0.40 ± 0.01 and for h at neutrally, 0.95 ± 0.04. The range of uncertainty indicated here is meant to include statistical uncertainty as well as the effect of possible systematic errors.Data for m and h for an extended stability range (1 > > – 3) are presented. Several formulas for m and h appearing in the literature have been used in a comparative study. But first all the formulas have been modified in accordance with the following assumptions: = 0.40 and ( h ) = 0 = 0.95; deviations from this result in the various studies are due to incomplete correction for flow distortion. After new corrections are introduced, the various formulas were compared with the present measurements and with each other. It is found that after this modification, the most generally used formulas for m and h for unstable conditions, i.e., those of Businger et al. (1971) and Dyer (1974) agree with each other to within ± 10% and with the present data. For stable conditions, the various formulas still disagree to some extent. The conclusion in relation to the present data is not as clear as for the unstable runs, because of increased scatter. It is, however, found that the modified curve of Businger et al. (1971) for h fits the data well, whereas for m , Dyer's (1974) curve appears to give slightly better agreement.  相似文献   

15.
E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c 2 /e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c 2 = 1.92, = 1.3, e = 1). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and e to 1.1 and 1.6, respectively, giving, with c 2 = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.  相似文献   

16.
The spectral equations of turbulent kinetic energy and temperature variance have been solved by using Onsager's energy cascade model and by extending Onsager's model to closure of terms that embody the interaction of turbulent and mean flow.The spectral model yields the following results: In a stably stratified shear flow, the peak wave numbers of the spectra of energy and temperature variance shift toward larger wave numbers as stability increases. In an unstably stratified flow, the peak wave numbers of energy spectra move toward smaller wave numbers as instability increases, whereas the opposite trend is observed for the peak wave numbers of temperature variance spectra. Hence, the peak wave numbers of temperature spectra show a discontinuity at the transition from stable to unstable stratification. At near neutral stratification, both spectra reveal a bimodal structure.The universal functions of the Monin-Obukhov similarity theory are predicted to behave as m ~ H ~ (- Z/L)-1/3 in an extremely unstable stratification and as m ~ H ~ z/L in an extremely stable stratification. For a stably stratified flow, a constant turbulent Prandtl number is expected.  相似文献   

17.
  相似文献   

18.
A two-dimensional mesoscale model has been developed to simulate the air flow over the Gulf Stream area where typically large gradients in surface temperature exist in the winter. Numerical simulations show that the magnitude and the maximum height of the mesoscale circulation that develops downwind of the Gulf Stream depends on both the initial geostrophic wind and the large-scale moisture. As expected, a highly convective Planetary Boundary Layer (PBL) develops over this area and it was found that the Gulf Stream plays an important role in generating the strong upward heat fluxes causing a farther seaward penetration as cold air advection takes place. Numerical results agree well with the observed surface fluxes of momentum and heat and the mesoscale variation of vertical velocities obtained using Doppler Radars for a typical cold air outbreak. Precipitation pattern predicted by the numerical model is also in agreement with the observations during the Genesis of Atlantic Lows Experiment (GALE).List of Symbols u east-west velocity [m s–1] - v north-south velocity [m s–1] - vertical velocity in coordinate [m s–1] - w vertical velocity inz coordinate [m s–1] - gq potential temperature [K] - q moisture [kg kg–1] - scaled pressure [J kg–1 K–1] - U g the east-south component of geostrophic wind [m s–1] - V g the north-south component of geostrophic wind [m s–1] - vertical coordinate following terrain - x east-west spatial coordinate [m] - y north-south spatial coordinate [m] - z vertical spatial coordinate [m] - t time coordinate [s] - g gravity [m2 s–1] - E terrain height [m] - H total height considered in the model [m] - q s saturated moisture [kg kg–1] - p pressure [mb] - p 00 reference pressure [mb] - P precipitation [kg m–2] - vertical lapse rate for potential temperature [K km–1] - L latent heat of condensation [J kg–1] - C p specific heat at constant pressure [J kg–1 K–1] - R gas constant for dry air [J kg–1 K–1] - R v gas constant for water vapor [J kg–1 K–1] - f Coriolis parameter (2 sin ) [s–1] - angular velocity of the earth [s–1] - latitude [o] - K H horizontal eddy exchange coefficient [m2 s–1] - t integration time interval [s] - x grid interval distance inx coordinate [m] - y grid interval distance iny coordinate [m] - adjustable coefficient inK H - subgrid momentum flux [m2 s–2] - subgrid potential temperature flux [m K s–1] - subgrid moisture flux [m kg kg–1 s–1] - u * friction velocity [m s–1] - * subgrid flux temperature [K] - q * subgrid flux moisture [kg kg–1] - w * subgrid convective velocity [m s–1] - z 0 surface roughness [m] - L Monin stability length [m] - s surface potential temperature [K] - k von Karman's constant (0.4) - v air kinematic viscosity coefficient [m2 s–1] - K M subgrid vertical eddy exchange coefficient for momentum [m2 s–1] - K subgrid vertical eddy exchange coefficient for heat [m2 s–1] - K q subgrid vertical eddy exchange coefficient for moisture [m2 s–1] - z i the height of PBL [m] - h s the height of surface layer [m]  相似文献   

19.
The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (Fa) and atmospheric moisture diffusivity (Kq). In a large ensemble of spinup experiments, the values of Fa and Kq are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of Fa and Kq, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.  相似文献   

20.
Progress in understanding how terrestrial ice volume is linked to Earths orbital configuration has been impeded by the cost of simulating climate system processes relevant to glaciation over orbital time scales (103–105 years). A compromise is usually made to represent the climate system by models that are averaged over one or more spatial dimensions or by three-dimensional models that are limited to simulating particular snapshots in time. We take advantage of the short equilibration time (10 years) of a climate model consisting of a three-dimensional atmosphere coupled to a simple slab ocean to derive the equilibrium climate response to accelerated variations in Earths orbital configuration over the past 165,000 years. Prominent decreases in ice melt and increases in snowfall are simulated during three time intervals near 26, 73, and 117 thousand years ago (ka) when aphelion was in late spring and obliquity was low. There were also significant decreases in ice melt and increases in snowfall near 97 and 142 ka when eccentricity was relatively large, aphelion was in late spring, and obliquity was high or near its long term mean. These glaciation-friendly time intervals correspond to prominent and secondary phases of terrestrial ice growth seen within the marine 18O record. Both dynamical and thermal effects contribute to the increases in snowfall during these periods, through increases in storm activity and the fraction of precipitation falling as snow. The majority of the mid- to high latitude response to orbital forcing is organized by the properties of sea ice, through its influence on radiative feedbacks that nearly double the size of the orbital forcing as well as its influence on the seasonal evolution of the latitudinal temperature gradient.  相似文献   

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

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