The Generalized Gielis Geometric Equation and Its Application

Comparison of Two Simplified Versions of the Gielis Equation for Describing the Shape of Bamboo Leaves

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.

Quantifying the Variation in the Geometries of the Outer Rims of Corolla Tubes of Vinca major L

Many geometries of plant organs can be described by the Gielis equation, a polar coordinate equation extended from the superellipse equation, r=a|cosm4φ|n2+|1ksinm4φ|n3−1/n1. Here, r is the polar radius corresponding to the polar angle φ; m is a positive integer that determines the number of angles of the Gielis curve when φ ∈ [0 to 2π); and the rest of the symbols are parameters to be estimated. The pentagonal radial symmetry of calyxes and corolla tubes in top view is a common feature in the flowers of many eudicots. However, prior studies have not tested whether the Gielis equation can depict the shapes of corolla tubes. We sampled randomly 366 flowers of Vinca major L., among which 360 had five petals and pentagonal corolla tubes, and six had four petals and quadrangular corolla tubes. We extracted the planar coordinates of the outer rims of corolla tubes (in top view) (ORCTs), and then fitted the data with two simplified versions of the Gielis equation with k = 1 and m = 5: r=acos54φn2+sin54φn3−1/n1 (Model 1), and r=acos54φn2+sin54φn2−1/n1 (Model 2). The adjusted root mean square error (RMSE_adj) was used to evaluate the goodness of fit of each model. In addition, to test whether ORCTs are radially symmetrical, we correlated the estimates of n₂ and n₃ in Model 1 on a log-log scale. The results validated the two simplified Gielis equations. The RMSE_adj values for all corolla tubes were smaller than 0.05 for both models. The numerical values of n₂ and n₃ were demonstrated to be statistically equal based on the regression analysis, which suggested that the ORCTs of V. major are radially symmetrical. It suggests that Model 1 can be replaced by the simpler Model 2 for fitting the ORCT in this species. This work indicates that the pentagonal or quadrangular corolla tubes (in top view) can both be modeled by the Gielis equation and demonstrates that the pentagonal or quadrangular corolla tubes of plants tend to form radial symmetrical geometries during their development and growth.

'biogeom': An R package for simulating and fitting natural shapes

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.

Comparison of a universal (but complex) model for avian egg shape with a simpler model

Recently, a universal equation by Narushin, Romanov, and Griffin (hereafter, the NRGE) was proposed to describe the shape of avian eggs. While NRGE can simulate the shape of spherical, ellipsoidal, ovoidal, and pyriform eggs, its predictions were not tested against actual data. Here, we tested the validity of the NRGE by fitting actual data of egg shapes and compared this with the predictions of our simpler model for egg shape (hereafter, the SGE). The eggs of nine bird species were sampled for this purpose. NRGE was found to fit the empirical data of egg shape well, but it did not define the egg length axis (i.e., the rotational symmetric axis), which significantly affected the prediction accuracy. The egg length axis under the NRGE is defined as the maximum distance between two points on the scanned perimeter of the egg's shape. In contrast, the SGE fitted the empirical data better, and had a smaller root-mean-square error than the NRGE for each of the nine eggs. Based on its mathematical simplicity and goodness-of-fit, the SGE appears to be a reliable and useful model for describing egg shape.

Evidence That Supertriangles Exist in Nature from the Vertical Projections of Koelreuteria paniculata Fruit

Many natural radial symmetrical shapes (e.g., sea stars) follow the Gielis equation (GE) or its twin equation (TGE). A supertriangle (three triangles arranged around a central polygon) represents such a shape, but no study has tested whether natural shapes can be represented as/are supertriangles or whether the GE or TGE can describe their shape. We collected 100 pieces of Koelreuteria paniculata fruit, which have a supertriangular shape, extracted the boundary coordinates for their vertical projections, and then fitted them with the GE and TGE. The adjusted root mean square errors (RMSEadj) of the two equations were always less than 0.08, and >70% were less than 0.05. For 57/100 fruit projections, the GE had a lower RMSEadj than the TGE, although overall differences in the goodness of fit were non-significant. However, the TGE produces more symmetrical shapes than the GE as the two parameters controlling the extent of symmetry in it are approximately equal. This work demonstrates that natural supertriangles exist, validates the use of the GE and TGE to model their shapes, and suggests that different complex radially symmetrical shapes can be generated by the same equation, implying that different types of biological symmetry may result from the same biophysical mechanisms.

