方法文章

利用数字叶形态学和单位叶面积质量通过化石叶片重建陆地古气候与古生态

3.6K 次观看

⸱

DOI：

10.3791/66838

⸱

2024年10月25日

 ,  ,  ,  ,  , 
1Department of Biology, University of Washington, 2Department of Paleobiology, National Museum of Natural History, Smithsonian Institution, 3Department of Geology and Geophysics, Program in Ecology, University of Wyoming, 4Department of Ecology and Evolutionary Biology, University of Michigan Herbarium, 5Department of Earth and Environmental Sciences, Wesleyan University, 6Department of Geosciences, Baylor University

本文内容

摘要

本方案展示了利用数字叶形态学和单位叶面积质量重建方法，对化石叶片的连续叶形态特征进行数字化测量与分析，以重建古气候和古生态环境。

摘要

气候和环境强烈影响植物叶片的大小、形状和锯齿特征（叶形态）。这些关系，特别是在木本双子叶被子植物中，已被用于开发基于叶片的古气候和古生态代用指标，应用于重建地球历史上过去约1.2亿年来的陆地古生态系统。此外，由于这些关系已在现存植物中得到记录，因此对于理解植物进化特征以及植物如何响应气候和环境变化具有重要意义。为了对现代和化石植物进行此类分析，必须采用可重复的方法准确测量叶片形态特征。本方案描述了一种基于计算机的方法，用于测量和分析现代及化石叶片中的多种叶形态变量。该方法可测量叶形态性状，特别是与叶缘锯齿、叶面积、叶裂程度和狭长性相关的变量，这些变量被用于数字叶形态代用法重建古气候；同时也可测量叶柄宽度和叶面积，用于重建单位叶面积质量，这是一种古生态代用指标。由于这种数字化叶片性状测量方法可应用于化石和现存植物，因此其应用不仅限于古气候和古生态重建，还可用于探索可能有助于理解叶片形态功能、叶片发育、叶性状系统发育关系以及植物进化的叶片性状。

引言

叶片是基本的生产单位，能够促进植物与其周围环境之间能量（如光、热）和物质（如二氧化碳、水蒸气）的交换1,2。为了实现这些功能，叶片必须在静止和有风的空气中依靠机械结构支撑自身重量以抵抗重力3,4。由于这些内在联系，叶片的大小、形状和叶缘锯齿特征（即形态学）在多个方面反映了其功能和生物力学特性，并为了解其生存环境与生态提供了重要线索。以往研究已量化了现代全球范围内叶片形态学与气候及生态之间的关系，建立了可用于化石叶片组合的代用指标5,6。这些代用指标为重建古气候和古生态提供了重要机会，有助于更深入理解地球历史上各系统之间复杂的相互作用。本文详细介绍了两种代用指标的应用方法：1）基于单位叶面积质量的古生态重建方法；2）数字叶片形态学用于古气候重建。

单位叶面积干重（MA）是新近纪植物学和古植物学中常被测定的植物性状。MA 的主要价值，尤其是对于化石重建而言，在于它是叶片经济谱的一部分；叶片经济谱是一组高度相关叶片性状的综合轴线，包括叶片光合速率、叶片寿命以及单位质量叶片的养分含量7。通过化石重建 MA 的能力，为我们了解这些原本无法获取的代谢和化学过程提供了窗口，并最终揭示有关植物生态策略和生态系统功能的有用信息。

Royer 等人5提出了一种基于叶片面积和叶柄宽度来估算木本双子叶被子植物（非单子叶植物）化石叶片MA的方法。理论上，叶柄起到悬臂的作用，将叶片重量维持在最佳位置3,4。因此，构成梁强度最主要部分的叶柄横截面积应与叶片质量高度相关。通过将叶柄形状简化为圆柱形管状结构，叶柄横截面积可用叶柄宽度的平方来表示，从而能够从二维化石中估算叶片质量（更多细节见 Royer 等人5）。叶片面积可直接测量。叶柄宽度的平方除以叶片面积（即叶柄度量指标；表1）共同构成了对化石MA的良好代用指标，使古植物学家得以进入现代基于性状的生态学研究领域。MA重建方法还已扩展至阔叶具柄裸子植物5,8、草本被子植物8以及蕨类植物9，这些类群所得到的关系不同于木本双子叶被子植物之间的关系，也彼此不同。扩展后的木本双子叶植物数据集以及用于重建样地水平MA均值与方差的新回归方程，使得推断化石植物群中木本双子叶被子植物叶片经济策略的多样性及其主要策略成为可能10。

一个世纪以来，人们已经注意到植物叶片形态特征与其生长气候之间的关系11,12。具体而言，木本双子叶被子植物叶片的形态特征与温度和湿度密切相关13。这一关系已成为众多单变量14,15,16,17 和多变量6,18,19,20,21,22 叶片形态古气候代用指标的基础。在过去约1.2亿年（白垩纪至今）的地球历史中，单变量和多变量叶片形态古气候方法已被广泛应用于各大洲以被子植物为主的化石植物群研究23。

叶片形态学古气候代用指标中利用的两个基本观测现象是：1）叶片大小与年均降水量（MAP）之间的关系，以及2）叶片锯齿（即叶片边缘向外突出的部分）与年均温度（MAT）之间的关系。具体而言，某一地区所有木本双子叶被子植物物种的平均叶片大小与年均降水量呈正相关，而该地区具有锯齿状叶片的木本双子叶被子植物物种所占比例，以及锯齿的大小和数量，则与年均温度呈负相关6,12,13,14,15,16,24。

这些叶片形态特征与气候之间关系的功能性联系，在理论和观测上均得到了有力支持1,2,25。例如，尽管较大的叶片可提供更大的光合作用表面积，但其需要更强的结构支撑，通过蒸腾作用会损失更多的水分，并由于更厚的边界层而保留更多的感热1,26,27。因此，在更湿润、更炎热的环境中，较大叶片更为常见，因为蒸腾作用增强所导致的水分流失能够有效冷却叶片，且在此类环境中水分补充相对充足，问题较小。相比之下，在干燥炎热的气候中，较小的叶片通过减少水分流失并增加感热散失来避免过热28,29。然而，对于其他叶片性状而言，究竟是哪些因素或因素组合对解释其功能性联系起主导作用，目前仍不明确。例如，针对叶片锯齿与年均温（MAT）关系的解释，已有多种假说被提出，包括叶片冷却、芽体高效包装、对较薄叶片的增强支撑与养分供给、通过排水器进行吐水作用，以及提高生长季早期的生产力等30,31,32,33。

大多数叶片形态古气候代用指标依赖于对叶片性状进行类别划分，而非对连续变量进行定量测量，这导致了若干潜在的缺陷。类别化方法排除了纳入与气候高度相关的连续性测量所捕获的更详细信息（例如锯齿数量、叶片狭长程度），从而可能降低古气候估算的准确性6,20,34。此外，在某些叶片性状评分方法中，被类别化评分的性状可能具有模糊性，导致可重复性问题，且部分性状与气候之间功能联系的实证支持有限6,15,16,35,36。

为解决这些不足，Huff 等人20提出了一种称为数字叶形态学（digital leaf physiognomy, DiLP）的方法，用于数字化测量连续的叶片性状。与以往方法相比，DiLP 的一个关键优势在于其依赖的性状具有以下特点：1）可在不同使用者之间可靠测量；2）在本质上是连续的；3）在功能上与气候相关联；4）在不同生长季节之间表现出表型可塑性6,37。这使得 DiLP 对年均温（MAT）和年均降水量（MAP）的估算比以往基于叶片形态学的古气候方法更为准确6。此外，该方法通过提供处理受损和不完整叶片的步骤，适应了化石记录不完整的特性。DiLP 方法已成功应用于来自多个大陆、跨越广泛地质年代的多种化石植物群6,38,39,40,41,42。

