A simple method using climatic variables to estimate canopy temperature,sensible and latent heat fluxes in a winter wheat field on the North China Plain

Estimation of evapotranspiration from a crop field is of great importance for detecting crop water status and proper irrigation scheduling. The Penman–Monteith equation is widely viewed as the best method to estimate evapotranspiration but it requires canopy resistance, which is very difficult to determine in practice. This paper presents a simple method simplified from the Penman–Monteith equation for estimating canopy temperature (T_c). The proposed method is a biophysically‐sound extended version of that proposed by Todorovic. The estimated canopy temperature is used to calculate sensible heat flux, and then latent heat flux is calculated as the residual of the surface energy balance. An eddy covariance (EC) system and an infrared thermometer (IRT) were installed in an irrigated winter wheat field on the North China Plain in 2004 and 2005, to measure T_c, and sensible and latent heat fluxes were used to test the modified Todorovic model (MTD). The results indicate that the original Todorovic model (TD) severely underestimates T_c and sensible heat flux, and hence severely overestimates the latent heat flux. However, the MTD model has good capability for estimating T_c, and gives acceptable results for latent heat flux at both half‐hourly and daily scales. The MTD model results also agreed well with the evapotranspiration calculated from the measured T_c. Copyright © 2008 John Wiley & Sons, Ltd.     

Assessment of detachment and sediment transport capacity of runoff by field experiments on a silt loam soil

The detachment capacity (D_c) and transport capacity (T_c) of overland flow are important variables in the assessment of soil erosion. They determine respectively the lower and upper limit of sediment transport by runoff and therefore control detachment and deposition pro‐cesses. In this study, the detachment and transport capacity of runoff was investigated by rainfall simulations and overland flow experiments on small field plots. On the bare field plots, it was found that T_c was strongly related to total runoff discharge. This was also observed for the plots covered by maize residues, but T_c was less due to the lower runoff velocity. A simple regression equation was derived to estimate T_c for both bare and covered soil. Comparing our observations with T_c equations mentioned in the literature revealed that T_c equations based on laboratory experiments overestimated, on average, our measurements. Although T_c can be assessed more easily in laboratory experiments, the applicability of the results to field conditions remains questionable. Detachment by runoff was also related to total runoff discharge. The D_c values were, however, 4–50 times smaller than the T_c at corresponding high and low runoff discharge. This indicates that detachment by runoff constitutes only part of the transported sediment. Interrill erosion supplies an important additional amount of sediment. In this study, however, only sealed soils were considered. In the case of freshly tilled, loose soils, the D_c of runoff may be larger, resulting in a larger contribution to the total soil loss. Copyright © 2008 John Wiley & Sons, Ltd.     

A note on tsunami amplitudes above submarine slides and slumps

Tsunami generated by submarine slumps and slides are investigated in the near-field, using simple source models, which consider the effects of source finiteness and directivity. Five simple two-dimensional kinematic models of submarine slumps and slides are described mathematically as combinations of spreading constant or slopping uplift functions. Tsunami waveforms for these models are computed using linearized shallow water theory for constant water depth and transform method of solution (Laplace in time and Fourier in space). Results for tsunami waveforms and tsunami peak amplitudes are presented for selected model parameters, for a time window of the order of the source duration.The results show that, at the time when the source process is completed, for slides that spread rapidly (c_R/c_T≥20, where c_R is the velocity of predominant spreading), the displacement of the free water surface above the source resembles the displacement of the ocean floor. As the velocity of spreading approaches the long wavelength tsunami velocity the tsunami waveform has progressively larger amplitude, and higher frequency content, in the direction of slide spreading. These large amplitudes are caused by wave focusing. For velocities of spreading smaller than the tsunami long wavelength velocity, the tsunami amplitudes in the direction of source propagation become small, but the high frequency (short) waves continue to be present. The large amplification for c_R/c_T1 is a near-field phenomenon, and at distances greater than several times the source dimension, the large amplitude and short wavelength pulse becomes dispersed.A comparison of peak tsunami amplitudes for five models plotted versus L/h (where L is characteristic length of the slide and h is the water depth) shows that for similar slide dimensions the peak tsunami amplitude is essentially model independent.     

