Climate, canopy conductance and leaf area development controls on evapotranspiration in a boreal coniferous forest over a 10-year period: A united model assessment

The aim of this work was to test a process-based model (hydrological model combined with forest growth model) on the simulation of seasonal variability of evapotranspiration (ET) in an even-aged boreal Scots pine (Pinus sylvestris L.) stand over a 10 year period (1999-2008). The water flux components (including canopy transpiration (E_t) and evaporation from canopy (E_c) and ground surface (E_g) were estimated in order to output the long-term stand water budget considering the interaction between climate variations and stand development. For validation, half-hourly data on eddy water vapor fluxes were measured during the 10 growing seasons (May-September). The model predicted well the seasonal course of ET compared to the measured values, but slightly underestimated the water fluxes both in non-drought and drought (2000, 2003 and 2006) years. The prediction accuracy was, on average, higher in drought years. The simulated ET over the 10 years explained, on average, 58% of the daily variations and 84% of the monthly amount of ET. Water amount from E_t contributed most to the ET, with the fractions of E_t, E_c and E_g being, on average, 67, 11 and 23% over the 10-year period, respectively. Regardless of weather conditions, the daily ET was strongly dependent on air temperature (T_a) and vapor pressure deficit (D_a), but less dependent on soil moisture (W_s). On cloudy and rainy days, there was a non-linear relationship between the ET and solar radiation (R_o). During drought years, the model predicted lower daily canopy stomatal conductance (g_cs) compared with non-drought years, leading to a lower level of E_t. The modeled daily g_cs responded well to D_a and W_s. In the model simulation, the annual LAI increased by 35% between 1999 and 2008. The ratio of E_c: ET correlated strongly with LAI. Furthermore, LAI reduced the proportion of E_g as a result of the increased share of E_c and E_t and radiation interception. Although the increase of LAI affected positively E_t, the contribution of E_t in ET was not significantly correlated with LAI. To conclude, although the model predicted reasonably well the seasonal course of ET, the calculation time steps of different processes in the model should be homogenized in the future to increase the prediction accuracy.     

Flows across high aspect ratio street canyons: Reynolds number independence revisited

The Reynolds number for flow in a street canyon, Re?=?U_refH/ν (where U_ref is a reference velocity, H the street canyon height, and ν the kinematic viscosity), cannot be matched between reduced-scale experiments and full-scale field measurements. This mismatch is often circumvented by satisfying the Re independence criterion, which states that above a critical Re (Re_c), the flow field remains invariant with Re. Re_c?=?11,000 is often adopted in reduced-scale experiments. In deep street canyons with height-to-width aspect ratio ≥?1.5, reduced-scale experiments have shown two recirculation vortices induced by the mean flows, but full-scale field measurements have observed only one vortex. We investigated this discrepancy by conducting water channel experiments with Re between 10⁴ and 10⁵ at three aspect ratios. The canyons with aspect ratio 1.0 have Re_c?=?11,000, the canyons with aspect ratio 1.5 have Re_c between 31,000 and 58,000, while the canyons with aspect ratio 2.0 have Re_c between 57,000 and 87,000. Therefore, the widely adopted Re_c?=?11,000 is not applicable for canyons with aspect ratio greater than 1.5. Our results also confirm that there is only one vortex in deep canyons at high Re. This single-vortex flow regime could change our fundamental understanding of deep canyons, which are often assumed to exhibit multiple-vortex flow regimes. Applications such as numerical model validation based on the multiple-vortex regime should be revisited. Our experimental data with Re up to 10⁵ could be used to validate numerical models at high Re.     

Identification of important factors for water vapor flux and CO2 exchange in a cropland

Water vapor flux and carbon dioxide (CO₂) exchange in croplands are crucial to water and carbon cycle research as well as to global warming evaluation. In this study, a standard three-layer feed-forward back propagation neural network technique associated with the Bayesian technique of automatic relevance determination (ARD) was employed to investigate water vapor and CO₂ exchange between the canopy of summer maize and atmosphere in responses to variations of environmental and physiological factors. These factors, namely the photosynthetically active radiation (PAR), air temperature (T), vapor pressure deficient (VPD), leaf-area index (LAI), soil water content in root zone (W), and friction velocity (U*), were used as inputs in neural network analysis. Results showed that PAR, VPD, T and LAI were the primary factors regulating both water vapor and CO₂ fluxes with VPD and W more critical to water vapor flux and PAR and T more crucial to CO₂ exchange. Furthermore, two time variables “day of the year (DOY)” and “time of the day (TOD)” could also improve the simulation results of neural network analysis. The important factors identified by the neural network technique used in this study were in the order of PAR > T > VPD > LAI > U* > TOD for water vapor flux and in the order of VPD > W > LAI > T > PAR > DOY for CO₂ exchange. This study suggests that neural network technique associated with ARD could be a useful tool for identifying important factors regulating water vapor and CO₂ fluxes in terrestrial ecosystem.     

Front velocity and structure of bottom gravity currents with a low volume of release propagating in a porous medium

The paper reports results of large eddy simulations of lock exchange compositional gravity currents with a low volume of release advancing in a horizontal, long channel. The channel contains an array of spanwise-oriented square cylinders. The cylinders are uniformly distributed within the whole channel. The flow past the individual cylinders is resolved by the numerical simulation. The paper discusses how the structure and evolution of the current change with the main geometrical parameters of the flow (e.g., solid volume fraction, ratio between the initial height of the region containing lock fluid and the channel depth, ratio between the initial length and height of the region containing lock fluid) and the Reynolds number. Though in all cases with a sufficiently large solid volume fraction the current transitions to a drag-dominated regime, the value of the power law coefficient, α, describing the front position’s variation with time (x _f  ~ t ^α , where t is the time measured from the removal of the lock gate) is different between full depth cases and partial depth cases. The paper also discusses how large eddy simulation (LES) results compare with findings based on shallow-water equations. In particular, LES results show that the values of α are not always equal to values predicted by shallow water theory for the limiting cases where the current height is comparable, or much smaller, than the channel depth.     

A peaked function for modeling temperature dependence of plant productivity

We suggests that temperature response of plant productivity can be modeled by the Arrhenius function modified to describe the effect of temperature on enzyme activity: G_A(T) = 2f(T)/(1 + f²(T)), where f(T) = exp(E_a/RT_opt − E_a/RT), R the universal gas constant, E_a the activation energy and T_opt is the optimal temperature. In common with other functions used for modeling the temperature response of plant productivity, the curve of function G is almost symmetrical and bell-shaped. The special convenience of G_A is that it relates the width of the “bell” to thermodynamic concepts, such as activation energy of chemical reactions converting carbon dioxide and water to carbohydrates.     

Parameter estimation for distributed delay based population models from laboratory data: egg hatching of Oulema duftschmidi Redthenbacher (Coleoptera,Chrysomelidae) as an example

The cohort development of poikilotherms under favorable temperature conditions may be described by the time-invariant distributed delay models that require three parameters, i.e. the thermal threshold T₀, the thermal constant W and the variability parameter H representing distributed maturation times, also known as time distributed or stochastic development. Here, a parsimonious method is developed that uses stage-frequency matrices obtained under controlled conditions. The analysis of these matrices permits the estimation of the three parameters, while the incorporation of both developmental threshold and thermal constant into the time-invariant distributed delay model permits the representation of stochastic cohort development under constant temperatures. The evaluated parameters can also be used in time-varying distributed delay models that are particularly useful for a wide range of fluctuating temperature conditions.We consider stage-frequency matrices obtained from cohorts of eggs and larvae of Oulema duftschmidi Redthenbacher at different constant temperatures. A visual inspection of the matrices clearly shows the temperature dependent as well as the time distributed passage from the egg to the larval stage. Under the current range of temperatures, a linear model for developmental rate satisfactorily describes egg development, and estimates the thermal threshold (T₀=11.2 °C), the thermal constant (W=81.3 day-degrees) and the delay order parameter (H=70). The resulting model can be used to represent the development of cohorts of O. duftschmidi eggs in a favorable temperature range.     

Theoretical analysis of change in forest carbon use efficiency with stand development: A case study on Hinoki Cypress (Chamaecyparis obtusa (Sieb. et Zucc.) Endl.) plantation from the seedling stage

Changes in carbon use efficiency (CUE), which is defined as the ratio of net primary production (NPP) to gross primary production (GPP), were estimated for the aerial parts of the Hinoki Cypress (Chamaecyparis obtusa (Sieb. et Zucc.) Endl.) with respect to stand development. The analysis incorporated previously published data from the early stages of stand development, namely the seedling stages of the cypress. For this analysis, a simple mathematical model to assess the changes in CUE was developed by incorporating data on physiological variables and mass of woody species. The CUE tended to increase with increases in the aboveground biomass of the stand, and then decreased gradually despite increases in the aboveground biomass. The CUE-value (0.28, 0.39) of the seedling stage was lower than that (0.33-0.58) of the young or mature trees. To examine the effect of physiological variables and mass on CUE, the ratios of the specific respiration rate to the specific photosynthetic rate (r/a) and the leaf biomass to the aboveground biomass or leaf mass ratio (y_L/y_T) were calculated. The low value of CUE at the seedling stage was due to the high ratio of specific respiration rate to specific photosynthetic rate r/a, but was not due to the high value of the leaf mass ratio y_L/y_T. In addition, the decline in CUE associated with older stages of stand development was due to the decreasing changes in y_L/y_T, and the r/a ratio did not influence the change in CUE.     

Estimation of some parameters of large fish populations by capture-recapture experiments

A certain number of tagged fish is liberated and assumed to be distributed randomly among a natural fish population. The fish are subjected to a number of fishing experiments within relatively short periods, and lie between equal intervening periods of durationT. Untagged fish are retained, while tagged fish are released during the fishing experiments. Denoting the catchability of untagged fish byq _u and that for tagged fish byq _t , it is assumed that they are related by the equation “q _u =cq _t ” wherec is a constant. Denoting the survival rates of tagged fish and the effective fishing effort of commerical fisheries per unit time from the (k-1)^th to thek ^th experiments by _t S _k andf _k , respectively, it is assumed that they vary from period to period. Assuming that during thek ^th experiment, the number of untagged fish captured and the experimental fishing rate of tagged fish are denoted by _u X _k and _t P _k , respectively, then $$\begin{array}{*{20}c} {\frac{{(_u X_k )^2 }}{{[_u X_{(k - 1)} ][_u X_{(k + 1)} ]}} = \frac{{_t S_k }}{{_t S_{(k + 1)} }} \cdot \frac{{e^{ - (1 - c)q_t f_{(k + 1)} T} }}{{e^{ - (1 - c)q_t f_k T} }} \cdot } \\ {\frac{{1 - c_t P_{(k - 1)} }}{{1 - c_t P_k }} \cdot \frac{{(_t P_k )^2 }}{{[_t P_{(k - 1)} ][_t P_{(k + 1)} ]}}.} \\ \end{array}$$ The above equation containsc as a single unknown, while all other terms are supplied by the capture-recapture experiments, exceptf _k andf _(k+1) which may be obtained from fisheries statistics. A number of the above equations are obtained from several experiments and can be combined into a single equation to obtain an overall estimate forc which can be used to derive estimates for experimental fishing rates, abundance, and instantaneous natural and fishing mortality rates for natural fish populations. These estimates are free from type (A) tagging errors, and have the advantage of taking into consideration the probable different behaviour of tagged and untagged fish.     

Effects of ocean acidification on the early developmental stages of the horned turban,Turbo cornutus

To estimate the impact of CO₂-driven ocean acidification on the early life stages of gastropods, the effects of increased partial pressure of seawater carbon dioxide (pCO₂) (800–2,000 μatm) on the early developmental stages and larval shell length of the commercially important gastropod, the horned turban snail, Turbo cornutus were investigated. Increase in experimental seawater pCO₂ had an increasingly negative impact on the early developmental rate; the proportion of embryos or larvae displaying retarded development increased at higher pCO₂. The proportion of embryos that developed to the 4-cell stage at 2 h after fertilization decreased linearly with increasing pCO₂. At ~1,000 μatm pCO₂, retarded development was observed in ~50 % of larvae. No embryos developed to the 4-cell stage at 2,000 μatm pCO₂ within 2 h of fertilization. A similar trend continued until 24–26 h after fertilization; the proportion of larvae attaining veliger stage by 24–26 h also decreased with increasing pCO₂. The shell length of T. cornutus veligers decreased gradually as seawater pCO₂ increased, but markedly decreased in seawater under nearly unsaturated and unsaturated conditions (≤1.04) of the aragonite saturation state (Ω _aragonite). The results indicate that increased pCO₂ seawater has a progressive and acute effect on embryonic and larval T. cornutus, and imply that the extended early developmental period and/or the downsized larval shell produced by ocean acidification will have a negative impact on survival, settlement and recruitment well into the future.     

Calculation and use of selectivity coefficients of feeding: Zooplankton grazing

A straightforward method for calculating selectivity coefficients (W_ij) of predation from raw data, mortality rates of prey, filtering rates, feeding rates and electivity indices is derived. Results from a comparison of selectivity coefficients for the copepod Diaptomus oregonensis grazing under a number of experimental conditions suggest that W_ij's for size-selective feeding are invariant, a conclusion also supported by the leaky-sieve model. Recommendations are made on how to use W_ij's in linear and nonlinear feeding constructs for zooplankton and other animals.     

Effects of changing climate on water and nitrogen availability with implications on the productivity of Norway spruce stands in Southern Finland

An integrated process-based model was used to study how the changing climate affects the availability of water and nitrogen, and consequently the dynamics of productivity of Norway spruce (Picea abies) on sites with different initial soil water conditions in southern Finland over a 100-year period. The sensitivity of the total stem volume growth in relation to short-term availability of water and nitrogen was also analyzed. We found that a high proportion (about 88–92%) of the total precipitation was lost in total evapotranspiration (incl. canopy evaporation (E_c), transpiration (E_t) and ground surface evaporation (E_g)), under both current and changing climate. Furthermore, under the changing climate the cumulative amount of E_c and E_g were significantly higher, while E_t was largely lower than under the current climate. Additionally, the elevated temperature and increased expansion of needle area index (L) enhanced E_c. Under the changing climate, the increasing soil water deficit (W_d) reduced the canopy stomatal conductance (g_cs), the E_t, humus yield (H, available nitrogen source) and nitrogen uptake (N_up) of the trees. During the latter phases of the simulation period, the canopy net photosynthesis (P_nc) was lower due to the reduced N_up and soil water availability. This also reduced the total stem volume production (V_s) on the site with the lower initial soil moisture content. The growth was slightly more sensitive to the change in precipitation than to the change in nitrogen content of the needles, when the elevated temperature was assumed. According to our findings, drought stress episodes may become more frequent under the changing climate. Thus, adaptive management strategies should be developed to sustain the productivity of Norway spruce in these conditions, and thus, to mitigate the adverse impacts of climate change.     

Experimental study on acute effects of the combined exposure to temperature increase and chlorination upon the marine copepod Acartia omorii

The acute effects of combined exposure to temperature increase and chlorination on a neritic marine copepod, Acartia omorii Bradford (collected offshore of Onjuku, Japan in 1982), were investigated in the laboratory. Continuous flow exposure and batch exposure modes were compared. Based on the results of continuous flow experiments, the 24-h median lethal concentration (24-h LC₅₀, in mg l^-1) of total residual chlorine was estimated using the multiple regression equation below, with a multiple correlation coefficient of 0.955: 24-h LC₅₀=2.988-0.034 dT-1.611 log₁₀ t where dT is temperature rise (°C) and t is exposure duration (min). In batch experiments, the predictive power of the multiple regression equation was reduced, probably due to variations in chlorine concentration during exposure duration.     

A new growth model for fishes and the estimation of optimum age of fish populations

Evidence is presented that a simple power equation of the formX _t –X _m =±b(|tࢤt _m|)^B can describe growth in length and weight of fishes, whereX _t denotes fish length or weight at aget, X _m is length or weight (L _m orW _m ) at a reference aget _m , andb andB are parameters to be estimated by the least squares. The optimum age of fish populations (t _y ) may be estimated by the equationMW _m /b=±y ^B (B/y-M) whereM denotes natural mortality and wherey=t _y -t _m .     

Energy conservation in hibernating endotherms: Why “suboptimal” temperatures are optimal

Many endotherms use facultative heterothermic responses of torpor or hibernation to conserve energy during periods of low energy availability. A common assumption when estimating winter energy budgets is that endotherms should hibernate at the ambient temperature (T_a) that minimizes torpid metabolic rate (TMR) and maximizes the duration of torpor bouts. However, previous studies of the energetic benefits of hibernation have assumed constant T_a within hibernacula. Here we use an individual-based energetic model to estimate overwinter energy expenditure of mammals hibernating at T_as that vary temporally. We show that, in accordance with the principles of Jenson's inequality, hibernators can conserve energy by selecting microclimates warmer than the single T_a value that minimizes TMR (T_min). As temporal variation in T_a increases, endotherms should choose microclimates with mean T_as progressively warmer than T_min. Further, as thermal conductance decreases, as it does with increasing body mass and use of social thermoregulation, the mean T_a that minimizes overwinter energy expenditure approaches, but never equals, T_min. We suggest that the commonly held assumption of stable microclimates in hibernacula has skewed the interpretation of the optimal expression of hibernation for energy conservation. Our results contradict much of the accepted understanding of hibernation energetics and add to a growing body of literature proposing that hibernating at a T_a warmer than T_min is optimal.     

Community metabolism of adenylates by microheterotrophs from the Los Angeles Harbor and Southern California coastal waters

The natural concentration (S _n) of dissolved total adenylates TA(=AMP+ADP+ATP) in coastal seawater from a depth of 1 m at 5 stations (California, USA) sampled periodically for 1 yr had a mean value ±1 SD of 2.8±1.7 nmol TA1^-1. The specific uptake rates of TA by microheterotrophs at a station inside the Los Angeles Harbor and at a station 1.5 km offshore in the San Pedro Channel were studied by simple uptake and saturation-type kinetic analysis using ³H-AMP as a tracer. Within the harbor, the specific uptake rate (nmol TA 10⁹ cell^-1 h^-1) at S _n ranged 10-fold from 0.028 in December to 0.28 in August. K _t (half-saturation constant) values always exceeded theS _n concentrations in any given month, and were greater in the harbor than in the channel. Generally, over 80% of biological uptake of ³H-AMP was associated with organisms <1.0 m, a size class accounting for about 20% of the total particulate adenylate concentration in the 0.2 to 203 m size fraction. Assuming steady-state conditions for the dissolved adenylate pool, we propose a model in which losses from this pool are balanced by inputs to the pool through inefficient feeding, lysis and decomposition of particulate adenylates.     

Fitting a parabola to growth data of fishes and some applications to fisheries

Denoting a fish length or weight at age t by X _t , a reference age by t _m , and the corresponding fish length or weight by X _m , the relation between age and length or weight may be described by a parabola as follows: $$\left| {X_t } \right. - X_m \left| = \right.a + b(\left| {t - t_m } \right.\left| ) \right. + c(\left| t \right. - t_m \left| ) \right.^2$$ or $$X_t = A + b(\left| {t - t_m } \right.\left| ) \right. + c(\left| t \right. - t_m \left| ) \right.^2$$ where a, b and c are constants. Each of the above Eqs. describes one curve at ages older than t _m and another one at younger ages, which is made possible by means of the transformation of t to (|t-t _m |). In 2 cases out of 10, the parabola takes the form of a cubic equation. Evidence is given that, as the growth data become fewer, the better fit of the parabola or cubic equation will probably be less in comparison to the von Bertalanffy equation (1938, 1949) as developed by Beverton and Holt (1957) and the power-growth equation (Rafail, 1971), and vice versa. This growth equation is used to derive models for estimating the optimum age and yield for fish populations.     

Diffusional mass transfer coefficient at the water–sediment interface for wind-induced flow in very shallow lagoons

Very shallow lagoons that are a few centimeters deep are common in the arid Andes of Northern Chile, Argentina, Bolivia and Perú. The dynamics of these lagoons are dominated by the water–sediment interface (WSI) and strong afternoon winds. Although many studies have examined the diffusional mass transfer coefficients (k _t ) of open channel flows, estimates for wind-induced flows are still unknown. The aim of this article is to propose and validate an analytical expression for computing k _t at the WSI for wind-induced flow. The laboratory measurements were conducted in a wind tunnel with a water tank of variable depth located at its downwind end. Natural muddy sediments were placed in the middle of the tank so that the dissolved oxygen (DO) was consumed in the sediments. The diffusional mass transfer coefficient that characterizes the DO uptake in the sediment was obtained from DO micro-profiles measured with an OX-25 Unisense microelectrode. Water velocity profiles were measured with a 2D side-view Sontek acoustic doppler velocimetry (ADV), and the wind shear velocity was computed based on wind velocity profiles that were measured with an Extech hot-wire anemometer. A total of 16 experiments were conducted with different water depths and wind shear stresses. The constants required by the model were determined from these experiments, and the analytical expression was successfully validated by the laboratory observations. The analytical expression obtained for computing k _t was also validated with field observations that were conducted in October, 2012, in Salar del Huasco, Northern Chile (20.274° S, 68.883° W, 3800 m above sea level). The comparison between the observed and predicted values of k _t provides a determination coefficient of r ² = 0.48 and a p value < 0.01. The results show that the value of k _t for wind-induced flow is proportional to the wind shear velocity and the inverse of the Reynolds number of the wind-induced current.