以下方案扩展了先前研究中描述的方法5,6,20,34，旨在说明利用木本双子叶被子植物化石叶片，结合DiLP和MA重建方法，重建古气候与古生态所需的各项操作步骤（有关本方案所测量与计算变量的解释，请参见表1）。此外，本方案还提供了用于记录和计算叶性状的步骤，这些性状虽未包含在DiLP或MA分析中，但易于实施，并可对叶片形态提供有价值的表征信息（表1）。本方案按以下顺序进行：1）化石叶片成像；2）叶片数字化处理，分为五种可能的处理情景；3）叶片数字化测量，同样按这五种可能的处理情景组织；4）使用R语言程序包dilp进行DiLP和MA分析10。

M的实验方案A 重建被嵌入在DiLP方案中，因为两者便于同时准备和测量。如果用户对M感兴趣A 仅用于分析时，无论叶缘是否具锯齿，均应遵循DiLP制备方案2中所述的制备步骤，以及仅包含叶柄宽度、叶柄面积和叶面积测量的测量步骤。随后，用户可在dilp R软件包中运行相应的函数以执行MA 重建

方案

1. 化石叶片成像

  1. 将叶片化石置于相机下方，并使用沙盘或橡皮泥等材料在化石下方进行固定，确保其尽可能平整放置。
    注意：当在同一岩块上拍摄多个标本时，最好分别对每个标本进行近距离特写拍摄，以确保化石细节清晰锐利。同时建议将化石放置在纯色哑光的深色背景上，例如黑色毛毡或绒布。
  2. 将比例尺水平放置，与叶片处于同一垂直平面内，靠近化石但不得遮挡化石的任何部分。若化石周围基质极少或没有基质，应将比例尺置于照片画面内，并保证其处于焦点之中。
  3. 使用相机三脚架或复制架，将相机正对化石叶片上方，使镜头与岩石表面平行。为确保叶片细节清晰成像，应在保持镜头/相机焦距范围的前提下，尽可能靠近化石拍摄，同时确保整个化石均位于照片取景范围内。
    注意：如条件允许，建议使用高分辨率数码相机，并配备手动对焦的微距镜头以及足够的景深，以便对后续处理所需的叶片结构实现清晰对焦。
  4. 采用间接照明方式，根据需要调整光照，使标本的整个轮廓清晰可见。通常情况下，每件化石均需重新调整照明条件。
  5. 拍摄化石叶片，并对图像文件进行适当命名标注。

2. 数字化制备