The response of CO2 flux to rain pulses at a saline desert

As the substantial component of the ecosystem respiration, soil CO₂ flux is strongly influenced by infrequent and unpredictable precipitation in arid region. In the current study, we investigated the response of soil CO₂ flux to rain pulses at a saline desert in western China. Soil CO₂ flux was measured continuously during the whole growing season of 2009 at six sites. We found that there were remarkable changes in amplitude or diurnal patterns of soil CO₂ flux induced by rainfall events: from bimodal before rain to a single peak after that. Further analysis indicated that there is a significant linear relationship (P < 0.001) between soil CO₂ flux and soil temperature (T_soil). However, a hysteresis between the waveform of diurnal course of CO₂ flux and T_soil was observed: with soil CO₂ flux always peaked earlier than T_soil. Furthermore, a double exponential decay function was fitted to the soil CO₂ flux after rainfall, and total carbon (C) releases were estimated by numerical integration for rainfall events. The relative enhancement and total C release, in association with the rain pulses, was linearly related to the amount of precipitation. According to the size and frequency of rainfall events, the total amount of C release induced by rain pulses was computed as much as 7.88 g C·m^–2 in 2009, equivalent to 10.25% of gross primary production. These results indicated that rain pulses played a significant role in the carbon budget of this saline desert ecosystem, and the size of them was a good indicator of rain‐induced flux enhancement. Copyright © 2012 John Wiley & Sons, Ltd.     

Notes on the variation of magnetization within basalt lava flows and dikes

Summary The magnetic properties of basaltic rocks are dominated by the contained primary Fe–Ti oxides. At solidus temperature (1000°C) the composition of these primary oxides is restricted to titanomagnetite (Fe_3-xTi_xO₄) and hemoilmenites (Fe_2-yTi_yO₃). The examination of 269 chemical analyses of the primary Fe–Ti oxides in basalts (in sensu lato) gives an average ofx=0.61 (T _c=168°C) for the titanomagnetites andy=0.89 (T _c=–121°C) for the hemoilmenites. If distinction is made between tholeiites, alkali basalts and andesites, a clear difference for thex-values is observed: the average for tholeiitesx=0.64 (T _c=144°C), for alkali basaltsx=0.52 (T _c=253°C), for andesitesx=0.38 (T _c=341°C).Environment of crystallization and cooling rate are major interrelated factors influencing subsequent changes in the mineralogy of the primary Fe–Ti oxides and resulting magnetic properties. This has been tested by studying the variation of magnetization and some of its parameters in three different basalt rock units: a dike, 180 cm, and two lava flows, 3 m and 33 m thick, respectively. Grain size and oxidation state of the titanomagnetites control the variation of magnetization in these basalt units.     

Nonlocal period parameters of frequency content characterization for near-fault ground motions

This paper focuses on the frequency property analysis of near-fault ground motions with and without distinct pulses, separately from the Chi-Chi and Northridge earthquakes. Ten scalar period parameters of ground motions, especially several nonlocal period parameters, are considered. Two new nonlocal parameters, namely the mean period of Hilbert marginal spectrum (T_mh) and the improved characteristic period (T_gi), are suggested. Moreover, comprehensive comparison and analysis indicate that T_mh, T_gi and T_avg (average spectral period) can distinguish the low-frequency components of near-fault ground motions; T_m (mean period of Fourier amplitude spectrum) and T_o (smoothed spectral predominant period) represent the moderate- and high-frequency components, respectively. The variance coefficient of predominant instantaneous frequency of Hilbert spectrum (H_cov) can be regarded as an alternative index to measure the non-stationary degree of near-fault ground motions. Finally, the velocity pulses and earthquake magnitude remarkably affect the frequency parameters of near-fault ground motions. Copyright © 2012 John Wiley & Sons, Ltd.     

Variation in canopy conductance of Cunninghamia lanceolata forest and its response to environmental factors after an extreme rainfall event

