Modeling the effect of stand and site characteristics on the probability of mistletoe infestation in Scots pine stands using remote sensing data

Over the past decade,the presence of mistletoe(Viscum album ssp.austriacum)in Scots pine stands has increased in many European countries.Understanding the factors that influence the occurrence of mistletoe in stands is key to making appropriate forest management decisions to limit damage and prevent the spread of mistletoe in the future.Therefore,the main objective of this study was to determine the probability of mistletoe occurrence in Scots pine stands in relation to stand-related endogenous factors such as age,top height,and stand density,as well as topographic and edaphic factors.We used unmanned aerial vehicle(UAV)imagery from 2,247 stands to detect mistletoe in Scots pine stands,while majority stand and site characteristics were calculated from airborne laser scanning(ALS)data.Information on stand age and site type from the State Forest database were also used.We found that mistletoe infestation in Scots pine stands is influenced by stand and site characteristics.We documented that the densest,tallest,and oldest stands were more susceptible to mistletoe infestation.Site type and specific microsite conditions associated with topography were also important factors driving mistletoe occurrence.In addition,climatic water balance was a significant factor in increasing the probability of mistletoe occurrence,which is important in the context of predicted temperature increases associated with climate change.Our results are important for better understanding patterns of mistletoe infestation and ecosystem functioning under climate change.In an era of climate change and technological development,the use of remote sensing methods to determine the risk of mistletoe infestation can be a very useful tool for managing forest ecosystems to maintain forest sustainability and prevent forest disturbance.

Predictive Vegetation Mapping Approach Based on Spectral Data, DEM and Generalized Additive Models

This study aims to provide a predictive vegetation mapping approach based on the spectral data, DEM and Generalized Additive Models (GAMs). GAMs were used as a prediction tool to describe the relationship between vegetation and environmental variables, as well as spectral variables. Based on the fitted GAMs model, probability map of species occurrence was generated and then vegetation type of each grid was defined according to the probability of species occurrence. Deviance analysis was employed to test the goodness of curve fitting and drop contribution calculation was used to evaluate the contribution of each predictor in the fitted GAMs models. Area under curve (AUC) of Receiver Operating Characteristic (ROC) curve was employed to assess the results maps of probability. The results showed that: 1) AUC values of the fitted GAMs models are very high which proves that integrating spectral data and environmental variables based on the GAMs is a feasible way to map the vegetation. 2) Prediction accuracy varies with plant community, and community with dense cover is better predicted than sparse plant community. 3) Both spectral variables and environmental variables play an important role in mapping the vegetation. However, the contribution of the same predictor in the GAMs models for different plant communities is different. 4) Insufficient resolution of spectral data, environmental data and confounding effects of land use and other variables which are not closely related to the environmental conditions are the major causes of imprecision.

Modeling hot strip rolling process under framework of generalized additive model

This research develops a new mathematical modeling method by combining industrial big data and process mechanism analysis under the framework of generalized additive models(GAM)to generate a practical model with generalization and precision.Specifically,the proposed modeling method includes the following steps.Firstly,the influence factors are screened using mechanism knowledge and data-mining methods.Secondly,the unary GAM without interactions including cleaning the data,building the sub-models,and verifying the sub-models.Subsequently,the interactions between the various factors are explored,and the binary GAM with interactions is constructed.The relationships among the sub-models are analyzed,and the integrated model is built.Finally,based on the proposed modeling method,two prediction models of mechanical property and deformation resistance for hot-rolled strips are established.Industrial actual data verification demonstrates that the new models have good prediction precision,and the mean absolute percentage errors of tensile strength,yield strength and deformation resistance are 2.54%,3.34%and 6.53%,respectively.And experimental results suggest that the proposed method offers a new approach to industrial process modeling.

Modeling and Fault Monitoring of Bioprocess Using Generalized Additive Models （GAMs） and Bootstrap