A General Model for Describing the Ovate Leaf Shape

Many plant species produce ovate leaves, but there is no general parametric model for describing this shape. Here, we used two empirical nonlinear equations, the beta and Lobry–Rosso–Flandrois (LRF) equations, and their modified forms (referred to as the Mbeta and MLRF equations for convenience), to generate bilaterally symmetrical curves along the x-axis to form ovate leaf shapes. In order to evaluate which of these four equations best describes the ovate leaf shape, we used 14 leaves from 7 Neocinnamomum species (Lauraceae) and 72 leaves from Chimonanthus praecox (Calycanthaceae). Using the AIC and adjusted root mean square error to compare the fitted results, the modified equations fitted the leaf shapes better than the unmodified equations. However, the MLRF equation provided the best overall fit. As the parameters of the MLRF equation represent leaf length, maximum leaf width, and the distance from leaf apex to the point associated with the maximum leaf width along the leaf length axis, these findings are potentially valuable for studying the influence of environmental factors on leaf shape, differences in leaf shape among closely related plant species with ovate leaf shapes, and the extent to which leaves are bilaterally symmetrical. This is the first work in which temperature-dependent developmental equations to describe the ovate leaf shape have been employed, as previous studies lacked similar leaf shape models. In addition, prior work seldom attempted to describe real ovate leaf shapes. Our work bridges the gap between theoretical leaf shape models and empirical leaf shape indices that cannot predict leaf shape profiles.

Application of an Ovate Leaf Shape Model to Evaluate Leaf Bilateral Asymmetry and Calculate Lamina Centroid Location

Leaf shape is an important leaf trait, with ovate leaves common in many floras. Recently, a new leaf shape model (referred to as the MLRF equation) derived from temperature-dependent bacterial growth was proposed and demonstrated to be valid in describing leaf boundaries of many species with ovate leaf shape. The MLRF model's parameters can provide valuable information of leaf shape, including the ratio of lamina width to length and the lamina centroid location on the lamina length axis. However, the model wasn't tested on a large sample of a single species, thereby limiting its overall evaluation for describing leaf boundaries, for evaluating lamina bilateral asymmetry and for calculating lamina centroid location. In this study, we further test the model using data from two Lauraceae species, Cinnamomum camphora and Machilus leptophylla, with >290 leaves for each species. The equation was found to be credible for describing those shapes, with all adjusted root-mean-square errors (RMSE) smaller than 0.05, indicating that the mean absolute deviation is smaller than 5% of the radius of an assumed circle whose area equals lamina area. It was also found that the larger the extent of lamina asymmetry, the larger the adjusted RMSE, with approximately 50% of unexplained variation by the model accounted for by the lamina asymmetry, implying that this model can help to quantify the leaf bilateral asymmetry in future studies. In addition, there was a significant difference between the two species in their centroid ratio, i.e., the distance from leaf petiole to the point on the lamina length axis associated with leaf maximum width to the leaf maximum length. It was found that a higher centroid ratio does not necessarily lead to a greater investment of mass to leaf petiole relative to lamina, which might depend on the petiole pattern.

A Superellipse with Deformation and Its Application in Describing the Cross-Sectional Shapes of a Square Bamboo

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.

Rate of Entropy Production in Evolving Interfaces and Membranes under Astigmatic Kinematics: Shape Evolution in Geometric-Dissipation Landscapes