Canopy conductance (g_c) is a crucial parameter in simulating evapotranspiration and modulating water exchange, but its variation mechanism has regional uncertainties and complex environmental co-controls. In addition, the effect of extreme rainfall on g_c cannot be ignored under the changing climate. Here, we investigated the variation and environmental controls on g_c and the effect of extreme rainfall events in a Cunninghamia lanceolata forest across the subtropical area of Southern China. In July 2020, an extreme rainstorm hit the source area of the Xin'an River, with the cumulative rainfall on July 7 and 8 reaching 216.6 mm. The thermal diffusion probe method was used to measure the density of sap flow, and the environmental factors such as air temperature (T_a), net solar radiation (R_n) and soil water content (SWC) were monitored during the growing seasons of 2020 and 2021. Ultimately, g_c obtained by the Penman-Monteith equation was adopted since the result from the Köstner equation was overestimated. g_c showed a unimodal curve on the diurnal scale, and this characteristic was more obvious after the extreme rainfall. Daily g_c appeared a fluctuating pattern with a maximum in summer. g_c was simultaneously affected by T_a, R_n, water vapour pressure difference (VPD), SWC, among which T_a was the most significant driving factor at both the diurnal and daily timescales. The regulation of T_a, VPD and SWC on g_c had obvious thresholds, and the most definite response mode was VPD (2020: 1.25 kPa; 2021: 0.95 kPa). SWC and T_a were the dominant factors after the rainfall period, and the promotion effect of VPD on post-rainy days turned to inhibition effect on typical sunny days. These findings will further reveal the water exchange mechanism between atmosphere and vegetation and impacts of environmental factors in subtropical coniferous forests, especially after the extreme rainfall events.     

A PHYSICAL SCALE MODEL STUDY OF THE COMPARATIVE PERFORMANCE OF TWO MODES OF OPERATION FOR FIXED-LOOP TURAM-TYPE EM SYSTEMS1

The results of a physical scale model study of the conventional mode of operation of fixed loop electromagnetic systems and an alternative mode called the tx-parallel mode in which traverses are run parallel to the long axis of the rectangular transmitter loop are presented. The results show that over thick or dipping conductors, the tx-parallel configuration provides coupling with the target which is comparable with that provided by the conventional configuration. In addition, the tx-parallel configuration is shown to provide more consistent indications of the direction and magnitude of conductor dip. Over wide conductors, where separate conventional surveys are needed to define the opposite edges of such conductors, it is shown that only a single tx-parallel survey is needed to locate both edges of the conductor. The tx-parallel results were found to allow better resolution of the individual anomalies due to closely spaced parallel conductors. The tx-parallel response of identical sheet conductors of opposite dip indicated that the response of the separate sheets could be recognized even when the two sheets were placed at zero separation. This was found to be due in part to spatial displacement of the individual current vortices within each conductor owing to their mutual repulsion.     

Relations of sand trapping efficiency and migration speed of transverse dunes to wind velocity

A self‐consistent model which describes transverse dune migration in equilibrium is introduced. It shows that an equilibrium expression for dune migration speed (c _d) must take into account sand trapping efficiency (T _E), and that T _E is strongly related to the wind speedup over the windward surface. An expression for sand trapping efficiency (T _E) is analytically derived from a microscale analysis of sand grain deposition on the slip face. Sand trapping efficiency (T _E) is mainly determined by shear velocity on a level surface (u_*(−∞)), and rapidly decreases as u_*(−∞) increases. For each dune height (H), dune migration speed (c _d) first increases, and then decreases monotonically after reaching the maximum, as the shear velocity on a level surface (u_*(−∞)) increases. Dune migration speed (c _d) is not inversely proportional to dune height (H). For low dunes, small sand trapping efficiency (T _E) suppresses c _d, whereas for high dunes, wind speedup and large T _E resist the decrease of c _d. Some field data show the same tendency. The dune‐to‐plane‐bed transition observed in subaqueous and venusian bedforms could be associated with the decrease of sand trapping efficiency (T _E). Copyright © 2000 John Wiley & Sons, Ltd.     