Fault monitoring of bioprocess is important to ensure safety of a reactor and maintain high quality of products. It is difficult to build an accurate mechanistic model for a bioprocess, so fault monitoring based on rich historical or online database is an effective way. A group of data based on bootstrap method could be resampling stochastically, improving generalization capability of model. In this paper, online fault monitoring of generalized additive models （GAMs） combining with bootstrap is proposed for glutamate fermentation process. GAMs and bootstrap are first used to decide confidence interval based on the online and off-line normal sampled data from glutamate fermentation experiments. Then GAMs are used to online fault monitoring for time, dissolved oxygen, oxygen uptake rate, and carbon dioxide evolution rate. The method can provide accurate fault alarm online and is helpful to provide useful information for removing fault and abnormal phenomena in the fermentation.

Simulating Potential Distribution of Tamarix chinensis in Yellow River Delta by Generalized Additive Models

There are typical ecosystems of littoral wetlands in the Yellow River Delta.In order to study the relationships between Tamarix chinensis and environmental variables and to predict T.chinensis potential distribution in the Yellow River Delta,641 vegetation samples and 964 soil samples were collected in the area in October of 2004,2005,2006 and 2007.The contents of soil organic matter,total phosphorus,salt,and soluble potassium were determined.Then,the analyzed data were interpolated into spatial raster data by Kriging interpolation method.Meanwhile,the digital elevation model,soil type map and landform unit map of the Yellow River Delta were also collected.Generalized Additive Models(GAMs) were employed to build species-environment model and then simulate the potential distribution of T.chinensis.The results indicated that the distribution of T.chinensis was mainly limited by soil salt content,total soil phosphorus content,soluble potassium content,soil type,landform unit,and elevation.The distribution probability of T.chinensis was produced with a lookup table generated by Grasp Module(based on GAMs) in software ArcView GIS 3.2.The AUC(Area Under Curve) value of validation and cross-validation of ROC(Receive Operating Characteristic) were both higher than 0.8,which suggested that the established model had a high precision for predicting species distribution.

Inference Procedures on the Generalized Poisson Distribution from Multiple Samples: Comparisons with Nonparametric Models for Analysis of Covariance (ANCOVA) of Count Data

Count data that exhibit over dispersion (variance of counts is larger than its mean) are commonly analyzed using discrete distributions such as negative binomial, Poisson inverse Gaussian and other models. The Poisson is characterized by the equality of mean and variance whereas the Negative Binomial and the Poisson inverse Gaussian have variance larger than the mean and therefore are more appropriate to model over-dispersed count data. As an alternative to these two models, we shall use the generalized Poisson distribution for group comparisons in the presence of multiple covariates. This problem is known as the ANCOVA and is solved for continuous data. Our objectives were to develop ANCOVA using the generalized Poisson distribution, and compare its goodness of fit to that of the nonparametric Generalized Additive Models. We used real life data to show that the model performs quite satisfactorily when compared to the nonparametric Generalized Additive Models.

Fitting Generalized Additive Logistic Regression Model with GAM Procedure

In dealing with nonparametric regression the GAM procedure is the most versatile of several new procedures. The terminology behind this procedure is more flexible than traditional parametric modeling tools. It relaxes the usual assumptions of parametric model and enables us to uncover structure to establish the relationship between independent variables and dependent variable in exponential family that may not be obvious otherwise. In this paper, we discussed two methods of fitting generalized additive logistic regression model, one based on Newton Raphson method and another based on iterative weighted least square method for first and second order Taylor series expansion. The use of the GAM procedure with the specified set of weights, using local scoring algorithm, was applied to real life data sets. The cubic spline smoother is applied to the independent variables. Based on nonparametric regression and smoothing techniques, this procedure provides powerful tools for data analysis.

From furnace up to freezer:Elevational patterns of plant diversity in Mount Palvar,a semi-arid Irano-Turanian mountain range of southwest Asia