This paper presents theory and simulation of viscous dissipation in evolving interfaces and membranes under kinematic conditions, known as astigmatic flow, ubiquitous during growth processes in nature. The essential aim is to characterize and explain the underlying connections between curvedness and shape evolution and the rate of entropy production due to viscous bending and torsion rates. The membrane dissipation model used here is known as the Boussinesq-Scriven fluid model. Since the standard approaches in morphological evolution are based on the average, Gaussian and deviatoric curvatures, which comingle shape with curvedness, this paper introduces a novel decoupled approach whereby shape is independent of curvedness. In this curvedness-shape landscape, the entropy production surface under constant homogeneous normal velocity decays with growth but oscillates with shape changes. Saddles and spheres are minima while cylindrical patches are maxima. The astigmatic flow trajectories on the entropy production surface, show that only cylinders and spheres grow under the constant shape. Small deviations from cylindrical shapes evolve towards spheres or saddles depending on the initial condition, where dissipation rates decrease. Taken together the results and analysis provide novel and significant relations between shape evolution and viscous dissipation in deforming viscous membrane and surfaces.

Increase in Absolute Leaf Water Content Tends to Keep Pace with That of Leaf Dry Mass—Evidence from Bamboo Plants

Leaves, as the most important photosynthetic organ of plants, are intimately associated with plant function and adaptation to environmental changes. The scaling relationship of the leaf dry mass (or the fresh mass) vs. leaf surface area has been referred to as “diminishing returns”, suggesting that the leaf area fails to increase in proportion to leaf dry mass (or fresh mass). However, previous studies used materials across different families, and there is lack of studies testing whether leaf fresh mass is proportional to the leaf dry mass for the species in the same family, and examining the influence of the scaling of leaf dry mass vs. fresh mass on two kinds of diminishing returns based on leaf dry mass and fresh mass. Bamboo plants (Poaceae: Bambusoideae) are good materials for doing such a study, which have astonishingly similar leaf shapes across species. Bamboo leaves have a typical parallel venation pattern. In general, a parallel venation pattern tends to produce a more stable symmetrical leaf shape than the pinnate and palmate venation patterns. The symmetrical parallel veins enable leaves to more regularly hold water, which is more likely to result in a proportional relationship between the leaf dry mass and absolute water content, which consequently determines whether the scaling exponent of the leaf dry mass vs. area is significantly different from (or the same as) that of the leaf fresh mass vs. area. In the present study, we used the data of 101 bamboo species, cultivars, forms and varieties (referred to as 101 (bamboo) taxa below for convenience) to analyze the scaling relationships between the leaf dry mass and area, and between leaf fresh mass and area. We found that the confidence intervals of the scaling exponents of the leaf fresh mass vs. dry mass of 68 out of the 101 taxa included unity, which indicates that for most bamboo species (67.3%), the increase in leaf water mass keeps pace with that of leaf dry mass. There was a significant scaling relationship between either leaf dry mass or fresh mass, and the leaf surface area for each studied species. We found that there was no significant difference between the scaling exponent of the leaf dry mass vs. leaf area and that of the leaf fresh mass vs. leaf area when the leaf dry mass was proportional to the leaf fresh mass. The goodness of fit to the linearized scaling relationship of the leaf fresh mass vs. area was better than that of the leaf dry mass vs. area for each of the 101 bamboo taxa. In addition, there were significant differences in the normalized constants of the leaf dry mass vs. fresh mass among the taxa (i.e., the differences in leaf water content), which implies the difference in the adaptabilities to different environments across the taxa.

Cross-Sectional Performance of Hollow Square Prisms with Rounded Edges

Hollow-section columns are one of the mechanically superior structures with high buckling resistance and high bending stiffness. The mechanical properties of the column are strongly influenced by the cross-sectional shape. Therefore, when evaluating the stability of a column against external forces, it is necessary to reproduce the cross-sectional shape accurately. In this study, we propose a mathematical method to describe a polygonal section with rounded edges and vertices. This mathematical model would be quite useful for analyzing the mechanical properties of plants and designing plant-mimicking functional structures, since the cross-sections of the actual plant culms and stems often show rounded polygons.