注意：本实验方案中使用的叶片结构学术语示意图见图1。请结合决策树（图2）及所提供的示例（图3），判断待测化石叶片适用的制备情形，并转至相应的章节进行操作。有关制备步骤的其他注意事项，请参考表2。若叶片属于情形1或情形5，则无法进行定量叶片形态测量的制备。

  1. 情景2：具有完整边缘的叶片，其面积或一半面积得以保存或可被重建。
    1. 在图像处理软件中打开文件（例如 Adobe Photoshop 或 GIMP）。如有必要，裁剪图像，这有助于减小最终文件大小，但需确保比例尺仍包含在内。
    2. 通过点击 图像 > 画布（Photoshop）；图像 > 画布大小（GIMP），将工作区域的宽度加倍。建议在当前画布的左侧或右侧添加新的画布区域。
    3. 如果叶片边缘需要部分重建，需判断叶片面积和形状是否能更可靠地从半片叶或整片叶进行测量（图3）。
    4. 将叶片从岩石基质中复制出来。使用套索工具描出整片叶或半片叶的轮廓，若存在叶柄也一并描出（见 表2）。复制并粘贴所选区域，并将其放置在画布的空白区域。建议粘贴两份该选区，其中一份保持未编辑状态，以便在需要重新开始制备过程时使用。
    5. 使用适当颜色的线条（通常在白色背景下为黑色）修复边缘的受损部分。绘制一条跨越受损边缘的线，以可靠地重建边缘，例如使用画笔工具或直线工具。确保线条足够粗以便可见（约1–2点粗细），且能连接受损区域两侧的边缘。
    6. 如果存在叶柄，使用套索工具将其从叶片上移除。
      1. 沿叶片基部的边缘目视追踪，直至其与叶柄接触的位置，该位置通常颜色较深且无明显叶脉。在此处放置一个套索点。在叶片的另一半也执行相同操作，并放置第二个点。
        注：如果叶基对称，该连线将大致垂直于叶柄；若不对称，则连线呈一定角度。
      2. 围绕整个叶柄完成选区。剪切并粘贴，或使用移动工具，将叶柄放置在叶片叶片旁但不接触叶片。
        注：对于心形或裂片状叶基，即叶基延伸至叶柄与叶片连接点以下的情况，叶柄可能位于叶基之上。应小心在叶柄实际连接处将其切出，紧密追踪叶柄边缘，并修复由此产生的边缘损伤。需注意，在大多数化石中这可能难以辨识。
    7. 如有必要，裁剪图像的最终区域以减小文件大小。参见 图3 了解完成制备后的图像应呈现的示例。
  2. 情景3：具锯齿的叶片，其面积或一半面积无法重建，但至少保留两个连续的齿以及至少25%的叶片
    注：此类叶片仅能测量齿部特征，因此仅针对这些测量进行制备。
    1. 在图像处理软件中打开文件（例如 Adobe Photoshop 或 GIMP）。如有必要，裁剪图像，这有助于减小最终文件大小，但需确保比例尺仍包含在内。
    2. 通过点击 图像 > 画布（Photoshop）；图像 > 画布大小 （GIMP），将工作区域的宽度扩大至三到四倍。建议在当前画布的左侧或右侧添加新的画布区域。
    3. 将叶片从岩石基质中复制出来。使用套索工具描出保存部分的叶片范围，包括存在的叶柄。无需精确描出边缘受损部分，因为这些部分将被移除。复制并粘贴所选区域，并将其放置在画布的空白区域。建议粘贴两份该选区，其中一份保持未编辑状态，以便在需要重新开始制备过程时使用。
    4. 如果存在叶柄，使用套索工具将其从叶片上移除。
      1. 沿叶片基部边缘目视追踪，直至其与叶柄接触的位置，该位置通常颜色较深且无明显叶脉。在此处放置一个套索点。在叶片的另一半也执行相同操作，并放置第二个点。
        注：如果叶基对称，该连线将大致垂直于叶柄；若不对称，则连线呈一定角度。
      2. 围绕整个叶柄完成选区。剪切并粘贴，或使用移动工具，将叶柄放置在叶片叶片旁但不接触叶片。
        注：对于心形或裂片状叶基，即叶基延伸至叶柄与叶片连接点以下的情况，叶柄可能位于叶基之上。应小心在叶柄实际连接处将其切出，紧密追踪叶柄边缘，并修复由此产生的边缘损伤。在大多数化石中这可能难以辨识。
    5. 使用套索工具移除靠近边缘受损部分的区域。
      1. 从沿边缘界定受损部分的一点开始，绘制一条从该点至主脉且垂直于主脉的直线（图4）。从最靠近损伤的保存完好的主齿缺刻处开始选区。这可确保在后续测量中不会将齿的侧面计入内部周长，也不会将次级齿误测为主齿。如果齿间距较远，此方法可能不适用，因为可能导致过多的保存边缘被移除（图4）。
      2. 沿主脉继续选区，直至与受损边缘另一边界齐平，然后绘制一条垂直于主脉至边缘的直线（图4）。
        注：对于羽状叶脉叶片（图1A），无论是否存在生长抑制性叶脉（参见 Ellis 等人43 的定义和示例），主脉即为主叶脉（即中脉）。对于掌状叶脉叶片（图1B、D），主脉为最近的主叶脉（例如 图4B）。对于羽状裂叶（图1C），若损伤位于羽状裂片上，则主脉为供给该裂片的叶脉（通常为二级叶脉）。
      3. 完成选区并删除该部分叶片。对叶片所有受损部分重复此操作。
    6. 复制并粘贴已制备的叶片，将其放置在画布的空白区域。
    7. 使用套索工具移除齿部。
      1. 从叶片顶端、某一裂片顶端或叶片碎片中最顶端的齿开始，沿该叶片、裂片或碎片上的每个主齿缺刻处进行选区（图5；参见 补充图1、补充图2 了解如何区分主齿与次级齿以及如何区分齿与裂片的技巧）。在切除齿部时务必遵循相应规则（表2；补充图3）。
        注：主齿通常向基部和顶端方向逐渐变小。
      2. 在选中最基部齿的顶端缺刻后，应用延伸规则（表2；补充图4）切除序列中的最后一个齿。
      3. 通过剪切并将齿粘贴至已去除齿的叶片旁（但不接触叶片）来移除齿部。如果制备过程中包含其他需要去除齿部的叶片裂片或碎片，请重复上述步骤，直至所有齿部均被移除。
    8. 如果创建了原始剪切叶片的额外副本，请删除该副本。如有必要，裁剪图像的最终区域以减小文件大小。参见 图3 了解完成制备后的图像应呈现的示例。
  3. 情景4：具锯齿的叶片，其面积或叶片面积得以保存或可被重建
    1. 在图像处理软件中打开文件（例如 Adobe Photoshop 或 GIMP）。如有必要，裁剪图像，这有助于减小最终文件大小，但需确保比例尺仍包含在内。
    2. 通过点击 图像 > 画布（Photoshop）；图像 > 画布大小 （GIMP），将工作区域的宽度扩大至三到四倍。建议在当前画布的左侧或右侧添加新的画布区域。
    3. 确定叶片的制备方式。叶片面积/形状测量需在整片叶或半片叶上进行，需判断哪种方式可获得更准确的测量结果。齿部测量应在所有保存的边缘段上进行。在某些情况下，叶片面积/形状测量所用的叶片子集可能与齿部变量测量所用的子集不同。
      注：在提供的示例（图6）中，判断半片叶比整片叶更可靠地重建。右下角保存的边缘（>1个保存的齿）被用于齿部测量。以下情景4的流程大致遵循该示例（图6），但在不同制备情境下细节可能略有差异。
    4. 将叶片从岩石基质中复制出来，确保包含所有保存的边缘。
      1. 使用套索工具描出叶片边缘，包括存在的叶柄。无需精确描出那些不会被纳入面积/形状测量的受损边缘部分，因为这些部分将被移除（例如 图6 中叶片的右半部分）。
      2. 复制并粘贴所选区域，并将其放置在画布的空白区域。建议粘贴两份该选区，其中一份保持未编辑状态，以便在需要重新开始制备过程时使用。
    5. 如果存在叶柄，使用套索工具将其从叶片上移除。
      1. 沿叶片基部边缘目视追踪，直至其与叶柄接触的位置，该位置通常颜色较深且无明显叶脉。在此处放置一个套索点。在叶片的另一半也执行相同操作，并放置第二个点。
        注：如果叶基对称，该连线将大致垂直于叶柄；若不对称，则连线呈一定角度。围绕整个叶柄完成选区。
      2. 剪切并粘贴，或使用移动工具，将叶柄放置在叶片旁但不接触叶片。
        注：对于心形或裂片状叶基，即叶基延伸至叶柄与叶片连接点以下的情况，叶柄可能位于叶基之上。应小心在叶柄连接处将其切出，紧密追踪叶柄边缘，并修复由此产生的边缘损伤（补充图5）。在大多数化石中这可能难以辨识。
    6. 复制并粘贴已移除叶柄的孤立叶片，创建第二份副本用于齿部测量，并将其放置在画布的空白区域。
    7. 制备一份用于叶片面积和形状测量的叶片版本。
      1. 若制备半片叶，使用套索工具修剪掉多余的叶片部分，仅保留完整的半片叶。若制备整片叶，则不移除任何叶片材料。
      2. 如有必要，使用适当颜色的线条（使用直线或画笔工具，通常在白色背景下为黑色线条）修复边缘的任何受损区域。确保线条足够粗以便可见（约1–2点粗细），且能连接受损区域两侧的边缘。
    8. 制备一份用于齿部测量的叶片版本。
      1. 使用套索工具移除靠近边缘受损部分的区域。
        1. 从沿边缘界定受损部分的一点开始，绘制一条从该点至主脉且垂直于主脉的直线（图4）。从最靠近损伤的保存完好的主齿缺刻处开始选区。这可确保在后续测量中不会将齿的侧面计入内部周长，也不会将次级齿误测为主齿。如果齿间距较远，此方法可能不适用，因为可能导致过多的保存边缘被移除（图4）。
        2. 沿主脉继续选区，直至与受损边缘另一边界齐平，然后绘制一条垂直于主脉至边缘的直线（图4）。
          注：对于羽状叶脉叶片（图1A），无论是否存在生长抑制性叶脉（参见 Ellis 等人43 的定义和示例），主脉即为主叶脉（即中脉）。对于掌状叶脉叶片（图1B、D），主脉为最近的主叶脉（例如 图4B）。对于羽状裂叶（图1C），若损伤位于羽状裂片上，则主脉为供给该裂片的叶脉（通常为二级叶脉）。
        3. 删除叶片的受损部分。对叶片所有受损部分重复此操作。
    9. 复制并粘贴已去除受损部分、为齿部测量准备的版本，并将其放置在画布的空白区域。
    10. 使用套索工具移除齿部。
      1. 从叶片顶端、某一裂片顶端或叶片碎片中最顶端的齿开始，沿叶片、裂片或碎片上的每个主齿缺刻处进行选区（图5；参见 补充图2 了解如何区分主齿与次级齿以及如何区分齿与裂片的技巧）。在切除齿部时务必遵循相应规则（表2；补充图2）。
        注：主齿通常向基部和顶端方向逐渐变小。
      2. 在选中最基部齿的顶端缺刻后，应用延伸规则（表2；补充图4）切除序列中的最后一个齿。
      3. 通过剪切并将齿粘贴至已去除齿的叶片旁（但不接触叶片）来移除齿部。如果制备过程中包含其他需要去除齿部的叶片裂片或碎片，请重复上述步骤，直至所有齿部均被移除。
    11. 如果创建了原始剪切叶片的额外副本，请删除该副本。如有必要，裁剪图像的最终区域以减小文件大小。参见 图3 了解完成制备后的图像应呈现的示例。

3. 数字化测量