Much of the world's biodiversity lies in heterogeneous mountain areas with their diverse environments.As an example,Iranian montane ranges are highly diverse,particularly in the Irano-Turanian phytogeographical region.Understanding plant diversity patterns with increasing elevation is of high significance,not least for conservation planning.We studied the pattern of species richness,Shannon diversity,endemic richness,endemics ratio,and richness of life forms along a 3900 m elevational transect in Mount Palvar,overlooking the Lut Desert in Southeast Iran.We also analyzed the effect of environmental variables on species turnover along the vertical gradient.A total of 120 vegetation plots(10 m×10 m)were sampled along the elevational transect containing species and environmental data.To discover plant diversity pattern along the elevational gradient,generalized additive model(GAM)was used.Non-metric multidimensional scaling(NMDS)was applied for illustrating the correlation between species composition and environmental variables.We found hump-shaped pattern for species richness,Shannon diversity,endemic richness,and species richness of different life forms,but a monotonic increasing pattern for ratio of endemic species from low to high elevations.Our study confirms the humped pattern of species richness peaking at intermediate elevations along a complete elevational gradient in a semi-arid mountain.The monotonic increase of endemics ratio with elevation in our area as a case study is consistent with global increase of endemism with elevation.According to our results,temperature and precipitation are two important climatic variables that drive elevational plant diversity,particularly in seasonally dry areas.Our study suggests that effective conservation and management are needed for this low latitude mountain area along with calling for long-term monitoring for species redistribution.

Grouping tree species to estimate basal area increment in temperate multispecies forests in Durango,Mexico

Multispecies forests have received increased scientific attention,driven by the hypothesis that biodiversity improves ecological resilience.However,a greater species diversity presents challenges for forest management and research.Our study aims to develop basal area growth models for tree species cohorts.The analysis is based on a dataset of 423 permanent plots(2,500 m^(2))located in temperate forests in Durango,Mexico.First,we define tree species cohorts based on individual and neighborhood-based variables using a combination of principal component and cluster analyses.Then,we estimate the basal area increment of each cohort through the generalized additive model to describe the effect of tree size,competition,stand density and site quality.The principal component and cluster analyses assign a total of 37 tree species to eight cohorts that differed primarily with regard to the distribution of tree size and vertical position within the community.The generalized additive models provide satisfactory estimates of tree growth for the species cohorts,explaining between 19 and 53 percent of the total variation of basal area increment,and highlight the following results:i)most cohorts show a"rise-and-fall"effect of tree size on tree growth;ii)surprisingly,the competition index"basal area of larger trees"had showed a positive effect in four of the eight cohorts;iii)stand density had a negative effect on basal area increment,though the effect was minor in medium-and high-density stands,and iv)basal area growth was positively correlated with site quality except for an oak cohort.The developed species cohorts and growth models provide insight into their particular ecological features and growth patterns that may support the development of sustainable management strategies for temperate multispecies forests.

Spatio-temporal variation of depth to groundwater level and its driving factors in arid and semi-arid regions of India

Climate change and increasing anthropogenic activities,such as over-exploitation of groundwater,are exerting unavoidable stress on groundwater resources.This study investigated the spatio-temporal variation of depth to groundwater level(DGWL)and the impacts of climatic(precipitation,maximum temperature,and minimum temperature)and anthropogenic(gross district product(GDP),population,and net irrigated area(NIA))variables on DGWL during 1994-2020.The study considered DGWL in 113 observation wells and piezometers located in arid western plains(Barmer and Jodhpur districts)and semi-arid eastern plains(Jaipur,Ajmer,Dausa,and Tonk districts)of Rajasthan State,India.Statistical methods were employed to examine the annual and seasonal patterns of DGWL,and the generalized additive model(GAM)was used to determine the impacts of climatic and anthropogenic variables on DGWL.During 1994-2020,except for Barmer District,where the mean annual DGWL was almost constant(around 26.50 m),all other districts exhibited increase in DGWL,with Ajmer District experiencing the most increase.The results also revealed that 36 observation wells and piezometers showed a statistically significant annual increasing trend in DGWL and 34 observation wells and piezometers exhibited a statistically significant decreasing trend in DGWL.Similarly,32 observation wells and piezometers showed an statistically significant increasing trend and 37 observation wells and piezometers showed a statistically significant decreasing trend in winter;33 observation wells and piezometers indicated a statistically significant increasing trend and 34 had a statistically significant decreasing trend in post-monsoon;35 observation wells and piezometers exhibited a statistically significant increasing trend and 32 observation wells and piezometers showed a statistically significant decreasing trend in pre-monsoon;and 36 observation wells and piezometers reflected a statistically significant increasing trend and 30 observation wells and piezometers reflected a statistically significant decreasing trend in monsoon.Interestingly,most of the observation wells and piezometers with increasing trends of DGWL were located in Dausa and Jaipur districts.Furthermore,the GAM analysis revealed that climatic variables,such as precipitation,significantly affected DGWL in Barmer District,and DGWL in all other districts was influenced by anthropogenic variables,including GDP,NIA,and population.As a result,stringent regulations should be implemented to curb excessive groundwater extraction,manage agricultural water demand,initiate proactive aquifer recharge programs,and strengthen sustainable management in these water-scarce regions.

