Abstract: The number and composition of species in a community can be quantified with alpha-diversity indices, including species richness (R), Simpson's index (D), and the Shannon-Wiener index (HGREEK TONOS). In forest communities, there are large variations in tree size among species and individuals of the same species, which result in differences in ecological processes and ecosystem functions. However, tree size inequality (TSI) has been largely neglected in studies using the available diversity indices. The TSI in the diameter at breast height (DBH) data for each of 999 20 m x 20 m forest census quadrats was quantified using the Gini index (GI), a measure of the inequality of size distribution. The generalized performance equation was used to describe the rotated and right-shifted Lorenz curve of the cumulative proportion of DBH and the cumulative proportion of number of trees per quadrat. We also examined the relationships of alpha-diversity indices with the GI using correlation tests. The generalized performance equation effectively described the rotated and right-shifted Lorenz curve of DBH distributions, with most root-mean-square errors (990 out of 999 quadrats) being < 0.0030. There were significant positive correlations between each of three alpha-diversity indices (i.e., R, D, and H') and the GI. Nevertheless, the total abundance of trees in each quadrat did not significantly influence the GI. This means that the TSI increased with increasing species diversity. Thus, two new indices are proposed that can balance alpha-diversity against the extent of TSI in the community: (1 – GI) x D, and (1 – GI) x H'. These new indices were significantly correlated with the original D and HGREEK TONOS, and did not increase the extent of variation within each group of indices. This study presents a useful tool for quantifying both species diversity and the variation in tree sizes in forest communities, especially in the face of cumulative species loss under global climate change.

Abstract: The goal of this study is to provide a rigorous tool to quantify the inequality of the leaf size distribution of an individual plant, thereby serving as a reference trait for quantifying plant adaptations to local environmental conditions. The tool to be presented and tested employs three components: (1) a performance equation (PE), which can produce flexible asymmetrical and symmetrical bell-shaped curves, (2) the Lorenz curve (i.e., the cumulative proportion of leaf size vs. the cumulative proportion of number of leaves), which is the basis for calculating, and (3) the Gini index, which measures the inequality of leaf size distribution. We sampled 12 individual plants of a dwarf bamboo and measured the area and dry mass of each leaf of each plant. We then developed a generalized performance equation (GPE) of which the PE is a special case and fitted the Lorenz curve to leaf size distribution using the GPE and PE. The GPE performed better than the PE in fitting the Lorenz curve. We compared the Gini index of leaf area distribution with that of leaf dry mass distribution and found that there was a significant difference between the two indices that might emerge from the scaling relationship between leaf dry mass and area. Nevertheless, there was a strong correlation between the two Gini indices (r2 = 0.9846). This study provides a promising tool based on the GPE for quantifying the inequality of leaf size distributions across individual plants and can be used to quantify plant adaptations to local environmental conditions.

Abstract: Stomata are essential for the exchange of water vapour and atmospheric gases between vascular plants and their external environments. The stomatal geometries of many plants appear to be elliptical. However, prior studies have not tested whether this is a mathematical reality, particularly since many natural shapes that appear to be ellipses are superellipses with greater or smaller edge curvature than predicted for an ellipse. Compared with the ellipse equation, the superellipse equation includes an additional parameter that allows generation of a larger range of shapes. We randomly selected 240 stomata from each of four Magnoliaceae species to test whether the stomatal geometries are superellipses or ellipses. The stomatal geometries for most stomata (943/960) were found to be described better using the superellipse equation. The traditional “elliptical stomata hypothesis” resulted in an underestimation of the area of stomata, whereas the superellipse equation accurately predicted stomatal area. This finding has important implications for the estimation of stomatal area in studies looking at stomatal shape, geometry, and function.