注意：数据录入模板表格作为补充文件1提供。有关测量步骤的其他注意事项，请参考表3。在情景1和情景5中，唯一需要的步骤是在数据录入表格中记录叶缘状态（步骤3.5）。

  1. 打开 ImageJ 软件44。设置将自动进行的测量项目（在程序安装后仅需执行一次）。
  2. 点击 分析 > 设置测量值，仅选择“面积”、“周长”和“Feret 直径”。确保小数位数设置为 3 位。
  3. 通过点击 文件 > 打开 或直接将图像拖放到已打开的 ImageJ 工具栏中，打开准备好的化石叶片图像。
  4. 为每一张新的叶片图像设置比例尺。
    注意：这是关键步骤，必须对每一张新的叶片图像执行此操作，以确保测量结果准确。
    1. 选择直线工具。放大比例尺条，尽可能沿比例尺画一条最长的直线。
    2. 点击 分析 > 设置比例尺。在“已知距离”栏中输入以厘米为单位的测量长度（以与现代校准数据集所用单位保持一致）。无需更改长度单位。点击 确定。
  5. 在数据录入电子表格中标记叶片为具齿（0）或全缘（1）。
  6. 如果存在叶柄，则测量叶柄宽度。测量应在仍保留在岩基质中的原始叶片标本上进行，因为这能提供更完整的背景信息。
    注意：若叶柄缺失，在某些情况下可测量中脉最基部位置的宽度以代替叶柄宽度。但仅当整个中脉宽度均被保存时才可进行此操作（即没有叶片膜质部分压缩在叶脉之上，或化石保存了叶片的远轴面），且该物种或形态型的其他标本需表明基部叶脉宽度与叶柄宽度相当。
    1. 在叶柄与叶片连接处画一条垂直于叶柄的直线；若连接点不对称，则在最基部的连接点处画一条垂直于叶柄的直线。此直线必须仔细绘制，建议放大该区域以便更精确地绘线。
      注意：在某些特殊情况下需调整此步骤，例如在最基部连接点处存在损伤、毛状体、蜜腺、刺或其他影响叶柄宽度准确测量的特征。此时应在该特征下方第一个可可靠测量的位置测量叶柄宽度。
    2. 点击 分析 > 测量，或使用快捷键测量所画直线的长度。再次绘制同一条线以记录测量位置，可通过点击 编辑 > 绘图 或使用快捷键实现。如有必要，使用主工具栏上的颜色选择器更改线条颜色。
    3. 完成绘线后，保存图像，建议使用修改后的文件名。
    4. 在数据录入电子表格的“叶柄宽度”项下记录该直线的长度。
  7. 通过将图像转换为黑白图像来准备叶片以进行后续测量。操作方法为点击 图像 > 类型 > 8 位。
  8. 点击 图像 > 调整 > 阈值 或使用快捷键对图像进行阈值化处理。将弹出标题为“阈值”的对话框，并使图像部分区域变为红色。若叶片颜色较浅而背景较暗，请点击 暗色背景。
  9. 使用滑动条调整阈值，直到叶片内部完全呈红色且与背景清晰区分。这是关键步骤，容易因操作不当导致数据不准确。务必通过放大叶片边缘部分确认红色区域恰好覆盖整个叶片（即叶片周长全部为红色且无多余区域）。叶片内部出现红色间隙可接受，不影响测量结果。
    注意：若叶片轮廓不清晰，首先尝试在放大状态下调整阈值，确认叶片轮廓完全呈红色。若化石与背景对比度过低导致无法可靠应用阈值，可使用画笔工具在对比度不足的区域为叶片边缘添加实线轮廓。或者，将图像返回图像处理软件（如 Adobe Photoshop 或 GIMP），调整孤立叶片图层的对比度或背景颜色，以增强区分度。
  10. 对按情景 2（步骤 2.1）和情景 4（步骤 2.4）准备的叶片测量叶面积和叶形。
    注意：参考 图 2 和 图 6B，确定在图像的哪些部分测量相应变量。若叶片按情景 3（步骤 2.2）准备，则跳过此步骤，直接进入步骤 3.11。
    1. 测量为叶面积和叶形分析准备的叶片（应仅去除叶柄，若存在叶柄）。选择 魔棒 工具，点击叶片内部，整个叶片应被黄色轮廓线包围——请确认轮廓正确。
    2. 点击 分析 > 测量 或使用快捷键进行测量。
    3. 若测量区域为完整叶片，在数据录入电子表格中记录面积、周长、Feret 和最小 Feret 值；若测量区域为半片叶片，则仅记录 Feret 值并进入下一步。
    4. 若测量区域为半片叶片，测量叶片的人工中线周长，即叶片被切成两半后形成的人工边缘长度（图 6B）。若测量区域为完整叶片，则跳过此步骤及步骤 3.10.5。
      注意：测量人工中线周长可用于从半片叶片计算叶片周长（见下文步骤 3.10.5）。叶片周长不用于 DiLP 和 MA 分析中的变量，但可用于其他有助于形态特征描述的变量（如形状因子、紧凑度；表 1）。
      1. 右键点击 直线 工具，选择“分段线”工具，沿人工中线周长完整描迹。
      2. 点击 分析 > 测量 或使用快捷键测量其长度。该测量值将用于下文公式计算叶片周长（步骤 3.10.5）。
    5. 若测量区域为半片叶片，在录入数据时对测量值进行修正：将面积乘以 2，最小 Feret 值乘以 2，并使用以下公式计算叶片周长：先从半片叶片周长中减去人工中线周长，再将结果乘以 2：
      ​叶片周长 = （周长 - 人工中线周长）× 2
    6. 若存在切除的叶柄，测量其面积；若不存在，则情景 2 的测量已完成，但情景 4 需继续进入下一步。
      1. 使用魔棒工具点击切除的叶柄，其应被黄色轮廓线包围。点击 分析 > 测量 或使用快捷键进行测量。在数据录入电子表格的“叶柄面积”项下记录面积。
    7. 对于情景 2（步骤 2.1），测量已完成；对于情景 4（步骤 2.3），进入下一步。
  11. 对按情景 3（步骤 2.2）和情景 4（步骤 2.3）准备的叶片测量齿相关变量。
    1. 测量原始叶片。使用魔棒工具选择原始叶片内部（即仍保留叶齿、用于齿测量的叶片；图 6B），其应被黄色轮廓线包围。点击 分析 > 测量 或使用快捷键进行测量。
      注意：根据叶片准备方式，可能需要测量多个不连续区域，并将面积相加（如 图 6B）。或者，可按住 Shift 键并使用魔棒工具同时选择多个区域。
    2. 在数据录入电子表格的“原始叶片面积”和“原始叶片周长”项下记录面积和周长。
    3. 测量内部原始叶片。选择内部原始叶片内部（即已去除叶齿、用于齿测量的叶片；图 6B），其应被黄色轮廓线包围。点击 分析 > 测量 或使用快捷键进行测量。
      注意：根据叶片准备方式，可能需要测量多个不连续区域，并将面积相加（如 图 6B）。或者，可按住 Shift 键并使用魔棒工具同时选择多个区域。
    4. 在数据录入电子表格的“内部原始叶片面积”和“内部原始叶片周长”项下记录面积和周长。
    5. 测量切割周长。为清晰观察叶片，需先移除阈值：点击阈值对话框中的 重置，或点击 编辑 > 撤销——后者通常也会移除黑白 8 位转换。选择分段线工具，沿原始叶片上的切割周长完整描迹。
    6. 点击 分析 > 测量 或使用快捷键进行测量。若存在多个部分，重复前述步骤测量各部分的切割周长。在数据录入电子表格的“切割周长长度”项下记录长度或总长度。
      注意：切割周长是通过去除损伤组织准备叶片时引入的。大多数情况下，其不同于人工中线周长（图 6B）。
    7. 统计主齿和次级齿的数量（若存在）。
      注意：参见 补充图 2，了解如何区分主齿和次级齿。
      1. 若尚未移除阈值，现在移除。点击阈值对话框中的 重置，或点击 编辑 > 撤销——后者通常也会移除黑白 8 位转换。
      2. 统计主齿数量（参见 补充图 2 了解区分主齿和次级齿的技巧）。选择 多点工具，可能需要先右键点击 点工具 以选择 多点工具。点击每个主齿进行编号。
      3. 若误选某点，按住 Alt 键（Windows 系统）或 Command/cmd 或 option 键（Mac 系统）并同时点击该点以删除。在数据录入电子表格的“主齿数量”项下记录最终数值。
      4. 点击 编辑 > 选择 > 取消选择 或使用键盘快捷键，清除多点工具的计数和标注。
      5. 统计总齿数（即叶片上所有主齿和次级齿的总数）。选择 多点工具，可能需要先右键点击 点工具 以选择 多点工具。点击每个齿（包括主齿和次级齿）进行编号。
        注意：统计总齿数而非次级齿数，可避免重复计数。通过从总齿数中减去主齿数即可得到次级齿数（见步骤 3.11.7.6）。
      6. 若误选某点，按住 Alt 键（Windows 系统）或 Command/cmd 或 option 键（Mac 系统）并同时点击该点以删除。用总齿数减去主齿数，得到次级齿数。在数据录入电子表格的“次级齿数量”项下记录该数值。
        注意：部分用户倾向于在准备叶片图像时进行齿数统计，而非在测量阶段进行。