Optimization of environmental variables in habitat suitability modeling for mantis shrimp Oratosquilla oratoria in the Haizhou Bay and adjacent waters

Habitat suitability index(HSI)models have been widely used to analyze the relationship between species abundance and environmental factors,and ultimately inform management of marine species.The response of species abundance to each environmental variable is different and habitat requirements may change over life history stages and seasons.Therefore,it is necessary to determine the optimal combination of environmental variables in HSI modelling.In this study,generalized additive models(GAMs)were used to determine which environmental variables to be included in the HSI models.Significant variables were retained and weighted in the HSI model according to their relative contribution(%)to the total deviation explained by the boosted regression tree(BRT).The HSI models were applied to evaluate the habitat suitability of mantis shrimp Oratosquilla oratoria in the Haizhou Bay and adjacent areas in 2011 and 2013–2017.Ontogenetic and seasonal variations in HSI models of mantis shrimp were also examined.Among the four models(non-optimized model,BRT informed HSI model,GAM informed HSI model,and both BRT and GAM informed HSI model),both BRT and GAM informed HSI model showed the best performance.Four environmental variables(bottom temperature,depth,distance offshore and sediment type)were selected in the HSI models for four groups(spring-juvenile,spring-adult,falljuvenile and fall-adult)of mantis shrimp.The distribution of habitat suitability showed similar patterns between juveniles and adults,but obvious seasonal variations were observed.This study suggests that the process of optimizing environmental variables in HSI models improves the performance of HSI models,and this optimization strategy could be extended to other marine organisms to enhance the understanding of the habitat suitability of target species.

Examining spatiotemporal distribution and CPUE-environment relationships for the jumbo flying squid Dosidicus gigas offshore Peru based on spatial autoregressive model

The spatiotemporal distribution and relationship between nominal catch-per-unit-ef fort(CPUE) and environment for the jumbo flying squid( Dosidicus gigas) were examined in of fshore Peruvian waters during 2009–2013. Three typical oceanographic factors aff ecting the squid habitat were investigated in this research, including sea surface temperature(SST), sea surface salinity(SSS) and sea surface height(SSH). We studied the CPUE-environment relationships for D. gigas using a spatially-lagged version of spatial autoregressive(SAR) model and a generalized additive model(GAM), with the latter for auxiliary and comparative purposes. The annual fishery centroids were distributed broadly in an area bounded by 79.5°–82.7°W and 11.9°–17.1°S, while the monthly fishery centroids were spatially close and lay in a smaller area bounded by 81.0°–81.2°W and 14.3°–15.4°S. Our results show that the preferred environmental ranges for D. gigas offshore Peru were 20.9°–21.9°C for SST, 35.16–35.32 for SSS and 27.2–31.5 cm for SSH in the areas bounded by 78°–80°W/82–84°W and 15°–18°S. Monthly spatial distributions during October to December were predicted using the calibrated GAM and SAR models and general similarities were found between the observed and predicted patterns for the nominal CPUE of D. gigas. The overall accuracies for the hotspots generated by the SAR model were much higher than those produced by the GAM model for all three months. Our results contribute to a better understanding of the spatiotemporal distributions of D. gigas off shore Peru, and off er a new SAR modeling method for advancing fishery science.