Estimating time parameters of overland flow on impervious surfaces by the particle tracking method

Abstract

Travel time and time of concentration T_c are important time parameters in hydrological designs. Although T_c is the time for the runoff to travel to the outlet from the most remote part of the catchment, most researchers have used an indirect method such as hydrograph analysis to estimate T_c. A quasi two-dimensional diffusion wave model with particle tracking for overland flow was developed to determine the travel time, and validated for runoff discharges, velocities, and depths. Travel times for 85%, 95% and 100% of particles arrival at the outlet of impervious surfaces (i.e. T_t85, T_t95, and T_t100) were determined for 530 model runs. The correlations between these travel times and T_c estimated from hydrograph analysis showed a significant agreement between T_c and T_t85. All the travel times showed nonlinear relationships with the input variables (plot length, slope, roughness coefficient, and effective rainfall intensity) but showed linear relationships with each other.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

Performance of methods for estimating the time of concentration in a watershed of a tropical region

The time of concentration (T_c) is a fundamental parameter in the design of hydrological projects for watersheds. In this study a graphical methodology is described for estimating T_c in a watershed, and this is applied to 17 rainfall–runoff events from a rural watershed located near the capital city of Mato Grosso do Sul State, in the Brazilian Cerrado. The T_c values obtained through the graphical method were compared to T_c values estimated using 20 equations from various references. The equations were selected by considering those that were not developed using data for urban watersheds, and the results of the graphical method were compared to those derived by applying the equations to sub-basin data. The graphical method was reliable in determining T_c, and Ventura’s equation was found to present the best performance for a rural watershed in a tropical climate region.     

Convection in a rotating cylinder with its sidewall having a finite thermal conductance

Abstract

An investigation is made of steady thermal convection of a Boussinesq fluid confined in a vertically-mounted rotating cylinder. The top and bottom endwall disks are thermal conductors at temperatures T_t and T_b with δT = T_t ? T_b >0. The vertical sidewall has a finite thermal conductance. A Newtonian heat flux condition is adopted at the sidewall. The Rayleigh number of the fluid system is large to render a boundary layer-type flow. Finite-difference numerical solutions to the full Navier-Stokes equations are obtained. The vertical motions within the buoyancy layer along the sidewall induce weak meridional flows in the interior. Because of the Coriolis acceleration, the meridional flows give rise to azimuthal flows relative to the rotating container. Strong vertical gradients of azimuthal flows exist in the regions near the endwalls. As the stratification effect increases, concentration of flow gradients in thin endwall boundary layers becomes more pronounced. The azimuthal flow field exhibits considerable horizontal gradients. The temperature field develops horizontal variations superposed on the dominant vertical distribution. As either the sidewall thermal conductance or the stratification effect decreases, the temperature distribution tends to the profile varying linearly with height. Comparisons of the sizes of the dynamic effects demonstrate that, in the bulk of flow field, the vertical shear of azimuthal velocity is supported by the horizontal temperature gradient, resulting in a thermal-wind relation.     

Turbulent transport of magnetic fields. V. Distribution of magnetic energy in a simple α2-dynamo

Abstract

In this paper we analyse the stationary mean energy density tensor T_ij = B_iB_j for the x ²-sphere. This model is one of the simplest possible turbulent dynamos, originally due to Krause and Steenbeck (1967): a conducting sphere of radius R with homogeneous, isotropic and stationary turbulent convection, no differential rotation and negligible resistivity. The stationary solution of the (linear) equation for T_ij is found analytically. Only T_rr , T _θθ and T _φφ are unequal to zero, and we present their dependence on the radial distance r.

The stationary solution depends on two coefficients describing the turbulent state: the diffusion coefficient β≈?u²?τ_c/3 and the vorticity coefficient γ ≈ ?|?×u|²?τ_c/3 where u(r, t) is the turbulent velocity and _c its correlation time. But the solution is independent of the dynamo coefficient α≈??u·?×u?τ_c/3 although α does occur in the equation for T_ij . This result confirms earlier conclusions that helicity is not required for magnetic field generation. In the stationary state, magnetic energy is generated by the vorticity and transported to the boundary, where it escapes at the same rate. The solution presented contains one free parameter that is connected with the distribution of B over spatial scales at the boundary, about which T_ij gives no information. We regard this investigation as a first step towards the analysis of more complicated, solar-type dynamos.     