4. 在 R 软件中运行分析

注意：以下步骤需要使用 R 软件包 dilp11。数据录入电子表格将被读入 R 并由该软件包使用。请参考数据录入电子表格中的“附加说明”标签页（补充文件 2）。该 R 脚本可同时分析多个位点，也可用于单个位点的分析。

  1. 使用您偏好的环境打开 R（推荐使用 R Studio）。有关 R 的入门介绍，可参考例如 https://cran.r-project.org/doc/manuals/r-release/R-intro.pdf
  2. 在 R 会话中安装 dilp 包. 有关如何安装该软件包及其相关函数的更多信息，请参见以下网站： https://cran.r-project.org/package=dilp
  3. 读取包含化石木本双子叶被子植物叶片性状数据的 .csv 文件（即在数据录入电子表格中记录的数据）。
  4. 运行该函数 dilp() 用于年均温（MAT）和年均降水（MAP）的重建及其相关误差。MAT 和 MAP 的结果来自多元线性回归模型（MLR；即 DiLP）以及两个单变量线性回归模型（SLR；即叶面积和叶缘分析）。运行该函数 lma() 单位叶面积质量（MA在形态类型和位点水平上的重建。
  5. 运行后 dilp() 建议通过检查返回结果中的离群值和错误对象，来排查潜在的数据采集问题并确认数据质量 dilp() 结果。或者，使用该功能 dilp加工() 随后 dilp_outliers() 和 dilp_errors()参考已制备的标本，对任何被标记的问题进行核查，必要时可重新测量。建议编辑原始数据文件，然后将其重新读入 R 中。
  6. 使用该函数判断化石地点是否位于校准数据集的叶片形态多变量空间范围内 dilp_cca().

结果

本研究采用先前发表的来自不列颠哥伦比亚省中南部早始新世麦卡比（McAbee）化石遗址的叶片形态测量数据集，以展示数字叶片形态学（DiLP）和单位面积叶片质量（MA）两种重建方法的代表性结果（Lowe 等人38；数据见补充文件2）。该遗址为在高地火山地貌背景下重建新生代最温暖时期（即早始新世气候最适宜期）的古气候与古生态提供了重要机会38,45,46,47。化石组合采自湖相沉积序列中两个独立的层位，分别命名为H1（厚28 cm）和H2（厚27 cm），通过普查法在狭窄的地层范围内进行整合取样，即收集或统计所有可归入特定形态类型的标本38,48。

McAbee 叶片形态数据通过了 dilp_errors() 标记的误差检查，经 dilp_outliers() 标记的七个离群值也进行了复核，以确保这些数值反映的是数据中真实的变异，而非方法学上的错误。随后，数据通过 dilp() 函数用于古气候重建，并通过 lma() 函数进行单位面积叶片质量的重建。

MA 重建结果及其95%预测区间的下限和上限已报告于 表4 在物种和位点水平上，使用Royer等人提出的方法进行计算。5 以及 Butrim 等人10重建值在 M 的范围内A 现代陆生物种的典型值（30-330 g/m）2)49. 使用 Royer 等人所述的阈值5，大多数物种具有重建的 MA 与叶片寿命相一致的 <1年（≤87克/米）2），约1年（88–128 克/平方米）2），而没有一种是典型的 >1年（≥129克/米）2）。位点M的重建A 均值和方差在McAbee反映了该地点叶片经济策略的普遍性与多样性10,50H1与H2时期在位点均值和方差之间无显著差异，因此没有证据表明这两个时间段内叶片经济策略的组成和多样性存在变化。此外，使用Royer等人提出的方程进行的位点均值重建结果也支持这一结论。5 以及 Butrim 等人10 非常相似。

使用Peppe等人6中提出的多元线性回归（DiLP）以及单一线性回归（叶缘和叶面积分析）方程对年均温（MAT）和年均降水量（MAP）进行重建的结果见表5。若化石叶片组合的叶片形态特征落在校准数据集的形态空间范围内，则古气候估算结果最为可靠。这一评估通过函数dilp_cca()执行的典范对应分析（CCA）步骤完成。McAbee H1和H2均位于校准数据集中观察到的叶片形态范围之内（图7A）。如果某遗址的重建值超出校准空间范围，则应谨慎解释其古气候重建结果（例如，需结合独立的证据进行比对；详见Peppe等人6的进一步讨论）。H1和H2的MAT与MAP重建结果均表明其属于温带季节性生物群系（图7B、C），这与独立证据高度一致，包括基于McAbee化石植物群和昆虫群最近现生亲缘种所推断的结果45。

化石叶片形态：羽状脉和掌状脉，具裂片，无裂片示意图。
图 1：本文中使用的叶片形态及结构术语。（A）羽状脉、无裂片、全缘叶；（B）掌状脉、无裂片、具锯齿叶；（C）羽状脉、具裂片、全缘叶；（D）掌状脉、具裂片、具锯齿叶。请点击此处查看此图的放大版本。

叶片形态分析流程图；完整边缘与锯齿边缘情景；叶片测量方法。
图 2：方法流程图。 该流程图展示了不同的叶片保存状态和叶片类型如何决定可可靠测量的叶片性状类型（黄色框）。这将决定本实验方案中采用的制备情景，以及数据录入表格中应填写数据的列（项目符号所示）。请点击此处查看此图的放大版本。

化石叶片复原情景；示意图展示各种制备方法及结果。
图3：不同的制备情景。 不同的制备情景，展示已完成数字化制备、可用于测量阶段的图像示例。(A) 情景1，全缘叶，其整体面积或半叶面积无法重建；(B) 情景5，具锯齿叶，其整体面积或半叶面积无法重建，且不具备≥2个连续锯齿和/或保存叶片比例不足25%；(C) 情景2，全缘叶，其整体面积或半叶面积得以保存或可被重建；(D) 情景3，具锯齿叶，其整体面积或半叶面积无法重建，但具有≥2个连续锯齿且保存叶片比例≥25%；(E) 情景4，具锯齿叶，其整体面积或半叶面积得以保存或可被重建。请点击此处查看此图的放大版本。

化石叶片重建示意图；叶片形态与模式分析，含比例尺参考。
图4：损伤去除示意图。 展示如何切除受损边缘以及邻近该受损边缘的叶片区域。红色虚线表示使用套索工具进行选择的方式。请注意，损伤边界的起始位置有意选在主齿的缺刻处（参见补充图2以帮助区分主齿与次级齿）。(A) 一种羽状脉叶片，其选择范围延伸至中脉。(B) 一种掌状脉叶片，其选择范围延伸至最近的主脉。请点击此处查看该图的放大版本。

叶脉模式分析示意图，展示分形几何结构及放大后的叶脉细节。
图5：展示如何裁剪出叶齿的示例。（A）红色虚线表示如何使用套索工具进行选择。请注意，本例中叶齿为复合型，因此仅在主级缺刻之间进行选择（参见补充图2以帮助区分主级叶齿与次级叶齿）；（B）叶齿选择方式的放大视图，红色圆点表示选择过程中鼠标点击的位置；（C）去除叶齿后叶片的复制图像。请点击此处查看该图的放大版本。

