Abstract: 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.

Abstract: Egg geometry can be described using Preston's equation, which has seldom been used to calculate egg volume (V) and surface area (S) to explore S versus V scaling relationships. Herein, we provide an explicit re-expression of Preston's equation (designated as EPE) to calculate V and S, assuming that an egg is a solid of revolution. The side (longitudinal) profiles of 2221 eggs of six avian species were digitized, and the EPE was used to describe each egg profile. The volumes of 486 eggs from two avian species predicted by the EPE were compared with those obtained using water displacement in graduated cylinders. There was no significant difference in V using the two methods, which verified the utility of the EPE and the hypothesis that eggs are solids of revolution. The data also indicated that V is proportional to the product of egg length (L) and maximum width (W) squared. A 2/3-power scaling relationship between S and V for each species was observed, that is, S is proportional to (LW2)(2/3). These results can be extended to describe the shapes of the eggs of other species to study the evolution of avian (and perhaps reptilian) eggs.

Abstract: The Dirichlet problem for the Laplace equation in normal-polar annuli is addressed by using a suitable Fourier-like technique. Attention is in particular focused on the wide class of domains whose boundaries are defined by the so-called superformula introduced by Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica© is developed in order to validate the proposed methodology. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained.

Abstract: Gielis curves and surfaces can describe a wide range of natural shapes and they have been used in various studies in biology and physics as descriptive tool. This has stimulated the generalization of widely used computational methods. Here we show that proper normalization of the Levenberg-Marquardt algorithm allows for efficient and robust reconstruction of Gielis curves, including self-intersecting and asymmetric curves, without increasing the overall complexity of the algorithm. Then, we show how complex curves of k-type can be constructed and how solutions to the Dirichlet problem for the Laplace equation on these complex domains can be derived using a semi-Fourier method. In all three methods, descriptive and computational power and efficiency is obtained in a surprisingly simple way.

Abstract: Simon Stevin (1548-1620) is mainly known for the decimal system and his Clootkrans proof. His influence is also profound in infinitesimal calculus, mechanics, and even in abstract algebra and today’s conception of polynomials, algorithms, and generic programming. Here we review his influence as assessed in generic programming. According to Dr. Stepanov, one of the most influential researchers in generic programming, Stevin’s work on polynomials can be regarded as the essence of generic programming: an algorithm from one domain can be applied in another similar domain.

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: Möbius bands have been studied extensively, mainly in topology. Generalized Möbius-Listing surfaces and bodies providing a full geometrical generalization, is a quite new field, motivated originally by solutions of boundary value problems. Analogous to cutting of the original Möbius band, for this class of surfaces and bodies, results have been obtained when cutting such bodies or surfaces. In general, cutting leads to interlinked and intertwined different surfaces or bodies, resulting in very complex systems. However, under certain conditions, the result of cutting can be a single surface or body, which reduces complexity considerably. Our research is motivated by this reduction of complexity. In the study of cutting Generalized Möbius-Listing bodies with polygons as cross section, the conditions under which a single body results, displaying the Möbius phenomenon of a one-sided body, have been determined for even and odd polygons. These conditions are based on congruence and rotational symmetry of the resulting cross sections after cutting, and on the knife cutting the origin. The Möbius phenomenon is important, since the process of cutting (or separation of zones in a GML body in general) then results in a single body, not in different, intertwined domains. In all previous works it was assumed that the cross section of the GML bodies is constant, but the main result of this paper is that it is sufficient that only one cross section on the whole GML structure meets the conditions for the Möbius phenomenon to occur. Several examples are given to illustrate this.