Linearity extensions of the market model:a case of the top 10 cryptocurrency prices during the pre‑COVID‑19 and COVID‑19 periods

This research investigates the appropriateness of the linear specification of the market model for modeling and forecasting the cryptocurrency prices during the pre-COVID-19 and COVID-19 periods.Two extensions are offered to compare the performance of the linear specification of the market model(LMM),which allows for the measurement of the cryptocurrency price beta risk.The first is the generalized additive model,which permits flexibility in the rigid shape of the linearity of the LMM.The second is the time-varying linearity specification of the LMM(Tv-LMM),which is based on the state space model form via the Kalman filter,allowing for the measurement of the time-varying beta risk of the cryptocurrency price.The analysis is performed using daily data from both time periods on the top 10 cryptocurrencies by adjusted market capitalization,using the Crypto Currency Index 30(CCI30)as a market proxy and 1-day and 7-day forward predictions.Such a comparison of cryptocurrency prices has yet to be undertaken in the literature.The empirical findings favor the Tv-LMM,which outperforms the others in terms of modeling and forecasting performance.This result suggests that the relationship between each cryptocurrency price and the CCI30 index should be locally instead of globally linear,especially during the COVID-19 period.

Model-based estimation of above-ground biomass in the miombo ecoregion of Zambia

Background：Information on above-ground biomass（AGB） is important for managing forest resource use at local levels,land management planning at regional levels,and carbon emissions reporting at national and international levels.In many tropical developing countries,this information may be unreliable or at a scale too coarse for use at local levels.There is a vital need to provide estimates of AGB with quantifiable uncertainty that can facilitate land use management and policy development improvements.Model-based methods provide an efficient framework to estimate AGB.Methods：Using National Forest Inventory（NFI） data for a^1,000,000 ha study area in the miombo ecoregion,Zambia,we estimated AGB using predicted canopy cover,environmental data,disturbance data,and Landsat 8 OLI satellite imagery.We assessed different combinations of these datasets using three models,a semiparametric generalized additive model（GAM） and two nonlinear models（sigmoidal and exponential）,employing a genetic algorithm for variable selection that minimized root mean square prediction error（RMSPE）,calculated through cross-validation.We compared model fit statistics to a null model as a baseline estimation method.Using bootstrap resampling methods,we calculated 95% confidence intervals for each model and compared results to a simple estimate of mean AGB from the NFI ground plot data.Results：Canopy cover,soil moisture,and vegetation indices were consistently selected as predictor variables.The sigmoidal model and the GAM performed similarly;for both models the RMSPE was -36.8 tonnes per hectare（i.e.,57% of the mean）.However,the sigmoidal model was approximately 30% more efficient than the GAM,assessed using bootstrapped variance estimates relative to a null model.After selecting the sigmoidal model,we estimated total AGB for the study area at 64,526,209 tonnes（＋/- 477,730）,with a confidence interval 20 times more precise than a simple designbased estimate.Conclusions：Our findings demonstrate that NFI data may be combined with freely available satellite imagery and soils data to estimate total AGB with quantifiable uncertainty,while also providing spatially explicit AGB maps useful for management,planning,and reporting purposes.

A spatially-explicit count data regression for modeling the density of forest cockchafer(Melolontha hippocastani) larvae in the Hessian Ried(Germany)