Abstract: Preston's equation is a general model describing the egg shape of birds. The parameters of Preston's equation are usually estimated after re-expressing it as the Todd-Smart equation and scaling the egg's actual length to two. This method assumes that the straight line through the two points on an egg's profile separated by the maximum distance (i.e., the longest axis of an egg's profile) is the mid-line. It hypothesizes that the photographed egg's profile is perfectly bilaterally symmetrical, which seldom holds true because of photographic errors and placement errors. The existing parameter estimation method for Preston's equation considers an angle of deviation for the longest axis of an egg's profile from the mid-line, which decreases prediction errors to a certain degree. Nevertheless, this method cannot provide an accurate estimate of the coordinates of the egg's center, and it leads to sub-optimal parameter estimation. Thus, it is better to account for the possible asymmetry between the two sides of an egg's profile along its mid-line when fitting egg-shape data. In this paper, we propose a method based on the optimization algorithm (optimPE) to fit egg-shape data and better estimate the parameters of Preston's equation by automatically searching for the optimal mid-line of an egg's profile and testing its validity using profiles of 59 bird eggs spanning a wide range of existing egg shapes. We further compared this method with the existing one based on multiple linear regression (lmPE). This study demonstrated the ability of the optimPE method to estimate numerical values of the parameters of Preston's equation and provide the theoretical egg length (i.e., the distance between two ends of the mid-line of an egg's profile) and the egg's maximum breadth. This provides a valuable approach for comparing egg shapes among conspecifics or across different species, or even different classes (e.g., birds and reptiles), in future investigations.

Abstract: Bamboo is an important component in subtropical and tropical forest communities. The plant has characteristic long lanceolate leaves with parallel venation. Prior studies have shown that the leaf shapes of this plant group can be well described by a simplified version (referred to as SGE-1) of the Gielis equation, a polar coordinate equation extended from the superellipse equation. SGE-1 with only two model parameters is less complex than the original Gielis equation with six parameters. Previous studies have seldom tested whether other simplified versions of the Gielis equation are superior to SGE-1 in fitting empirical leaf shape data. In the present study, we compared a three-parameter Gielis equation (referred to as SGE-2) with the two-parameter SGE-1 using the leaf boundary coordinate data of six bamboo species within the same genus that have representative long lanceolate leaves, with >300 leaves for each species. We sampled 2000 data points at approximately equidistant locations on the boundary of each leaf, and estimated the parameters for the two models. The root–mean–square error (RMSE) between the observed and predicted radii from the polar point to data points on the boundary of each leaf was used as a measure of the model goodness of fit, and the mean percent error between the RMSEs from fitting SGE-1 and SGE-2 was used to examine whether the introduction of an additional parameter in SGE-1 remarkably improves the model’s fitting. We found that the RMSE value of SGE-2 was always smaller than that of SGE-1. The mean percent errors among the two models ranged from 7.5% to 20% across the six species. These results indicate that SGE-2 is superior to SGE-1 and should be used in fitting leaf shapes. We argue that the results of the current study can be potentially extended to other lanceolate leaf shapes.

Abstract: Many natural objects exhibit radial or axial symmetry in a single plane. However, a universal tool for simulating and fitting the shapes of such objects is lacking. Herein, we present an R package called 'biogeom' that simulates and fits many shapes found in nature. The package incorporates novel universal parametric equations that generate the profiles of bird eggs, flowers, linear and lanceolate leaves, seeds, starfish, and tree-rings, and three growth-rate equations that generate the profiles of ovate leaves and the ontogenetic growth curves of animals and plants. 'biogeom' includes several empirical datasets comprising the boundary coordinates of bird eggs, fruits, lanceolate and ovate leaves, tree rings, seeds, and sea stars. The package can also be applied to other kinds of natural shapes similar to those in the datasets. In addition, the package includes sigmoid curves derived from the three growth-rate equations, which can be used to model animal and plant growth trajectories and predict the times associated with maximum growth rate. 'biogeom' can quantify the intra- or interspecific similarity of natural outlines, and it provides quantitative information of shape and ontogenetic modification of shape with important ecological and evolutionary implications for the growth and form of the living world.