化石叶片测量示意图，记录面积、形状和叶齿参数，用于古植物学分析。
图6：制备方案4的示意图。 以按方案4制备的一个叶片为例，展示其制备决策和测量步骤。(A) 一种制备方案，其中决定采用半片叶片以获得最可靠的叶形和叶面积测量结果，并保留两侧中段的叶缘用于叶齿参数测量。(B) 一个示例，展示在制备后的叶片各部分所测量的不同变量。加粗文字标示了用于DiLP和MA分析所需的测量参数，而非加粗文字（叶片周长、最小费雷特直径和人工中周长）则标示了虽非必需但对其他形态特征描述（例如形状因子和紧密度；表1）有帮助的测量参数。请点击此处查看该图的放大版本。

气候数据分析；年均降水量与温度关系；多色生物群落带图。
图7：代表性结果。 来自Lowe等人38在早始新世McAbee化石层采集的两个化石层位（H1和H2）的结果。（A）典型对应分析，显示校准数据集中多变量叶片形态特征的表现。校准数据来自Peppe等人6。两个McAbee层位的叶片形态特征叠加显示，位于校准空间之内。（B和C）利用Peppe等人6提出的方程估算的温度和降水量及其相关不确定性（模型标准误差），并将两个McAbee层位的结果叠加于Whittaker生物群落图上。（B）使用数字叶片形态学（DiLP）多元线性回归模型（MLR）重建的估算值，（C）使用叶片面积分析（LAA）和叶片边缘分析（LMA）单一线性回归（SLR）方程重建的估算值，并将两个McAbee层位的结果叠加于Whittaker生物群落图上。 请点击此处查看该图的放大版本。

表1：叶片形态学变量。 本方案中测量和/或计算并用于预测模型中的变量，以重建叶片单位面积干重（MA）、年均温（MAT）和年均降水量（MAP）。MAT和MAP采用Peppe等6提出的方程进行重建，其中数字叶片形态学（DiLP）采用多元分析方法，叶缘分析（LMA）和叶面积分析（LAA）采用单变量分析方法。标记为“其他”的变量未用于MA、DiLP、LMA和LAA分析中，但仍根据本方案进行测量和计算，因其易于实施并可提供对叶片形态特征的有用描述。请点击此处下载该表格。

表2： 准备步骤的其他注意事项和说明。请点击此处下载该表格。

表3： 测量步骤的其他注意事项和说明。 请点击此处下载该表格。

表4：基于Lowe等人38的研究对McAbee化石层位的叶片单位面积干重（MA）及其95%预测区间上下限的重建结果。重建分别基于形态型均值5、样地均值5,10以及样地方差10进行。请点击此处下载该表格。

表5：利用Peppe等人6提出的数字叶形态学（DiLP）的多元线性回归（MLR）以及叶缘分析（LMA）和叶面积分析（LAA）的单一线性回归（SLR），对中国台湾早始新世McAbee化石层位H1和H2层的年均温度（MAT）与年均降水量（MAP）进行的重建结果。请点击此处下载该表格。

补充图 1：哈佛森林的红栎（Quercus rubra）叶片，展示裂片与锯齿的判别规则。 文中已定义线段 p 和 d。比例尺 = 1 cm。请点击此处下载该文件。

补充图 2：哈佛森林的Betula lutea叶片，展示区分次级齿与主齿的规则。 图中孤立的叶片部分已放大 2 倍。蓝色线条连接切入程度最大的叶凹（即主级叶凹），与这些叶凹相关的齿视为主齿（蓝色箭头）。红色圆点标记的齿可被识别为次级齿，因其顶端叶凹的切入程度较小。红色箭头所指的齿与主齿的切入程度相似，但可通过其相对较细的主脉来判断为次级齿。比例尺 = 1 cm。请点击此处下载该文件。

补充图3：牙齿状突起选择、羽状裂片规则及裂片优先规则示意图。（A）来自Huyck保护区的Hamamelis virginiana叶片的牙齿状突起选择。深色区域代表被纳入总体牙齿状突起统计的叶组织，因为次级齿与主齿已作区分。（B）来自IES的Quercus alba叶片展示了裂片优先规则。深色区域被测量为裂片，未着色区域被测量为牙齿状突起，但根据裂片优先规则，所有突出部分均被视为裂片。比例尺 = 1 cm。请点击此处下载该文件。

补充图4： Acer saccharum 来自阿勒格尼国家森林的叶片，用于说明延伸规则和单齿规则。 虚线表示所选牙齿，实线表示相应牙齿的对称轴。黑色区域为用于压平叶片以便拍照的配重。比例尺 = 1 cm。 请点击此处下载此文件。

补充图5：展示如何从心形基部上方切割叶柄的正确方法。 请点击此处下载该文件。

补充文件 1：所有测量的数字叶形态变量的数据录入模板。此文件不应被修改，因为它将用作 R 软件包的输入文件。请点击此处下载该文件。

补充文件 2：来自 Lowe 等人38研究中 McAbee 化石层的示例数据。 该数据用于生成图 7，并用于代表性结果的讨论。请点击此处下载该文件。

补充文件3：化石数字化叶片形态学规则文档。 请点击此处下载该文件。

讨论

本文介绍了如何在木本双子叶被子植物化石叶片上测量叶片形态的连续性状，并将其应用于基于现代校准数据建立的代用指标，以重建古气候和古生态。为此，必须确保方法学步骤与代用指标校准数据集中的方法保持一致5,6,10。这一考虑应从化石叶片采集阶段即开始，尤其涉及样本量问题。建议尽可能在较窄的地层范围内合并化石叶片组合，以获得足够数量的可测量标本和形态类型，从而最大限度地减少时间平均效应。同时建议，将古气候重建限定在具有至少350个可鉴定标本以及至少15–20种木本双子叶被子植物形态类型的遗址上19,51,52。此外，在选择用于分析的叶片时，建议尽可能多地测量每种形态类型的叶片数量，并至少选取能够代表该形态类型内叶片形态变异性的标本。

在实施制备和测量步骤时，必须格外注意与校准数据集保持一致。制备阶段的操作最有可能因操作者不同而引入主观性并导致结果差异。然而，若仔细遵循本方案，并经常参考附加注意事项表格（表2、表3）及规则文档（补充文件3），该方法可获得客观且可重复的叶片形态学测量结果。对于初次使用该方法的用户，建议由经验更丰富的人士协助确认叶片是否已正确制备。在测量叶柄宽度以用于MA重建时需特别谨慎。由于这些数值会被平方，测量中的任何不准确都会被放大。不完整的保存状态或样品损伤可能改变叶柄的尺寸，应仔细避免。

这些方法存在一些值得注意的局限性。最重要的是，dilp R 软件包中包含的代用指标仅适用于木本双子叶被子植物，因此可能排除了古代群落中其他重要的植物类群。然而，已有针对具叶柄的广义裸子植物5,8、草本被子植物8以及蕨类植物9的基于叶片叶柄的物种级MA代用指标发表，用户可根据需要单独引入这些方法。在木本双子叶被子植物以外的植物类群缺失的情况下，对样地水平MA均值和方差的重建影响可能最为显著，因为这将导致对整个群落内经济策略的评估不完整。系统发育历史会影响叶片锯齿的出现23，因此分析具有新颖分类组成的化石群落时，可能会给估算结果带来不确定性，尽管这种潜在影响尚未经过检验和证实。