Background： In this paper, a regression model for predicting the spatial distribution of forest cockchafer larvae in the Hessian Ried region （Germany） is presented. The forest cockchafer, a native biotic pest, is a major cause of damage in forests in this region particularly during the regeneration phase. The model developed in this study is based on a systematic sample inventory of forest cockchafer larvae by excavation across the Hessian Ried. These forest cockchafer larvae data were characterized by excess zeros and overdispersion. Methods： Using specific generalized additive regression models, different discrete distributions, including the Poisson, negative binomial and zero-inflated Poisson distributions, were compared. The methodology employed allowed the simultaneous estimation of non-linear model effects of causal covariates and, to account for spatial autocorrelation, of a 2-dimensional spatial trend function. In the validation of the models, both the Akaike information criterion （AIC） and more detailed graphical procedures based on randomized quantile residuals were used. Results： The negative binomial distribution was superior to the Poisson and the zero-inflated Poisson distributions, providing a near perfect fit to the data, which was proven in an extensive validation process. The causal predictors found to affect the density of larvae significantly were distance to water table and percentage of pure clay layer in the soil to a depth of I m. Model predictions showed that larva density increased with an increase in distance to the water table up to almost 4 m, after which it remained constant, and with a reduction in the percentage of pure clay layer. However this latter correlation was weak and requires further investigation. The 2-dimensional trend function indicated a strong spatial effect, and thus explained by far the highest proportion of variation in larva density. Conclusions： As such the model can be used to support forest practitioners in their decision making for regeneration and forest protection planning in the Hessian predicting future spatial patterns of the larva density is still comparatively weak. Ried. However, the application of the model for somewhat limited because the causal effects are

Standardizing CPUE of Ommastrephes bartramii for Chinese squid-jigging fishery in Northwest Pacific Ocean

Generalized linear models （GLM） and generalized additive models （GAM） were used to standardize catch per unit fishing effort （CPUE） of Ommastrephes bartramii for Chinese squid-jigging fishery in Northwest Pacific Ocean. Three groups of variables were considered in the standardization： spatial variables （longitude and latitude）, temporal variables （year and month） and environmental variables, including sea surface temperature （SST）, sea surface salinity （SSS） and sea level height （SLH）. CPUE was treated as the dependent variable and its error distribution was assumed to be log-normal in each model. The model selections of GLM and GAM were based on the finite sample-corrected Akaike information criterion （AICC） and pseudo-coefficient （Pcf） combined P-value, respectively. Both GAM and GLM analysis showed that the month was the most important variable affecting CPUE and could explain 21.3% of variability in CPUE while other variables only explained 8.66%. The interaction of spatial and temporal variables weakly influenced the CPUE. Moreover, spatio-temporal factors may be more important in influencing the CPUE of this squid than environmental variables. The standardized and nominal CPUEs were similar and had the same trends in spatio-temporal distribution, but the standardized CPUE values tended to be smaller than the nominal CPUE. The CPUE tended to have much higher monthly variation than annual variations and their values increased with month. The CPUE became higher with increasing latitude-high CPUE usually occurred in 145°E-148°E and 149°E-162°E. The CPUE was higher when SST was 14-21℃ and the SLH from -22 cm to -18 cm. In this study, GAM tended to be more suitable than GLM in analysis of CPUE.

The study on fishing ground of neon flying squid, Ommastrephes bartrami,and ocean environment based on remote sensing data in the Northwest Pacific Ocean

The relationships between the neon flying squid, Ommastrephes bartrami, and the relative ocean environmental factors are analyzed. The environmental factors collected are sea surface temperature （SST）, chlorophyll concentration （Chl-α） and sea surface height （SSH） from NASA, as well as the yields of neon flying squid in the North Pacific Ocean. The results show that the favorable temperature for neon flying squid living is 10℃-22℃ and the favorite temperature is between 15℃-17℃. The Chl-α concentration is 0.1-0.6 mg/m^3. When Chl-α concentration changes to 0.12-0.14 mg/m^3, the probability of forming fishing ground becomes very high. In most fishing grounds, the SSH is higher than the mean SSH. The generalized additive model （GAM） was applied to analyze the correlations between neon flying squid and ocean environmental factors. Every year, squids migrate northward from June to August and return southward during October-November, and the characteristics of the both migrations are very different. When squids migrate to the north, most relationships between the yields and SST are positive. The relationships are negative when squids move to southward. The relationships between the yields and Chl-a concentrations are negative from June to October, and insignificant in November. There is no obvious correlation between the catches of squid and longitude, but good with latitude.

Growth and form of Quercus robur and Fraxinus excelsior respond distinctly different to initial growing space: results from 24-year-old Nelder experiments