Rundle模式地震滞后环和突变现象研究

基于耗散结构理论详细研究了地震活动的主要特征。研究结果表明，地震活动呈现双稳态及尖拐突变现象，主要表现是：发生在应力临界值区间[-pc,pc]内的滞后现象和发生在临界值pc和-pc处的突变现象，前对应系统变能积累过程，后（即强震过程）地该系统能量的主要耗散机制。     

Lapse time and frequency-dependent attenuation characteristics of coda waves in the Northwestern Himalayas

We analyzed the local earthquakes waveform recorded on a broadband seismic network in the northwestern Himalayan Region to compute lapse time and frequency dependence of coda Q (Q _c). The observed Q _c values increase with increasing lapse time at all frequency bands. The increase in Q _c values with lapse time is attributed to an increase in Q _c with depth. This implies that attenuation decreases with increasing depth. The approximate radius of medium contributing to coda generation varies from 55 to 130 km. By comparing the Q _c values with those from other regions of the world, we find that they are similar to those obtained from tectonically active regions. The estimated Q _c values show a frequency-dependent relationship, Q _c = Q ₀ f ⁿ , where Q ₀ is Q _c at 1 Hz and n represents degree of frequency dependence. They represent the level of heterogeneity and tectonic activity in an area. Our results show that northwest Himalayas are highly heterogeneous and tectonically very active. Q ₀ increases from 113 ± 7 to 243 ± 10 and n decreases from 1.01 ± 0.05 to 0.85 ± 0.03 when lapse time increases from 30 to 70 s. As larger time window sees the effect of deeper part of the Earth, it is concluded that Q ₀ increases and n decreases with increasing depth; i.e., heterogeneity decreases with depth in the study area.     

TRANSIENT ELECTROMAGNETIC FIELD OF A POLARIZABLE HALF-SPACE DUE TO VARIOUS CURRENT PULSES *

In this paper, the expressions for the induced voltage function V(T) and mutual impedance function Z_m(T) have been derived for specific range of time T and specific values of induced polarization parameters a for ramp and saw-tooth type of current pulses. The computed results for various cases are also presented. The low values of induced polarization parameter a represent the medium possessing membrane polarization, whereas high values exhibit electrode polarization medium. The method has practical applicability and is best suited for the interpretation of transient electromagnetic fields over a polarizable half-space.     

APPROXIMATE INVERSION OF AIRBORNE EM DATA FROM A MULTILAYERED GROUND1

A complex transfer function c (or generalized skin depth) can be derived from data for the secondary magnetic field measured by a dipole system with small coil spacing at height h above the ground. This function has a useful property: For a uniform or layered ground, the real part of c yields the‘ centroid depth’z* of the in-phase current system as a function of frequency. This parameter can be combined with the apparent resistivity ρ_a derived by conventional methods. The function ρ_a(z*), if known over a broad frequency range, yields a smoothed approximation of the true distribution ρ(z) without an initial model. The relations between ρ_a(z*) and ρ(z) are studied for a number of multilayer models. An example of the application of the ρ_a*) algorithm to data from a groundwater survey is given.     

MEASURED UNDERWATER NEAR-FIELD E-PATTERNS OF A PULSED,HORIZONTAL DIPOLE ANTENNA IN AIR: COMPARISON WITH THE THEORY OF THE CONTINUOUS WAVE,INFINITESIMAL ELECTRIC DIPOLE1