化石叶片还需要得到充分保存，才能进行超出叶缘状态的叶片形态计量学定量测量。对于DiLP方法而言，这一点对于全缘叶尤其重要，因为只有当整个叶片或半个叶片被保存下来或可被重建时，它们才能提供除叶缘状态之外的其他信息。同样，叶片只有在满足以下条件时才能用于MA重建分析：（1）叶柄在叶片基部的着生部位被保存下来，或在特定情况下，叶片基部及中脉最基部区域被保存（见步骤3.6中的注释）；（2）能够估算叶片的大小，可通过完整叶片测量或半叶重建实现。这意味着某些形态类型可能完全无法纳入站点级别的MA分析。最后，时间是本实验方案的一个限制因素，因为用于古气候重建的单变量替代方法所需时间相对较短。

尽管存在这些局限性，DiLP 和 M 的使用仍然具有重要意义A 重建方法相较于其他方法仍具有若干优势。MA 重建是化石记录中复原叶片经济策略的少数可行方法之一，利用二维叶柄宽度和叶面积测量值，可基于常见的印痕/压型叶化石进行重建。对于DiLP方法而言，引入多个与气候功能相关联的连续性测量指标，有助于提高测量的可重复性以及气候重建结果的准确性6,13本方案旨在通过利用叶片碎片进行叶缘锯齿度测量，以应对化石记录的不完整性。尽管叶面积的连续测量可提供更多关于叶片大小的信息，但为增加样本量，可结合使用叶片大小分级的方法来补充基于连续测量的DiLP MAP估算结果。16,53 或通过结合叶片面积的叶脉比例估算42,54,55与大多数复杂方法一样，随着使用者在操作中积累经验并增强信心，该方案的时间效率将会提高，尤其是在准备步骤方面。按照本方案进行位点水平DiLP测量的事实已表明 >150 现代6,10,56 迄今已有至少22个化石组合证实了其可行性6,38,39,40,41,42最后，叶片形态的全面测量具有超出本文讨论范围的应用价值，可有助于描述植物生态学、生理学、进化及发育的其他方面，并适用于现代56 以及古生态学研究40.

总之，实施本文详述的方法可使用户利用可靠且可重复的方法重建古气候与古生态。这些方法为展示气候和生态系统对环境扰动的过往响应实例提供了重要机会，并有助于进一步揭示地球自然系统之间复杂的相互作用。

致谢

AJL 感谢华盛顿大学 2020-2022 年本科 Team Leaf 团队为开发 DiLP 有效培训材料所提供的动力和建议。AGF、AB、DJP 和 DLR 感谢卫斯理大学和贝勒大学的众多本科生，他们测量了现生和化石叶片，其反馈对于修改和更新本实验方案至关重要。作者感谢 PBot 定量性状工作组及 PBOT 团队鼓励开展本协议的规范化工作，使其更广泛地惠及科研群体。本研究得到了美国国家科学基金会（DLR 获资助项目 EAR-0742363，DJP 获资助项目 EAR-132552）和贝勒大学（青年研究者发展计划资助 DJP）的支持。我们感谢两位匿名审稿人及评审编辑所提出的反馈意见，这些意见有助于提升本实验方案的清晰性和全面性。

材料

本文使用的材料清单
姓名公司目录编号评论
复制架或三脚架用于化石摄影
数码相机用于化石摄影，建议使用高分辨率相机
图像编辑软件 用于数字处理。示例包括 Adobe Photoshop 和 GIMP，后者为免费软件（https://www.gimp.org/）
ImageJ 软件IJ1.46pr用于进行数字测量，免费软件（https://imagej.net/ij/index.html）
Microsoft ExcelMicrosoft或类似的数据录入软件
R 软件The R Foundation用于运行提供的 R 脚本（https://www.r-project.org/）。RStudio 提供用户友好的 R 环境（https://posit.co/download/rstudio-desktop/）。两者均为免费软件。
dilp R 程序包可按照此处说明进行安装：https://github.com/mjbutrim/dilp