Initial growing space is of critical importance to growth and quality development of individual trees. We investigated how mortality, growth （diameter at breast height, total height）, natural pruning （height to first dead and first live branch and branchiness） and stem and crown form of 24-year-old pedunculate oak （Quercus robur [L.]） and European ash （Fraxinus excelsior [L.]） were affected by initial spacing. Data were recorded from two replicate single-species Nelder wheels located in southern Germany with eight initial stocking regimes varying from 1,020 to 30,780 seedlings·ha？1. Mortality substantially decreased with increasing initial growing space but significantly differed among the two species, averaging 59% and 15% for oak and ash plots, respectively. In contrast to oak, the low self-thinning rate found in the ash plots over the investigated study period resulted in a high number of smaller intermediate or suppressed trees, eventually retarding individual tree as well as overall stand development. As a result, oak gained greater stem dimensions throughout all initial spacing regimes and the average height of ash significantly increased with initial growing space. The survival of lower crown class ashes also appeared to accelerate self-pruning dynamics. In comparison to oak, we observed less dead and live primary branches as well as a smaller number of epicormic shoots along the first 6 m of the lower stem of dominant and co-dominant ashes in all spacing regimes. Whereas stem form of both species was hardly affected by initial growing space, the percentage of brushy crowns significantly increased with initial spacing in oak and ash. Our findings suggest that initial stockings of ca. 12,000 seedlings per hectare in oak and 2,500 seedlings per hectare in ash will guarantee a sufficient number of at least 300 potential crop trees per hectare in pure oak and ash plantations at the end of the self-thinning phase, respectively. If the problem of epicormic shoots and inadequate self-pruning can be controlled with trainer species, the initial stocking may be reduced significantly in oak.

Probabilistic Precipitation Forecasting Based on Ensemble Output Using Generalized Additive Models and Bayesian Model Averaging

A probabilistic precipitation forecasting model using generalized additive models （GAMs） and Bayesian model averaging （BMA） was proposed in this paper. GAMs were used to fit the spatial-temporal precipitation models to individual ensemble member forecasts. The distributions of the precipitation occurrence and the cumulative precipitation amount were represented simultaneously by a single Tweedie distribution. BMA was then used as a post-processing method to combine the individual models to form a more skillful probabilistic forecasting model. The mixing weights were estimated using the expectation-maximization algorithm. The residual diagnostics was used to examine if the fitted BMA forecasting model had fully captured the spatial and temporal variations of precipitation. The proposed method was applied to daily observations at the Yishusi River basin for July 2007 using the National Centers for Environmental Prediction ensemble forecasts. By applying scoring rules, the BMA forecasts were verified and showed better performances compared with the empirical probabilistic ensemble forecasts, particularly for extreme precipitation. Finally, possible improvements and a^plication of this method to the downscaling of climate change scenarios were discussed.

Spatio-seasonal patterns of fish diversity,Haizhou Bay,China

Spatial-seasonal patterns in fish diversity in Haizhou Bay were studied based on stratified random surveys conducted in 2011.Principal component analysis was conducted to distinguish different diversity components,and the relationships among 11 diversity indices were explored.Generalized additive models were constructed to examine the environmental effects on diversity indices.Eleven diversity indices were grouped into four components:(1) species numbers and richness,(2) heterogeneous indices,(3) evenness,and(4) taxonomic relatedness.The results show that diversity indices among different components are complementary.Spatial patterns show that fish diversity was higher in coastal areas,which was affected by complex bottom topography and spatial variations of water mass and currents.Seasonal trends could be best explained by the seasonal migration of dominant fish species.Fish diversity generally declined with increasing depth except for taxonomic distinctness,which increased with latitude.In addition,bottom temperature had a significant effect on diversity index of richness.These results indicate that substrate complexity and environmental gradients had important influences on fish diversity patterns,and these factors should be considered in fishery resource management and conservation.Furthermore,diversity in two functional groups(demersal/pelagic fishes) was influenced by different environmental factors.Therefore,the distribution of individual species or new indicators in diversity should be applied to examine spatio-seasonal variations in fish diversity.