At Delft Geotechnics the technique of ground-penetrating radar is in use for the detection of buried objects such as pipes. To enable us to give our ‘measurements in the field’ a more quantitative interpretation than can be deduced from these alone, a series of experiments has been started under well-defined conditions. A cylindrical vessel containing water simulates wet soil. Mounted horizontally above the water surface is a pulsed triangular half-wave dipole which is used as a transmitting antenna (TA). It has a carrier-frequency of about 160 MHz and a pulse repetition-frequency of about 50 kHz. A movable receiving dipole (‘probe’) in the water measures the transverse, mutually orthogonal E_φ,- and E_θ-components of the pulses as a function of probe-position (r, θ, φ) and of the height h of the TA above the water surface. When these pulses are Fourier-transformed, the transverse electric fields E_φ and E_θ at 200 MHz are obtained. The resulting field patterns are compared with computational results on the basis of the theory of the continuous wave, infinitesimal electric dipole (‘point dipole’). It can be concluded that:

• 1 Far-field conditions have not fully developed at a depth of about 2.50 m, the largest value of the radius r at which field patterns were measured, although it represents a distance of about 15 wavelengths.
• 2 The attenuation constant of the tapwater used, as deduced from E-field measurements for θ= 0, 2.50 m < r < 2.75 m, is slightly less than the value measured using a network analyser and air line combination, in agreement with (1).
• 3 E _φ field patterns calculated using the value of the conductivity σ corresponding to the former value of the attenuation constant agree reasonably well with the measured patterns for r≤ 2.50 m and for θ < 20° at all antenna heights considered. Calculated Eφ patterns do not agree so well with the measured patterns when h is close to zero. With increasing height the agreement inproves.
• 4 In accordance with the theory of the point-dipole, the angular distribution of the radiation patterns of the TA becomes wider as the frequency decreases.
• 5 The normalized underwater pulse-spectra shift to lower frequencies with increasing r. This can be explained since the attenuation constant of the water rises with rising frequency.

     

Comparison of Magnetic Properties and Shadow Area of Leading and Trailing Spots with Different Asymmetries

The magnetic properties of the shadow of magnetic-related leading and trailing spots (those connected by forces lines of magnetic field, which are calculated from a field in potential approximation) are studied in this work. The correlations are established between individual characteristics of the field in the spot shadow and these characteristics from the shadow area S for spot pairs, for which the minimum angle between the measured vector of magnetic induction B in the shadow of the leading (L) spot and positive normal to the solar surface is lower than in the trailing (F) spot (α_min-L < α_min-F) and, vice versa, when α_min-L > α_min-F. It is shown that the α_min-L(S_L), α_min-F(S_F), B_max-L(S_L) and B_max-F(S_F) correlations are similar behaviorally and quantitatively for two groups of spots with different asymmetries of a magnetically connected field (B_{max-L, F} is the maximum of magnetic induction in the shadow of leading and trailing spots). The correlation between the average angles within the spot shadow 〈α_{L, F}〉 and the area of the spot shadow S_{L, F} and between the average value of magnetic induction in the spot shadow 〈B_{L, F}〉 differ in two cases. In most studied spot pairs, the leading spot is closer to the dividing line of polarity between the spots rather than the trailing one.     

Characterization of volcanic thermal anomalies by means of sub-pixel temperature distribution analysis

The simultaneous solution of the Planck equation (involving the widely used “dual-band” technique) using two shortwave infrared (SWIR) bands allows for an estimate of the fractional area of the hottest part of an active lava flow (f _h) and the background temperature of the cooler crust (T _c). The use of a high spectral and spatial resolution imaging spectrometer with a wide dynamic range of 15 bits (DAIS 7915) in the wavelength range from 0.501 to 12.67 μm resulted in the identification of crustal temperature and fractional areas for an intra-crater hot spot at Mount Etna, Italy. This study indicates the existence of a relationship between these T _c and f _h extracted from DAIS and Landsat TM data. When the dual band equation system is performed on a lava flow, a logarithmic distribution is obtained from a plot of the fractional area of the hottest temperature vs. the temperature of the cooler crust. An entirely different distribution is obtained over active degassing vents, where increases in T _c occur without any increase in f _h. This result indicates that we can use scatter plots of T _c vs. fh to discriminate between different types of volcanic activity, in this case between degassing vents and lava flows, using satellite thermal data.