参考文献

  1. Gates, D. M. Transpiration and leaf temperature. Ann Rev Plant Physiol. 19 (1), 211-238 (1968).
  2. Givnish, T. On the adaptive significance of leaf form. Topics Plant Pop Biol. , 375-407 (1979).
  3. Niklas, K. J. Plant biomechanics: an engineering approach to plant form and function. , University of Chicago Press. (1992).
  4. Niklas, K. J. Plant Allometry: The Scaling of Form and Process. , University of Chicago Press. (1994).
  5. Royer, D. L., et al. Fossil leaf economics quantified: Calibration, Eocene case study, and implications. Paleobiology. 33 (4), 574-589 (2007).
  6. Peppe, D. J., et al. Sensitivity of leaf size and shape to climate: global patterns and paleoclimatic applications. New Phytol. 190 (3), 724-739 (2011).
  7. Wright, I. J., et al. The worldwide leaf economics spectrum. Nature. 428 (6985), 821-827 (2004).
  8. Royer, D. L., Miller, I. M., Peppe, D. J., Hickey, L. J. Leaf economic traits from fossils support a weedy habit for early angiosperms. Am J Botany. 97 (3), 438-445 (2010).
  9. Peppe, D. J., et al. Biomechanical and leaf-climate relationships: a comparison of ferns and seed plants. Am J Botany. 101 (2), 338-347 (2014).
  10. Butrim, M., et al. dilp: Reconstruct paleoclimate and paleoecology with leaf physiognomy. R package version 1.1.0. , https://cran.r-project.org/web/packages/dilp/index.html (2024).
  11. Bailey, I. W., Sinnott, E. W. A botanical index of Cretaceous and Tertiary climates. Science. 41 (1066), 831-834 (1915).
  12. Bailey, I. W., Sinnott, E. W. The climatic distribution of certain types of angiosperm leaves. Am J Botany. 3 (1), 24-39 (1916).
  13. Peppe, D. J., Baumgartner, A., Flynn, A. G., Blonder, B. Reconstructing paleoclimate and paleoecology using fossil leaves. Meth Paleoecol: Reconstructing Cenozoic Terre Environ Ecol Comm. , 289-317 (2018).
  14. Wolfe, J. A. Temperature parameters of humid to mesic forests of eastern Asia and relation to forests of other regions of the Northern Hemisphere and Australasia. , U.S. Govt. Print. Off. (1979).
  15. Wilf, P. When are leaves good thermometers? A new case for Leaf Margin Analysis. Paleobiology. 23 (3), 373-390 (1997).
  16. Wilf, P., Wing, S. L., Greenwood, D. R., Greenwood, C. L. Using fossil leaves as paleoprecipitation indicators: An Eocene example. Geology. 26 (3), 203-206 (1998).
  17. Miller, I. M., Brandon, M. T., Hickey, L. J. Using leaf margin analysis to estimate the mid-Cretaceous (Albian) paleolatitude of the Baja BC block. Earth Planetary Sci Lett. 245 (1), 95-114 (2006).
  18. Wing, S. L., Greenwood, D. R. Fossils and fossil climate: The case for equable continental interiors in the Eocene. Philosophical Transact: Biol Sci. 341 (1297), 243-252 (1993).
  19. Wolfe, J. A. A method of obtaining climatic parameters from leaf assemblages. U S Geol Sur. (2040), 1-71 (1993).
  20. Huff, P. M., Wilf, P., Azumah, E. J. Digital future for paleoclimate estimation from fossil leaves? Preliminary results. Palaios. 18 (3), 266-274 (2003).
  21. Spicer, R. A. Recent and future developments of CLAMP: Building on the legacy of Jack A. Wolfe. Adv Angiosperm Paleobotany Paleoclimatic Reconstruct. 258, 109-118 (2007).
  22. Yang, J., Spicer, R. A., Spicer, T. E. V., Li, C. S. CLAMP Online': a new web-based palaeoclimate tool and its application to the terrestrial Paleogene and Neogene of North America. Palaeobiodivers Palaeoenviron. 91 (3), 163-183 (2011).
  23. Little, S. A., Kembel, S. W., Wilf, P. Paleotemperature proxies from leaf fossils reinterpreted in light of evolutionary history. PLoS One. 5 (12), e15161(2010).
  24. Webb, L. J. A physiognomic classification of Australian rain forests. J Ecol. 47 (3), 551-570 (1959).
  25. Schmerler, S. B., et al. Evolution of leaf form correlates with tropical-temperate transitions in Viburnum (Adoxaceae). Proc Royal Soci B: Biol Sci. 279 (1744), 3905-3913 (2012).
  26. Nicotra, A. B., et al. The evolution and functional significance of leaf shape in the angiosperms. Funct Plant Biol. 38 (7), 535-552 (2011).
  27. Leigh, A., Sevanto, S., Close, J. D., Nicotra, A. B. The influence of leaf size and shape on leaf thermal dynamics: does theory hold up under natural conditions. Plant Cell Environ. 40 (2), 237-248 (2017).
  28. Wright, I. J., et al. Global climatic drivers of leaf size. Science. 357 (6354), 917-921 (2017).
  29. Givnish, T. J. Leaf and canopy adaptations in tropical forests. Physiol Ecol Plants Wet Tropics. 12, 51-84 (1984).
  30. Feild, T. S., Sage, T. L., Czerniak, C., Iles, W. J. Hydathodal leaf teeth of Chloranthus japonicus (Chloranthaceae) prevent guttation-induced flooding of the mesophyll. Plant, Cell Environ. 28 (9), 1179-1190 (2005).
  31. Royer, D. L., Wilf, P. Why do toothed leaves correlate with cold climates? Gas exchange at leaf margins provides new insights into a classic paleotemperature proxy. Int J Plant Sci. 167 (1), 11-18 (2006).
  32. Edwards, E. J., Spriggs, E. L., Chatelet, D. S., Donoghue, M. J. Unpacking a century-old mystery: Winter buds and the latitudinal gradient in leaf form. Botanical Soc Am. 103 (6), 975-978 (2016).
  33. Givnish, T. J., Kriebel, R. Causes of ecological gradients in leaf margin entirety: Evaluating the roles of biomechanics, hydraulics, vein geometry, and bud packing. Am J Botany. 104 (3), 354-366 (2017).
  34. Royer, D. L., Wilf, P., Janesko, D. A., Kowalski, E. A., Dilcher, D. L. Correlations of climate and plant ecology to leaf size and shape: potential proxies for the fossil record. Am J Botany. 92 (7), 1141-1151 (2005).
  35. Green, W. A. Loosening the CLAMP: an exploratory graphical approach to the climate leaf analysis multivariate program. Palaeontol Electro. 9 (2), 19(2006).
  36. Peppe, D. J., Royer, D. L., Wilf, P., Kowalski, E. A. Quantification of large uncertainties in fossil leaf paleoaltimetry. Tectonics. 29 (3), 002549(2010).
  37. Royer, D. L., Meyerson, L. A., Robertson, K. M., Adams, J. M. Phenotypic plasticity of leaf shape along a temperature gradient in Acer rubrum. PLoS One. 4 (10), e7653(2009).
  38. Lowe, A. J., et al. Plant community ecology and climate on an upland volcanic landscape during the early Eocene climatic optimum: McAbee fossil beds, British Columbia, Canada. Palaeogeograp, Palaeoclimatol, Palaeoecol. 511, 433-448 (2018).
  39. Flynn, A. G., Peppe, D. J. Early Paleocene tropical forest from the Ojo Alamo Sandstone, San Juan Basin, New Mexico, USA. Paleobiology. 45 (4), 612-635 (2019).
  40. Allen, S. E., Lowe, A. J., Peppe, D. J., Meyer, H. W. Paleoclimate and paleoecology of the latest Eocene Florissant flora of central Colorado, U.S.A. Palaeogeograp, Palaeoclimatol, Palaeoecol. 551, 109678(2020).
  41. Baumgartner, A., Peppe, D. J. Paleoenvironmental changes in the Hiwegi Formation (lower Miocene) of Rusinga Island, Lake Victoria, Kenya. Palaeogeograp, Palaeoclimatol, Palaeoecol. 574, 110458(2021).
  42. Wagner, J. D., Peppe, D. J., O'keefe, J. M., Denison, C. N. Plant community change across the Paleocene-Eocene boundary in the Gulf coastal plain, Central Texas. Palaios. 38 (10), 436-451 (2023).
  43. Ellis, B., et al. Manual of leaf architecture. , Cornell University Press. Ithaca, NY. (2009).
  44. Schneider, C. A., Rasband, W. S., Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat Meth. 9, 671-675 (2012).
  45. Greenwood, D. R., Pigg, K. B., Basinger, J. F., DeVore, M. L. A review of paleobotanical studies of the Early Eocene Okanagan (Okanogan) Highlands floras of British Columbia, Canada, and Washington, USA. Canadian J Earth Sci. 53 (6), 548-564 (2016).
  46. Lowe, A. J., West, C. K., Greenwood, D. R. Volcaniclastic lithostratigraphy and paleoenvironment of the lower Eocene McAbee fossil beds, Kamloops Group, British Columbia, Canada. Canadian J Earth Sci. 55 (8), 923-934 (2018).
  47. Lowe, A. J., et al. Dynamics of deposition and fossil preservation at the early Eocene Okanagan Highlands of British Columbia, Canada: insights from organic geochemistry. Palaios. 37 (5), 185-200 (2022).
  48. Gushulak, C. A. C., West, C. K., Greenwood, D. R. Paleoclimate and precipitation seasonality of the Early Eocene McAbee megaflora, Kamloops Group, British Columbia. Canadian J Earth Sci. 53 (6), 591-604 (2016).
  49. Poorter, H., Niinemets, Ü, Poorter, L., Wright, I. J., Villar, R. Causes and consequences of variation in leaf mass per area (LMA): a meta-analysis. New Phytol. 182 (3), 565-588 (2009).
  50. Enquist, B. J., et al. Scaling from traits to ecosystems. Adv Ecol Res. 52, 249-318 (2015).
  51. Burnham, R. J., Wing, S. L., Parker, G. G. The reflection of deciduous forest communities in leaf litter: implications for autochthonous litter assemblages from the fossil record. Paleobiology. 18 (1), 30-49 (1992).
  52. Burnham, R. J., Ellis, B., Johnson, K. R. Modern tropical forest taphonomy: does high biodiversity affect paleoclimatic interpretations. PALAIOS. 20 (5), 439-451 (2005).
  53. Jacobs, B. F. Estimation of low-latitude paleoclimates using fossil angiosperm leaves: examples from the Miocene Tugen Hills, Kenya. Paleobiology. 28 (3), 399-421 (2002).
  54. Sack, L., et al. Developmentally based scaling of leaf venation architecture explains global ecological patterns. Nat Comm. 3 (1), 837(2012).
  55. Merkhofer, L., et al. Resolving Australian analogs for an Eocene Patagonian paleorainforest using leaf size and floristics. Am J Botany. 102 (7), 1160-1173 (2015).
  56. Baumgartner, A., Donahoo, M., Chitwood, D. H., Peppe, D. J. The influences of environmental change and development on leaf shape in Vitis. Am J Botany. 107 (4), 676-688 (2020).

重印与许可

标签

ImageJ
利用数字叶形态学和单位叶面积质量通过化石叶片重建陆地古气候与古生态
联系我们
向图书馆推荐

隐私政策

使用条款

政策

版权所有 © 2026 MyJoVE Corporation。保留所有权利。