Abstract: Plants have diverse leaf shapes that have evolved to adapt to the environments they have experienced over their evolutionary history. Leaf shape and leaf size can greatly influence the growth rate, competitive ability, and productivity of plants. However, researchers have long struggled to decide how to properly quantify the complexity of leaf shape. Prior studies recommended the leaf roundness index (RI = 4πA/P2) or dissection index (DI = ), where P is leaf perimeter and A is leaf area. However, these two indices merely measure the extent of the deviation of leaf shape from a circle, which is usually invalid as leaves are seldom circular. In this study, we proposed a simple measure, named the ellipticalness index (EI), for quantifying the complexity of leaf shape based on the hypothesis that the shape of any oval leaf can be regarded as a variation from a standard ellipse. 2220 leaves from nine species of Magnoliaceae were sampled to check the validity of the EI. We also tested the validity of the Montgomery equation (ME), which assumes a proportional relationship between leaf area and the product of leaf length and width, because the EI actually comes from the proportionality coefficient of the ME. We also compared the ME with five other models of leaf area. The ME was found to be the best model for calculating leaf area based on consideration of the trade-off between model fit vs. complexity, which strongly supported the robustness of the EI for describing oval leaf shape. The new index can account for both leaf shape and size, and we conclude that it is a promising method for quantifying and comparing oval leaf shapes across species in future studies.

Abstract: The shape of leaf laminae exhibits considerable diversity and complexity that reflects adaptations to environmental factors such as ambient light and precipitation as well as phyletic legacy. Many leaves appear to be elliptical which may represent a ‘default’ developmental condition. However, whether their geometry truly conforms to the ellipse equation (EE), i.e., (x/a)2 + (y/b)2 = 1, remains conjectural. One alternative is described by the superellipse equation (SE), a generalized version of EE, i.e., |x/a|n +|y/b|n = 1. To test the efficacy of EE versus SE to describe leaf geometry, the leaf shapes of two Michelia species (i.e., M. cavaleriei var. platypetala, and M. maudiae), were investigated using 60 leaves from each species. Analysis shows that the majority of leaves (118 out of 120) had adjusted root-mean-square errors of < 0.05 for the nonlinear fitting of SE to leaf geometry, i.e., the mean absolute deviation from the polar point to leaf marginal points was smaller than 5% of the radius of a hypothesized circle with its area equaling leaf area. The estimates of n for the two species were ˂ 2, indicating that all sampled leaves conformed to SE and not to EE. This study confirms the existence of SE in leaves, linking this to its potential functional advantages, particularly the possible influence of leaf shape on hydraulic conductance.

Abstract: Leaf shape is closely related to economics of leaf support and leaf functions, including light interception, water use, and CO2 uptake, so correct quantification of leaf shape is helpful for studies of leaf structure/function relationships. There are some extant indices for quantifying leaf shape, including the leaf width/length ratio (W/L), leaf shape fractal dimension (FD), leaf dissection index, leaf roundness index, standardized bilateral symmetrical index, etc. W/L ratio is the simplest to calculate, and recent studies have shown the importance of the W/L ratio in explaining the scaling exponent of leaf dry mass vs. leaf surface area and that of leaf surface area vs. leaf length. Nevertheless, whether the W/L ratio could reflect sufficient geometrical information of leaf shape has been not tested. The FD might be the most accurate measure for the complexity of leaf shape because it can characterize the extent of the self-similarity and other planar geometrical features of leaf shape. However, it is unknown how strongly different indices of leaf shape complexity correlate with each other, especially whether W/L ratio and FD are highly correlated. In this study, the leaves of nine Magnoliaceae species (>140 leaves for each species) were chosen for the study. We calculated the FD value for each leaf using the box-counting approach, and measured leaf fresh mass, surface area, perimeter, length, and width. We found that FD is significantly correlated to the W/L ratio and leaf length. However, the correlation between FD and the W/L ratio was far stronger than that between FD and leaf length for each of the nine species. There were no strong correlations between FD and other leaf characteristics, including leaf area, ratio of leaf perimeter to area, fresh mass, ratio of leaf fresh mass to area, and leaf roundness index. Given the strong correlation between FD and W/L, we suggest that the simpler index, W/L ratio, can provide sufficient information of leaf shape for similarly-shaped leaves. Future studies are needed to characterize the relationships among FD and W/L in leaves with strongly varying shape, e.g., in highly dissected leaves.