Abstract: The principle of similarity (Thompson, 1917) states that the weight of an organism follows the 3/2-power law of its surface area and is proportional to its volume on the condition that the density is constant. However, the allometric relationship between leaf weight and leaf area has been reported to greatly deviate from the 3/2-power law, with the irregularity of leaf density largely ignored for explaining this deviation. Here, we choose 11 bamboo species to explore the allometric relationships among leaf area (A), density (ρ), length (L), thickness (T), and weight (W). Because the edge of a bamboo leaf follows a simplified two-parameter Gielis equation, we could show that A ∝ L2 and that A ∝ T2. This then allowed us to derive the density-thickness allometry ρ ∝ Tb and the weight-area allometry W ∝ A(b+3)/2 ≈ A9/8, where b approximates −3/4. Leaf density is strikingly negatively associated with leaf thickness, and it is this inverse relationship that results in the weight-area allometry to deviate from the 3/2-power law. In conclusion, although plants are prone to invest less dry mass and thus produce thinner leaves when the leaf area is sufficient for photosynthesis, such leaf thinning needs to be accompanied with elevated density to ensure structural stability. The findings provide the insights on the evolutionary clue about the biomass investment and output of photosynthetic organs of plants. Because of the importance of leaves, plants could have enhanced the ratio of dry material per unit area of leaf in order to increase the efficiency of photosynthesis, relative the other parts of plants. Although the conclusion is drawn only based on 11 bamboo species, it should also be applicable to the other plants, especially considering previous works on the exponent of the weight-area relationship being less than 3/2 in plants.

Abstract: Delve into the diverse aspects of plant morphology, their responses to global climate change, and the spatiotemporal dynamics of forest productivity. Join us on a journey through the intricate web of plant characteristics and their impact on the environment.

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.

Abstract: Plant leaves exhibit diverse shapes that enable them to utilize a light resource maximally. If there were a general parametric model that could be used to calculate leaf area for different leaf shapes, it would help to elucidate the adaptive evolutional link among plants with the same or similar leaf shapes. We propose a simplified version of the original Gielis equation (SGE), which was developed to describe a variety of object shapes ranging from a droplet to an arbitrary polygon. We used this equation to fit the leaf profiles of 53 species (among which, 48 bamboo plants, 5 woody plants, and 10 geographical populations of a woody plant), totaling 3310 leaves. A third parameter (namely, the floating ratio c in leaf length) was introduced to account for the case when the theoretical leaf length deviates from the observed leaf length. For most datasets, the estimates of c were greater than zero but less than 10%, indicating that the leaf length predicted by the SGE was usually smaller than the actual length. However, the predicted leaf areas approximated their actual values after considering the floating ratios in leaf length. For most datasets, the mean percent errors of leaf areas were lower than 6%, except for a pooled dataset with 42 bamboo species. For the elliptical, lanceolate, linear, obovate, and ovate shapes, although the SGE did not fit the leaf edge perfectly, after adjusting the parameter c, there were small deviations of the predicted leaf areas from the actual values. This illustrates that leaves with different shapes might have similar functional features for photosynthesis, since the leaf areas can be described by the same equation. The anisotropy expressed as a difference in leaf shape for some plants might be an adaptive response to enable them to adapt to different habitats.

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 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: The Lamé curve is an extension of an ellipse, the latter being a special case. Dr. Johan Gielis further extended the Lamé curve in the polar coordinate system by introducing additional parameters (n1, n2, n3; m): rφ=1Acosm4φn2+1Bsinm4φn3−1/n1, which can be applied to model natural geometries. Here, r is the polar radius corresponding to the polar angle φ; A, B, n1, n2 and n3 are parameters to be estimated; m is the positive real number that determines the number of angles of the Gielis curve. Most prior studies on the Gielis equation focused mainly on its applications. However, the Gielis equation can also generate a large number of shapes that are rotationally symmetric and axisymmetric when A = B and n2 = n3, interrelated with the parameter m, with the parameters n1 and n2 determining the shapes of the curves. In this paper, we prove the relationship between m and the rotational symmetry and axial symmetry of the Gielis curve from a theoretical point of view with the condition A = B, n2 = n3. We also set n1 and n2 to take negative real numbers rather than only taking positive real numbers, then classify the curves based on extremal properties of r(φ) at φ = 0, π/m when n1 and n2 are in different intervals, and analyze how n1, n2 precisely affect the shapes of Gielis curves.