Abstract: Many cross-sectional shapes of plants have been found to approximate a superellipse rather than an ellipse. Square bamboos, belonging to the genus Chimonobambusa (Poaceae), are a group of plants with round-edged square-like culm cross sections. The initial application of superellipses to model these culm cross sections has focused on Chimonobambusa quadrangularis (Franceschi) Makino. However, there is a need for large scale empirical data to confirm this hypothesis. In this study, approximately 750 cross sections from 30 culms of C. utilis were scanned to obtain cross-sectional boundary coordinates. A superellipse exhibits a centrosymmetry, but in nature the cross sections of culms usually deviate from a standard circle, ellipse, or superellipse because of the influences of the environment and terrain, resulting in different bending and torsion forces during growth. Thus, more natural cross-sectional shapes appear to have the form of a deformed superellipse. The superellipse equation with a deformation parameter (SEDP) was used to fit boundary data. We find that the cross-sectional shapes (including outer and inner rings) of C. utilis can be well described by SEDP. The adjusted root-mean-square error of SEDP is smaller than that of the superellipse equation without a deformation parameter. A major finding is that the cross-sectional shapes can be divided into two types of superellipse curves: hyperellipses and hypoellipses, even for cross sections from the same culm. There are two proportional relationships between ring area and the product of ring length and width for both the outer and inner rings. The proportionality coefficients are significantly different, as a consequence of the two different superellipse types (i.e., hyperellipses and hypoellipses). The difference in the proportionality coefficients between hyperellipses and hypoellipses for outer rings is greater than that for inner rings. This work informs our understanding and quantifying of the longitudinal deformation of plant stems for future studies to assess the influences of the environment on stem development. This work is also informative for understanding the deviation of natural shapes from a strict rotational symmetry.

Abstract: Starting from logarithmic, sinusoidal and power spirals, it is shown how these spirals are connected directly with Chebyshev polynomials, Lamé curves, with allometry and Antonelli-metrics in Finsler geometry. Curvature is a crucial concept in geometry both for closed curves and equiangular spirals, and allowed Dillen to give a general definition of spirals. Many natural shapes can be described as a combination of one of two basic shapes in nature—circle and spiral—with Gielis transformations. Using this idea, shape description itself is used to develop a novel approach to anisotropic curvature in nature. Various examples are discussed, including fusion in flowers and its connection to the recently described pseudo-Chebyshev functions.

Abstract: There is convincing evidence for a scaling relationship between leaf dry weight (DW) and leaf surface area (A) for broad-leaved plants, and most estimates of the scaling exponent of DW vs. A are greater than unity. However, the scaling relationship of leaf fresh weight (FW) vs. A has been largely neglected. In the present study, we examined whether there is a statistically strong scaling relationship between FW and A and compared the goodness of fit to that of DW vs. A. Between 250 and 520 leaves from each of 12 bamboo species within 2 genera (Phyllostachys and Pleioblastus) were investigated. The reduced major axis regression protocols were used to determine scaling relationships. The fit for the linearized scaling relationship of FW vs. A was compared with that of DW vs. A using the coefficient of determination (i.e., r2). A stronger scaling relationship between FW and A than that between DW and A was observed for each of the 12 bamboo species investigated. Among the 12 species examined, five had significantly smaller scaling exponents of FW vs. A compared to those of DW vs. A; only one species had a scaling exponent of FW vs. A greater than that of DW vs. A. No significant difference between the two scaling exponents was observed for the remaining 6 species. Researchers conducting future studies might be well advised to consider the influence of leaf fresh weight when exploring the scaling relationships of foliar biomass allocation patterns.

Abstract: The leaf area, as an important leaf functional trait, is thought to be related to leaf length and width. Our recent study showed that the Montgomery equation, which assumes that leaf area is proportional to the product of leaf length and width, applied to different leaf shapes, and the coefficient of proportionality (namely the Montgomery parameter) range from 1/2 to π/4. However, no relevant geometrical evidence has previously been provided to support the above findings. Here, four types of representative leaf shapes (the elliptical, sectorial, linear, and triangular shapes) were studied. We derived the range of the estimate of the Montgomery parameter for every type. For the elliptical and triangular leaf shapes, the estimates are π/4 and 1/2, respectively; for the linear leaf shape, especially for the plants of Poaceae that can be described by the simplified Gielis equation, the estimate ranges from 0.6795 to π/4; for the sectorial leaf shape, the estimate ranges from 1/2 to π/4. The estimates based on the observations of actual leaves support the above theoretical results. The results obtained here show that the coefficient of proportionality of leaf area versus the product of leaf length and width only varies in a small range, maintaining the allometric relationship for leaf area and thereby suggesting that the proportional relationship between leaf area and the product of leaf length and width broadly remains stable during leaf evolution.