هل هو بخيل في الأيض أم مسرف؟ رؤى حول جلد الكسلان الأرضي وعلم وظائف الحرارة الذي كشفت عنه النمذجة الفيزيائية الحيوية وقياسات درجة الحرارة القديمة باستخدام نظائر الكتلة. Metabolic skinflint or spendthrift? Insights into ground sloth integument and thermophysiology revealed by biophysical modeling and clumped isotope paleothermometry

المجلة: Journal of Mammalian Evolution، المجلد: 32، العدد: 1
DOI: https://doi.org/10.1007/s10914-024-09743-2
PMID: https://pubmed.ncbi.nlm.nih.gov/39822851
تاريخ النشر: 2025-01-14

هل هو بخيل في الأيض أم مسرف؟ رؤى حول جلد الكسلان الأرضي وعلم وظائف الحرارة الذي كشفت عنه النمذجة الفيزيائية الحيوية وقياسات درجة الحرارة القديمة باستخدام نظائر الكتلة.

مايكل د. ديك • وارن ب. بورتر • بول د. ماثيوسون • ديفيد م. لوفليس • راندون ج. فلوريس أرادنا ك. تريباتي روبرت أ. إيجل دارين م. شوارترز ميخائيل ت. بوتشر

تم القبول: 17 ديسمبر 2024 / نُشر على الإنترنت: 14 يناير 2025
© المؤلفون 2024، نشر مصحح 2025

الملخص

تم التعرف على بقايا الميجاثير منذ القرن الثامن عشر وكانت من بين أول الفقاريات الضخمة التي تم دراستها. بينما تظهر عدة أمثلة من الأنسجة المحفوظة تغطية كثيفة من الفراء للكسالى الأرضية الأصغر التي تعيش في المناخات الباردة مثل الميلودون والنوتروثيريوبس، فإن القليل جداً معروف عن جلد الميجاثير. بافتراض أن الميجاثير كان لديه أيض مشابه لأثداء الثدييات المشيمية، تم الافتراض سابقاً أنه سيكون لديه القليل من الفراء أو لا يوجد على الإطلاق لأنه حقق أحجاماً جسمية عملاقة. هنا يتم اختبار “نموذج الجلد الخالي من الشعر” باستخدام التحليلات الجيوكيميائية لتقدير درجة حرارة الجسم لإنشاء نماذج جديدة لأيض الكسالى الأرضية، وتغطية الفراء، والمناخ القديم باستخدام برنامج Niche Mapper. تشير المحاكاة التي تفترض نشاطاً أيضياً مشابهاً لتلك الموجودة في الزنارثرا الحديثة إلى أن تغطية الفراء القليلة ستؤدي إلى إجهاد بارد عبر معظم النطاقات العرضية التي كانت تسكنها الكسالى الأرضية المنقرضة. على وجه التحديد، كان الإريموتيريوم يحتاج بشكل أساسي إلى فراء كثيف بعمق 10 مم مع تداعيات لتغيرات موسمية في عمق الفراء في أقصى خطوط العرض الشمالية وفراء قليل في المناطق الاستوائية؛ بينما كان الميجاثير يحتاج إلى فراء كثيف بعمق 30 مم على مدار السنة في نطاقه الحصري من المناخات الأكثر برودة وجفافاً؛ وكان الميلودون والنوتروثيريوبس يحتاجان إلى فراء كثيف. الفرو لتجنب الضغط الحراري، متطابقًا مع بقايا الجلد لكلا الجنسين، مما يشير أيضًا إلى استخدام التنظيم الحراري السلوكي. علاوة على ذلك، أظهرت بيانات قياس درجة الحرارة باستخدام النظائر المتجمعة من الأسنان المحفوظة لأربعة أجناس من الكسلان الأرضي درجات حرارة جسم معاد بناؤها أقل من تلك التي تم الإبلاغ عنها سابقًا للثدييات الأرضية الكبيرة. سمح هذا المزيج من الأيض المنخفض والفرو الكثيف لكسالى الأرض بالعيش في بيئات متنوعة.

الكلمات الرئيسية درجة حرارة الجسم ⋅ الفراء ⋅ الجيولوجيا الكيميائية ⋅ الأيض ⋅ المناخ القديم ⋅ علم الحفريات ⋅ الجلد ⋅ زينارثرا
△ مايكل ت. بوتشر
mtbutcher@ysu.edu
1 قسم العلوم الكيميائية والبيولوجية، جامعة يونغستاون ستيت، يونغستاون، أوهايو، الولايات المتحدة الأمريكية
2 قسم البيولوجيا التكاملي، جامعة ويسكونسن ماديسون، ماديسون، ويسكونسن، الولايات المتحدة الأمريكية
3 قسم علوم الأرض، جامعة ويسكونسن ماديسون، ماديسون، ويسكونسن، الولايات المتحدة الأمريكية
4 قسم علوم الأرض والفضاء، جامعة كاليفورنيا – لوس أنجلوس، لوس أنجلوس، كاليفورنيا، الولايات المتحدة الأمريكية
5 قسم العلوم الجوية والمحيطية، معهد البيئة والاستدامة، مركز القيادة المتنوعة في العلوم، جامعة كاليفورنيا – لوس أنجلوس، لوس أنجلوس، كاليفورنيا، الولايات المتحدة الأمريكية
قسم علوم الأرض، جامعة بويز ستيت، بويز، أيداهو، الولايات المتحدة الأمريكية

مقدمة

تمتلك الزواحف ذات الأطراف الغريبة تاريخًا غنيًا ومتعدد الجوانب محفوظًا في السجل الأحفوري يعود إلى 60 مليون سنة مضت (MYA) إلى عصر الإيوسين. تتكون الزواحف ذات الأطراف الغريبة من الفصائل Cingulata، التي تضم المدرعات الحية والبانباثير والغيلبتودون المنقرضة، وPilosa، التي تحتوي على آكلات النمل والكسالى. بينما تعتبر الأجناس Choloepus (الكسالى ذات الأصبعين) وBradypus (الكسالى ذات الثلاثة أصابع) هي السلالات الحية الوحيدة من الكسالى (الشجرية)، كانت هناك العديد من الأجناس المنقرضة من الكسالى (أي الشجرية، الأرضية، والمائية) التي احتلت في السابق موائل متنوعة خلال العصر الأوليغوسيني واستمرت حتى أواخر العصر الجليدي. ازدهرت الكسالى الأرضية في جميع أنحاء أمريكا الجنوبية وانتشرت شمالًا لتستعمر في النهاية أمريكا الشمالية خلال التبادل البيولوجي العظيم الأمريكي في أواخر العصر السينوزوي (GABI) (فيلدهامر وآخرون 2007). عدة كسالى أرضية منقرضة
أصبحت أنواع الكسلان مثل إيريموثيريوم وميغاثيريوم أكبر الكائنات ذات المفاصل الغريبة التي عاشت على الإطلاق، حيث وصلت إلى أحجام تعادل أو تتجاوز أحجام الفيلة الأفريقية الحديثة في السافانا.
حتى الآن، لم يتم العثور على جلد محفوظ يمكن أن يُنسب بشكل قاطع إلى أكبر الكسلان الأرضية إيريموثيريوم وميغاثيريوم. تم فقدان العينة الوحيدة من ميغاثيريوم التي تحتوي على عينة من الجلد قبل أن يتم جردها (بوجوس وسالاس 2004). ومع ذلك، تم إعادة بناء كلا النوعين تاريخياً على أنهما كانا يمتلكان فروًا كثيفًا وطويلًا يشبه ذلك الخاص بالكسلان الأرضية الأصغر (مثل ميلودون ونوثروثيريوبس) والكسلان الشجري الحالي. ومع ذلك، تم تحدي هذا المظهر بناءً على الديناميات الحرارية التي تشير إلى أن أكبر الكسلان الأرضية كانت تمتلك جلدًا خاليًا إلى حد كبير من الشعر (ماكناب 1985؛ فاريña 2002). استخدمت هذه التحليلات المعادلة التالية للتوصيل الحراري: ، حيث هو كثافة تدفق الحرارة في الجلد، هو الموصلية الحرارية للجلد، و s هو سمك الجلد، ودرجة حرارة الجسم الأساسية ودرجة حرارة البيئة يتم تمثيلها بـ و ، على التوالي. افترض فاريña (2002) أنه نظرًا لأحجامها الضخمة، فإن الكسلان الأرضي الكبير لن يكون محايدًا حراريًا إذا كان مغطى بفرو سميك. ومع ذلك، إذا كانت غلاف ميغاثيريوم تشبه بشكل وثيق غلاف أكبر الثدييات الأرضية الحالية (أي، الفيلة، فرس النهر، ووحيد القرن)، فسيكون في منطقة محايدة حراريًا حتى لو كانت درجات الحرارة المحيطة قريبة من في المتوسط. تُعرف منطقة الحرارة المحايدة بأنها نطاق من درجات الحرارة المحيطة حيث يمكن للكائنات الداخلية أن تبقى في توازن حراري دون إنتاج حرارة أيضية أو الحاجة إلى فقدان الحرارة بالتبخر. بينما قد يكون من المغري إجراء مقارنات مع حيوانات الميجافونا ذات الفراء الكثيف مثل الماموث الصوفي، يجب ملاحظة أن هذه المقارنات غير قابلة للتطبيق على الميجاثير. كان لدى الماموث الصوفي مجموعة من التكيفات للعيش في التندرا المتجمدة في سيبيريا مثل تقليل حجم الأذنين والذيل ووجود طبقة عازلة من الدهون (Boeskorov et al. 2016)، بالإضافة إلى الهيموغلوبين الذي منع دمائهم من التجمد (Campbell et al. 2010).
تم استحضار “نموذج الجلد الخالي من الشعر” مرة أخرى مؤخرًا من قبل ليندسي وآخرين (2020) فيما يتعلق بتجمع ضخم من بقايا إيريموثيريوم التي تم استردادها من تانكي لوما في الإكوادور. كان يُعتقد أن سبب وفاة هذه الحيوانات هو المرض (أو الجفاف) نتيجة تلوث فضلاتها للموحلة. افترض ليندسي وآخرون (2020) أن سلوك الاستحمام سيكون مفيدًا إذا كان لدى E. laurillardi جلد خالي من الشعر، حيث سيوفر الحماية من الشمس والحشرات الاستوائية، وبالتالي سيكون مشابهًا لاستراتيجية الحماية المستخدمة من قبل أكبر الثدييات الأرضية الحية.
على الرغم من الأدلة التي تدعم نموذج خالٍ من الشعر (أو شبه خالٍ من الشعر) في أكبر الكسلان الأرضية التي عاشت على الإطلاق
عاشت، كانت هناك افتراضات حاسمة مطبقة على المنهجيات السابقة التي قد لا تكون قابلة للتطبيق. أولاً، معادلة الموصلية الحرارية المستخدمة من قبل فاريña (2002) تأخذ في الاعتبار فقط درجات حرارة الجسم الأساسية وتأثيرات درجات الحرارة المحيطة وحدها على التوازن الحراري دون الأخذ في الاعتبار عوامل بيئية أخرى. بينما تعتبر هذه الحسابات مفيدة في تحديد تدفق التوصيل الحراري، كان ينبغي أن تؤخذ في الاعتبار عوامل بيئية أخرى قد تؤثر على التوازن الحراري مثل الرطوبة النسبية وسرعة الرياح، بالإضافة إلى التغيرات في درجة الحرارة المحيطة عبر التوزيع الجغرافي الواسع لأجناس الكسلان الأرضي. ثانياً، لم يكن النموذج السابق (فاريña 2002) دقيقًا في حساب معدل الأيض الأساسي للزواحف الحالية في حسابات تدفق الحرارة، أو جلد الكسلان الحي السميك الذي سيؤثر على فقدان الحرارة. بدلاً من ذلك، تم تكبير معدل الأيض الأساسي لإنسان عارٍ بشكل متساوي إلى حجم ميغاثيريوم مع افتراض درجة حرارة الجسم الأساسية وسمك الجلد مشابهين لتلك الخاصة بالبشر (فاريña 2002). ثالثاً، تم تحديد معدل الأيض الأساسي لميغاثيريوم الأمريكي عن طريق تقليل معدل الأيض المقدر إلى النصف، مما أدى إلى الاستنتاجات بأنه لا يزال يمكن أن يكون بلا شعر وأن درجات الحرارة المنخفضة تصل إلى ستكون ضمن منطقة الحرارة المحايدة لها. رابعًا، تقديرات منطقة الحرارة المحايدة المستمدة من النموذج السابق لفصائل الكسلان الأرضي الأصغر غير متسقة مع بيانات المناخ القديم المعروفة والجلد المحفوظ. درجات الحرارة لـ و تم حساب C، على التوالي، على أنها ضمن منطقة الحرارة المحايدة لميلودون داروينى بلا شعر تمامًا مقابل نظيره الذي يمتلك فروًا بطول 4 سم، ولكن مع كل منهما له كتلة مكافئة لـ . (فاريña 2002).
تشير الأدلة الأحفورية المعروفة عن الميلودون، مع ذلك، إلى أنه كان مغطى بشكل واسع بالفرو، الذي كان طوله يصل إلى 6.5-10 سم في بعض مناطق جسمه (كولينز 1933؛ سانتوس وآخرون 2024). تم العثور على نظائر مستقرة في عينات معاصرة من مينا نوتيوماتودون تعود إلى قبل سنوات أيضًا أنتجت تقديرات لدرجة الحرارة السنوية المتوسطة تتراوح بين و (جونزاليس-غواردا وآخرون 2017). علاوة على ذلك، فإن العلاقة بين كتلة الجسم ودرجة حرارة الجسم في الثدييات معقدة ويُزعم أنها تتوسط بشكل أساسي من خلال علم البيئة (كلارك وروثيري 2008). تشير البيانات الحديثة حول طاقة الثدييات الكبيرة الحالية إلى أن المنطقة الحرارية المحايدة لا تمثل سوى نطاقات درجات الحرارة البيئية المواتية حيث لا يحتاج الجسم إلى إنفاق الكثير من الطاقة بخلاف معدل الأيض الأساسي، وليس النطاق الكامل للوظائف الطبيعية المرتبطة بالنشاط اليومي مثل البحث عن الطعام. وبالتالي، تم اعتماد المناطق التوجيهية، ومناطق التحمل، ومناطق البقاء من قبل مؤلفين آخرين (ميتشل وآخرون 2018؛ بويلز وآخرون 2019؛ ليفيسك ومارشال 2021)، والتي تشمل درجات الحرارة التي يمكن أن يعمل فيها الكائن الحي ذو الدم الحار بشكل طبيعي، ودرجات الحرارة حيث
يحدث إجهاد حراري، ودرجات الحرارة التي يحدث فيها الموت بسبب الإجهاد الحراري، على التوالي.
مؤخراً، حاول فاريلا وآخرون (2024) تقدير معدلات الأيض لعدة أنواع من الزينارثرا المنقرضة من خلال قياس قطر الثقوب المغذية في عظمة الفخذ. أدت التقديرات الخاصة بتدفق الدم المستمدة من معادلة طورها سيمور وآخرون (2019) إلى معدلات أيض في الزينارثرا الأحفورية قابلة للمقارنة مع الثدييات المشيمية غير الزينارثرا من نفس الحجم. بينما تدعم هذه النتائج التكهنات بأن الكسلان الأرضي كان يمتلك أيض مشيمي نموذجي بسبب زيادة كتلة العضلات الهيكلية مقارنة بالكسالى الشجرية الحية (فاريña 2002)، وواجه ظروفاً بيئية أقل من الحيادية الحرارية المقترحة (تيخادا وآخرون 2021)، تشير بيانات أخرى إلى أيض أكثر شبهاً بالزينارثرا. على سبيل المثال، أظهرت تقديرات درجة حرارة الجسم في نوتروثيريوبس (استناداً إلى الأحماض الأمينية المحفوظة في الكولاجين العظمي) أن لديها نطاقاً من قيم درجة حرارة الجسم مشابهة للزينارثرا الحية (هو 1967). تشير تقديرات مماثلة لمعدلات الأيض والموصلية الحرارية لكل من ميلودون ونوتروثيريوبس إلى أن كلاً منهما كان لديه معدل أيض أساسي منخفض وعزل نسبي ضعيف مقارنةً بالثدييات الأخرى من نفس الحجم (ماكناب 1985). من المحتمل أن يكون الكسلان الأرضي قد امتلك معدلات أيض أقل من الثدييات الأخرى من حجمه، على الرغم من أن هذا الاستنتاج المؤقت لا يعني أنهم كانوا بطيئي الحركة مثل الكسلان الشجري الحي. باستخدام طول الخطوة من مسار أحفوري في بيهوان كوه، الأرجنتين، المنسوب إلى ميغاثيريوم، توقع بلانكو وتزيرفونوغورا (2003) أن متوسط سرعة مشيه كان ، وهو أسرع بعشر مرات من المشي المعلق الحالي لبراديبوس (Gorvet et al. 2020; McKamy et al. 2023). سرعة قصوى تصل إلى ما يقرب من ، وهو لا يزال أبطأ بكثير من السرعة القصوى لـ تم قياسها للفيلة (هوتشينسون وآخرون 2006)، وقد تم تقديرها أيضًا لميجاثيروم (بيليه وآخرون 2013). على الرغم من أن الكسلان الأرضي والفيلة والكسالى الشجرية لديهم أنماط مختلفة من الحركة، فإن السرعات المبلغ عنها تؤكد المبدأ العام الذي ينص على أن الحيوانات ذات معدلات الأيض المنخفضة تميل إلى التحرك بشكل أبطأ من تلك ذات معدلات الأيض الأعلى.
هدف هذا البحث هو اختبار نموذج الجلد الخالي من الشعر لأكبر حيوانات الكسلان الأرضية، إيريموثيريوم وميغاثيريوم، ووضع قيود جديدة على تنظيم درجة حرارتها. لاختبار النموذج بدقة من خلال وسائل كمية، تم إدخال العديد من العوامل المناخية والفسيولوجية، مثل درجة حرارة الجسم الأساسية، ومعدل الأيض، وتقديرات المناخ القديم. في الوقت نفسه، سيوفر دمج هذه المتغيرات مزيدًا من الفهم لإعادة بناء جلد حيوانات الكسلان الأرضية الكبيرة حتى يتم اكتشاف عينات من الجلد و/أو الفراء. بالإضافة إلى القيود الفسيولوجية والبيئية على الجلد، سيساهم ذلك في فهم أفضل لتطور درجة حرارة الجسم الأساسية و/أو معدل الأيض عبر مختلف الأنواع المنقرضة والحالية.
فروع النشوء والتطور للـ Folivora هي نتيجة متوقعة لهذه الدراسة.

المواد والأساليب

عينات الحفريات وتحضير العينات

عينات الأسنان من إريموتيريوم ميغاثيريوم بروماجاثيريوم ) و نوتروثيريوبس ( تم استخدام ( ) في سلسلة من التحليلات الجيوكيميائية (الجدول 1). تم استخدام عينات من العظام القشرية من نفس العينة أو من عينة أخرى في نفس المنطقة العامة لتحديد درجة التحول الجيولوجي من خلال تحليل امتصاص العناصر الأرضية النادرة (REE) في كل من عظام الأسنان والعينات القشرية. كانت العينات المستخرجة من مجموعات متحف التاريخ الطبيعي في شيكاغو، إلينوي (FMNH) وجامعة فلوريدا في غينزفيل، فلوريدا (UF). لم تكن هناك حاجة للحصول على تصاريح للدراسة الحالية، التي امتثلت لجميع اللوائح القانونية ذات الصلة. بينما لم يتم الحصول على عينات من Mylodon لهذه الدراسة، تم استنتاج تقديرات درجة حرارة الجسم لهذا الجنس من القيم المستخلصة من Promegatherium نظرًا لأن كلا الجنسين كانا متشابهين في حجم الجسم. هذه المقارنة المستندة إلى حجم الجسم بدلاً من النسب التطورية مبررة نظرًا لأن Bradypus و Choloepus متشابهان في الحجم والقدرات الأيضية على الرغم من كونهما بعيدين عن بعضهما البعض (McNab 1985؛ Cliffe et al. 2023). استراتيجيات أخذ العينات أخذت في الاعتبار توزيع جغرافي جماعي واسع لهذه الأجناس من أقصى الجنوب في نطاقها في الأرجنتين إلى أقصى الشمال في نطاقها في فلوريدا.
نظرًا لأن العاج الخارجي لأسنان الزينارثران الأحفورية عرضة لامتصاص كميات كبيرة من عناصر الأرض النادرة (MacFadden et al. 2010)، تم أخذ عينة من العاج الداخلي بشكل تفضيلي، حيث أظهرت دراسة حديثة (Larmon et al. 2019) باستخدام الكاثودولومينسنس أن هذه المادة السنية هي الأكثر مقاومة للتغيرات الجيولوجية وتظهر أقل كمية من التمعدن مقارنةً بطبقات العاج الأخرى. تقريبًا تم جمع مسحوق العاج الداخلي (وعظم القشرة) باستخدام مثقاب يدوي كهربائي مزود برأس مثقاب ماسي بقطر 0.7 مم. تم أخذ عينات مثقوبة من العاج الداخلي من عدة مواقع عشوائية على طول السن بدءًا من القمة إلى الجذر (الشكل 1أ). لإزالة أي ملوثات عضوية، تم غسل العينات بـ من لمدة 3 ساعات في أنابيب سعة 15 مل. ظلت الأنابيب مفتوحة حتى يمكن أن يتم إطلاقه أثناء التنقية. ثم تم ترسيب العينات في جهاز الطرد المركزي بسرعة 10 آلاف دورة في الدقيقة (نصف قطر الطرد المركزي القوة الطرد المركزية ) لمدة 10 دقائق. على وجه التحديد، تم طرد جميع العينات في جهاز الطرد المركزي عند درجة حرارة ثابتة من ، تليها 3 شطفات في درجة حرارة الغرفة الماء. تم تصريف أي سائل متبقي قبل
الجدول 1 تقديرات درجة حرارة الجسم من تحليلات النظائر المتجمعة ونسب التحول من تحليلات العناصر الأرضية النادرة (REE). الاختصارات: SD، الانحراف المعياري؛ SE، الخطأ المعياري
النوع/ معرف العينة محلية بيئة الترسيب REE (%o، (%o, V-SMOW) ± الانحراف المعياري
إيريموثيريوم UF 95869 إنغليس 1 أ، مقاطعة سيتروس، فلوريدا، الولايات المتحدة الأمريكية حفرة انهيارية ٥ 0.42
نوتروثيريوبس UF 87131 ليسي 1 أ، مقاطعة هيلزبورو، فلوريدا، الولايات المتحدة الأمريكية مصب ٤ 0.86
ميغاثيريوم FMNH P13725 وادي تارِيخا، الأرجنتين نهري ٤ 0.12
إيريموثيريوم UF 312730 هالي 7G، مقاطعة ألاچوا، فلوريدا، الولايات المتحدة الأمريكية بحيري/ كارس ٤ 0.02
بروماجاثيريوم FMNH 14404 بوابة كورال كيمادو، الأرجنتين أرضي غير محدد ٣ ٣٦.٦
عدد قياسات النظائر المستقرة المكررة من نفس العينة
نسبة عناصر الأرض النادرة بين العاج والعظم القشري (انظر الطرق)
درجة الحرارة (أندرسون وآخرون 2021)
(حيث “w” تشير إلى مياه تكوين المعادن؛ ليكوييه وآخرون 2010)
تجفيفه طوال الليل في فرن عند كان من الضروري زيادة وقت التجفيف في الفرن لبعض العينات لضمان التجفيف الكامل للملوثات بحيث تم تجفيف العينات إلى كتلة مركزة.

تحليل التغيرات الجيوكيميائية

لتقييم كمية التحول الجيولوجي الموجودة في العينات، تم إجراء تحليلات العناصر الأرضية النادرة على العينات الفرعية. ) من العاج الداخلي المحضر والعظم القشري (MacFadden et al. 2010). للتحليل، تم وزن العينات بواسطة ميزان دقيق وتم هضمها في خليط من 0.1 مل من 29 م HF و 0.3 على صفيحة ساخنة في تم تجفيف العينات وإعادة حلها في 0.5 مل من ثلاث مرات قبل التخفيف بالوزن إلى لأغراض التحليل. تم إعداد العينات جنبًا إلى جنب مع مجموعة من المساحيق القياسية الجيولوجية والعظمية وعمليات الفراغ الكلي وفقًا لنفس إجراء الذوبان، لاستخدامها كمواد معايرة ومواد مرجعية خارجية. تم إجراء التحليلات على جهاز مطياف الكتلة أحادي الرباعي Thermo Scientific iCAP-RQ باستخدام نازع رذاذ ميكرومست PFA؛ غرفة رذاذ كوانتية مبردة بتقنية بيلتييه، وحقن كوانتي 2.5 مم وشعلة كوانتية؛ مخاريط نيكل، إدخال مصفاة عالية المصفوفة (3.5 مم). تم سحب العينات باستخدام جهاز أخذ العينات التلقائي ESI SC 4DX FAST باستخدام طريقة FAST 2.0 مل، مع وقت استقرار قدره 30 ثانية وشطف مزدوج (4 ثوانٍ لكل منهما) في تم تتبع HF بعد كل تحليل. تم خلط العينات في الوقت الحقيقي مع محلول 2 ppb In باستخدام جهاز أخذ العينات التلقائي FAST، الذي تم استخدامه كمعيار داخلي للتجربة. التحليلات المبلغ عنها هي متوسط 5 قياسات، والتي
تتكون من 25 تمريرة على ، ، ، لمدة 10 مللي ثانية لكل منها. تم تطبيع العدادات الخام إلى كثافة موحدة على مدار التشغيل من خلال التطبيع مقابل In. تم تحويل التركيزات من العدادات المعيارية باستخدام برنامج Thermo Scientific Qtegra (eQuant) من خلال الانحدار الخطي لمعايير المعايرة (BHVO-2، BCR-2، BIR-1، AGV-2؛ Jochum et al. 2016)، مع ربطها بالفراغ الإجرائي الكلي. تم مراقبة الدقة والموثوقية باستخدام مسحوق الصخور الكاملة W2 (Jochum et al. 2016)، رماد العظام NIST1400، وجبة العظام NIST1486 ومواد العظام الأحفورية NOCS (Chavagnac et al. 2007) وعمليات تخفيف متعددة لمحلول معيار مختلط بتركيز 100 جزء في المليون من العناصر الأرضية النادرة (BDH82025-926) مخفف إلى تركيزات مشابهة لتلك الخاصة بالعينات غير المعروفة.
تم تطبيع جميع القياسات مقابل الصخر الزيتي الأسترالي بعد العصر الأركي (PAAS؛ نانس وتايلور 1976) لتعويض تأثيرات تأثير وفرة الأعداد الزوجية والفردية لأودو-هاركنز. بينما أبلغ ماكفادن وآخرون (2010) عن زيادة كبيرة في امتصاص العناصر الأرضية النادرة في الأحافير مقارنة بالعظام الحديثة وأنسجة الأسنان، تم إنشاء مؤشر للعناصر الأرضية النادرة لتقييم درجة التحول الجيولوجي في أسنان الأحافير من نوع زينارثران مقارنة بعينات المينا من مواقع مشابهة. مجموع جميع العناصر الأرضية النادرة المطبوعة ( تم تقسيم العناصر الأرضية النادرة الموجودة في العاج العادي على مجموع جميع العناصر الأرضية النادرة المعيارية الموجودة في العظم القشري والتعبير عنها كنسبة. عينة/ العظم القشري). النسب تشير التغيرات الجيوكيميائية الكبيرة إلى أنه يجب التعامل مع بيانات الجيوكيمياء النظيرية بحذر، في حين أن نسب العينات بين 0.5 و 1.0 تشير إلى تغيرات خفيفة إلى متوسطة.
الشكل 1 أ. سن Megatherium FMNH P13725 مع المناطق المأخوذة عينة منها موضحة بالأسهم الصفراء؛ مقياس الرسم يساوي 30 مم. ب. نسبة الفوسفات في أسنان الكسلان الأرضي المأخوذة من العينة مقارنة بالنطاقات المبلغ عنها لبيوآباتيت الثدييات المشيمية الحالية؛ تشير المناطق المظللة باللون الأزرق والرمادي إلى النسبة المئوية 95 و 99 على التوالي؛
بيانات عظام الثدييات والمينا من Eagle وآخرون (2011). ج. مخطط تشتت للعلاقات بين ، ودرجات الحرارة المقدرة للأسنان المأخوذة؛ الرموز ذات الحدود الخفيفة تمثل عينات فردية، والرموز ذات الحدود السميكة وأشرطة الخطأ تمثل المتوسط والانحرافات المعيارية للعينات
عينات ونسب (مماثل لمعظم عينات المينا الأحفوري من غير الزينارثران) يشير إلى عينة بها القليل من التغير الجيوكيميائي أو لا توجد تغيرات على الإطلاق (ماكفادن وآخرون 2010).
تم إعداد عينات للفوسفات تحليل يتبع بروتوكول ترسيب الفوسفات الفضي السريع الموضح في Mine et al. (2017) والمفصل في Conte et al. (2024). باختصار، تم وزن حوالي 1.0 ملغ من مادة العينة في أنابيب ميكروسنترفيج سعة 2 مل وذوبانها طوال الليل في حمض النيتريك 2.0 م. فلوريد الكالسيوم ثم تم ترسيبه بإضافة 2.9 م من حمض الهيدروفلوريك (HF) و 2.0 م من هيدروكسيد الصوديوم تمت طرد العينات لتكوين راسب وتم نقل الراسب المحتوي على أيون الفوسفات إلى أنبوب ميكروسنترفيوج جديد ونظيف. تم شطف الحبيبة بمحلول فلوريد الصوديوم بتركيز 0.1 م لاستعادة الفوسفات المتبقي، وتم دمج العينة الناتجة من السائل العلوي
مع العينة الأولى. تم ترسيب الفوسفات العيني كفوسفات فضي مع محلول أمين فضي يحتوي على 1.09 م من هيدروكسيد الأمونيوم و 0.37 م نترات الفضة في نطاق pH من 5.5 إلى 7.5. تم إجراء التعديلات على pH من خلال إضافة كميات صغيرة من أو 2 م NaOH وتم السماح للعينات بالتفاعل لمدة 10 دقائق. تم طرد العينات لتكوين راسب. البلورات وبعد إزالة السائل العلوي تم شطفها خمس مرات بالماء المقطر وتجفيفها طوال الليل في الفرن. تم وزن العينات المجففة إلى في كبسولات فضية محكمة التحضير وتم تحليلها ثلاث مرات باستخدام محلل العناصر الحرارية من ThermoFisher (TC/EA) مقترن بمقياس الكتلة لنسبة النظائر بتدفق مستمر Delta V plus في مختبر النظام البيئي للنظائر المستقرة بجامعة UC Merced (SIELO). تم ترسيب وتحليل المادة المرجعية NIST SRM 120 C.
جنبًا إلى جنب مع عينات غير معروفة كمعيار تحضير وأسفر عن قيمة لـ مواد مرجعية USGS و USGS 81 ( ) تم استخدامها لتصحيح الخطية والانحراف والتوحيد.

كيمياء النظائر المستقرة للبيوآباتيت وتقديرات درجة حرارة الجسم

تعتبر قياسات نظائر الكربونات المتجمعة طريقة لتحديد درجات حرارة تكوين المعادن استنادًا إلى القياسات (التي يتم الإبلاغ عنها على أنها قيم) وفرة الروابط في وحدات الكربونات (غوش وآخرون 2006). القيم الموجودة في البيوآباتيت السني تعكس درجة حرارة جسم الكائن الحي الذي شكلها (Eagle et al. 2010). علاوة على ذلك، تستطيع قياسات المعادن الحيوية الفقارية التمييز بين الفسيولوجيا الداخلية والخارجية (Eagle et al. 2011, 2015; Griffiths et al. 2023)، وكذلك بين ذوات الدم الحار ذات معدلات الأيض المختلفة. على سبيل المثال، لعبت هذه التقنيات دورًا حاسمًا في حل الفروق في درجة حرارة الجسم بين أسماك القرش ذات الدم الحار إقليميًا، والثدييات الكبيرة البرية والبحرية، وبين الثدييات والطيور. كما أن القياسات المستمدة من الأحافير أسفرت أيضًا عن درجات حرارة جسم معقولة للفقرات المنقرضة عبر مجموعة من العصور (Eagle et al. 2010, 2011, 2015; Griffiths et al. 2023).
تم تحليل جميع العينات باستخدام مطيافية الكتلة لنسبة النظائر المستقرة. على وجه التحديد، تم هضم عينات من العاج الداخلي لكسلان الأرض في حمام حمضي مشترك يحتوي على تفاعل حمض الفوسفوريك عند تليها خطوات متعددة من التنقية بالتبريد العميق والمرور عبر إما عمود كروماتوغرافيا الغاز (مطياف الكتلة MAT 253، ثيرمو فيشر ساينتيفيك، الولايات المتحدة الأمريكية) أو فخ الامتصاص (مطياف الكتلة Nu Perspective، Nu Instruments-AMETEK، المملكة المتحدة) المحشو بـ Porapak Type-QTM 50/80. بعد تنقية العينة، ثم تم إدخال الغاز إلى أنظمة مطياف الكتلة لنسبة النظائر (IRMS) عبر كتلة تغيير تلقائية، مما سمح بإجراء قياسات متناوبة مستمرة للعينة والمرجع. غاز على الكواشف المهيأة لقياس النظائر الجزيئية ذات الكتل من 44 إلى 49. تم إجراء القياسات على مطياف الكتلة Thermo MAT 253. من مادة العينة) تتكون من تسعة كتل من 10 دورات بإجمالي 720 ثانية من وقت التكامل عند إشارة على مدار كتلة مع تعديلات الضغط بين كل كتلة للعودة إلى شعاع بقوة 16 فولت في بداية كل كتلة جديدة. تم أخذ القياسات باستخدام مطياف الكتلة Nu Perspective ( من مادة العينة) تتكون من ثلاث كتل من 20 دورة بإجمالي لدمج إشارة تراوحت من 80 إلى 30 نانو أمبير على الكتلة 44 على مدار كل عملية اكتساب. تم استخدام معايير الكربونات ETH-1 و ETH-2 لإجراء تصحيح غير خطي على العينة.
غير معروفة. تم بعد ذلك إسقاط هذه القيم المصححة غير الخطية في إطار مرجعي I-CDES (Bernasconi et al. 2021) باستخدام ETH-1-2 و ETH-3 جنبًا إلى جنب مع ثلاثة معايير داخلية. لم يتم تضمين ETH-4 في التصحيحات، بل تم استخدامه كتحقق من جودة البيانات. تم تطبيق التصحيحات على عينات الأسنان وبيانات المعايير على متوسط متحرك من 10 معايير على كل جانب من العينة المعطاة وتم حسابها باستخدام برنامج Easotope (جون وبوين 2016). تم اشتقاق درجات الحرارة المعاد بناؤها من -طريقة معايرة درجة الحرارة (أندرسون وآخرون 2021) التي أظهرت مؤخرًا أن استخدام إطار المرجع I-CDES الأحدث كانت معايرة البيوأباتيت غير قابلة للتمييز إحصائيًا عن المعايرة المستخدمة في هذه التحليلات (أندرسون وآخرون 2021؛ غريفيثس وآخرون 2023). المعادلات المستخدمة لتقدير تكوين المياه أثناء التمعدن منشورة في مكان آخر (ليكوير وآخرون 2010).
تعتبر قيود مقبولة في تطبيق قياس درجة الحرارة باستخدام نظائر الكتل لتقييم فسيولوجيا الفقاريات هي الحفاظ على توقيعات النظائر الأصلية وإمكانية حدوث تغييرات بعد الإيداع، وكيف تؤثر هذه العوامل على دقة التقديرات. بينما يمكن أن يكون تحديد التغيرات بعد الإيداع معقدًا، تم استخدام النهجين التاليين: (1) تحديد تركيز العناصر الأرضية النادرة في أسنان الكسلان الأحفوري و(2) مقارنة القياسات لـ من كل من مجموعة الكربونات ومجموعة الفوسفات في البيوآباتيت الأحفوري. نظرًا لأن مجموعة الكربونات قد تكون أكثر عرضة للتغيير الكيميائي من مجموعة الفوسفات، فقد تم استخدام الفارق بين نطاق الكربونات والفوسفات الذي يُرى في الكائنات الحية الحالية للإشارة إلى تفضيل الدياجينيسيس لمجموعة الكربونات (Eagle et al. 2011). يمكن أن يتم امتصاص العناصر الأرضية النادرة بواسطة البيوآباتيت خلال عملية التحجر، وقد يشير محتوى العناصر الأرضية النادرة المرتفع جدًا إلى إعادة بلورة الفوسفات الأحفوري (McFadden et al. 2010).

نمذجة الحيوانات والمناخ القديم

تم تقييم تقديرات درجة حرارة الجسم المستمدة من الأدبيات للزواحف الحية (McNab 1985؛ Cliffe et al. 2018؛ Seitz و Puig 2018) من خلال الانحدار الخطي لتحديد العلاقة بين درجة حرارة الجسم الأساسية وحجم الجسم لجميع المحاكاة الأولية. درجة الحرارة (بالـ تم تحويل البيانات إلى لوغاريتم ورسمها مقابل لوغاريتم الكتلة (بالكيلوغرام) لتحديد العلاقة المتوقعة. لم تُظهر هذه البيانات أي علاقة بين الأنواع الحالية من الزينارثران (الشكل 2)؛ ومع ذلك، تتراوح كتل أجسام الزينارثران الحالية بين 0.2 و 45.2 كجم، وبالتالي لا تقارب أي من الأنواع المستخدمة في هذه العلاقة الحجم الجسدي للكسالى الأرضية المنقرضة. ومع ذلك، كانت محاكاة النمذجة تهدف إلى تمثيل حيوان بالغ متوسط الحجم في أي نقطة زمنية جيولوجية معينة، لكنها ليست مبنية على عينة محددة مثل تلك المستخدمة في تحليلات النظائر المتجمعة. الجسم
الشكل 2: انحدار كتلة الجسم الحالية للزواحف الغريبة مقابل درجة حرارة الجسم ; ). تشمل البيانات ثلاثة أنواع من الكسلان ( ثمانية أنواع من المدرع ) وثلاثة أنواع من آكلات النمل ( ) وتم الحصول عليها من الأدبيات (ماكناب 1985؛ كليف وآخرون 2018؛ سايتز وبويغ 2018)
تم إدخال تقديرات درجة الحرارة، المستمدة من تحليلات مقياس الحرارة القديمة باستخدام نظائر متجمعة بين مجموعة متنوعة من المتغيرات البيئية والفسيولوجية، في برنامج نمذجة البيولوجيا الفيزيائية الذي تم اختباره على نطاق واسع، Niche Mapper. يقوم النموذج الفرعي للمناخ الدقيق في Niche Mapper بتحويل مدخلات قيم المناخ الكلي إلى أكثر المناخات الدقيقة حرارة (أقل ظل) وأبردها (أقصى ظل) في نفس الوقت من 2 متر فوق سطح الأرض إلى 2 متر تحت سطح الأرض. ثم يستخدم النموذج الفرعي للحيوانات النتائج الساعية للعثور على بيئة دقيقة مناسبة للراحة الحرارية المثلى، مع السماح للمستخدم بتعريف مجموعة متنوعة من المدخلات الشكلية والسلوكية والفسيولوجية، بما في ذلك هندسة جسم الحيوان، ومعدل الأيض المحدد من قبل المستخدم، وسلوكيات التعشيش والحفر، وما إلى ذلك. يحل كلا النموذجين الفرعيين توازنات الحرارة والكتلة للركائز المحلية والحيوان، على التوالي، بالإضافة إلى العثور على أفضل مناخ دقيق متاح لكل ساعة من يوم محاكى للحفاظ على الحيوان المحاكى ضمن أقرب نطاق ممكن من منطقة التوازن الحراري.
تم استخدام العمل الأخير مع Niche Mapper للتنبؤ بالقيود البيئية والفسيولوجية على أنواع منقرضة مختلفة، بما في ذلك الماموث (Wang et al. 2018, 2020) والديناصورات الترياسية (Lovelace et al. 2020). تم تقدير معدل الأيض الأساسي باستخدام علاقة قياس الأيض المنشورة المستمدة من 15 نوعًا حاليًا من الزينارثرا. ؛ وايت وسيمور 2004)، حيث M هو كتلة الجسم بالجرامات، ومعدل الأيض الأساسي الناتج في تم تحويله إلى وحدة الواط. تم استخدام مضاعف النشاط 2 لتقدير قيم معدل الأيض الميداني. العلاقة التناسبية ( ؛ فاريña 2002)، حيث M هو كتلة الجسم بالكيلوغرام، تم اشتقاقه في الأصل من كليبر (1932) واستخدم لمقارنة معدلات الأيض المستمدة من ثدييات ذات حجم جسم مشابه مع الأيض المشيمي النموذجي.
تم تجهيز نماذج الكمبيوتر ثلاثية الأبعاد لكل من Eremotherium وMylodon باستخدام منحنيات القاعدة غير المنتظمة المعقولة (NURBS) واستخدمت لتقدير أحجام الجسم، ومساحات السطح، ومعدل فقدان الحرارة، وكانت هذه البيانات بمثابة مؤشرات للمعلمات المحسوبة نفسها في النماذج الأيضية لـ Megatherium وNothrotheriops، على التوالي. تم بناء النماذج في ZBrush بواسطة فناني الباليون، وتم استخدامها في هذه الدراسة بإذن. على الرغم من العلاقة البعيدة من الناحية التطورية، تم استخدام النموذج ثلاثي الأبعاد لـ Mylodon بشكل خاص كمؤشر لـ Nothrotheriops نظرًا لأن كلا النوعين كانا أكثر تشابهًا في حجم الجسم مقارنةً بأي منهما مع الميجاثير. تقديرات الكتلة لـ Eremotherium ( ; باربوسا وآخرون 2023)، ميغاثيريم ( ; براسي وغاردينر 2015)، ميلودون ( ؛ تم إدخال تقديرات الكتلة لـ Nothrotheriops (463 كجم؛ مكدوغال 2005) في هذه المحاكاة. بينما تم الإبلاغ عن تقديرات كتلة أعلى في الأدبيات لكل من Eremotherium و Megatherium ( ; دانتاس وآخرون ؛ فاريña وآخرون 1998)، تم اعتبار هذه تقديرات عالية هنا، وتم استخدام كتل جسمية أكثر تحفظًا بدلاً من ذلك. تم اختبار نماذج تغطية الجسم بالكامل بسمك فراء يبلغ 10 و30 و50 مم لجميع الأنواع الأربعة. كانت كثافات الفراء التي تم اختبارها تتكون بشكل خاص من فراء “كثيف” يحتوي على 2000 شعرة. (Lovelace وآخرون 2020) و 1,300 شعرة (على سبيل المثال، كثافة شعر جانب الحصان؛ تريجير 1965)، بالإضافة إلى “شعر نادر” يتكون من 8.5 شعرة (على سبيل المثال، كثافة شعر جانب الخنزير؛ تريجير 1965) و 0.07 شعرة/ (أقصى كثافة لفراء الفيل؛ مايرفولد وآخرون 2012).
تم استخدام محاكاة غرف التمثيل الغذائي لكل نوع لتقدير نطاق درجة الحرارة المحايدة تقريبًا، بالإضافة إلى درجة الحرارة الحرجة الدنيا ودرجات الحرارة الحرجة العليا. تم اشتقاق نطاقات درجة حرارة الجسم الأساسية النهائية المستخدمة من تلك التي تم تحديدها من خلال تحليل الباليوثيرمومترية لل isotopes المجمعة، بالإضافة إلى درجات حرارة الجسم الأساسية المنشورة من الزينارثران الحالية (مك ناب 1985؛ كليف).
تم تحديد مناطق الحرارة المحايدة الناتجة ضمن نطاق مستهدف من معدل الأيض الأساسي مع الحيوان الذي تم نمذجته وهو واقف بلا تعرض لأشعة الشمس. سرعات الرياح الثابتة تبلغ 0.1 تم افتراضها أيضًا لمحاكاة غرفة الأيض. تم تحديد درجة الحرارة الحرجة العليا عندما انخفض معدل الأيض إلى أقل من 5% من نطاق معدل الأيض الأساسي المستهدف، بينما تم تحديد درجة الحرارة الحرجة السفلى عندما كان معدل الأيض أعلى من نطاق معدل الأيض الأساسي المستهدف وفقًا للتحليلات السابقة (Lovelace et al. 2020). تم تكرار هذه العملية لكل نموذج جلد مقترح.
لأخذ الظروف البيئية لمواقع الحفريات (بخلاف درجات الحرارة المحيطة) خلال العصر الجليدي في الاعتبار، فإن قاعدة بيانات علم الحفريات (https://paleobiodb.org) تم استخدامه لتحديد النطاقات الشمالية والجنوبية الأكثر بعدًا، بالإضافة إلى نطاق الأنواع في خطوط العرض المتوسطة لاختبار مجموعة متنوعة من الظروف المناخية لكل نوع. تم الحصول على متوسط درجات الحرارة السنوية المقابلة للأعمار الطبقية التقريبية للمواقع المحاكاة من محاكاة PALEO-PGEM v1.0 المنشورة سابقًا لدرجات الحرارة العالمية المتوسطة السنوية المختلفة على مدى 5 ملايين سنة مع إجمالي 5001 قياس في فترات زمنية مدتها 1000 عام (Holden et al. 2019). ثم تم استخدام مستخرج المناخ العالمي في Niche Mapper لتحديد الفرق المعامل من قيم المناخ الحديثة لدرجات الحرارة القصوى والدنيا الشهرية لعمر الموقع الطبقي المقدر. تم تقريب الرطوبة النسبية، وتغطية السحب، وسرعة الرياح من القيم الحديثة المقاسة بواسطة قاعدة بيانات مستخرج المناخ العالمي في Niche Mapper لمناطق مماثلة. يقوم نموذج المناخ الجزئي في Niche Mapper بتحويل بيانات المناخ الشهرية القصوى/الدنيا المحددة من قبل المستخدم إلى قيم ساعة باستخدام دوال جيبية تفترض نفس ظروف الطقس لكل يوم من أيام شهر معين. ثم يقوم النموذج بحل توازن الحرارة وكتلة الماء لسطح الأرض ودرجات الحرارة وسرعات الرياح والرطوبة (2 م فوق وتحت مستوى سطح الأرض). يتم حساب الإشعاع الشمسي من خوارزمية إشعاع شمسي مضمنة (McCullough and Porter 1971). تم حساب الإشعاع الحراري للسماء من تحليل انحدار تجريبي بناءً على المدخلات لدرجة حرارة الهواء على ارتفاع 2 م فوق سطح الأرض (Idso and Jackson 1969).
أخيرًا، بالنسبة لمحاكاة Niche Mapper لكسالى الأرض في بيئاتها المعنية، تم تحديد منطقة الراحة الحرارية لمتطلبات معدل الأيض اليومي المقدر لكل نوع بناءً على الحدود العليا والسفلى المحتملة للنشاط الأيضي اليومي. تم حساب إجمالي الإنفاق الأيضي اليومي باستخدام تكامل شبه منحرف لتوقعات معدل الأيض الساعي. تم تحديد الحد الأدنى لمنطقة الهدف من خلال ضرب معدل الأيض الأساسي المقدر في بواسطة 24، بينما تم تحديد الحد الأقصى من خلال ضرب معدل الأيض الميداني المقدر
(أي، معدل الأيض الأساسي مضروبًا في عامل 2 محول إلى ) بنفس العامل. ضمن هذه المنطقة، يمكن للحيوانات المحاكية أن تعمل بمستويات أيض طبيعية طالما كانت الموارد متاحة بسهولة. من أجل البساطة، تُشار إلى مناطق التحمل والبقاء فوق الحدود العليا للحدود العليا لمنطقة الراحة الحرارية المحاكية هنا باسم “ضغط البرد” حيث إن الحيوان المحاكي سيكون مضطرًا لإنفاق طاقة أيض إضافية لمنع انخفاض حرارة الجسم. وبالمثل، تُشار إلى مناطق التحمل والبقاء تحت الحدود السفلى لمنطقة الراحة الحرارية باسم “ضغط الحرارة” حيث سيكون على الحيوان المحاكي خفض معدل أيضه تحت معدل الأيض الأساسي المقدر لمنع ارتفاع حرارة الجسم. قدمت نتائج Fariña (2002) تحليل حساسية لمقارنة تحملات الحرارة لكسالى الأرض بمعدلات أيض أكثر تشابهًا مع الثدييات المشيمية الكبيرة، بدلاً من أيض يشبه الزين-أرترا.

تحليل إحصائي

تنبؤات معدل الأيض اليومي (بالـ ) التي تم إنشاؤها بواسطة محاكاة المناخ لكل نموذج من نماذج الغلاف (تغطية الفراء) تم رسمها كسلسلة من المنحنيات لتحديد أي نموذج يناسب أفضل تقديرات الحد الأدنى والحد الأقصى لإنفاق معدل الأيض اليومي، أو الساعات في اليوم حيث يمكن الحفاظ على معدل الأيض الأساسي دون الحاجة للبحث عن الظل، لكل نوع. بمجرد تحديد نماذج الغلاف التي تحتوي على أقل قدر من الضغط الحراري، تم مقارنة تقديرات الحد الأقصى لساعات النشاط الأيضي فوق مستويات معدل الأيض الأساسي عبر التوزيع الجغرافي لكل نوع جنبًا إلى جنب مع تلك الخاصة بالثدييات المشيمية الحالية ذات أحجام مشابهة. تم تعديل متوسطات استهلاك الطعام اليومي في المستخرجة من نماذج الغلاف التي تحتوي على أقل قدر من الضغط الحراري لكل بيئة محاكية باستخدام منحنى انحدار متعدد الحدود من الدرجة الثانية للإشارة إلى الاختلافات المحتملة في تناول الطعام بين الأنواع. تم اشتقاق القيمة الغذائية للأوراق، والأعشاب، والنباتات الخشبية من مكتبة تكوين علف DairyOne (https://apps. dairyone.com/feedcomposition/؛ انظر المورد عبر الإنترنت 1 للحصول على القائمة الكاملة من معلمات المناخ والفسيولوجيا)، وتم تقدير الكتلة الرطبة بالكيلوغرام من العلف المستهلك أيضًا من هذه المحاكاة. تم أيضًا تعديل بيانات الطعام المستهلكة باستخدام انحدارات متعددة الحدود من الدرجة الثانية لتقييم أنماط استهلاك الطعام الشهري إحصائيًا في المواطن المتوقعة للأنواع الأحفورية المأخوذة.

نتائج

تحليل نظائر الكتلة المتجمعة وREER

تراجعت نسب مؤشر REE المحسوبة لمعظم العينات التي تم تحليلها هنا تحت قيمة 1.0 مما يشير إلى القليل من التغيير الدياجيني في العينات (الجدول 1). كانت الاستثناء هي عينة Promegatherium (FMNH P14511)، التي أعطت قيمة مؤشر REE قدرها 36.6 واعتبرت على الأرجح أنها تغيرت دياجينيًا أو عينة غير نقية. كانت عينة Eremotherium (UF 312730) لديها أدنى نسبة REE (0.02) التي تتماشى مع قيم مؤشر REE من العاج إلى العظم القشري المحددة لأفراد آخرين تم العثور عليهم في Halie 7G (النطاق؛ 0.004-0.01؛ MacFadden et al. 2010؛ المورد عبر الإنترنت 2). كل من عينة Eremotherium (نسبة REE؛ 0.42) وعينة Nothrotheriops (نسبة REE؛ 0.86) المستردة من Inglis 1 A وLeisey 1 A، على التوالي، أيضًا كانت ضمن التباين المعروف لقيم مؤشر REE المسجلة للأحفوريات من الأفراد في كلا الموقعين (النطاق؛ 0.33-0.72؛ MacFadden et al. 2010). أدت بيانات النظائر المتجمعة المقابلة من عينة Promegatherium FMNH P14511 إلى تقدير متوسط درجة حرارة قدرها وتم تحديدها على أنها منخفضة جدًا لدرجة حرارة جسم واقعية لثديي كبير، حيث كانت أكثر تشابهًا مع الحد الأدنى من تقديرات درجة حرارة الجسم لل ectotherms التي تم تحديدها من خلال تحليلات سابقة مماثلة (Eagle et al. 2011، 2015؛ Griffiths et al. 2023). يُلاحظ أن الانزياح الكربوني-الفوسفات كان ضمن حدود معظم الفقاريات الحالية، مما يشير إلى أن تغيير الجزء الكربوني قد يكون طفيفًا بدلاً من كونه واسع النطاق (الشكل 1ب).
الشكل 3 مخطط صندوق وشوكة لتقديرات درجة الحرارة المحيطة للبيئات التي تسكنها الأنواع الأربعة من كسالى الأرض المنقرضة المأخوذة
على الرغم من انخفاض مؤشر REE الخاص به، أظهرت عينة Eremotherium UF 312730 تكرارًا ضعيفًا بين التحليلات (خطأ قياسي مرتفع؛ الجدول 1) وأنتجت تقديرات درجة حرارة نظائر متجمعة والتي من غير المحتمل أن تكون في نطاق فسيولوجي. كان هذا على الأرجح نتيجة للتلوث العضوي الذي لم يتم إزالته بنجاح من خلال التنظيف المسبق لعينة السن و/أو مراحل التنظيف غير المكتملة من تنقية العينة، حيث لم يشير أي من الفوسفات وREE إلى الدياجينيس (الجدول 1؛ الشكل 1ب). تدعم هذه الفرضية القيم المرتفعة غير الطبيعية الملاحظة (الشكل 1ج؛ المورد عبر الإنترنت 2) في هذه العينة والتي تشير إلى وجود تلوث يخلق تداخلات إيزوبارية لقياسات النظائر المتجمعة. كانت العينات الثلاث المتبقية من الأسنان لعينة Eremotherium وNothrotheriops وMegatherium التي كانت تحتوي على REE منخفض وcarbonate-phosphate انزياحات في نطاق الثدييات الحالية (الجدول 1؛ الشكل 1ب) قد yielded درجات حرارة متجمعة قدرها , و , و ، على التوالي، مع توافق جيد بين التحليلات المكررة و القيم المنخفضة نسبيًا (الشكل 1ج). وفقًا للمعايير المعمول بها، لا يوجد سبب للاشتباه في أن هذه العينات قد تعرضت للتأثيرات الدياجينية أو العيوب التحليلية، وبالتالي فإن التقديرات الناتجة من المحتمل أن تمثل درجات حرارة الجسم الأساسية لكسالى الأرض المنقرضة.

تقديرات المناخ القديم ومعلمات الفسيولوجيا

تم نمذجة مجموعة واسعة من درجات الحرارة المحيطة لأربعة أجناس من كسالى الأرض التي تم تحليلها (الشكل 3). تم التنبؤ بأن أكبرها، Eremotherium، قد عانت من درجات حرارة محيطة متغيرة (النطاق: ) في نطاقه الجغرافي الواسع، على الرغم من أن القيم الأدنى تم تحديدها بواسطة القيم الشاذة. كانت الفروق المتوقعة في درجات الحرارة المحيطة لـ Nothrotheriops مشابهة في القيم المطلقة لـ Eremotherium، ولكن تم تقدير درجات حرارة محيطة دنيا أكثر برودة لهذا النوع (النطاق: ج). ومع ذلك، تم التنبؤ بأن كل من ميغاثيريوم وميلودون كانا يسكنان بشكل أساسي بيئات باردة إلى معتدلة، حيث تتداخل قيم درجات الحرارة الدنيا والقصوى مع نطاقات درجات الحرارة لإيريموثيريوم ونوثروثيريوبس (الشكل 4). لم تتجاوز قيم معدل الأيض الأساسي المقدرة لجميع الأنواع الأربعة المأخوذة عينات منها (الشكل 4)؛ كانت جميع تقديرات معدل الأيض أقل من تلك التي قدرتها التحليلات السابقة لفاريña (2002).

محاكاة غرف الأيض ونطاقات المنطقة الحرارية المحايدة

كشفت تحليلات غرفة الأيض عن مناطق حرارية محايدة واسعة للفرو الكثيف (1,300-2,000 شعرة/ عبر جميع الأجناس الأربعة التي انخفضت مع زيادة سمك الفراء
الشكل 4: رسم بياني عمودي لمعدل الأيض الأساسي المقدر (BMR) ومعدل الأيض الميداني (FMR) في الكسلان الأرضي. القيم بوحدات الواط (W) المحسوبة من معادلة قياس الأيض في الزينارثرا (وايت وسيمور 2004) مقارنة بتلك (على اليمين) المستندة إلى ‘الأيض المشيمي النموذجي’ (كلايبر 1932؛ فاريña 2002)
(الشكل 5). كان لكل من إيريموثيريوم وميغاثيريوم مناطق حرارية متعادلة مشابهة بسمك فراء يبلغ 10 مم (النطاق؛ ، المجمعة لكلا النوعين)، والتي تتداخل مع القيم المقدرة لدرجات الحرارة المحيطة التي كان من الممكن أن تتعرض لها كلا الجنسين طوال توزيعهما الجغرافي (الشكل 3). من المدهش أن الميجاثيريوم، الذي كان لديه فرو كثيف وسمك معطف يبلغ 50 مم، كان لديه أدنى درجة حرارة حرجة دنيا بين الأربعة أجناس. ) التي تعادل التقديرات السابقة لميلودون ذو الفراء الكثيف ضمن منطقة الحرارة المحايدة (فاريña، 2002)، ولكنها أقل من تقديرات درجة الحرارة الحرجة الدنيا لإيريموثيريوم مع نفس ميزات الغلاف الجلدي المطبقة هنا ( ظروف الفراء الكثيف لميلودون ونوثروثيريوبس مع سماكات الفراء كما أن لديها مناطق حرارية محايدة تتداخل مع نطاقات درجات الحرارة المحيطة المحتملة التي من المحتمل أن يكون كل نوع قد تعرض لها، بالإضافة إلى درجات حرارة محيطة تتجاوز بكثير أدنى التقديرات التي تم نمذجتها هنا. و ، على التوالي). اتجاهات متسقة لمناطق حرارية محايدة ضيقة ناتجة عن فراء رقيق ( شعر ) مع تقديرات درجة حرارة حرجة دنيا أعلى من ودرجة الحرارة الحرجة العليا
تقديرات تتراوح من تم تحديدها لجميع الأنواع الأربعة التي تم تحليلها.

محاكاة المناخ الصغير وتحمل درجات الحرارة للجلد

تمت ملاحظة نتائج مماثلة تفضل تغطية الفراء الكثيف بشكل عام عندما تم إدخال عدة متغيرات مناخية، بما في ذلك الرطوبة النسبية، والهطول، وسرعة الرياح، ودرجة الحرارة المحيطة في نماذج المناخ القديم المحسنة. ومن الجدير بالذكر أن Eremotherium تم تقديره بأنه لا يعاني من ضغط حراري في جميع البيئات التي تحتوي على عزل من الفراء الكثيف بسمك 10 مم (1,300-2,000 شعرة). )، ولكن كان من الممكن أن يتعرض لضغط بارد بارز في أقصى شمال نطاق توزيعه، على الرغم من عدم حدوث ضغط حراري في أقصى الجنوب خلال أبرد أشهر السنة في كلا الموقعين المأخوذين (الشكل 6). تتناقض هذه النتيجة مع الاتجاهات المحددة التي أشار إليها تغطية الفراء النادرة في أقصى الشمال وأقصى الجنوب، والتي أظهرت باستمرار عدم وجود ضغط حراري خلال أشهر السنة الأكثر دفئًا فقط. كثيف تم تحديد أن عزل الفراء لـ Eremotherium كان ضارًا في الغالب بسبب النتائج
الشكل 5 المناطق الحرارية المحايدة لـ إيريموثيريوم؛ ب. ميغاثيريوم؛ ج. ميلودون؛ و د. نوتروثيريوبس. تم اشتقاق مناطق الحرارة المحايدة المقدرة من محاكاة غرف التمثيل الغذائي مع كثافات فرو مختلفة وسماكات معطف.
في إجمالي إجهاد الحرارة، باستثناء أشهر الشتاء فقط في أقصى خطوط العرض الشمالية وعلى مدار السنة في أقصى خطوط العرض الجنوبية. وعلى العكس من ذلك، تم نمذجة إيريموثيريوم ليكون بلا إجهاد حراري مع فراء رقيق يبلغ 50 مم، ولكن مع تغطية فراء أكثر كثافة مما يؤدي إلى إجهاد حراري مستمر في المناطق النيوتropicalية (الشكل 6).
أظهرت نماذج المناخ القديم المتطورة أن الميجاثيروم كان سيتجنب الضغط الحراري على مدار السنة مع عزل كامل للجسم من الفرو الكثيف بسمك 30 مم في جميع البيئات الثلاث المحاكية، بينما أدى الفرو الخفيف الكامل (بكلا العمقين والكثافات) إلى الضغط البارد على مدار السنة (الشكل 6). كما كان من الممكن أن يتجنب كل من الميلودون والنوتروثيريوبس أي ضغط حراري مع عزل من الفرو الكثيف (الكثافة؛ شعر ). على وجه الخصوص، الميلودون
تم تصميمه ليكون بلا إجهاد حراري مع فرو بسمك 10 مم خلال أشهر السنة الأكثر حرارة، وأيضًا مع فرو بسمك 30 مم خلال أشهر السنة الأكثر برودة في المناطق الشمالية والمتوسطة (الشكل 6). تم تحديد أن كل من الميلودون في أقصى جنوب نطاقه ونوثروثيريوبس عبر نطاقه الجغرافي بالكامل لا يعانيان من إجهاد حراري مع فرو كثيف بسمك 50 مم.

ساعات النشاط الأيضي النشطة اليومية

تقديرات الساعات في اليوم التي يمكن فيها الحفاظ على معدل الأيض الأساسي دون البحث عن الظل للظروف الجلدية التي تؤدي إلى أقل إجهاد حراري (حسب تغطية الفراء، العمق، والكثافة) عبر الجغرافيا المعروفة
الشكل 6 الاتجاهات النموذجية لنشاط الأيض اليومي لكسلان الأرض (كيلوجول في اليوم) -المحور) استجابةً للعوامل المناخية على مدار السنة (شهر؛ -المحور) عبر أقصى الشمال، والوسط، وأقصى الجنوب لكل جنس. معدل الأيض المستهدف لإريموتيريوم
توزيعات كل نوع تم تحليلها كانت بمتوسط حوالي 15 ساعة في اليوم (الشكل 7). بينما تم تصميم الميجاثير لتكون لديها نفقات أيض كبيرة في المناخات الباردة (مثل النطاقات الجغرافية الشمالية لـ Eremotherium؛
يوم ; معدل الأيض المستهدف لميغاثيروم = 51,944 يوم ; معدل الأيض المستهدف لميلودون يوم ; معدل الأيض المستهدف لنوتروثيريوبس
نطاقات جغرافية متوسطة العرض وجنوبية لميجاثيريوم)، أكبر عدد من ساعات النشاط الأيضي المستهلكة يوميًا ( تم التنبؤ بها في النطاقات العرضية الجنوبية لكل من ميغاثيريوم وميلودون.
الشكل 7: رسم بياني عمودي لعدد الساعات المتوسطة في اليوم التي يمكن أن يتجاوز فيها معدل الأيض المستويات القاعدية في أقصى الشمال والوسط وأقصى الجنوب لكل جنس من أجناس الكسلان الأرضي مقارنةً بالثدييات المشيمية الحالية ذات الأحجام الجسمانية المماثلة وكسلان الأشجار الحاليين. تظهر كظلال نظائر حجم الجسم من اليسار إلى اليمين: الفيل الأفريقي في السافانا (Loxodonta africana؛ لجنس Eremotherium)، الفيل الآسيوي (Elephas maximus؛ لجنس Megatherium)، وحيد القرن الهندي.
(وحيد القرن الهندي؛ بالنسبة لـ Mylodon)، الدب البني (Ursus arctos؛ بالنسبة لـ Nothrotheriops)، Bradypus و Choloepus. تم الحصول على الظلال من PhyloPic.org (Michaud M، Taylor J، Keesey M، Morrow C، Jenkins X) وهي تحت الملكية العامة (ساعات النشاط الإجمالية للأنماط الحالية كانت مستندة إلى ساعات النوم اليومية المبلغ عنها في الأدبيات (Stelmock و Dean 1986؛ Deka و Sarma 2015؛ Gravett et al. 2017؛ Cliffe et al. 2023)
أظهر نوتروثيريوبس قيمة ثابتة للنشاط الأيضي الأقصى تبلغ 15 ساعة في اليوم. في جميع أنحاء نطاقها الجغرافي (الشكل 7). تشير البيانات المقارنة حول النفقات الطاقية اليومية إلى أن نظائر الحجم الجسدي الحالية لـ Eremotherium و Megatherium و Mylodon (مثل الفيل الأفريقي، الفيل الهندي، ووحيد القرن الهندي، على التوالي) جميعها لها قيم تبلغ 22 ساعة في اليوم. (ديكا وسارما 2015؛ غرافيت وآخرون 2017) بينما تم قياس ذلك بالنسبة لنوتروثيريوبس الأصغر بكثير، وهو نظير لدب بني (أورسوس أركتوس)، عند 20 ساعة في اليوم (ستيلموك ودين، 1986). باختصار، تم التنبؤ بأن جميع أنواع الكسلان الأرضي لديها نشاط أيض أقل من نظائر الثدييات المشيمية الحديثة ذات الكتل الجسمانية المماثلة، لكنها أكثر نشاطًا بكثير من الكسلان الشجري الحالي (على سبيل المثال، ; كليف وآخرون 2023).

الاستهلاك الغذائي الشهري

تشير التوزيعات المقدرة لاستهلاك الطاقة عبر الأجناس الأربعة إلى عدم وجود اتجاهات ذات دلالة إحصائية في استهلاك الطعام الشهري مع تغير البيئة (الشكل 8). بينما تم إثبات زيادة عامة في استهلاك الطعام بالنسبة للميغاثير الكبيرة الجسم، تم نمذجة إيريموثيريوم ليشهد زيادة قليلة أو معدومة في
استهلاك الطعام ( ) في جميع أنحاء توزيعها الجغرافي الواسع. تم تصميم جنس ميغاثيريوم ليشهد زيادة معتدلة في تناول الطعام الشهري في أقصى خطوط العرض الجنوبية , بينما توقعت محاكاة ميلودون زيادة مماثلة في تناول الطعام في خطوط العرض الجنوبية , مع تداخل ملحوظ مع تلك الخاصة بميغاثيريوم. أخيرًا، أظهرت اتجاهات البيانات لنوتروثيريوبس انخفاضًا معتدلًا في تناول الطعام الشهري في خطوط العرض الجنوبية خلال أشهر الصيف.

نقاش

الكيمياء الجيولوجية للأسنان الأحفورية وتفسيرات الفسيولوجيا

تشير التشابهات في درجات حرارة النظائر المتجمعة من الأسنان المستعادة من مواقع مختلفة، إلى جانب قيم مؤشر REE المنخفضة، إلى أن هذه التقديرات قد تعكس درجات حرارة الجسم الأساسية الفعلية للكسالى الأرضية المنقرضة. هذه الدرجات حرارة أقل بكثير من درجات الحرارة المستمدة من النظائر المتجمعة لغيرها من الثدييات الكبيرة
الشكل 8 مخطط انتشار لكتل الطعام الرطب المستهلكة يوميًا على مدار السنة لكل نموذج من نماذج الجلد التي تظهر الحد الأدنى من الضغط الحراري للمدى الشمالي (الأزرق)، والوسط (الأخضر)، والمدى الجنوبي (الأحمر). تم ملاءمة بيانات كل جنس من الكسالى الأرضية مع انحدارات متعددة الحدود من الدرجة الثانية: إيريموثيريوم (y , ); ميغاثيريوم ( , ); ميلودون ( ; نوتروثيريوبس . 0.669)
الثدييات بما في ذلك قيم تقارب لفيل هندي وماموثات منقرضة (إيجل وآخرون 2010)، و لحنش بالين (غريفيث وآخرون 2023). تم الحصول على مادة الأسنان من عينات الكسالى الأرضية التي تظهر القليل من التغيرات الدياجينية من ثلاثة بيئات ترسيب مختلفة وأنتجت تقديرات درجات حرارة ضمن نطاق درجات حرارة الجسم للزواحف الحالية. نظرًا لأن التقديرات لنوتروثيريوبس تم اشتقاقها من عينة سن واحدة ذات نسبة REE مرتفعة نسبيًا (0.86)، يجب توخي بعض الحذر في تفسير النتائج من هذه العينة. القياسات التي تم الحصول عليها من عينة إيريموثيريوم UF 312730 تعطي قيمًا أعلى في كل من و مما قد يشير إلى تلوث العينات مما يؤدي إلى تداخل إيزوباري وتقدير مفرط لاحق لـ ودرجات حرارة المستمدة بشكل غير واقعي (الشكل 1c). بعيدًا عن ذلك، تشير هذه البيانات الجديدة إلى أن الكسالى الأرضية المنقرضة قد كانت غير متجانسة حراريًا، ومن المحتمل أنها تقلبت في درجات حرارة جسمها الأساسية كوسيلة رئيسية لتنظيم الحرارة. تقديرات درجة حرارة الجسم المستمدة من تحليلات النظائر المتجمعة هي بالفعل ضمن نطاق الخطأ لدرجات حرارة الجسم المعروفة للزواحف الحالية ( )، مما يؤكد مرة أخرى قيم درجة حرارة الجسم المحددة لمحاكاة نيش مابير ويقترح أن التشابهات الشاملة الملحوظة في الفسيولوجيا الحرارية هي سمة سلفية في هذه المجموعة.
تقديرات النظائر المتجمعة لدرجة حرارة الجسم، ومع ذلك، هي في الطرف الأدنى من محاكاة نيش مابير للتقدير المتوسط لـ المحدد لـ UF 95869. يمكن تفسير ذلك من خلال القياسات المأخوذة من مواقع الأسنان المحددة التي تعكس فترات محتملة من السكون مشابهة لما تم اقتراحه للطرودونتيدات (تاجليافينتور وآخرون 2023)، قد تعكس التقديرات من الأسنان أيضًا اختلافات طفيفة في درجة حرارة الجسم الأساسية ودرجة حرارة تكلس الأسنان كما تم استكشافه في تفسير قيم النظائر المتجمعة المستمدة من أسماك القرش الداخلية الإقليمية (غريفيث وآخرون 2023). يجب استكشاف هذه الفرضيات في التحليلات المستقبلية التي تركز على أخذ عينات من مناطق محددة من الأورثودنتين الداخلي لتحديد ما إذا كانت درجات حرارة الجسم الأساسية تتقلب طوال حياة الحيوانات. تم تطبيق مثل هذه الطرق سابقًا على الأسنان من إيريموثيريوم لاستنتاج الظروف البيئية خلال حياة فرد واحد، باستخدام تحديدًا قياسات و (لارمون وآخرون 2019). قد تكون تحليلات المزيد من العينات التي تغطي نطاقًا أوسع من ظروف المناخ مشابهة للنهج الذي تم استخدامه على بقايا الديناصورات غير الطائرة (فريكي ورودجرز 2000؛ أمويت وآخرون 2006) مفيدة أيضًا. ومع ذلك، تشير قياسات النظائر المتجمعة جنبًا إلى جنب مع نتائج تحليلات النمذجة الفسيولوجية هنا بشكل جماعي إلى أن الكسالى الأرضية المنقرضة كانت لديها درجات حرارة جسم أقل ومعدلات أيض أساسية أقل من معظم الثدييات المشيمية الكبيرة.

مقارنات مع التحليلات السابقة

تشير تقديرات درجة حرارة الجسم المنخفضة التي تم الحصول عليها من أسنان الكسالى الأرضية الأحفورية إلى معدلات أيض أقل مشابهة للزواحف الحالية وتتناقض مع القيم الأعلى لمعدلات الأيض بناءً على تقدير معدل تدفق الدم (فاريلّا وآخرون 2024). افترض فاريلّا وآخرون (2024) أن أقطار الثقوب الغذائية الصغيرة لعظم الفخذ، وبالتالي معدلات تدفق الدم المقيدة، إلى جانب عادات الزواحف الحالية الحفرية والشجرية كلها مرتبطة بمعدلات أيض منخفضة. ومع ذلك، تظهر البيانات الخام من تلك الدراسة تداخلًا بين أقطار الثقوب الغذائية لعظم الفخذ ومعدلات تدفق الدم للزواحف الحالية والثدييات المشيمية ذات معدلات الأيض المرتفعة. على سبيل المثال، الثقوب الغذائية لكل من النمس العملاق (Myrmecophaga tridactyla) وختم النمر (Hydrurga leptonyx) لها قطر مقاس 0.9 مم على الرغم من عاداتها الوظيفية المختلفة تمامًا وبيئتها. نظرًا لهذه الأنواع من التناقضات في تحليلهم، نحذر من أن استخدام قياس قطر الثقب الغذائي وحده ليس مؤشرًا موثوقًا لمعدل الأيض، وقد يعكس على الأرجح أثرًا لحجم الجسم المماثل. تقديرات درجة حرارة الجسم الأساسية التي تم الحصول عليها من تحليلات النظائر المتجمعة المعتمدة هنا قد تقدم تقديرًا أكثر مباشرة لدرجة حرارة الجسم الأساسية وتوسع فهمنا للقيود الفسيولوجية التي قد تواجهها الكسالى الأرضية المنقرضة.
كما أسفرت محاكاة المناخ القديم عن قيم معدلات الأيض التي تختلف اختلافًا جذريًا عن تلك الخاصة بدراسة فاريña (2002) حول طاقة الكسالى الأرضية الحرارية. بالإضافة إلى تقديرات محسنة لمعدلات الأيض الأساسية الخاصة بالمجموعة للكسالى الأرضية، تم إدخال العديد من العوامل البيئية الأخرى بخلاف درجات الحرارة المحيطة وتم نمذجتها عبر التوزيعات الجغرافية لكل جنس لتحديد حدود أكثر دقة للحياد الحراري. أنتجت تقديرات معدل الأيض الأساسي المستمدة من معادلة قياس الزواحف (وايت وسيمور 2004) معدلات أيض تمثل فقط من تلك الخاصة بالثدييات ذات الأيض المشيمي النموذجي كما تحددها علاقة كليبر (1932) (فاريña 2002). تم التنبؤ بمعدلات أيض أساسية أقل بكثير في الدراسة الحالية، مما أدى إلى نتائج مختلفة بشكل ملحوظ لنطاقات تحمل الحرارة عبر ظروف الجلد المختلفة المختبرة مقارنة بتلك التي تم الإبلاغ عنها بالفعل. اقترحت الدراسة السابقة الرئيسية (فاريña 2002) أن ميغاثيريوم ذو الفراء القليل وكتلة جسم وأيض كان مثل الثدييات المشيمية النموذجية سيكون ضمن منطقة حرارته المحايدة عند درجات حرارة محيطة تبلغ مقابل النتائج من محاكاة المناخ الحالية، التي تقترح أن ميغاثيريوم في ظل ظروف مماثلة كان سيتعرض لضغط برودة.
تشير التحديدات التقليدية لمناطق الحياد الحراري باستخدام أيض الزواحف كبديل إلى أن
إيريموثيريوم وميغاثيريوم ذوي الفراء الكثيف وعمق الفراء 10 مم سيكونان مشابهين لتلك النماذج ذات الفراء القليل (انظر الشكل 5). ومع ذلك، أظهرت نماذج المناخ الصغير المستخدمة هنا أن العوامل البيئية بخلاف درجة الحرارة المحيطة (أي، الرطوبة النسبية، سرعة الرياح، إلخ) والنشاط الأيضي (أي، الحركة، البحث عن الطعام، إلخ) كانت مسؤولة عن النتيجة المتوقعة بأن سماكات 10 مم الكثيفة ستفيد إيريموثيريوم فقط خلال أشهر السنة الأكثر دفئًا في أقصى شماله وأقصى جنوبه على مدار السنة. كما أظهرت نماذج المناخ الصغير الجديدة أن وجود فراء كثيف بسماكة 10 مم سيؤدي إلى ضغط حراري مستمر لإيريموثيريوم في المناطق الاستوائية، ومع ذلك ضغط برودة مستمر لميغاثيريوم في جميع أنحاء توزيعه الجغرافي (الشكل 6). على الرغم من أن الضغط الحراري تم نمذجته ليكون موجودًا في السكان الاستوائيين لإيريموثيريوم الذين لديهم تغطية كاملة للفراء بعمق ستسمح حالة عمق الفرو الكثيف بسمك 10 مم بعدم تعرضه لضغوط حرارية شديدة على مدار السنة، حتى خلال أشهر السنة الأكثر حرارة في نطاقه الجغرافي بالكامل. تشير النتائج للمنطقة الجنوبية الأكثر لميلودون والبيانات الخاصة بنوثروثيريوبس في نطاقه بالكامل مع عمق فرو كثيف بسمك 50 مم، إلى حماية سنوية شبه كاملة من الضغوط الحرارية في بيئاتهم المعنية، مع ضغوط برد خفيفة خلال الأشهر الثلاثة الأكثر برودة في السنة (درجات الحرارة المتوقعة؛ ). تتعارض هذه النتائج أيضًا مع التقديرات السابقة (فاريña 2002) التي تفيد بأن الحياد الحراري لميلودون، الذي يمتلك معدل الأيض الأساسي نصف معدل الأيض المشيمي النموذجي وغطاءً كثيفًا بعمق 40 مم، سيكون عند درجة حرارة محيطة من .

الحياد الحراري للميغاثير والتوزيع الجغرافي

أظهرت البيانات من تحليلات المناخ القديم رؤى حاسمة حول عتبات التنظيم الحراري للكسالى الأرضية المنقرضة وكيف كانت ستؤثر على مظهر حياتها من خلال حالة غلافها. على وجه الخصوص، تشير النتائج المتعلقة بـ Eremotherium إلى أن مجموعة واسعة من نماذج الغلاف قد تكون ممكنة اعتمادًا على درجات الحرارة المحيطة في بيئاتها السابقة. نظرًا لأن Eremotherium كان له توزيع جغرافي يمتد شمالًا حتى فرجينيا في أمريكا الشمالية وجنوبًا حتى البرازيل في أمريكا الجنوبية، فمن المحتمل أن يكون لهذا الجنس درجات متفاوتة من تغطية الفراء عبر نطاقه الجغرافي. كانت تغطية الفراء الكاملة على الجسم بشكل متقطع ستكون مفيدة في النطاقات المتوسطة مع عدم تعرضه لضغوط حرارية على مدار السنة. كما أن التغيرات في عمق الفراء كانت ستكون مفيدة خلال الأشهر الأكثر دفئًا لـ Eremotherium في أقصى شمال نطاقه، على الرغم من عدم وجود دليل على تساقط الفراء في الزواحف الحالية. ومع ذلك، فإن الفراء الكثيف بسمك 10 مم كان يمكن أن يمنع الضغوط الحرارية في درجات الحرارة المحيطة المتوسطة.
الشكل 9 أ. خريطة التوزيع الجغرافي لمواقع بقايا الحفريات لعناكب الأرض Eremotherium وMegatherium وMylodon وNothrotheriops من عصر البليستوسين التي تمتد من 2.5 مليون سنة إلى 11,000 سنة مضت. المناطق المظللة باللون الأزرق تمثل مدى الأنهار الجليدية خلال أقصى فترة جليدية قبل 21,000 سنة، والدوائر الحمراء والخضراء تشير إلى المواقع التي تم تصوير الفن القديم فيها. البيانات: تم إنشاء خريطة الحدود السياسية الأساسية باستخدام ‘خريطة العالم’ من QGIS3؛ تم تعديل وتحديث طبقة الجليد المرسومة يدويًا بعد البيانات المنشورة (Ehlers and Hughes 2018)؛ تم استرجاع البيانات المكانية لتوزيعات الأنواع من
نطاقات درجات الحرارة كان من المتوقع أن يتعرض إيريموثيريوم (الشكل 6). علاوة على ذلك، لأن إيريموثيريوم عبر برزخ بنما إلى أمريكا الشمالية ولكنه تراجع لاحقًا إلى أمريكا الجنوبية خلال أقصى فترة جليدية (ماكدونالد ولوندليوس 2009)، فمن الممكن أن يكون قد خضع لهجرات خلال أوقات الموسم الشديد. نظرًا لوجود نموذج لوجود ضغط بارد معتدل مع فرو كثيف بعمق 10 مم وعدم وجود حيادية حرارية مع عمق فرو يبلغ 50 مم في أقصى النطاقات الجنوبية، فإنه يُقترح أن إيريموثيريوم كان متكيفًا جيدًا مع الظروف النيوتروبية التي يُتوقع أن تكون سائدة في
قاعدة بيانات علم الأحياء القديمة (paleobiodb.org). ب. إعادة بناء بالغ ورضيع من إيريموثيريوم أثناء تبديل الفراء الموسمي وهما يتدحرجان في بركة طينية خارج عرين كائن نوتروثيريوبس الناشئ في أوائل الربيع في فلوريدا خلال العصر الجليدي (الموقع موضح بالدائرة الخضراء في اللوحة 9a). ج. إعادة بناء ميغاثيريوم وهو يتحمل برودة ورياح بامباس الأرجنتين مع فرو سميك (الموقع موضح بالدائرة الحمراء في اللوحة 9a). فن باليو من هايدر سيد جعفري
أقصى مدى جنوبي لنطاقه (الشكل 9أ). تتفق هذه الفرضية بشكل ملحوظ مع غياب إيريموثيريوم في أمريكا الشمالية خلال أقصى حد جليدي.
على الرغم من تشابه حجم الجسم مع إيريموثيريوم، من المتوقع أن يكون جنس ميغاثيريوم قد عانى من إجهاد بارد مستمر مع تغطية شعرية نادرة فقط في جميع أنحاء نطاقه الجغرافي. يمكن تفسير ذلك من خلال معرفة أن ميغاثيريوم كان يعيش في بيئات باردة وجافة في السهول في الأرجنتين وتشيلي. باستخدام أساليب نمذجة مشابهة للفيزيولوجيا كما هو موضح هنا، كشفت تحليلات الحساسية السابقة (لوفليس وآخرون 2020) أن الظروف المحيطة
كانت قيم درجة الحرارة وسرعة الرياح هما العاملان المناخيان اللذان يؤثران بشكل كبير على القيم المتوقعة لمعدل الأيض الميداني، في حين أن متغيرات الرطوبة النسبية وتغطية السحب كان لها تأثير ضئيل على الأيض المحاكى. وبالتالي، قد تكون الفراء الكثيف لميجاثيريوم ضرورية لمواجهة التعرض الأكبر للعوامل الجوية (مثل الرياح أو الأمطار أو الثلوج) في البيئات الأكثر انفتاحًا. تعكس حالة الجلد هذه نتائج محاكاة المناخ الحالية التي تشير إلى أنه مع تغطية كاملة للجسم بالفراء بعمق 10 مم، كان ميجاثيريوم سيعاني من إجهاد برودة مستمر بغض النظر عن كثافة الفراء، لكنه كان سيحقق حيادية حرارية كاملة في جميع التوزيعات العرضية مع الفراء الكثيف. شعر فرو بعمق 30 مم. وعلى العكس، من خلال وجود فرو كثيف بعمق 50 مم، كان من الممكن أن يتعرض الميجاثيريوم لضغط حراري مستمر بسبب العزل الزائد (الشكل 6). علاوة على ذلك، بالنظر إلى أن أعماق الفرو نتج عن ذلك إجهاد حراري متوقع لـ Eremotherium بحجم مماثل (باستثناء أبرد الشهور في أقصى شمال نطاقه)، ومن المعقول أن نفترض أن الظروف النيوترونية قد عملت كمتدرج حراري يقيّد Megatherium إلى خطوط العرض الجنوبية واحتلال المناطق الغربية من بيرو والإكوادور (الشكل 9أ؛ بوجوس وسالاس 2004). تدعم هذه الفرضية نماذج التوزيعات البيئية لمختلف الزينارثرانات من العصر الجليدي التي قدمها فاريلّا وآخرون (2018) والتي تظهر تداخلًا ضئيلًا أو معدومًا بين النطاقات الجغرافية لـ Eremotherium وMegatherium مما يشير إلى اختلافات واضحة في تفضيل الموائل.

تنظيم الحرارة في الميلودون والنوتروثيريوبس

تم تصميم كل من Mylodon و Nothrotheriops ليحتاجا إلى كثافة عمق الفراء عبر نطاقات التعرض لدرجات الحرارة المحيطة المتوقعة، والتي تم ملاحظتها مباشرة في الفراء المحفوظ من كلا الجنسين (ريدج وود 1901؛ هاوسمان 1929). ومع ذلك، نظرًا لأن الميلودون كان معروفًا بأنه يسكن جبال الأنديز في أمريكا الجنوبية، فإن النوثروثيريوبس كان من الممكن أن يتعرض لدرجات حرارة المناخ القديم التي تصل إلى (الشكل 3)، كان من المتوقع أن يكون هناك تنظيم حراري سلوكي في المناخات الباردة بالإضافة إلى وجود غطاء سميك من الفرو. إحدى الفرضيات هي أنه كان يمكن الحفاظ على الحياد الحراري من خلال الاحتماء في الكهوف خلال فترات الضغط الحراري. تكشف تقديرات متوسط درجة الحرارة السنوية للكهوف في الجنوب الغربي الأمريكي حيث تم توثيق بقايا Nothrotheriops عن نطاقات درجات حرارة مستقرة كانت ستوفر الحماية ضد فترات الضغط الحراري (6.6-17 “; ماكدونالد 2022؛ الشكل 9ب). المواقع الخاصة بـ Mylodon في تشيلي مثل كهف تشيكا تظهر نطاقًا مشابهًا من درجات الحرارة المستقرة (4.83-17.4 “; نهمه وآخرون 2022). فرضية ثانية هي أن كلا النوعين قد يكونا قد خضعا
فترات من السكون و/أو السبات كوسيلة سلوكية لتحمل البرد.
كشفت الدراسات النظيرية والتافونومية عن ‘النمر’ الباتاغوني المنقرض (Panthera onca mesembrina) أن علامات العض المنسوبة إلى هذا المفترس المنقرض قد وُجدت على عظام الأفراد اليافعين والبالغين من ميلودون في كهوف ميلودون (بريفوستي ومارتن 2013؛ مارتن 2018). بينما كان ميلودون أصغر من الميجاثير، إلا أن الفرد البالغ كان لا يزال أكبر بكثير من اليغور، مما يشير إلى أن الأفراد الذين وُجدت عليهم علامات العض في الكهوف قد تم افتراسهم أثناء خضوعهم للسبات، عندما كانوا الأكثر عرضة للهجوم. يمكن أيضًا ربط الدعم الإضافي للفرضيات المذكورة أعلاه باستراتيجيات التنظيم الحراري السلوكية التي لوحظت في الأرماديلوس الحديثة. على سبيل المثال، يُعرف الأرماديلو ذو التسعة أحزمة (Dasypus novemcinctus) بأنه يسعى للملجأ خلال فترات البرد الشديد ويستخدم مزيجًا من السبات مع البحث المتقطع عن الطعام (مكبي وباكر 1982؛ نايت وآخرون 2020). تم توثيق السبات المنتظم أيضًا في الأرماديلو بيشي (Zaedyus pichiy) (سوبرينا وبويلي 2007؛ سوبرينا وآبا 2014)، وهو النوع الوحيد المعروف من الزينارثرا الذي يخضع لسبات حقيقي.
من الممكن أيضًا أن يكون ميلودون ونوثروثيريوبس قد كان لديهما مستويات نشاط أعلى من تلك المقدرة في التحليل الحالي. ومع ذلك، فإن المتطلبات الأيضية الأعلى في الكسلان الأرضي الأصغر ستكون مفيدة فقط في درجات حرارة أقل بكثير من نطاق درجات الحرارة المحيطة المتوقعة لكلا النوعين. وبناءً عليه، فإن تفسيرًا محتملاً آخر لطاقة حرارتهم ينطوي على وجود نظام غذائي أغنى بالطاقة من نظائرهم من الميجاثير. اقترح تحليل حديث للأحماض الأمينية (تيخادا وآخرون 2021) أن نظام ميلودون الغذائي قد يتكون من بروتين حيواني بسبب توقيعات النيتروجين الموجودة في مواد عظامهم المتحجرة. على الرغم من أنه لم يكن من الممكن تحديد التوقيتات الدقيقة لتوقيعات النيتروجين التي تشير إلى التغذية الكالورية، فإن استهلاك البروتين الحيواني كان له فوائد تنظيم حراري لميلودون. يوفر النظام الغذائي الغني بالبروتين بعض المرونة في الحياد الحراري في المناخات الباردة، ومن الممكن أن يكون ميلودون قد استهلك المواد الحيوانية بشكل متقطع خلال الأشهر الأكثر برودة، بالإضافة إلى البحث عن الأعشاب (فاريل وآخرون 2023). على العكس من ذلك، لم تُظهر نظائر النيتروجين المستعادة من بقايا نوثروثيريوبس دليلًا على نظام غذائي آكل للحوم، بل على العكس، كان هذا الجنس يتناول نباتات صحراوية غنية بالكربوهيدرات مثل الجلوبمالو وجذور الكاسافا (بوانار وآخرون 1998).
كان لدى الميلودونتيدات أيضًا عظام جلدية بأشكال وأحجام متنوعة قد وفرت الحماية ضد المفترسات (كولينز 1933؛ مكدويل 2018؛ توليدو وآخرون 2021). نظرًا للطبيعة الوعائية لهذه الهياكل في الكسلان الأرضي بشكل عام، افترض كارتيل وبوهوركويز (1986) أن العظام الجلدية قد لعبت دورًا في التنظيم الحراري أيضًا. ومع ذلك، فإن الفحوصات النسيجية
لعظام ميلودون بواسطة توليدو وآخرون (2021) لا تظهر أي ميزات تشير إلى تخزين الأنسجة الدهنية على غرار تلك الموجودة في الأرماديلوس الحالية (كرمبوتش وآخرون 2009). بينما لم يتم إجراء قياسات التوصيل الحراري بعد، فقد أشار توليدو وآخرون (2021) إلى أن فرو الميلودونتيدات مثل ميلودون كان سيوفر كمية أكبر من العزل مقارنة بالعظام الوعائية.

العوامل الجلدية والسلوكية

تم إعادة بناء الكسلان الأرضي العملاق تاريخيًا على أنه نوع ذو فرو كثيف، والبيانات الجديدة المقدمة من محاكاة نيش مابير لا تتحدى كثيرًا هذا المظهر الجلدي. على العكس، تتعارض نتائج النمذجة مع إعادة بناء الجلد الشبيه بالفيل ذو الفرو القليل من إيريموثيريوم وميغاثيريوم (فارينا 2002). على الرغم من التناقضات الواضحة بين الدراسات حول الحياد الحراري وتأثيراته على جلد الكسلان الأرضي، يمكن الآن استنتاج العديد من التفضيلات السلوكية والبيئية الحرجة. على سبيل المثال، تشير بنية الأسنان الدقيقة وأنماط التآكل الدقيقة (غرين وكالتوف 2015)، بالإضافة إلى التحليلات النظيرية (دانتاس وآخرون 2021)، إلى نظام غذائي نباتي بالكامل لإيريموثيريوم وميغاثيريوم على الرغم من التكهنات السابقة بأنهما كانا ربما آكلي لحوم أو/و يمارسان سلوك النهب (فارينا وبلانكو 1996). في الواقع، كان من الممكن أن يتكون نظام كلا النوعين بشكل أساسي من الأوراق والنباتات العالية الأخرى (بارغو وآخرون 2006؛ غرين وكالتوف 2015؛ بوشيرينس وآخرون 2017)، والتي تتداخل مع تلك الخاصة بأقرب أقاربهم الحاليين من نوع براديبوس (كليف وآخرون 2015؛ باولي وآخرون 2016؛ بريسلي وآخرون 2019؛ ديلسوك وآخرون 2019؛ تيخادا وآخرون 2024). الأوراق فقيرة في القيمة الغذائية مما يتوافق مع معدلات الأيض المنخفضة للغاية التي لوحظت في براديبوس فارياجاتوس (كليف وآخرون 2015؛ باولي وآخرون 2016) وقد تكون عاملًا رئيسيًا في سبب انخفاض معدلات الأيض لدى الميجاثير مقارنة بالثدييات ذات الحجم الكبير.
واحدة من النتائج الرئيسية لدراستنا هي أن الميجاثير من المتوقع أن يكون لديها معدل أيض أساسي أكثر تشابهًا مع الزينارثرا الحديثة مقارنة بالثدييات المشيمية الأخرى. ومع ذلك، كان من الممكن أيضًا أن تظهر عددًا من السلوكيات المشابهة للثدييات الكبيرة الحالية ذات الأيض المشيمي النموذجي. من الجدير بالذكر أن سرير العظام لإيريموثيريوم في تانكي لوما في جنوب غرب الإكوادور وُجد أنه يحتوي على عدة أفراد ماتوا في مستنقع راكد. تم تفسير التجميع كدليل على سلوك الرعي في الميجاثير (ليندسي وآخرون 2020). كما تم اقتراح أن إيريموثيريوم شارك في سلوك الاستحمام لأغراض تنظيم الحرارة على غرار ما لوحظ في فرس النهر الحديث (ستيرز وآخرون 2018). من المحتمل أن يكون التباين المودل في كل من كثافة وعمق الفرو عبر النطاق الجغرافي لإيريموثيريوم مرتبطًا بشكل أوسع باستراتيجية تنظيم حراري مماثلة. الكائنات الكبيرة المهاجرة
مثل البيسون الأمريكي (Bison bison) تتخلص من فروها الشتوي الكثيف إلى فرو أخف ومتقطع يتزامن مع سلوك الاستحمام (مكميلان وآخرون 2000). كما كان من الممكن أن تشارك ميغاثيريوم الكبيرة في سلوك الاستحمام ولكنها لم تكن بحاجة إلى التخلص من الفرو بسبب حيادها الحراري على مدار السنة عندما يكون لديها فرو كثيف وعميق بسمك 30 مم (الأشكال 5 و6).
من المقترح أيضًا أن الديناميات السلوكية الاجتماعية للميغاثير المنقرضة والميودونت الكبيرة (مثل ليستودون) تتماشى مع تلك التي لوحظت في الثدييات المشيمية الكبيرة الحالية (تومسيني وآخرون 2020؛ ليندسي وآخرون 2020)، لكنها تختلف كثيرًا عن أنماط الحياة الانفرادية للكسلان الشجري الحالي (تاوب وآخرون 1999). قد تكون إحدى مزايا السلوك الاجتماعي الأكثر تعقيدًا هي الرعاية الأبوية المشتركة (التي تم استنتاجها أيضًا من سرير العظام لإيريموثيريوم في تانكي لوما؛ ليندسي وآخرون 2020) مع عدة بالغين يحميون اليافعين الأكثر عرضة للافتراس. على العكس من ذلك، تربي أنثى براديبوس الوحيدة نسلها الوحيد، وإذا سقط من الشجرة أثناء الرضاعة/الفطام، فإنها لا تستثمر الطاقة الأيضية لاستعادته (تاوب وآخرون 1999).
مرة أخرى، تتنبأ النتائج الحالية بوجود جلد كثيف الشعر لميلودون ونوتروثيريوبس كما يتضح من جلودهم المحنطة (ريدج وود 1901؛ هاوسمان 1929). نظرًا لعدم وجود طبقات عظمية أحفورية تتكون من عدة أفراد، يُقترح أن كلا الجنسين قد عاشا أنماط حياة انفرادية مشابهة لكسالى الأشجار والدرعيات الحالية (نابولي وآخرون 1990). كما أن الميلودونيدات تمتلك عظام أطراف قوية كانت مقاومة لقوى الانحناء مشابهة للدرعيات الحالية (كوبلو وآخرون 2015)، مما يشير إلى أن بعض الأجناس كانت قادرة على حفر أوكارها أو جحورها (بارغو وآخرون 2000). قد يكون سلوك العيش المنفرد كما هو الحال في الدببة البنية الحديثة (غونزاليس-برناردو وآخرون 2020)، وهو مشابه في حجم الجسم لكلا الجنسين الأصغر من كسالى الأرض، قد كان أيضًا نموذجيًا لميلودون ونوتروثيريوبس في المساكن الكهفية. علاوة على ذلك، قد يكون هناك تفضيل بيئي يمكن أن يتداخل مع الدببة وهو النظام الغذائي. نوع حديث من الكسالى ذو الأصابع الاثنين، تشولوبيوس هوفمانني، هو قريب وثيق لميلودون (بريسلي وآخرون 2019؛ ديلسوك وآخرون 2019؛ تيخادا وآخرون 2024)، وقد عُرف بأنه يستهلك أحيانًا اللحم في الأسر (هايسن 2011). كما أن أنواع جنس تشولوبيوس أكثر نشاطًا مقارنةً ببراديبوس فاريغاتوس (سونكويست ومونتغمري 1973)، الذي يتخصص في نظام غذائي من الأوراق على الأرجح مشابه لأقاربه من الميجاثير.
أخيرًا، بينما من المغري التركيز على القيود النشوء كسبب وحيد لمعدلات الأيض القاعدية المنخفضة المتوقعة للكسالى الأرضية المنقرضة، يجب أن نلاحظ أن الصفات يمكن أن تُفقد أو تُكتسب ثانويًا على مدار التاريخ التطوري لمجموعة معينة. قبل GABI، كانت مجموعة الحيوانات العاشبة الكبيرة جدًا في أمريكا الجنوبية تهيمن عليها الحيوانات غير النوتونية والأسترا بوثيرات التي كانت تتغذى على العشب والأوراق، والتي انقرضت بحلول نهاية العصر الميوسيني قبل حدوث أولى حالات الكائنات الكبيرة بشكل استثنائي.
الميغاثيرينات مثل إيريموثيريوم وميغاثيريوم (كروفت وآخرون 2020). كانت أكبر الأنواع من النوتونغولات أعضاء في عائلات توكودونتيداي، ليونتينيداي، وهومالودوثيريداي. بينما عاشت التوكودونتيدات والماكراتشينييدات جنبًا إلى جنب مع الميغاثيرينات خلال أواخر العصر الجليدي (ووصلت إلى كتل جسمية من ؛ كروفت وآخرون 2020؛ نيلسون وآخرون 2023)، الأستراوبثيرات، اللينتينيدات، والهومالودوثيرات التي تتراوح كتلها الجسمية من 375 إلى (كروفت وآخرون 2020؛ نيلسون وآخرون 2023) انقرضت في نهاية العصر الميوسيني الأوسط إلى المتأخر. بالمقابل، كانت معظم أجناس الكسلان الأرضي الميوسيني في أمريكا الجنوبية أصغر بكثير مقارنة بالنوتونغولات والأسترا بوثيرات المعاصرة، والكسلان الأرضي الأكبر الذي سيتطور في العصر البليوسيني والبلستوسيني. ; بانت وآخرون 2014: الجدول s1). كان من الممكن أن يؤدي انخفاض معدل الأيض في الكسلان الأرضي إلى تقليل متطلبات الموارد وقد يحد من المنافسة مع النوتونغولات والأستراپوثير. يمكن فهم ذلك بدوره كـ “خطوة بيئية” نحو احتلال الفجوة الغذائية للمتصفح دون الحاجة إلى الكثير من السعرات الحرارية مثل الثدييات ذات الحجم المماثل التي تتمتع بعملية أيض مشيمية نموذجية. تدعم الاتجاهات الملحوظة في كمية الطعام اليومي المستهلك في إريموتيريوم وميغاثيريوم هذه الفكرة من خلال الإشارة إلى عدم حدوث تغيير كبير في كمية الطعام المستهلكة عبر مناطقهم الجغرافية المختلفة (الشكل 7). يمكن أن يكون هذا العامل البيئي الرئيسي قد حد أيضًا من المنافسة داخل النوع على الموارد في الميغاثير، ولكنه قلل بشكل ملحوظ من ضغط المنافسة بين الأنواع مع قطعان من الفيلانيات (مثل الماموث والماستودون) التي دخلت أمريكا الجنوبية بعد GABI.

القيود

بينما أدت استخدام قيم مؤشر العناصر الأرضية النادرة في الدراسة الحالية إلى قيم ضمن الحدود المقبولة كما أبلغ عنها ماكفادن وآخرون (2010)، هناك تناقضات بين هذه الطريقة والفوسفات. البيانات. على وجه التحديد، يحتوي FMNH P14511 على مؤشر REE مرتفع بشكل غير طبيعي يبلغ 36.6 وفوسفات مقبول القياسات ضمن نطاق الثدييات المشيمية الموجودة. مؤشر REE هو التقنية الأولى المستخدمة لتقييم درجة التحول الجيولوجي في أسنان الزينارثران الأحفورية، ولم يتم اختبار هذه الطريقة بعد على عينات خارج المواقع من العصر البليوسيني والعصر الجليدي في فلوريدا. يمكن تفسير التركيز العالي لعناصر الأرض النادرة في العاج الداخلي لـ FMNH P14511 بوجوده في رواسب الميوسين الأقدم، على الرغم من أن هذا يتعارض مع الفوسفات. القياسات التي تشير إلى أنه محفوظ بشكل جيد. علاوة على ذلك، فإن عينة من Bradypus (UF 25980) المذكورة في المواد التكميلية لـ MacFadden وآخرون (2010) لديها قيمة مؤشر العاج/العظم تبلغ 30.6، وهو ما يقارن بقيمة مؤشر FMNH P14511. ومع ذلك، فإن المعايير التي وضعها MacFadden وآخرون
القيم المقبولة لمؤشر REE تستند إلى قياسات الثدييات الحاملة للمينا من نفس المواقع التي وُجدت فيها أحافير الزنجبيل بدلاً من العظام الحديثة، والعاج، والمينا. لم نتمكن من الحصول على عينات من المينا من الثدييات المشيمية غير الزنجبيلية من نفس المواقع، والتي كانت ستساعد في تحديد درجة التحول النسبي في أسنان الكسلان الأرضي المستخدمة في هذه الدراسة. كما أن درجات حرارة النظائر المتجمعة تستند أيضًا إلى يتم استبداله في ومواقع OH من الهيدروكسيباتيت، بينما تحدث تركيزات العناصر الأرضية النادرة على موقع الهيدروكسيباتيت (كون وآخرون 1999)، مما يشير إلى أن التغير الجيوكيميائي قد لا يكون هو نفسه في أي من مواقع الربط. استنادًا إلى هذه العوامل، وعدم احتساب مؤشر العناصر الأرضية النادرة لتلوث الكربونات في UF 312730 (الشكل 1c؛ المورد الإلكتروني 2)، نجادل بأن استخدام مؤشر العناصر الأرضية النادرة بمفرده لا يوفر دائمًا مؤشرًا قويًا على التحول الجيولوجي أو عدمه.
أظهر سواريز وباسي (2014) أن درجات حرارة النظائر المتجمعة الموجودة في العظام الأحفورية يمكن أن تكون متوسطة بين درجات الحرارة المحيطة خلال فترة التمعدن ودرجة حرارة الجسم. ومع ذلك، تم استخدام العظام بشكل أساسي في تحليلهم، والتي تكون أكثر عرضة للتغيرات الدياجينية مقارنة بالأنسجة السنية مثل المينا والعاج (ماكفادن وآخرون 2010). كما أن مجموعة بياناتهم لم تشمل العاج الداخلي من الزواحف الأحفورية، والذي يكون أقل عرضة للتغيرات الدياجينية من العاج الخارجي ولكنه ليس مقاومًا مثل المينا. أظهرت الكاثودولومينسنس التي أجراها لارمون وآخرون (2019) أن العاج الداخلي في سن إيريموثيريوم يظهر عددًا قليلاً جدًا من الشقوق حيث يمكن أن يحدث امتصاص المعادن، واستنتجوا أن القليل من التغيرات الدياجينية قد حدثت. علاوة على ذلك، لم يقارن سواريز وباسي (2014) تركيزات الفوسفات في عيناتهم مقارنة بعظام الثدييات الحالية، في حين أظهرت دراساتنا السابقة (إيجل وآخرون 2010، 2011) أن درجات الحرارة المستخلصة من أسنان الفقاريات المنقرضة باستخدام مطيافية الكتلة للنظائر المتجمعة تعكس درجة حرارة الجسم الأساسية للكائن، وليس متوسطًا بين درجة حرارة الجسم ودرجة حرارة البيئة.
البيانات الجديدة التي تم توليدها هنا تقدم عدة دلالات مثيرة للاهتمام حول جلد الكسلان الأرضي، والطاقة الحرارية، والبيئة، والسلوك. ومع ذلك، يُعترف بأن النتائج قد لا تكون حاسمة. إن وجود درجات حرارة جسم منخفضة تم تحديدها بواسطة قياس درجة الحرارة القديمة باستخدام نظائر متجمعة يشير إلى أن الكسلان الأرضي كان متغاير الحرارة، مشابهًا لأقرب أقاربه الأحياء. ومع ذلك، هناك حاجة إلى مزيد من التحليلات مع حجم عينة أكبر يتكون من عينات من بيئات ترسيب متنوعة للتحقق من هذه الفرضية. نظرًا لأن بعض عينات العظام القشرية لم تأتِ من نفس الأفراد كما هو الحال مع عينات الأسنان المعنية، فقد يؤثر ذلك على قيم نسب العناصر الأرضية النادرة بسبب الاختلافات المحتملة في امتصاص العناصر الأرضية النادرة.
خلال عملية التحجر. بينما تم أخذ عينات من أقسام أخرى من الأسنان الأحفورية المستخدمة في التحليلات الجيوكيميائية، لم يتم استخدامها بسبب قيود الوقت والتمويل لهذه الدراسة. يمكن أن تشير التحليلات الإضافية لهذه العينات بشكل أفضل إلى دليل على حالة غير متجانسة سابقة حيث يمكن أن تكون الطبقة الداخلية من الأورثودنتين قد سجلت إشارات نظيرية مختلفة، وبالتالي تقديرات مختلفة لدرجة حرارة الجسم.
التمثيلات النمذجة المقدمة هنا افترضت عمق وكثافة الفراء المتسقة عبر الجسم بالكامل، في حين أن هذه المعايير يمكن أن تختلف حول مساحة سطح الجسم ويجب أخذ هذا العامل في الاعتبار عند النظر في الآثار الكاملة للنتائج. يجب أيضًا أخذ عوامل أخرى مثل سمك الجلد، وسمك شعيرات الشعر، والفراء السفلي، ووجود/غياب شعيرات الشعر المجوفة في الاعتبار مع التحليلات المستقبلية. من المؤكد أننا اختبرنا فقط متوسط كتل الجسم للبالغين، وقد تختلف نتائجنا إذا تم اختبارها مع الأخذ في الاعتبار التباينات داخل الأنواع في كتل الجسم. تم استبعاد مناخ العصر البليوسيني من هذا التحليل بسبب أن متغيرات مناخ العصر الجليدي كانت أسهل للاستخدام لاستقراء قيم شهرية معقولة من بيانات المناخ الحديثة. إن تضمين بيانات مناخ عصر البليوسين سيعطي مزيدًا من الرؤية حول كيفية اختلاف الغلاف الجلدي والتغطية المرتبطة بالفراء في السكان السابقين لكل من إيريموثيريوم وميغاثيريوم. أخيرًا، تغطي الأعمار الجيولوجية للمواقع التي تم أخذ عينات منها لمحاكاة المناخ الصغير نطاقًا واسعًا من الزمن الجيولوجي يشمل كامل فترة العصر الجليدي. منذ سنوات، ولم تكن تواريخ أكثر دقة متاحة في قاعدة بيانات علم الأحياء القديمة (المورد الإلكتروني 1). كان من الممكن أن تحسن التقديرات الأكثر دقة للمناخ القديم وتاريخ مواقع الكسلان الأرضي بشكل كبير محاكاة المناخ الدقيق.

الاستنتاجات

تقدم قياسات النظائر المتجمعة في العاج الداخلي أدلة مقنعة على أن الكسلان الأرضي المنقرض كان لديه درجات حرارة جسم مشابهة لتلك الموجودة في الكسلان الحديث، وقد يكون لديه قدرة على التغير الحراري، بغض النظر عن الفروق المطلقة في أحجام الجسم. هذه النتائج، جنبًا إلى جنب مع نتائج النمذجة الفيزيائية الحيوية الجديدة في Niche Mapper، تظهر أن درجات حرارة الجسم المنخفضة والتمثيل الغذائي الأساسي، بالإضافة إلى المتغيرات البيئية، أثرت على توزيع الكسلان الأرضي الجغرافي، وقدرته على تنظيم الحرارة، وشكل الجلد. النموذج “عديم الشعر للجلد” الذي تم اقتراحه سابقًا لـ Megatherium غير قابل للدفاع عندما تكون تقديرات التمثيل الغذائي الأساسي للكسلان الحديث أكثر دقة، بينما تم اقتراح أن Eremotherium المرتبط ارتباطًا وثيقًا كان لديه شعر قليل في المناخات النيوتروبية. من المحتمل أن يكون للكسالى الأرضية العملاقة مجموعة متنوعة من مظهر الجلد من تغطية شعر كثيف وأرق في النيوتروبات إلى تغطية شعر كثيف وأسمك في البيئات الباردة. من المحتمل أن تكون الكسلان الأرضية الأصغر…
كانت تحتاج إلى فراء كثيف وسميك مدعوم بتنظيم حراري سلوكي لتحمل المناخات الباردة على مدار العام. تشير بيانات النمذجة أيضًا إلى قابلية التكيف المدهشة لفيزيولوجيا الزنارثرا في كل من أكبر وأصغر الكسلان الأرضي. علاوة على ذلك، يُعتقد أن نماذج الجلد المتوقعة لأكبر الأنواع موثوقة ما لم يتم دحضها بواسطة عينة مستقبلية تحتوي على جلد محفوظ بشكل واسع. يجب أن تفحص النماذج المستقبلية المحتملة أنواعًا أخرى (مثل Megalonyx وThalassocnus وLestodon، إلخ) وتضم دراسات نسيجية للمساعدة في تحديد مدى انخفاض معدلات الأيض الأساسية للكسلان الأرضي التي قد تكون مفيدة، وبالتالي توفير مزيد من الرؤية حول فيزيولوجيا هذه الأنواع وغيرها من الزنارثرا الأحفورية وكيف أثرت معدلات الأيض على مظهر جلدها ومعدلات نموها وبيئتها وسلوكها.
معلومات إضافية النسخة الإلكترونية تحتوي على مواد إضافية متاحة فيhttps://doi.org/10.1007/s10914-024-09743-2.
الشكر والتقدير نشكر بشكل كبير V. Sinkonen وY. Kutucu وجامعة الجمهورية لتزويدنا بنماذج ثلاثية الأبعاد من الكسلان الأرضي لهذا البحث وإضافة المنحنيات إليها. كما نشكر W. Simpson (متحف الحقل في شيكاغو) وB. MacFadden وR. Hulbert (جامعة فلوريدا) لتوفير عينات أحفورية للتحليلات الجيوكيميائية. تم دعم R.J. Flores من خلال منحة علوم الطاقة الأساسية DE-FG02-8361ER16402 وتم دعم A.K. Tripati من خلال منحة من مؤسسة Heising-Simons 2022-3314. تم توفير التمويل للبنية التحتية التحليلية لمختبر جيولوجيا النظائر في جامعة بويز ستايت المستخدم لقياس العناصر النادرة بواسطة LA-ICPMS في هذه الدراسة من خلال منح مؤسسة العلوم الوطنية EAR-1735889 وEAR-1920336. شكر خاص لـ M. Kohn (جامعة بويز ستايت) على ملاحظاته حول تفسير بيانات العناصر الأرضية النادرة. بالإضافة إلى ذلك، نشكر H.G. McDonald وJ. Tejada وR. McAfee وM. Curiel وS. Lee على رؤاهم وملاحظاتهم وتعليقاتهم النقدية.
مساهمات المؤلفين: قام م.د.د. بتطوير المفاهيم والأسلوب، وجمع وتحليل البيانات، وإعداد الأشكال، وصياغة وتنقيح المخطوطة؛ قام و.ب.ب. بتطوير الشيفرة وتسهيل تحليل البيانات؛ قام ب.د.م. بتطوير الشيفرة وتسهيل تحليل البيانات؛ قام د.م.ل. بأداء رسم البيانات؛ قام ر.ج.ف. بتحليل البيانات؛ قام أ.ك.ت. بتحليل البيانات؛ قام ر.أ.ع. بتطوير نهج التحليل وتفسير البيانات؛ قام د.م.س. بتطوير نهج التحليل وتفسير البيانات؛ قام م.ت.ب. بتطوير المفاهيم والأسلوب، وصياغة وتنقيح المخطوطة.
تمويل منح علوم الطاقة الأساسية DE-FG02-8361ER16402؛ منحة مؤسسة هايسينغ-سيمونز 2022-3314؛ منح المؤسسة الوطنية للعلوم EAR-1735889 وEAR-1920336.
توفر البيانات جميع البيانات التي تم جمعها وتحليلها لهذه الدراسة مدرجة في هذه المقالة المنشورة تماشيًا مع سياسة الإفصاح الخاصة بالمجلة. يشارك المؤلفون بيانات جديدة جنبًا إلى جنب مع نتائج هذا التقرير؛ تم إيداع مدخلات Niche Mapper في Zenodo (https://doi.org/10.5281/zenodo.10626431) وتم إيداع بيانات النظائر المتجمعة في EarthChem.https://doi.org/10.60520/IEDA/11312 7).

الإعلانات

المصالح المتنافسة يعلن المؤلفون عدم وجود مصالح متنافسة.
الوصول المفتوح هذه المقالة مرخصة بموجب رخصة المشاع الإبداعي النسب 4.0 الدولية، التي تسمح بالاستخدام والمشاركة والتكيف والتوزيع وإعادة الإنتاج بأي وسيلة أو صيغة، طالما أنك تعطي الائتمان المناسب للمؤلفين الأصليين والمصدر، وتوفر رابطًا لرخصة المشاع الإبداعي، وتوضح ما إذا كانت هناك تغييرات قد أُجريت. الصور أو المواد الأخرى من طرف ثالث في هذه المقالة مشمولة في رخصة المشاع الإبداعي الخاصة بالمقالة، ما لم يُشار إلى خلاف ذلك في سطر الائتمان للمواد. إذا لم تكن المادة مشمولة في رخصة المشاع الإبداعي الخاصة بالمقالة وكان استخدامك المقصود غير مسموح به بموجب اللوائح القانونية أو يتجاوز الاستخدام المسموح به، فسيتعين عليك الحصول على إذن مباشرة من صاحب حقوق الطبع والنشر. لعرض نسخة من هذه الرخصة، قم بزيارةhttp://creativecommons.org/licenses/by/4.0/.

References

Amiot R, Lécuyer C, Buffetaut E, Escarguel G, Fluteau F, Martineau F (2006) Oxygen isotopes from biogenic apatites suggest widespread endothermy in Cretaceous dinosaurs. Earth Planet Sci Lett 246(1-2):41-54. https://doi.org/10.1016/j.epsl.2006.04.018
Anderson NT, Kelson JR, Kele S, Daëron M, Bonifacie M, Horita J, et al (2021) A unified clumped isotope thermometer calibration ( ) using carbonate-based standardization. Geophy Res Lett 48(7):e2020GL092069. https://doi.org/10.1029/2020G L092069
Barbosa FHDS, Alves-Silva L, Liparini A, Porpino KDO (2023) Reviewing the body size of some extinct Brazilian Quaternary xenarthrans. J Quat Sci https://doi.org/10.1002/jqs. 3560
Bargo MS, Vizcaíno SF, Archuby FM, Blanco RE (2000) Limb bone proportions, strength and digging in some Lujanian (Late Pleis-tocene-Early Holocene) mylodontid ground sloths (Mammalia, Xenarthra). J Vertebr Paleontol 20(3):601-610. https://doi.org/10 .1671/0272-4634(2000)020[0601:LBPSAD]2.0.CO;2
Bargo MS, De Iuliis G, Vizcaíno SF (2006) Hypsodonty in Pleistocene ground sloths. Acta Paleontol Pol 51(1):53-61.
Bernasconi SM, Daëron M, Bergmann KD, Bonifacie M, Meckler AN, Affek HP et al (2021) InterCarb: A community effort to improve interlaboratory standardization of the carbonate clumped isotope thermometer using carbonate standards. Geochem Geophys Geosyst 22(5):e2020GC009588 https://doi.org/10.1029/2020GC009 588.
Billet G, Germain D, Ruf I, de Muizon C, Hautier L (2013) The inner ear of Megatherium and the evolution of the vestibular system in sloths. J Anat 223(6):557-567. https://doi.org/10.1111/joa. 12114
Blanco RE, Czerwonogora A (2003) The gait of Megatherium Cuvier 1796. Senckenb Biol 83: 61-68.
Bocherens H, Cotte M, Bonini RA, Straccia P, Scian D, Soibelzon L, Prevosti FJ (2017). Isotopic insight on paleodiet of extinct Pleistocene megafaunal xenarthrans from Argentina. Gondwana Res 48:7-14. https://doi.org/10.1016/j.gr.2017.04.003
Boeskorov GG, Mashchenko EN, Plotnikov VV, Shchelchkova MV, Protopopov AV, Solomonov NG (2016). Adaptation of the woolly mammoth Mammuthus primigenius (Blumenbach, 1799) to habitat conditions in the glacial period. Contemp Probl Ecol 9:54453. https://doi.org/10.1134/S1995425516050024
Boyles JG, Levesque DL, Nowack J, Wojciechowski MS, Stawski C, Fuller A, et al (2019) An oversimplification of physiological principles leads to flawed macroecological analyses. Ecol Evol 9(21):12020-12025. https://doi.org/10.1002/ece3.5721
Brassey CA, Gardiner JD (2015) An advanced shape-fitting algorithm applied to quadrupedal mammals: improving volumetric mass estimates. R Soc Open Sci 2(8):150302. https://doi.org/10.1098 /rsos. 150302
Campbell KL, Roberts JE, Watson LN, Stetefeld J, Sloan AM, Signore AV, Howatt JW, Tame JR, Rohland N, Shen TJ, Austin JJ (2010). Substitutions in woolly mammoth hemoglobin confer biochemical properties adaptive for cold tolerance. Nat Genet 42(6):53640. https://doi.org/10.1038/ng. 574
Cartelle C, Bohórquez GA (1986). Presença de ossículos dérmicos em Eremotherium laurillardi (Lund) Cartelle & Bohórquez 1982 (Edentanta, Megatheriidae). Iheringía Ser Geol 11:9-14.
Chavagnac V, Milton JA, Green DRH, Breuer J, Bruguier O, Jacob DE, et al. (2007). Towards the development of a fossil bone geochemical standard: An inter-laboratory study. Anal Chim Acta 599(2):177-190
Christiansen P, Fariña RA (2003) Mass estimation of two fossil ground sloths (Mammalia, Xenarthra, Mylodontidae). Senckenberg Biol 83(1):95-101
Clarke A, Rothery P (2008) Scaling of body temperature in mammals and birds. Funct Ecol 22(1):58-67. https://doi.org/10.1111/j. 136 5-2435.2007.01341.x
Cliffe RN, Haupt RJ, Avey-Arroyo JA, Wilson RP (2015) Sloths like it hot: ambient temperature modulates food intake in the brownthroated sloth (Bradypus variegatus). PeerJ 3:e875. https://doi.o rg/10.7717/peerj. 875
Cliffe RN, Scantlebury DM, Kennedy SJ, Avey-Arroyo J, Mindich D, Wilson RP (2018) The metabolic response of the Bradypus sloth to temperature. PeerJ 6:e5600. https://doi.org/10.7717/peerj. 5600
Cliffe RN, Haupt RJ, Kennedy S, Felton C, Williams HJ, Avey-Arroyo J, Wilson R (2023) The behaviour and activity budgets of two sympatric sloths; Bradypus variegatus and Choloepus hoffmanni. PeerJ 11:e15430. https://doi.org/10.7717/peerj. 15430
Collins RL (1933) Mylodont (ground sloth) dermal ossicles from Colombia, South America. J Wash Acad Sci 23(9):426-429.
Conte GL, Lopes Le, Mine AH, Trayler RB, Kim SL (2024) SPORA, a new silver phosphate precipitation protocol for oxygen isotope analysis of small, organic-rich bioappatite samples. Chem Geol 65:122000. https://doi.org/10.1016/.chemgeo.2024. 122000
Copploe JV II, Blob RW, Parrish JHA, Butcher MT (2015) In vivo strains in the femur of the nine-banded armadillo (Dasypus novemcinctus). J Morphol 276:889-899. https://doi.org/10.1002 /jmor. 20387
Croft DA, Gelfo JN, López GM (2020) Splendid innovation: the extinct South American native ungulates. Ann Rev Earth Planet Sci 48:259-290.
Dantas, MAT, Omena ÉC, da Silva JLL, Sial A (2021) Could Eremotherium laurillardi (Lund, 1842) (Megatheriidae, Xenarthra) be an omnivore species? Anu Inst Geocienc 44:1-5. https://doi.org/10. 11137/1982-3908_2021_44_36492
Dantas MAT, Cherkinsky A, Bocherens H, Drefahl M, Bernardes C, de Melo França L (2017) Isotopic paleoecology of the Pleistocene megamammals from the Brazilian Intertropical Region: Feeding ecology ( ), niche breadth and overlap. Quat Sci Rev 170:152-163. https://doi.org/10.1016/j.quascirev.2017.06.030
Deka RJ, Sarma NK (2015) Studies on feeding behaviour and daily activities of Rhinoceros unicornis in natural and captive condition of Assam. Indiana J Anim Res 49(4):542-545.
Delsuc, F, Kuch M, Gibb GC, Karpinski, E, Hackenberger D, Szpak P, Billet G (2019) Ancient mitogenomes reveal the evolutionary history and biogeography of sloths. Curr Biol 29(12):2031-2042. https://doi.org/10.1016/j.cub.2019.05.043
Eagle RA, Schauble EA, Tripati AK, Tütken T, Hulbert RC, Eiler JM (2010) Body temperatures of modern and extinct vertebrates from bond abundances in bioapatite. Proc Natl Acad Sci USA 107(23):10377-10382. https://doi.org/10.1073/pnas. 0911115107
Eagle RA, Tütken T, Martin TS, Tripati AK, Fricke HC, Connely M, et al (2011) Dinosaur body temperatures determined from isotopic ( ) ordering in fossil biominerals. Science 333(6041):443445. https://doi.org/10.1126/science. 1206196
Eagle RA, Enriquez M, Grellet-Tinner G, Pérez-Huerta A, Hu D, Tütken T, et al (2015) Isotopic ordering in eggshells reflects body temperatures and suggests differing thermophysiology in two Cretaceous dinosaurs. Nat Commun 6(1):8296. https://doi.org/1 0.1038/ncomms9296
Ehlers J, Gibbard PL, Hughes PD (2018) Quaternary glaciations and chronology. In: Menzies J, van der Meer JJM (eds) Past Glacial Environments. Elsevier, Amsterdam, pp 77-101.
Fariña RA (2002) Megatherium, the hairless: appearance of the great Quaternary sloths (Mammalia; Xenarthra). Ameghiniana 39:241-244.
Fariña RA, Blanco RE (1996) Megatherium, the stabber. Proc R Soc B263(1377):1725-1729.
Fariña RA, Vizcaíno SF, Bargo MS (1998) Body mass estimations in Lujanian (late Pleistocene-early Holocene of South America) mammal megafauna. Mastozool Neotrop 5(2):87-108.
Feldhamer GA, Drickamer LC, Vessey SH, Merritt JF, Krajewski C (2007) Mammalogy: Adaptation, Diversity, and Ecology. John Hopkins University Press, Baltimore.
Fricke HC, Rogers R (2000) Multiple taxa and multiple locality approach to providing oxygen isotope evidence for endothermic homeothermy in theropod dinosaurs. Geology 28(9):799-802.
Ghosh P, Adkins J, Affek H, Balta B, Guo W, Schauble EA, et al (2006) bonds in carbonate minerals: a new kind of paleothermometer. Geochim Cosmochim Acta 70(6):1439-1456. https://d oi.org/10.1016/j.gca.2005.11.014
González-Bernardo E, Russo LF, Valderrábano E, Fernández Á, Penteriani V (2020) Denning in brown bears. Ecol Evol 10(13):68446862. https://doi.org/10.1002/ece3.6372
González-Guarda E, Domingo L, Tornero C, Pino M, Fernández MH, Sevilla P, et al (2017). Late Pleistocene ecological, environmental and climatic reconstruction based on megafauna stable isotopes from northwestern Chilean Patagonia. Quat Sci Rev 170:188202. https://doi.org/10.1016/j.quascirev.2017.06.035
Gorvet MA, Wakeling JM, Morgan DM, Segura DH, Avey-Arroyo J, Butcher MT (2020) Keep calm and hang on: EMG activation in the forelimb musculature of three-toed sloths (Bradypus variegatus). J Exp Biol 223(14):jeb218370. https://doi.org/10.1242/j eb. 218370
Gravett N, Bhagwandin A, Sutcliffe R, Landen K, Chase MJ, Lyamin OI, et al (2017) Inactivity/sleep in two wild free-roaming African elephant matriarchs-Does large body size make elephants the shortest mammalian sleepers? PL oS ONE 12(3):e0171903. https ://doi.org/10.1371/journal.pone. 0171903
Green, JL, Kalthoff, DC (2015) Xenarthran dental microstructure and dental microwear analyses, with new data for Megatherium americanum (Megatheriidae). J Mammal 96(4):645-657. https://doi.o rg/10.1093/mammal/gyv045
Griffiths ML, Eagle RA, Kim SL, Flores R, Becker MA, Maisch IV HM, et al (2023) Thermal physiology of extinct megatooth sharks. Proc Natl Acad Sci USA 120(27):e2218153120. https://d oi.org/10.1073/pnas. 2218153120
Hausman LA (1929) The ovate bodies of the hair of Nothrotherium shastense. Am J Sci 106:331-333.
Hayssen V (2011) Choloepus hoffmanni (Pilosa: Megalonychidae). Mammal Species 43(873):37-55. https://doi.org/10.1644/873.1
Ho TY (1967) Relationship between amino acid contents of mammalian bone collagen and body temperature as a basis for estimation of body temperature of prehistoric mammals. Comp Biochem Physiol 22:113-119.
Holden PB, Edwards NR, Rangel TF, Pereira EB, Tran GT, Wilkinson RD (2019) PALEO-PGEM v1. 0: a statistical emulator of Plio-cene-Pleistocene climate. Geosci Model Dev 12(12):5137-5155. https://doi.org/10.5194/gmd-12-5137-2019
Hutchinson JR, Schwerda D, Famini DJ, Dale RH, Fischer MS, Kram R (2006) The locomotor kinematics of Asian and African elephants:
changes with speed and size. J Exp Biol 209(19):3812-3827. http s://doi.org/10.1242/jeb. 02443
Idso SB, Jackson RD (1969) Thermal radiation from the atmosphere. J Geophys Res 74(23):5397-5403.
Jochum KP, Weis U, Schwager B, Stoll B, Wilson SA, Haug GH, et al. (2016). Reference values following ISO guidelines for frequently requested rock reference materials. Geostands Geoanal Res 40(3):333-350.
John CM, Bowen D (2016) Community software for challenging isotope analysis: First applications of ‘Easotope’ to clumped isotopes. Rapid Commun Mass Spectrom 30:2285-2300. https://doi .org/10.1002/rcm. 7720
Kleiber M (1932) Body size and metabolism. Hilgardia 6(11):315-353.
Knight F, Connor C, Venkataramanan R, Asher R (2020) Body temperatures, life history, and skeletal morphology in the nine-banded armadillo (Dasypus novemcinctus). PCI Ecol. https://doi.org/10 .17863/CAM. 50971
Kohn MJ, Schoeninger MJ, Barker WW (1999) Altered states: effects of diagenesis on fossil tooth chemistry. Geochim Cosmochim Acta 63(18):2737-47. https://doi.org/10.1016/S0016-7037(99)0 0208-2
Krmpotic CM, Carlini AA, Galliari FC, Favaron P, Miglino MA, Scarano AC, Barbeito CG (2014). Ontogenetic variation in the stratum granulosum of the epidermis of Chaetophractus vellerosus (Xenarthra, Dasypodidae) in relation to the development of cornified scales. Zoology 117(6):392-397 https://doi.org/10.101 6/j.zool.2014.06.003.
Larmon JT, McDonald HG, Ambrose S, DeSantis LR, Lucero LJ (2019) A year in the life of a giant ground sloth during the Last Glacial Maximum in Belize. Sci Adv 5(2): 1-9. https://doi.org/1 0.1126/sciadv.aau1200
Lécuyer C, Balter V, Martineau F, Fourel F, Bernard A, Amiot R, Gardien V, Otero O, Legendre S, Panczer G, Simon L, Martini, R (2010) Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at . Geochim Cosmochim Acta 74(7):2072-2081. https://doi.org/10.1016/j.gca.2009.12.02 4
Levesque DL, Marshall KE (2021) Do endotherms have thermal performance curves? J Exp Biol 224(3):jeb141309. https://doi.org/1 0.1242/jeb. 141309
Lindsey EL, Reyes EXL, Matzke GE, Rice KA, McDonald HG (2020) A monodominant late-Pleistocene megafauna locality from Santa Elena, Ecuador: Insight on the biology and behavior of giant ground sloths. Paleogeogr Paleoclimatol Paleoecol 544:109599. https://doi.org/10.1016/j.palaeo.2020.109599
Lovelace DM, Hartman SA, Mathewson PD, Linzmeier BJ, Porter WP (2020) Modeling Dragons: Using linked mechanistic physiological and microclimate models to explore environmental, physiological, and morphological constraints on the early evolution of dinosaurs. PLS ONE 15(5):e0223872. https://doi.org/10.17026/d ans-238-9pjs
MacFadden BJ, DeSantis LR, Hochstein JL Kamenov GD (2010) Physical properties, geochemistry, and diagenesis of xenarthran teeth: prospects for interpreting the paleoecology of extinct species. Paleogeogr Paleoclimatol Paleoecol 291(3-4):180-189. http s://doi.org/10.1016/j.palaeo.2010.02.021
Martin FM (2018) Cueva del Milodón. The hunting grounds of the Patagonian panther. Quat Int 466:212-222. https://doi.org/10.101 6/j.quaint.2016.05.005
McBee K, Baker RJ (1982) Dasypus novemcinctus. Mammal Species 162:1-9.
McCullough EC, Porter WP (1971) Computing clear day solar radiation spectra for the terrestrial ecological environment. Ecology 52:1008-1015.
McDonald HG (2005) Paleoecology of extinct xenarthrans and the Great American Biotic Interchange. Bull Fla State Mus Nat Hist 45(4):313-333.
McDonald HG (2022) Paleoecology of the Eextinct Shasta Ground Sloth, Nothrotheriops shastensis (Xenarthra, Nothrotheridae): The Physical environment. N M Mus Nat Sci Bull 88:33-43
McDonald HG, Lundelius Jr EL (2009) The giant ground sloth Eremotherium laurillardi (Xenarthra. Megatheriidae) in Texas. Bull Mus North Ariz 65:407-421.
McKamy AJ, Young MW, Mossor AM, Young, JW, Avey-Arroyo JA, Granatosky MC, Butcher MT (2023) Pump the Brakes! The hindlimbs of three-toed sloths control suspensory locomotion. J Exp Biol 226(8):jeb245622. https://doi.org/10.1242/jeb. 245622
McMillan BR, Cottam MR, Kaufman DW (2000) Wallowing behavior of American bison (Bison bison) in tallgrass prairie: an examination of alternate explanations. Am Mid Nat 144(1):159-167.
McNab, BK (1985) Energetics, population biology, and distribution of xenarthrans, living and extinct. In: Montgomery GG (ed) The Evolution and Ecology of Armadillos, Sloths, and Vermilinguas. Smithsonian Institution Press, Washington, D.C. pp 219-232.
Mine AH, Waldeck A, Olack G, Hoerner ME, Alex S, Colman AS (2017) Microprecipitation and analysis of phosphate for paleoclimate and biogeochemistry research. Chem Geol 460:1-4 https://doi.org/10.1016/j.chemgeo.2017.03.032.
Mitchell D, Snelling EP, Hetem RS, Maloney SK, Strauss WM, Fuller A (2018) Revisiting concepts of thermal physiology: predicting responses of mammals to climate change. J Anim Ecol 87(4):956-973. https://doi.org/10.1111/1365-2656.12818
Myhrvold CL, Stone HA, Bou-Zeid E (2012) What is the use of elephant hair? PL oS ONE 7:e47018. https://doi.org/10.1371/journa 1.pone. 0047018
Nance WB, Taylor SR (1976) Rare earth element patterns and crustal evolution-I. Australian post-Archean sedimentary rocks. Geochim Cosmochim Acta 40:1539-1551.
Naples VL (1990) Morphological changes in the facial region and a model of dental growth and wear pattern development in Nothrotheriops shastensis. J Vert Paleontol 10(3):372-389.
Nehme C, Todisco D, Breitenbach SF, Couchoud I, Marchegiano M, Peral M, et al (2022) Holocene hydroclimate variability along the Southern Patagonian margin (Chile) reconstructed from Cueva Chica speleothems. Glob Planet Change 222:104050. https://do i.org/10.1016/j.gloplacha.2023.104050
Nelson A, Engelman RK, Croft DA (2023) How to weigh a fossil mammal? South American notoungulates as a case study for estimating body mass in extinct clades. J Mammal Evol 30(3):773-809 h ttps://doi.org/10.1007/s10914-023-09669-1
Pant SRA, Goswami A, Finarelli JA (2014) Complex body size trends in the evolution of sloths (Xenarthra: Pilosa). BMC Evol Biol 14:184. https://doi.org/10.1186/s12862-014-0184-1
Pauli JN, Peery MZ, Fountain ED, Karasov WH (2016) Arboreal folivores limit their energetic output, all the way to slothfulness. Am Nat 188(2):196-204.
Poinar HN, Hofreiter M, Spaulding WG, Martin PS, Stankiewicz BA, Bland H, Evershed RP, Possnert G, Paabo S (1998) Molecular coproscopy: dung and diet of the extinct ground sloth Nothrotheriops shastensis. Science 128(5375):402-406.
Presslee S, Slater GJ, Pujos F, Forasiepi AM, Fischer R, Molloy K, et al. (2019). Palaeoproteomics resolves sloth relationships. Nat Ecol Evol 3(7):1121-1130. https://doi.org/10.1038/s41559-019-0 909-z
Prevosti FJ, Martin FM (2013) Paleoecology of the mammalian predator guild of Southern Patagonia during the latest Pleistocene: ecomorphology, stable isotopes, and taphonomy. Quat Int 305:74-84. https://doi.org/10.1016/j.quaint.2012.12.039
Pujos F, Salas R (2004) A systematic reassessment and paleogeographic review of fossil Xenarthra from Peru. Bull Inst Fr Et Andines 33:331-377.
Ridewood WG (1901). On the structure of the hairs of Mylodon listai and other South American Edentata. J Cell Sci 2(175):393-411. h ttps://doi.org/10.1242/jcs.s2-44.175.393
Santos LE, Ajala-Batista L, Carlini AA, de Araujo Monteiro-Filho EL (2024) Mylodon darwinii (Owen, 1840): hair morphology of an extinct sloth. Zoomorphology 2024:1-9.
Seitz VP, Puig S (2018) Aboveground activity, reproduction, body temperature and weight of armadillos (Xenarthra, Chlamyphoridae) according to atmospheric conditions in the central Monte (Argentina). Mammal Biol 88:43-51.
Seymour R S, Hu Q, Snelling E P, White C R (2019). Interspecific scaling of blood flow rates and arterial sizes in mammals. J Exp Biol 222(7):jeb199554. https://doi.org/10.1242/Jen. 199554
Stears K, McCauley DJ, Finlay JC, Mpemba J, Warrington IT, Mutayoba BM, et al (2018) Effects of the hippopotamus on the chemistry and ecology of a changing watershed. Proc Nat Acad Sci USA 115(22):E5028-E5037. https://doi.org/10.1073/pnas. 1800407115
Stelmock J, Dean FC (1986) Brown bear activity and habitat use, Denali National Park-1980. Int Conf Bear Res Manag 6:155-167.
Suárez MB, Passey BH (2014) Assessment of the clumped isotope composition of fossil bone carbonate as a recorder of subsurface temperatures. Geochim Cosmochim Acta 140:142-159. https:/ dol.org/10.1016|)- gca.2014.05.026
Sunquist ME, Montgomery GG (1973) Activity patterns and rates of movement of two-toed and three-toed sloths (Choloepus hoffmanni and Bradypus infuscatus). J Mammal 54(4):946-954.
Superina M, Abba AM (2014) Zaedyus pichiy (Cingulata: Dasypodidae). Mammal Species 46(905):1-10. https://doi.org/10.1644/90 5.1
Superina M, Boily P (2007) Hibernation and daily torpor in an armadillo, the pichi (Zaedyus pichiy). Comp Biochem Physiol A Mol Int Physiol 148(4):893-898. https://doi.org/10.1016/j.cbpa. 2007 . 09.005
Tagliavento M, Davies AJ, Bernecker M, Staudigel PT, Dawson RR, Dietzel M, et al (2023) Evidence for heterothermic endothermy and reptile-like eggshell mineralization in Troodon, a non-avian maniraptoran theropod. Proc Natl Acad Sci 120(15):e2213987120. https://doi.org/10.1073/pnas. 2213987120
Taube E, Vié JC, Fournier P, Genty C, Duplantier JM (1999) Distribution of two sympatric species of sloths (Choloepus didactylus and Bradypus tridactylus) along the Sinnamary River, French Guiana 1. Biotropica 31(4):686-691.
Tejada JV, Flynn JJ, MacPhee R, O’Connell TC, Cerling TE, Bermudez L, et al (2021) Isotope data from amino acids indicate Darwin’s ground sloth was not an herbivore. Sci Reports 11(1):1-10. https://doi.org/10.1038/s41598-021-97996-9
Tejada JV, Antoine PO, Münch P, Billet G, Hautier L, Delsuc F, Condamine FL (2024) Bayesian total-evidence dating revisits sloth phylogeny and biogeography: a cautionary tale on morphological clock analyses. Syst Biol 73(1):125-139. https://doi.org/10.1093 /sysbio/syad069
Toledo N, Boscaini A, Pérez LM (2021) The dermal armor of mylodontid sloths (Mammalia, Xenarthra) from Cueva del Milodon (Ultima Esperanza, Chile). J Morphol 282(4):612-627. https://do i.org/10.1002/jmor. 21333
Tomassini RL, Montalvo C I, Garrone MC, Domingo L, Ferigolo J, Cruz LE, et al (2020) Gregariousness in the giant sloth Lestodon (Xenarthra): multi-proxy approach of a bonebed from the Last Maximum Glacial of Argentine Pampas. Sci Reports 10(1):19055. https://doi.org/10.1038/s41598-020-67863-0
Tregear RT (1965) Hair density, wind speed, and heat loss in mammals. J Appl Physiol 20(4):796-801.
Varela L, Tambusso PS, Patiño SJ, Di Giacomo M, & Fariña RA (2018) Potential distribution of fossil xenarthrans in South America during the late Pleistocene: co-occurrence and provincialism. J Mamm Evol 25:539-550. https://doi.org/10.1007/s10914-017-9 406-9
Varela L, Clavijo L, Tambusso PS, Fariña RA (2023). A window into a late Pleistocene megafauna community: Stable isotopes show niche partitioning among herbivorous taxa at the Arroyo del Vizcaíno site (Uruguay). Quat Sci Rev 317:108286. https://doi.org/1 0.1016/j.quascirev.2023.108286
Varela L, Tambusso S, Fariña R. (2024). Femora nutrient foramina and aerobic capacity in giant extinct xenarthrans. PeerJ 12:e17815. ht tps://doi.org/10.7717/peerj. 17815
Wang Y, Porter WP, Mathewson PD, Miller PA, Graham RW, Williams JW (2018) Mechanistic modeling of environmental drivers of
woolly mammoth carrying capacity declines on St. Paul Island. Ecology 99(12):2721-2730. https://doi.org/10.1002/ecy. 2524
Wang Y, Widga C, Graham RW, McGuire JL, Porter WP, D. Wårlind, D, Williams JW (2020) Caught in a bottleneck: Habitat loss for woolly mammoths in central North America and the ice-free corridor during the last deglaciation. Glob Ecol Biogeogr 30(2):527542. https://doi.org/10.1111/geb. 13238
White CR, Seymour RS (2004) Does basal metabolic rate contain a useful signal? Mammalian BMR allometry and correlations with a selection of physiological, ecological, and life-history variables. Physiol Biochem Zool 77(6):929-941.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Journal: Journal of Mammalian Evolution, Volume: 32, Issue: 1
DOI: https://doi.org/10.1007/s10914-024-09743-2
PMID: https://pubmed.ncbi.nlm.nih.gov/39822851
Publication Date: 2025-01-14

Metabolic skinflint or spendthrift? Insights into ground sloth integument and thermophysiology revealed by biophysical modeling and clumped isotope paleothermometry

Michael D. Deak • Warren P. Porter • Paul D. Mathewson • David M. Lovelace • Randon J. Flores . Aradhna K. Tripati ⋅ Robert A. Eagle ⋅ Darin M. Schwartz ⋅ Michael T. Butcher

Accepted: 17 December 2024 / Published online: 14 January 2025
© The Author(s) 2024, corrected publication 2025

Abstract

Remains of megatheres have been known since the 18th -century and were among the first megafaunal vertebrates to be studied. While several examples of preserved integument show a thick coverage of fur for smaller ground sloths living in cold climates such as Mylodon and Nothrotheriops, comparatively very little is known about megathere skin. Assuming a typical placental mammal metabolism, it was previously hypothesized that megatheres would have had little-to-no fur as they achieved giant body sizes. Here the “hairless model of integument” is tested using geochemical analyses to estimate body temperature to generate novel models of ground sloth metabolism, fur coverage, and paleoclimate with Niche Mapper software. The simulations assuming metabolic activity akin to those of modern xenarthrans suggest that sparse fur coverage would have resulted in cold stress across most latitudinal ranges inhabited by extinct ground sloths. Specifically, Eremotherium predominantly required dense 10 mm fur with implications for seasonal changes of coat depth in northernmost latitudes and sparse fur in the tropics; Megatherium required dense 30 mm fur year-round in its exclusive range of cooler, drier climates; Mylodon and Nothrotheriops required dense fur to avoid thermal stress, matching the integument remains of both genera, and further implying the use of behavioral thermoregulation. Moreover, clumped isotope paleothermometry data from the preserved teeth of four genera of ground sloth yielded reconstructed body temperatures lower than those previously reported for large terrestrial mammals ( ). This combination of low metabolisms and thick fur allowed ground sloths to inhabit various environments.

Keywords Body temperature ⋅ Fur ⋅ Geochemistry ⋅ Metabolism ⋅ Paleoclimate ⋅ Paleontology ⋅ Skin ⋅ Xenarthra
△ Michael T. Butcher
mtbutcher@ysu.edu
1 Department of Chemical and Biological Sciences, Youngstown State University, Youngstown, OH, USA
2 Department of Integrative Biology, University of Wisconsin Madison, Madison, WI, USA
3 Department of Geosciences, University of Wisconsin Madison, Madison, WI, USA
4 Department of Earth and Space Sciences, University of California – Los Angeles, Los Angeles, CA, USA
5 Department of Atmospheric and Oceanic Sciences, Institute of the Environment and Sustainability, Center for Diverse Leadership in Science, University of California – Los Angeles, Los Angeles, CA, USA
Department of Geosciences, Boise State University, Boisie, ID, USA

Introduction

Xenarthrans have a rich and diverse history preserved in the fossil record dating back 60 million years ago (MYA) to the Eocene Epoch. The Xenarthra consists of the clades Cingulata, comprising living armadillos and the extinct pampatheres and glyptodonts, and Pilosa, which contains anteaters and sloths. While the genera Choloepus (two-toed sloths) and Bradypus (three-toed sloths) are the only extant linages of (tree) sloths, numerous extinct genera of sloths (i.e., arboreal, terrestrial, and aquatic) once occupied various niches during the Oligocene and persisted until the Late Pleistocene. Ground sloths thrived throughout South America and dispersed north to eventually colonize North America during the late Cenozoic Great American Biotic Interchange (GABI) (Feldhamer et al. 2007). Several extinct ground
sloth taxa such as Eremotherium and Megatherium became the largest xenarthrans ever in existence, reaching sizes comparable to or exceeding that of modern African savanna elephants.
To date, preserved integument that is definitively attributed to the largest ground sloths Eremotherium and Megatherium has yet to be found. The only specimen of Megatherium having a sample of skin was lost before it could be inventoried (Pujos and Salas 2004). Nonetheless, both taxa have been historically reconstructed as having a thick coat of long, shaggy fur similar to that of smaller ground sloths (e.g., Mylodon and Nothrotheriops) and extant tree sloths. This appearance, however, was challenged on the basis of thermal energetics suggesting a largely hairless integument in the largest ground sloths (McNab 1985; Fariña 2002). These analyses utilized the following thermal conductivity equation: , where is the heat flow density of the skin, is the thermal conductivity of the skin, s is the skin thickness, and core body temperature and ambient temperature is represented by and , respectively. Fariña (2002) hypothesized that given their massive sizes, the large-bodied ground sloths would not have been thermally neutral if they were covered in thick fur. However, if the integument of Megatherium closely mirrored that of the largest extant terrestrial mammals (i.e., elephants, hippopotamuses, and rhinoceroses), then it would be in a thermoneutral zone even if the surrounding temperatures were nearly on average. A thermoneutral zone is defined as a range of ambient temperatures where an endotherm can remain in heat balance without metabolic heat production or requiring evaporative heat loss. While it is tempting to make comparisons to extensively furry megafauna such as the woolly mammoth, it should be noted that these comparisons are not applicable to megatheres. Woolly mammoths had a suite of adaptations for living in the frigid Siberia tundra such as reduced ear and tail sizes and an insulating layer of fat (Boeskorov et al. 2016), as well as hemoglobin that prevented their blood from freezing (Campbell et al. 2010).
The “hairless model of integument” was again recently evoked by Lindsey et al. (2020) regarding a mass assemblage of Eremotherium remains recovered from Tanque Loma in Ecuador. The cause of death of these animals was believed to be disease (or drought) from fecal matter contamination of the marsh. Lindsey et al. (2020) postulated that wallowing behavior would be advantageous if E. laurillardi had a hairless integument, as it would offer protection against the sun and tropical insects, and thus would be similar to a protective strategy used by the largest-bodied extant terrestrial mammals.
Despite evidence in favor of a hairless (or nearly hairless) model of integument in the largest ground sloths that ever
lived, there were critical assumptions applied to the previous methodologies that may not be tenable. First, the thermal conductivity equation used by Fariña (2002) only accounts for core body temperatures and the effects of ambient temperatures alone on thermoneutrality without accounting for other environmental factors. While this calculation is useful in determining heat conduction flux, other environmental factors that could impact thermoneutrality such as relative humidity and wind speed, in addition to variations in ambient temperature across the broad geographic distribution of the ground sloth genera should have been considered. Second, the previous model (Fariña 2002) did not accurately account for the basal metabolism of extant xenarthrans for calculations of heat flux, or the thick skin of living sloths that would affect heat loss. Instead, the basal metabolic rate of a naked human was isometrically scaled up to the size of Megatherium assuming similar core body temperature and skin thickness to that of hominids (Fariña 2002). Third, determinations of the basal metabolic rate of Megatherium americanum were made by halving the estimated metabolic rate, thus leading to the conclusions that it could still be furless and that temperatures as low as would be within its thermoneutral zone. Fourth, thermoneutral zone estimates derived from the previous model for smaller ground sloth taxa are inconsistent with known paleoclimate data and preserved integument. Temperatures of and C , respectively, have also been calculated as being within the thermoneutral zone for a completely furless Mylodon darwinii versus a conspecific with a fur coat 4 cm long, but with each having an equivalent mass of . (Fariña 2002).
The fossil evidence known for Mylodon, however, indicates that it was extensively covered in fur, which was as long as 6.5-10 cm in some regions of its body (Collins 1933; Santos et al. 2024). Stable isotopes found in contemporary samples of Notiomastodon enamel dating to years ago have also produced mean annual temperature estimates ranging between and (González-Guarda et al. 2017). Moreover, the relationship between body mass and body temperature in mammals is complex and purported to be mediated primarily through ecology (Clarke and Rothery 2008). Recent data on large extant mammal energetics indicate that the thermoneutral zone only accounts for favorable ambient temperature ranges where little energy needs to be expended beyond the basal metabolic rate, and not the full extent of the normal functions associated with daily activity such as foraging. Thus, prescriptive zones, tolerance zones, and survival zones have been adopted by other authors (Mitchell et al. 2018; Boyles et al. 2019; Levesque and Marshall 2021), which encompass temperatures where an endotherm can function normally, temperatures where
thermal stress occurs, and temperatures where death by thermal stress occurs, respectively.
Recently, Varela et al. (2024) attempted to estimate the metabolic rates of several extinct xenarthrans by measuring the diameter of the femoral nutrient foramina. Estimates of blood flow obtained from an equation developed by Seymour et al. (2019) resulted in metabolic rates in fossil xenarthrans comparable to non-xenarthran placental mammals of similar size. While these findings lend support to speculation that ground sloths possessed a typical placental metabolism due to increased skeletal muscle mass compared to living tree sloths (Fariña 2002), and faced environmental conditions below proposed thermal neutrality (Tejada et al. 2021), other data suggests a more xenarthran-like metabolism. For example, estimates of the body temperature in Nothrotheriops (based on amino acids preserved in bone collagen) showed that it had a range of body temperature values akin to extant xenarthrans (Ho 1967). Similar estimates of the metabolic rates and thermal conductivity for both Mylodon and Nothrotheriops suggest that each had low basal metabolic rate and relatively poor insulation compared with other mammals of similar size (McNab 1985). Ground sloths plausibly may have had lower metabolic rates than other mammals of their size, although this tentative conclusion does not imply that they were as slow-moving as living tree sloths. Using step length from a fossilized trackway at Pehuén Có, Argentina attributed to Megatherium, Blanco and Czerwonogora (2003) predicted that its average walking speed was , which is 10 times faster than extant Bradypus suspensory walking below-branch (Gorvet et al. 2020; McKamy et al. 2023). A maximum speed of nearly , which is still considerably slower than the top velocity of measured for elephants (Hutchinson et al. 2006), also has been estimated for Megatherium (Billet et al. 2013). Despite ground sloths, elephants, and tree sloths having different modes of locomotion, their reported speeds emphasize the general principle that animals with low metabolic rates tend to move more slowly than those with higher metabolic rates.
The objective of this research is to further test the hairless model of integument for the largest ground sloths, Eremotherium and Megatherium, and to place novel constraints on their thermoregulation. To rigorously test the model through quantitative means, numerous climatic and physiological factors were input, for example core body temperature, metabolic rate, and estimates of paleoclimate. Simultaneously, integrating these variables will provide more insight into the reconstruction of large ground sloth integument until skin and/or fur samples are discovered. In addition to physiological and ecological constraints on integument, a better understanding of the evolution of core body temperature and/or metabolic rate across various extinct and extant
phylogenetic branches of the Folivora is an expected outcome of this study.

Materials and methods

Fossil specimens and sample preparation

Tooth samples from Eremotherium ( ), Megatherium , Promegatherium ( ), and Nothrotheriops ( ) were used in a series of geochemical analyses (Table 1). Cortical bone samples from either the same specimen or from another in the same general locality were also used to determine the degree of diagenesis via an analysis of rare earth element (REE) uptake in both cortical bone and tooth samples. The specimens obtained were from the collections of the Field Museum of Natural History in Chicago, IL (FMNH) and the University of Florida in Gainesville, FL (UF). No permits were required for the present study, which complied with all relevant legal regulations. While no samples of Mylodon were obtained for this study, body temperature estimates for this genus were extrapolated from values obtained from Promegatherium given that both genera were similar in body size. This comparison based on body size rather than phylogeny is justifiable given that Bradypus and Choloepus are similar in size and metabolic capabilities despite being distantly related ( McNab 1985 ; Cliffe et al. 2023). The sampling strategy accounted for a wide collective geographic distribution of these genera from the southernmost extreme of their range in Argentina to the northernmost extreme of their range in Florida.
Given that the outer dentin of fossil xenarthran teeth are susceptible to large quantities of rare earth element uptake (MacFadden et al. 2010), inner orthodentine was preferentially sampled, as a recent study (Larmon et al. 2019) using cathodoluminescence revealed that this dental material is the most resistant tissue to diagenesis and exhibits the least amount of mineral impregnation compared to other dentine layers. Approximately of powderized internal orthodentine (and cortical bone) were harvested with a power hand drill fit with a 0.7 mm diamond-tipped drill bit. Drilled samples of inner orthodentine were taken from several random sites along the entire length of the tooth starting at the apex to the root (Fig. 1a). To remove any organic contaminants, the samples were washed with of for 3 h in 15 ml conicals. The tubes remained open so that could be released during purification. The samples were then pelleted in a centrifuge at 10 k rpm (centrifugation radius ; centrifugal force ) for 10 min . Specifically, all samples were centrifuged at a constant temperature of , followed by 3 rinses in room temperature water. Any remaining liquid was decanted before
Table 1 Body temperature estimates from clumped isotopic analyses and diagenesis ratios from rare earth element (REE) analyses. Abbreviations: SD, standard deviation; SE, standard error
Genus/ Sample ID Locality Depositional environment REE (%o, (%o, V-SMOW) ± SD
Eremotherium UF 95869 Inglis 1 A, Citrus County, Florida, USA Sinkhole 5 0.42
Nothrotheriops UF 87131 Leisey 1 A, Hillsborough County, Florida, USA Estuary 4 0.86
Megatherium FMNH P13725 Tarija Valley, Argentina Fluvial 4 0.12
Eremotherium UF 312730 Halie 7G, Alachua County, Florida, USA Lacustrine/ karst 4 0.02
Promegatherium FMNH 14404 Puerta de Corral Quemado, Argentina Terrestrial indet. 3 36.6
Number of replica stable isotope measurements from the same sample
Rare earth element ratio between dentin and cortical bone (see Methods)
temperature (Anderson et al. 2021)
(where “w” signifies mineral formation waters; Lécuyer et al. 2010)
drying overnight in an oven at . Extra drying time in the oven was necessary for some of the samples to ensure complete desiccation of contaminants so that the samples were dried to a concentrated mass.

Geochemical alteration analysis

To assess the amount of diagenesis present in the samples, REE analyses were performed on subsamples ( ) of prepared inner orthodentine and cortical bone (MacFadden et al. 2010). For analysis, samples were weighed by microbalance and digested in a mixture of 0.1 mL 29 M HF and 0.3 on a hotplate at . Samples were dried down and redissolved in 0.5 ml of three times prior to diluting gravimetrically to for analysis. Samples were prepared along-side a suite of geologic and bone standard powders and total procedural blank following the same dissolution procedure, for use as calibration and external reference materials. Analyses were performed on a Thermo Scientific, iCAP-RQ, single quadrupole mass spectrometer using a PFA micromist nebulizer; peltier-cooled quartz, cyclonic spray chamber, 2.5 mm quartz injector and quartz torch; Ni cones, high matrix ( 3.5 mm ) skimmer insert. Samples were uptaken using an ESI SC 4DX FAST autosampler using the FAST 2.0 mL method, with a stabilization time of 30 s and dual rinse ( 4 s each) in trace HF following each analysis. Samples were mixed in realtime with a 2 ppb In solution using the FAST autosampler, which was used as an internal standard for the experiment. Reported analyses are the mean of 5 measurements, which
consist of 25 sweeps on , , , for 10 ms each. Raw counts were normalized to a uniform intensity over the run by normalizing against In. Concentrations were converted from normalized counts using Thermo Scientific Qtegra software (eQuant) by linear regression of the calibration standards (BHVO-2, BCR-2, BIR-1, AGV-2; Jochum et al. 2016), anchored to the total procedural blank. Accuracy and precision were monitored using whole rock powder W2 (Jochum et al. 2016), NIST1400 bone ash, NIST1486 bone meal and NOCS fossil bone material (Chavagnac et al. 2007) and multiple dilutions of a mixed 100 ppm REE standard solution (BDH82025-926) diluted to similar concentrations as the unknowns.
All measurements were normalized against PostArchean Australian Shale (PAAS; Nance and Taylor 1976) to compensate for the effects of the Oddo-Harkins even-odd abundance effect. While MacFadden et al. (2010) reported a significant increase in REE uptake in fossils compared to modern bone and dental tissues, an REE index was established to assess the degree of diagenesis in fossil xenarthran teeth compared to enamel samples from similar localities. The sum of all normalized REE ( ) present in orthodentine were divided by the sum of all normalized REE found in cortical bone and expressed as a ratio ( sample/ cortical bone). Ratios indicate significant geochemical alteration signifying that isotopic geochemical data should be treated with caution, whereas sample ratios between 0.5 and 1.0 indicate mild to moderately altered
Fig. 1 a. The tooth of Megatherium FMNH P13725 with sampled regions indicated by yellow arrows; scale bar equals 30 mm . b. Percentage of phosphate in the sampled ground sloth teeth compared to the ranges reported for extant placental mammal bioapatite; blue and gray shaded areas indicate 95th and 99th percentile, respectively;
mammal bone and enamel data from Eagle et al. (2011). c. Scatter plot of the relationships between , and estimated temperatures of the sampled teeth; symbols with light outlines represent individual samples, and symbols with bold outlines and error bars represent the mean and standard deviations of the samples
samples, and ratios (akin to most samples of fossil enamel from non-xenarthrans) indicates a sample with little to no geochemical alteration (MacFadden et al. 2010).
Samples were prepared for phosphate analysis following the rapid silver phosphate precipitation protocol outlined in Mine et al. (2017) and detailed in Conte et al. (2024). In brief, approximately 1.0 mg of sample material was weighed out in 2 mL microcentrifuge tubes and dissolved overnight in 2.0 M nitric acid ( ). Calcium fluoride was then precipitated with the addition of 2.9 M hydrofluoric acid (HF) and 2.0 M sodium hydroxide . Samples were centrifuged to pellet and the phosphate ion containing supernatant was transferred to a new clean microcentrifuge tube. The pellet was rinsed with 0.1 M sodium fluoride to recover residual phosphate, and the resulting supernatant aliquot was combined
with the first aliquot. Sample phosphate was precipitated as silver phosphate with a silver amine solution containing 1.09 M ammonium hydroxide and 0.37 M silver nitrate at a pH range of 5.5 to 7.5. Adjustments to pH were made with the addition of small aliquots of or 2 M NaOH and samples were allowed to react for 10 min . Samples were centrifuged to pellet crystals and following the removal of supernatant fluid were rinsed five times with distilled water and dried overnight in a oven. Dried samples were weighed out to in packed silver capsules and analyzed in triplicate using a ThermoFisher Thermal Conversion Elemental Analyzer (TC/EA) paired with a Delta V plus Continuous Flow Isotope Ratio Mass Spectrometer in the Stable Isotope Ecosystem Laboratory of UC Merced (SIELO). Reference material NIST SRM 120 C was precipitated and analyzed
alongside sample unknowns as a preparation standard and yielded a value of . Reference materials USGS and USGS 81 ( ) were used for linearity and drift corrections and standardization.

Bioapatite stable isotope geochemistry and body temperature estimates

Carbonate clumped isotope thermometry is an approach to determine mineral formation temperatures based on measurements (which are reported as values) of the abundance of bonds in carbonate moieties (Ghosh et al. 2006). The values found in tooth bioapatite reflect the body temperature of the organism that formed them (Eagle et al. 2010). Furthermore, measurements of vertebrate biominerals are able to distinguish between endothermic and ectothermic physiology (Eagle et al. 2011, 2015; Griffiths et al. 2023), as well as between endotherms with different metabolic rates. For example, these techniques have played a crucial role in resolving body temperature differences between regionally endothermic sharks, large terrestrial and marine mammals, and between mammals and birds. Measurements from fossils also yielded plausible body temperatures for extinct vertebrates across of a range of epochs (Eagle et al. 2010, 2011, 2015; Griffiths et al. 2023).
All samples were analyzed via stable isotope ratio mass spectrometry. Specifically, samples of ground sloth inner orthodentine were digested in a common acid bath containing phosphoric acid reacted at followed by multiple steps of cryogenic purification and passage through either a gas chromatograph column (MAT 253 mass spectrometer, Thermo Fisher Scientific, USA) or adsorption trap (Nu Perspective mass spectrometer, Nu Instruments-AMETEK, UK) packed with Porapak Type-QTM 50/80. Following sample purification, gas was then introduced into the isotope ratio mass spectrometer (IRMS) systems via an automated changeover block, thus allowing for continuous alternating measurements of sample and reference gas on detectors configured to measure isotopologues of masses 44 to 49 . Measurements performed on the Thermo MAT 253 mass spectrometer ( of sample material) consisted of nine blocks of 10 cycles totaling 720 s of integration time at a signal of over the course of a block with pressure adjustments between each block to return to a beam of 16 V at the start of each new block. Measurements taken with the Nu Perspective mass spectrometer ( of sample material) consisted of three blocks of 20 cycles totaling in of integration at a signal which ranged from 80 to 30 nA on mass 44 over the course of each acquisition. Carbonate standards ETH-1 and ETH-2 were used to perform a nonlinearity correction on sample
unknowns. These nonlinearity corrected values were then projected into an I-CDES reference frame (Bernasconi et al. 2021) using ETH-1-2, and ETH-3 along with three in-house standards. ETH-4 was not included in the corrections and was instead used as a check for data quality. Corrections were applied to tooth samples and standard data over a moving average of 10 standards on either side of a given sample and were calculated using Easotope software (John and Bowen 2016). Reconstructed temperatures were derived from a -temperature calibration method (Anderson et al. 2021) that has recently demonstrated that using the latest I-CDES reference frame a bioapatite calibration was statistically indistinguishable from the calibration used in these analyses (Anderson et al. 2021; Griffiths et al. 2023). Equations used to estimate the formation of waters during mineralization are published elsewhere (Lécuyer et al. 2010).
An accepted limitation in the application of clumped isotope thermometry for assessing vertebrate physiology is the preservation of original isotope signatures and the possibility of post-depositional diagenesis, and how these factor influence the accuracy of the estimates. While precisely constraining diagenesis can be complex, the following two approaches were employed: (1) determining the REE concentration of fossil ground sloth teeth and (2) comparing measurements of of both the carbonate and phosphate moieties of fossil bioapatite. As the carbonate moiety may be more susceptible to chemical alteration than the phosphate moiety deviation, the offset between the range of carbonate and phosphate seen in extant organisms has been used to indicate preferential diagenesis of the carbonate moiety (Eagle et al. 2011). Rare earth elements can be taken up by bioapatite during the fossilization process and very high REE content may be indicative of recrystallization of the fossil apatite (McFadden et al. 2010).

Modeling of animals and paleoclimate

Body temperature estimates obtained from the literature for extant xenarthrans (McNab 1985; Cliffe et al. 2018; Seitz and Puig 2018) were evaluated via linear regression to determine the relationship between core body temperature and body size for all initial simulations. Temperature (in ) was log-transformed and plotted against the log of mass (in kg) to establish an expected relationship. These data showed no relationship among extant xenarthrans (Fig. 2); however, extant xenarthran body masses range between 0.2 and 45.2 kg , thus none of the taxa used for this scaling relationship approximate the body size of extinct ground sloths. Nonetheless, modeling simulations were meant to represent an average adult-sized animal at any given point in geologic time, but are not meant to be based on a specific specimen such as those used in the clumped isotope analyses. Body
Fig. 2 Regression of extant xenarthran body mass against body temperature ; ). The data include three sloth species ( ), eight armadillo species ( ), and three anteater species ( ) and were obtained from the literature (McNab 1985; Cliffe et al. 2018; Seitz and Puig 2018)
temperature estimates, derived from clumped isotope paleothermometer analyses among various other environmental and physiological variables were input into the broadly tested biophysiological modeling software, Niche Mapper. The microclimate sub-model of Niche Mapper converts inputs of macroclimate values into the warmest (minimum shade) and coolest (maximum shade) microclimates simultaneously from 2 m above ground to 2 m below ground. The animal sub-model then uses the hourly results in order to find a microenvironment suitable for optimal thermal comfort, while allowing the user to define various morphological, behavioral, and physiological inputs, including animal body geometry, user defined metabolic rate, denning and burrowing behaviors, etc. Both sub-models solve for heat and mass balances of the local substrates and the animal, respectively, in addition to finding the best available microclimate for each hour of a simulated day to keep the simulated animal within the closest possible range of the thermoneutral zone.
Recent work with Niche Mapper has been used to predict both environmental and physiological constraints on various extinct taxa, including mammoths (Wang et al. 2018, 2020) and Triassic dinosaurs (Lovelace et al. 2020). Basal metabolic rate was estimated using a published metabolic scaling relationship derived from 15 extant xenarthran species ( ; White and Seymour 2004), where M is body mass in grams, and the resulting basal metabolic rate in was converted into the unit of watts. An activity multiplier of 2 was used to estimate values of field metabolic rate. The scaling relationship ( ; Fariña 2002), where M is body mass in kg , was originally derived from Kleiber (1932) and it was used to compare metabolic rates derived from a mammal of similar body size with the typical placental metabolism.
Three-dimensional computer models for both Eremotherium and Mylodon were fitted with non-uniform rational basis spline (NURBS) curves and used to estimate body volumes, surface areas, and rate of heat dissipation, and these data further served as proxies for the same computed parameters in metabolic models for Megatherium and Nothrotheriops, respectively. The models were constructed in ZBrush by paleoartists, and we’re used in this study with permission. While distantly related phylogenetically, the 3D model of Mylodon was specifically used as a proxy for Nothrotheriops given that both taxa were more similar in body size than either were to the megatheres. Mass estimates for Eremotherium ( ; Barbosa et al. 2023), Megatherium ( ; Brassey and Gardiner 2015), Mylodon ( ; Christiansen and Fariña 2003), and Nothrotheriops ( 463 kg ; McDonald 2005) were input into these simulations. While higher mass estimates have been reported in the literature for both Eremotherium and Megatherium ( ; Dantas et al. ; Fariña et al. 1998), these were treated as high-end estimates herein, and more conservative body masses were used instead. Models of full body fur coverage were tested at coat thicknesses of 10, 30, and 50 mm for all four taxa. Fur densities that were tested specifically consisted of “dense” fur of 2,000 hairs (Lovelace et al. 2020) and 1,300 hairs (e.g., horse flank fur density; Tregear 1965), in addition to “sparse” fur of 8.5 hairs (e.g., pig flank fur density; Tregear 1965) and 0.07 hairs/ (maximum elephant fur density; Myhrvold et al. 2012).
Metabolic chamber simulations for each taxon were used to estimate approximate thermoneutral zone, as well as lower critical temperature and upper critical temperatures. The final core body temperature ranges used were derived from those determined by the clumped isotope paleothermometry analysis, as well as published core body temperatures from extant xenarthrans (McNab 1985; Cliffe
et al. 2018; Seitz and Puig 2018). Output thermoneutral zones were determined within a target range of basal metabolic rate with the animal modeled as standing still with no exposure to sunlight. Constant wind speeds of 0.1 were also assumed for the metabolic chamber simulations. The upper critical temperature was determined when the metabolic rate fell below 5% of the target basal metabolic rate range, while the lower critical temperature was determined when the metabolic rate was above of the target basal metabolic rate range as per previous analyses (Lovelace et al. 2020). This process was repeated for each proposed integument model.
To account for environmental conditions of the fossil localities (beyond ambient temperatures) during the Pleistocene, the Paleobiology Database (https://paleobiodb.org) was used to determine the northernmost and southernmost ranges, in addition to a mid-latitude species range to test for a variety of climatic conditions for each taxon. Average yearly temperatures corresponding to the approximate stratigraphical ages of the simulated localities were obtained from previously published PALEO-PGEM v1.0 simulations of various global mean annual temperatures over a span of 5 million years with a total of 5,001 measurements in 1,000year increments (Holden et al. 2019). The global climate extractor in Niche Mapper was then used to determine the factored difference from modern climate values of monthly maximum and minimum temperature for the estimated stratigraphical age of the locality. Relative humidity, cloud cover, and wind speed were approximated from modern values measured by global climate extractor database in Niche Mapper for similar regions. The microclimate model in Niche Mapper converts these user-specified monthly max/ min climate data into hourly values using sinusoidal functions that assume the same weather conditions for each day of a given month. The model then solves for heat and water mass balance for the ground surface and the temperatures, wind speeds, and humidity ( 2 m above and below ground level). Solar radiation is computed from an embedded solar radiation algorithm (McCullough and Porter 1971). Sky thermal radiation was computed from an empirical regression based on the input for the 2 m above ground air temperature (Idso and Jackson 1969).
Finally, for the Niche Mapper simulations of ground sloths in their respective environments, a thermal comfort zone was determined for the estimated daily metabolic rate requirements for each taxon based on likely upper and lower bounds of daily metabolic activity. Total daily metabolic expenditure was calculated using a trapezoidal integration of hourly metabolic rate predictions. The lower bound of the target zone was determined by multiplying estimated basal metabolic rate in by 24 , while the upper bound was determined by multiplying estimated field metabolic rate
(i.e., basal metabolic rate multiplied by a factor of 2 converted to ) by the same factor. Within this zone, the simulated animals could function at normal metabolic levels as long as resources were readily available. For simplicity, zones of tolerance and survival above the upper bounds of the upper limits of the modeled thermal comfort zone are referred to herein as “cold stress” since the simulated animal would be expending additional metabolic energy to prevent hypothermia. Likewise, zones of tolerance and survival below the lower bounds of the thermal comfort zone are referred to as “heat stress” as the simulated animal would have to lower its metabolic rate below its estimated basal metabolic rate to prevent hyperthermia. The results of Fariña (2002) provided a sensitivity analysis to compare the thermal tolerances of ground sloths with metabolic rates more akin to large placental mammals, as opposed to a xen-arthran-like metabolism.

Statistical analysis

Predictions of daily metabolic rate (in ) that were generated by the climate simulations of each integument (fur coverage) model were plotted as a series of curves to determine which model best fit the estimated minimum and maximum daily metabolic rate expenditures, or the hours per day where basal metabolic rate can be maintained without the need for seeking shade, for each taxon. Once the integument models with the least amount of thermal stress were established, estimates of maximum hours of metabolic activity beyond levels of basal metabolic rate were compared across the geographic distribution of each taxon along with those of extant placental mammals of similar size ranges. Daily food consumption averages in extrapolated from integument models with the least amount of thermal stress for each simulated environment were fitted with a second-order polynomial regression curve to indicate possible differences in food intake among taxa. Nutritional value of leaves, grasses, and woody plants were derived from the DairyOne feed composition library (https://apps. dairyone.com/feedcomposition/; see Online Resource 1 for full list of climate and physiology parameters), and the wet mass in kg of forage consumed was additionally estimated from these simulations. The consumed food data were also fitted with second order polynomial regressions to statistically evaluate patterns of monthly food consumption in the respective predicted habitats of the sampled fossil taxa.

Results

Clumped isotope paleothermometry and REE ratios

Calculated REE index ratios for a majority of specimens analyzed here fell below a value of 1.0 indicating little diagenetic alteration in the samples (Table 1). An exception was the Promegatherium specimen (FMNH P14511), which yielded a REE index value of 36.6 and was considered to be likely diagenetically altered or an impure sample. The Eremotherium specimen (UF 312730) had the lowest REE ratio (0.02) that is consistent with both dentine to cortical bone REE index values determined for other individuals found at Halie 7G (range; 0.004-0.01; MacFadden et al. 2010; Online Resource 2). Both the Eremotherium (REE ratio; 0.42 ) and Nothrotheriops (REE ratio; 0.86) tooth specimens recovered from Inglis 1 A and Leisey 1 A , respectively, also fell within the known variation of REE index values recorded for fossils from individuals at both localities (range; 0.33-0.72; MacFadden et al. 2010). The corresponding clumped isotope data from the Promegatherium specimen FMNH P14511 resulted in a mean temperature estimate of and was determined to be too low for a realistic body temperature of a large mammal, being more similar to the lower end of body temperature estimates of ectotherms determined by similar previous analyses (Eagle et al. 2011, 2015; Griffiths et al. 2023). It is noted that the carbonate-phosphate offset was within the bounds of most extant vertebrates, indicating that alteration of the carbonate moiety may be subtle rather than extensive (Fig. 1b).
Fig. 3 Box and whisker plot of ambient temperature estimates for the environments inhabited by the four extinct ground sloth taxa sampled
Despite its low REE index, the Eremotherium specimen UF 312730 exhibited poor replication between analyses (high standard error; Table 1) and yielded clumped isotopic temperature estimates which are implausible to be in a physiologic range. This was likely the consequence of organic contamination that had not successfully been removed by pre-cleaning of the tooth specimen and/or incomplete cleanup stages of sample purification, as neither phosphate and REE indicated diagenesis (Table 1; Fig. 1b). This hypothesis is supported by the observed abnormally high values (Fig. 1c; Online Resource 2) in this sample indicative of the presence contamination creating isobaric interferences for clumped isotope measurements. The three remaining teeth specimens of Eremotherium, Nothrotheriops, and Megatherium that had low REE and carbonatephosphate offsets in the range of extant mammals (Table 1; Fig. 1b) yielded clumped isotope temperatures of , and , respectively, with good agreement between replicate analyses and relatively low values (Fig. 1c). As per established criteria, there is no reason to suspect that these samples are compromised by diagenesis or analytical artifacts, thus the resulting estimates are likely to represent the core body temperatures of extinct ground sloths.

Estimates of paleoclimate and physiological parameters

A wide range of ambient temperatures were modeled for the four genera of ground sloths analyzed (Fig. 3). The largest, Eremotherium, was predicted to have experienced variable ambient temperatures (range: ) throughout its extensive geographic range, though the lowest values were determined by outliers. The predicted ambient temperature differences for Nothrotheriops were similar in absolute values to Eremotherium, but colder minimum ambient temperatures were estimated for this taxon (range: C). Both Megatherium and Mylodon, however, were predicted to primarily have inhabited cold-to-temperate environments with minimum and maximum temperature values overlapping with the temperature ranges for Eremotherium and Nothrotheriops (Fig. 4). Estimated basal metabolic rate values for all four taxa sampled did not exceed (Fig. 4); all metabolic rate estimates were less than those approximated by previous analysis of Fariña (2002).

Metabolic chamber simulations and thermoneutral zone ranges

Metabolic chamber analyses revealed broad thermoneutral zones for dense fur (1,300-2,000 hairs/ ) across all four genera that decreased with an increase in coat thickness
Fig. 4 Bar chart of estimated basal metabolic rate (BMR) and field metabolic rate (FMR) in ground sloths. Values in unit watts (W) calculated from the xenarthran metabolic scaling equation (White and Seymour 2004) are compared to those (right) based on a ‘typical placental’ metabolism (Kleiber 1932; Fariña 2002)
(Fig. 5). Both Eremotherium and Megatherium had similar thermoneutral zones with coat thicknesses of 10 mm (range; , combined for both taxa), which overlap with the estimated ambient temperature values that both genera would have experienced throughout their geographic distribution (Fig. 3). Surprisingly, Megatherium with a dense fur and a coat thickness of 50 mm had the lowest lower critical temperature among the four genera ( ) that is equivalent to previous estimates for a densely furred Mylodon within its thermoneutral zone (Fariña, 2002), but lower than lower critical temperature estimates of Eremotherium with the same integument features applied herein ( ). Dense fur conditions for Mylodon and Nothrotheriops with coat thicknesses of also have thermoneutral zones that overlap with likely ambient temperature ranges that each taxon would have experienced, as well as ambient temperatures well beyond the lowest estimates modeled here ( and , respectively). Consistent trends of narrow thermoneutral zones resulting from sparse fur ( hairs ) with lower critical temperature estimates higher than and upper critical temperature
estimates ranging from were determined for all four taxa analyzed.

Microclimate simulations and integument thermal tolerances

Similar outcomes favoring dense fur coverages were generally observed when multiple climatic variables, including relative humidity, precipitation, wind speed, and ambient temperature were input in the refined paleoclimate models. Notably, Eremotherium was estimated to have no thermal stress in all environments with dense, 10 mm fur insulation (1,300-2,000 hairs/ ), but would have experienced prominent cold stress in the northernmost range of its distribution, although no thermal stress occurring in the southernmost latitudes during the coldest months of the year in both localities sampled (Fig. 6). This result contrasts with the resolved trends that sparse fur coverage indicated in northernmost and southernmost latitudes, which consistently showed no thermal stress during only the warmest months of the year. Dense fur insulation for Eremotherium was determined to be mostly detrimental by resulting
Fig. 5 Thermoneutral zones of . Eremotherium; b. Megatherium; c. Mylodon; and d. Nothrotheriops. Estimated thermoneutral zones were derived from metabolic chamber simulations with different fur densities and coat thicknesses
in overall heat stress, except for during the winter months in only the northernmost latitudes and year-round in the southernmost latitudes. Conversely, Eremotherium additionally was modeled to have no thermal stress with sparse 50 mm fur, but with more dense fur coverage resulting in constant heat stress in neotropical regions (Fig. 6).
The refined paleoclimate models showed that Megatherium would have avoided thermal stress year-round with full body, dense 30 mm fur insulation in all three of the simulated environments, whereas full body, sparse fur (both depths and densities) resulted in cold stress throughout the year (Fig. 6). Both Mylodon and Nothrotheriops would also have avoided any thermal stress with dense fur insulation (density; hairs ). In particular, Mylodon
was modeled to have no thermal stress with 10 mm fur during the warmest months of the year and also with 30 mm fur during the coldest months at northern and mid latitudes (Fig. 6). Both Mylodon at the southernmost extent of its range and Nothrotheriops across its entire geographic range were determined to have no heat stress with dense fur with a coat thickness of 50 mm .

Daily hours of active metabolic expenditure

Estimates of hours per day where basal metabolic rate can be maintained without seeking shade for the integument conditions resulting in the least thermal stress (by fur coverage, depth, and density) across the known geographic
Fig. 6 Modeled trends of ground sloth daily metabolic activity (kJ day -axis) in response to climatic factors over the course of a year (month; -axis) across the northernmost, middle, and southernmost range of each genus. Eremotherium target metabolic rate
distributions of each taxon analyzed averaged approximately 15 h day (Fig. 7). While the megatheres were modeled to have large metabolic expenditure in cold climates (e.g., northern geographic ranges of Eremotherium;
day ; Megatherium target metabolic rate=51,944 day ; Mylodon target metabolic rate day ; Nothrotheriops target metabolic rate
mid-latitudinal and southern geographic ranges of Megatherium), the greatest number of hours of metabolic activity expended per day ( ) were predicted in the southern latitudinal ranges of both Megatherium and Mylodon.
Fig. 7 Bar chart of the average number of hours per day where metabolic rate can surpass basal levels in the northernmost middle, and southernmost range of each ground sloth genus compared to extant placental mammals of similar body sizes and extant tree sloths. Shown as silhouettes are the body size analogues in order left-to-right: African savanna elephant (Loxodonta africana; for Eremotherium), Asian elephant (Elephas maximus; for Megatherium), Indian rhinoceros
(Rhinoceros unicornis; for Mylodon), brown bear (Ursus arctos; for Nothrotheriops), Bradypus and Choloepus. Silhouettes were obtained from PhyloPic.org (Michaud M, Taylor J, Keesey M, Morrow C, Jenkins X) and are under public domain (Total hours of activity for extant analogues were based on daily hours of sleep reported in the literature (Stelmock and Dean 1986; Deka and Sarma 2015; Gravett et al. 2017; Cliffe et al. 2023)
Nothrotheriops showed a consistent maximum metabolic activity value of 15 h day throughout its geographic range (Fig. 7). Comparative data on daily energetic expenditures indicate that extant body size analogues of Eremotherium, Megatherium, and Mylodon (e.g., African elephant, Indian elephant, and Indian rhinoceros, respectively) all have values of 22 h day (Deka and Sarma 2015; Gravett et al. 2017) while that for much smaller Nothrotheriops, a brown bear (Ursus arctos) analog, was measured at 20 h day (Stelmock and Dean, 1986). In sum, all ground sloth taxa were predicted to have less metabolic activity than modern placental mammalian analogues of similar body masses, but are far more active than extant tree sloths (e.g., ; Cliffe et al. 2023).

Monthly dietary intake

The estimated distributions of energy intake across all four genera indicated no statistically significant trends in monthly food consumption with a change in environment (Fig. 8). While a general increase in food consumption was demonstrated for the large-bodied megatheres, Eremotherium was modeled to experience little-to-no increase in
food consumption ( ) throughout its wide geographic distribution. The genus Megatherium was modeled to experience a moderate increase in monthly food intake in the southernmost latitudes , whereas the simulations of Mylodon predicted a similar increase in food intake in southern latitudes , notably overlapping with that of Megatherium. Finally, the data trends for Nothrotheirops showed a moderately strong decrease in monthly food intake in the southern latitudes during the summer months.

Discussion

Geochemistry of fossil teeth and interpretations of physiology

The similarities in clumped isotope temperatures from teeth recovered from different locations, coupled with the low REE index values, provides confidence that these estimates may reflect the actual core body temperatures of extinct ground sloths. These temperatures are notably lower than the clumped isotope derived temperatures of other large
Fig. 8 Scatter plot of kilograms of wet food mass consumed per day throughout the year for each integument model exhibiting minimal thermal stress for northernmost (blue), middle (green), and southernmost ranges (red). Data for each genus of ground sloth are fit with second order polynomial regressions: Eremotherium (y , ); Megatherium ( , ); Mylodon ( ; Nothrotheriops . 0.669)
mammals including values of approximately for an Indian elephant and extinct mammoths (Eagle et al. 2010), and for a baleen whale (Griffiths et al. 2023). Tooth material from ground sloth specimens showing little diagenetic alteration were obtained from three different depositional environments and produced temperature estimates within the range of extant xenarthran body temperatures. Given that the estimates for Nothrotheriops were derived from a single tooth sample with a relatively high REE ratio (0.86), some degree of caution should be given to the interpretation of the results from this specimen. Measurements obtained from the Eremotherium sample UF 312730 yield higher values in both and which could indicate contamination of the samples resulting in isobaric interference and a subsequent overestimation of and unrealistically cold derived temperatures (Fig. 1c). This aside, these novel data suggest that extinct ground sloths may have been heterothermic, and likely fluctuated their core body temperatures as a main means of thermoregulation. The core body temperature estimates derived from the clumped isotope analyses are indeed within the range of error of known body temperatures for extant xenarthrans ( ), further verifying body temperature values determined for the Niche Mapper simulations and suggesting that the observed ubiquitous similarities in thermophysiology are an ancestral trait in this clade.
Clumped isotope estimates of body temperature, however, are on the lower end of Niche Mapper simulations for the mean estimate of determined for UF 95869. This could be explained by the measurements taken from specific tooth loci reflecting potential periods of torpor akin to what has been proposed for troodontids (Tagliaventoor et al. 2023), estimates from teeth might may also reflect slight differences in core body temperature and the temperature of tooth mineralization as has been explored in interpreting clumped isotope values derived from regionally endothermic sharks (Griffiths et al. 2023). These hypotheses should be explored in future analyses focusing on taking samples from specific regions of internal orthodentine to determine if core body temperatures were fluctuating throughout the life of the animals. Such methods were previously applied to teeth from Eremotherium to infer the environmental conditions during the life of a single individual, specifically using and measurements (Larmon et al. 2019). Analyses of more specimens spanning a broader range of climate conditions similar to the approaches that have been used on non-avian dinosaur remains (Fricke and Rodgers 2000 ; Amoit et al. 2006) may also be informative. Nevertheless, the clumped isotope measurements along with the results of the physiological modeling analyses herein collectively indicate that extinct ground sloths had lower body temperatures and basal metabolic rates than most large terrestrial placental mammals.

Comparisons with previous analyses

The low body temperature estimates obtained from fossil ground sloth teeth suggest lower metabolic rates akin to extant xenarthrans and contradict the higher values of metabolic rates based on estimated blood flow rate (Varela et al. 2024). Varela et al. (2024) hypothesized that small nutrient foramen diameters of the femur, and thus restricted blood flow rates, along with the fossorial and arboreal habits of extant xenarthrans all correlated with low metabolic rates. However, raw data from that study shows overlap between femoral nutrient foramen diameters and blood flow rates of extant xenarthrans and placental mammals with elevated metabolic rates. For example, the nutrient foramina of both the giant anteater (Myrmecophaga tridactyla) and the leopard seal (Hydrurga leptonyx) have a measured diameter of 0.9 mm despite their extremely different functional habits and ecology. Given these types of discrepancies in their analysis, we caution that using measured nutrient foramen diameter alone is not a reliable predictor of metabolic rate, and it could likely reflect an artifact of similar body size. The core body temperature estimates obtained from the established clumped isotope analyses used herein potentially offer a more direct estimation of core body temperature and expand our understanding of the physiological limitations potentially experienced by extinct ground sloths.
The paleoclimate simulations also yielded metabolic rate values drastically different values from those of the Fariña (2002) study on ground sloth thermal energetics. In addition to refined estimates of clade specific basal metabolic rates for ground sloths, numerous other environmental factors beyond ambient temperatures were input and modeled across the geographic distributions of each genus to determine more exact boundaries of thermal neutrality. The basal metabolic rate estimates derived from the xenarthran scaling equation (White and Seymour 2004) produced metabolic rates representing only of those for mammals with typical placental metabolism as determined by the relationship of Kleiber (1932) (Fariña 2002). Substantially lower basal metabolic rates were predicted in the present study, thus producing markedly different results for thermal tolerance ranges across the various integument conditions tested compared with those already reported. The main previous study (Fariña 2002) proposed that a sparsely furred Megatherium with a body mass of and metabolism that was that of a typical placental mammal would be within its thermoneutral zone at ambient temperatures of versus the results from the present climate simulations, which suggest that Megatherium under similar conditions would have been cold stressed.
Traditional determinations of thermoneutral zones using the xenarthran metabolism as a proxy suggest
Eremotherium and Megatherium with dense fur and a coat depth of 10 mm would be akin to those of sparsely furred models (see Fig. 5). However, the microclimate models employed herein demonstrated that environmental factors aside from ambient temperature (i.e., relative humidity, wind speed, etc.) and metabolic activity (i.e., locomotion, foraging, etc.) accounted for the predicted outcome that dense 10 mm thicknesses would benefit Eremotherium only during the warmest months of the year at its northernmost and southernmost ranges year-round. The novel microclimate models also showed that having dense fur at a coat thickness of 10 mm would result in constant heat stress for Eremotherium in the neotropics, yet constant cold stress for Megatherium throughout its geographic distribution (Fig. 6). Despite heat stress modeled to be present in the tropical populations of Eremotherium having a full body fur coverage at a depth of , the condition of dense 10 mm fur depths would allow it to exhibit no intense thermal stress year-round, even during the hottest months of the year throughout its entire geographic range. The findings for southernmost region for Mylodon and those data for Nothrotheriops throughout its entire range with dense 50 mm fur depths, suggest nearly complete yearly protection from thermal stress in their respective environments, with mild cold stress for the three coldest months of the year (predicted temperatures; ). These results also contrast with previous estimates (Fariña 2002) that thermal neutrality for Mylodon having a basal metabolic rate half of the typical placental metabolism and dense fur 40 mm deep would be at an ambient temperature of .

Megathere thermal neutrality and geographic distribution

The data from the paleoclimate analyses revealed critical insight into the thermoregulatory thresholds of extinct ground sloths and how these would have impacted their life appearance by condition of their integument. In particular, the results for Eremotherium suggest that a wide range of integument models could be plausible depending on ambient temperatures of their past environments. Given that Eremotherium had a geographic distribution spanning as far north as Virginia in North America and as far south as Brazil in South America, it is plausible that this genus had varying degrees of fur coverage throughout its geographic range. Sparse, full body fur coverage would have been beneficial at mid-ranges with no thermal stress experienced year-round. Changes in fur depth would also have been beneficial during warmer months for Eremotherium at the most northern extent of its range, although there is no evidence of coat shedding in extant xenarthrans. Nonetheless, dense 10 mm fur could have prevented thermal stress at average ambient
Fig. 9 a. Map of geographic distribution for Pleistocene-aged sites of fossil remains of ground sloths Eremotherium, Megatherium, Mylodon, and Nothrotheriops spanning from 2.5 MYA to 11,000 years ago. Areas shaded blue represent the extent of glaciers during the last glacial maximum 21,000 years ago and red and green circles indicate locations where the paleo art is depicted. Data: political boundary base map was generated using QGIS3’s ‘World Map’; hand-drawn ice sheet overlay informed and modified after published data (Ehlers and Hughes 2018); spatial data for species distributions was retrieved from
temperature ranges of that Eremotherium was expected to have experienced (Fig. 6). Moreover, because Eremotherium crossed the Isthmus of Panama into North America but subsequently retreated to South America during the last glacial maximum (McDonald and Lundelius 2009) it is possible for it to have undergone migrations during times of extreme seasonality. Due to the modeled presence of mild cold stress with dense 10 mm fur and no thermal neutrality with a 50 mm fur depth in the southernmost ranges, it is therefore suggested that Eremotherium was well-adapted to the neotropical conditions expected to be prevalent in the
the Paleobiology Database (paleobiodb.org). b. Reconstruction of an adult and infant Eremotherium going through seasonal fur shedding while wallowing in a mudhole outside the cave den of an emerging Nothrotheriops in the early spring of Florida during the Pleistocene (location indicated by green circle in panel 9a). c. Reconstruction of Megatherium enduring the cool, windy Pampas of Argentina with a thick fur coat (location indicated by red circle in panel 9a). Paleoart by Haider Syed Jaffri
southernmost extent of its range (Fig. 9a). This hypothesis notably agrees with the absence of Eremotherium in North America during the last glacial maximum.
Despite being similar in body size to Eremotherium, the genus Megatherium is expected to have exhibited constant cold stress with only sparse fur coverage throughout its geographic range. This can be explained by Megatherium being known to have inhabited cooler, dry environments in the pampas of Argentina and Chile. Using similar modeling approaches to physiology as presented herein, previous sensitivity analyses (Lovelace et al. 2020) revealed that ambient
temperature and wind speed values were the two climatic factors that most strongly influence predicted values of field metabolic rate, whereas the variables of relative humidity and cloud coverage had little effect on simulated metabolism. Thus, dense fur for Megatherium might have been needed to counteract greater exposure to the elements (e.g., wind, rain, or snow) in more open environments (Fig. 9c). This integument condition is reflected in the present climate simulation results which indicate that with a full body coverage of fur at depths of 10 mm , Megatherium would have been constantly cold-stressed regardless of fur density, but it would have complete thermal neutrality in all latitudinal distributions with dense ( hairs ) fur with a coat depth of 30 mm . Conversely, by having dense 50 mm fur, Megatherium would have been constantly heat-stressed by over insulation (Fig. 6). Furthermore, considering that fur depths of resulted in predicted heat stress for similarly-sized Eremotherium (with the exception of the coldest months at northern extents of its range), it is reasonable to speculate that neotropical conditions may have acted as a thermal gradient restricting Megatherium to southern latitudes and occupation of the westernmost regions of Peru and Ecuador (Fig. 9a; Pujos and Salas 2004). This hypothesis is supported by models of environmental distributions of various Pleistocene xenarthrans by Varela et al. (2018) which demonstrate little-to-no overlap between the geographic ranges of Eremotherium and Megatherium suggesting distinct differences in habitat preference.

Thermoregulation in Mylodon and Nothrotheriops

Both Mylodon and Nothrotheriops were modeled to require dense fur depths across their respective ranges of predicted ambient temperature exposure, which is directly observed in the preserved pelts from both genera (Ridgewood 1901; Hausman 1929). Nevertheless, given that Mylodon was known to inhabit the Andes Mountains in South America, and Nothrotheriops would have experienced paleoclimate temperatures as low as (Fig. 3), behavioral thermoregulation would have been expected in cold climates in addition to having thick fur coverage. One supposition is that thermal neutrality could have been maintained by taking shelter in caves during periods of thermal stress. Mean annual temperature estimates of caves in the American southwest where remains of Nothrotheriops have been documented reveal stable temperature ranges that would have provided protection against periods of thermal stress (6.6-17 ; McDonald 2022; Fig. 9b). Localities for Mylodon in Chile such as Chica cave exhibit a similar range of stable temperatures (4.83-17.4 ; Nehme et al. 2022). A second hypothesis is that both taxa may have undergone
periods of torpor and/or hibernation as a behavioral means of cold tolerance.
Isotopic and taphonomic studies on the extinct Patagonian ‘panther’ (Panthera onca mesembrina) revealed that bite marks attributed to this extinct predator have been found on the bones of juvenile and adult individuals of Mylodon in Cueva del Milodón (Prevosti and Martin 2013; Martin 2018). While Mylodon was smaller than the megatheres, an adult individual was still much larger than a jaguar, which suggests that individuals found with bite marks in caves were preyed upon while they were undergoing torpor, when they would have been most vulnerable to attack. Additional support for the above hypotheses proposed herein can also be linked to the behavioral thermoregulation strategies observed in modern armadillos. For example, the ninebanded armadillo (Dasypus novemcinctus) is known to seek shelter during periods of extreme cold and employ a combination of torpor with episodic foraging (McBee and Baker 1982; Knight et al. 2020). Regular torpor also has been documented in the pichi armadillo (Zaedyus pichiy) (Superina and Boily 2007; Superina and Abba 2014), and it is the only known extant xenarthran to undergo true hibernation.
It is also possible that Mylodon and Nothrotheriops may have had higher activity levels than those estimated in the present analysis. However, higher metabolic demands in the smaller ground sloths would only be useful in temperatures well below the range of ambient temperatures predicted for both taxa. Accordingly, another potential explanation for their thermal energetics involves having a diet richer in energy than their megathere conspecifics. A recent analysis of amino acids (Tejada et al. 2021) suggested that the diet of Mylodon could have consisted of animal protein due to the nitrogen signatures found in their fossilized bone material. Though the exact timings of nitrogen signatures indicative of omnivory could not be determined, animal protein consumption would have had thermoregulatory benefits for Mylodon. A protein-rich diet provides some flexibility to thermal neutrality in cold climates, and it is possible that Mylodon consumed animal matter sporadically during the coldest months, in addition to foraging on grasses (Varela et al. 2023). Conversely, nitrogen isotopes recovered from Nothrotheriops fossil remains definitively did not show evidence for an omnivorous diet, but instead, this genus was eating carbohydrate-rich desert plants such as globemallow and casava root (Poinar et al. 1998).
Mylodontids also possessed dermal ossicles of varying shapes and sizes which may have provided protection against predators (Collins 1933; McDonald 2018; Toledo et al. 2021). Given the vascularized nature of these structures in ground sloths as a whole, Cartelle and Bohórquez (1986) hypothesized that the dermal ossicles may have played a role in thermoregulation as well. Yet, histological examinations
of Mylodon ossicles by Toledo et al. (2021) show no features indicating the storage of adipose tissue akin to those in extant armadillos (Krmpotic et al. 2009). While measurements of the thermal conductance have yet to be performed, Toledo et al. (2021) implied that the fur of mylodontids such as Mylodon would have provided a greater amount of insulation compared to vascularized osteoderms.

Integument and behavioral factors

Giant ground sloths have historically been reconstructed as being densely furred taxa and the novel data provided by the Niche Mapper simulations do little to challenge this integument appearance. In contrast, the modeling results disagree with more sparsely furred, elephant-like integument reconstructions of Eremotherium and Megatherium (Fariña 2002). Despite the obvious discrepancies between studies on thermal neutrality and its effects on ground sloth integument, several critical behavioral and ecological preferences can now be inferred. For example, dental microstructure and microware patterns (Green and Kalthoff 2015), in addition to isotopic analyses (Dantas et al. 2021), suggest a completely herbivorous diet for Eremotherium and Megatherium despite previous speculation that they were possibly omnivorous and/or engaged in scavenging behavior (Fariña and Blanco 1996). Indeed, the diets of both taxa would have primarily consisted of leaves and other high-growing vegetation (Bargo et al. 2006; Green and Kalthoff 2015; Bocherens et al. 2017), which overlap with those of their closest extant relatives Bradypus spp. (Cliffe et al. 2015; Pauli et al. 2016; Preslee et al. 2019; Delsuc et al. 2019; Tejada et al. 2024). Leaves are poor in nutritional value correlating with the extremely low metabolic rates observed in Bradypus variegatus (Cliffe et al. 2015; Pauli et al. 2016) and may be a primary factor for why megatheres would have had low metabolic rates for mammals of their massive body size.
One of the main outcomes of our is that megatheres are predicted to have a basal metabolic rate more similar to modern xenarthrans than other placental mammals. Yet, they also could have exhibited a number of behaviors similar to large extant mammals with the typical placental metabolism. Notably, the bone bed of Eremotherium at Tanque Loma in southwest Ecuador was found to contain multiple individuals that died in a stagnant marsh or swamp. The assemblage was interpreted as evidence of herding behavior in megatheres (Lindsey et al. 2020). It was further suggested that Eremotherium engaged in wallowing behavior for thermoregulatory purposes akin to that observed in a modern hippopotamus (Stears et al. 2018). The modeled variation in both the density and depth of fur across the geographic range of Eremotherium is arguably more broadly correlated with a similar thermoregulatory strategy. Large migratory
ungulates such as the American bison (Bison bison) shed their dense winter coat to a lighter, patchy coat that coincides with wallowing behavior (McMillan et al. 2000). Large Megatherium also could have engaged in wallowing but would not have required shedding due to its year-round thermal neutrality when having a dense, deep 30 mm fur coat (Figs. 5 and 6).
Social behavior dynamics of the extinct megatheres and large mylodonts (e.g., Lestodon) are further proposed to be consistent with those observed in large extant placental mammals (Tomassini et al. 2020; Lindsey et al. 2020), yet differ greatly from the solitary lifestyles of extant tree sloths (Taube et al. 1999). One advantage to more complex social behavior could have been shared parental care (also inferred from the bonebed of Eremotherium at Tonque Loma; Lindsey et al. 2020) with multiple adults protecting juveniles that are most susceptible to predation. Conversely, a sole female Bradypus spp. raises her single offspring, and if it falls from the tree during lactation/weaning, then she does not invest the metabolic energy to retrieve it (Taube et al. 1999).
Again, the present findings predict densely furred integument for Mylodon and Nothrotheriops as demonstrated by their mummified pelts (Ridgewood 1901; Hausman 1929). Given that there is a lack of fossil bonebeds consisting of multiple individuals, it is suggested that both genera potentially lived solitary lifestyles similar to extant tree sloths and armadillos (Naples et al. 1990). Mylodontids also possess robust limb bones that were resistant to bending forces akin to extant armadillos (Copploe et al. 2015), suggesting that some genera were capable of digging their own dens or burrows (Bargo et al. 2000). Solitary denning behavior as seen in modern brown bears (González-Bernardo et al. 2020), a body size analog of both smaller ground sloth taxa, also may have been typical of Mylodon and Nothrotheriops in cave dwellings. Moreover, an ecological preference that could have overlapped with ursids is diet. A modern species of two-toed sloth, Choloepus hoffmanni, is a close relative to Mylodon (Preslee etc. al. 2019; Delsuc et al. 2019; Tejada et al. 2024), and has been known to occasionally consume meat in captivity (Hayssen 2011). Species of the genus Choloepus are also more active compared to Bradypus variegatus (Sunquist and Montgomery 1973), which specializes in a diet of leaves presumably akin to its megathere relatives.
Last, while it is tempting to focus on phylogenetic constraint as the singular reason for the predicted low basal metabolic rates of extinct ground sloths, it should be noted that traits can be secondarily lost or acquired throughout the evolutionary history of a clade. Prior to GABI, the very large herbivore guild of South America was dominated by notoungulates and astrapotheres that were mixed grazers and browsers that went extinct by the end of the Miocene prior to the earliest occurrences of the exceptionally large
megatherines such as Eremotherium and Megatherium (Croft et al. 2020). The largest of the notoungulates were members of the families Toxodontidae, Leontiniidae, and Homalodotheriidae. While toxodontids and macraucheniids lived alongside megatherines through the late Pleistocene (and reached body masses of ; Croft et al. 2020; Nelson et al. 2023), astraoptheres, leontinids, and homalodotheres ranging in body masses from 375 to (Croft et al. 2020; Nelson et al. 2023) went extinct by the end of the middle to late Miocene. By contrast, most South American Miocene ground sloth genera were much smaller in size compared to contemporary notoungulates and astrapotheres, and larger ground sloths that would evolve in the Pliocene and Pleistocene ( ; Pant et al. 2014: table s1). A lower metabolic rate in ground sloths would have reduced resource demands and potentially limited competition with notoungulates and astrapotheres. This in turn could be understood as an “ecological steppingstone” towards occupying the vacant browser niche without requiring as many calories as a similarly-sized mammal with a typical placental metabolism. The observed trends in daily food intake in Eremotherium and Megatherium support this notion by suggesting little change in the amount of food consumed across their various geographic regions (Fig. 7). This major ecological factor also could have limited intraspecific competition for resources in megatheres, but notably reduced the pressure of interspecific competition with herds of proboscideans (e.g., mammoths and mastodons) with much greater energetic demands that entered South America following GABI.

Limitations

While use of REE index values in the present study produced values within acceptable limits as reported by MacFadden et al. (2010), there are inconsistencies between this method and the phosphate data. Namely, FMNH P14511 has an abnormally high REE index of 36.6 and acceptable phosphate measurements within the range of extant placental mammals. The REE index is the first technique used to assess the degree of diagenesis in fossil xenarthran teeth, and this method has yet to be tested on specimens outside of the Pliocene and Pleistocene localities of Florida. The high REE concentration in the inner orthodentine of FMNH P14511 could be explained by it being found in older Miocene deposits, although this is inconsistent with the phosphate measurements which indicate that it is well-preserved. Moreover, a specimen of Bradypus (UF 25980) reported in the supplementary materials of MacFadden et al. (2010) has a dentine/bone index value of 30.6, which is comparable to the index value of FMNH P14511. That said, the criteria set by MacFadden et
al. (2010) of acceptable REE index values are based on measurements of enamel-bearing mammals from the same localities as fossil xenarthrans rather than from modern bone, dentine, and enamel. We were unable to procure enamel samples from non-xenarthran placental mammals from the same localities, which would have aided in determining the relative degree of diagenesis in the ground sloth teeth used in this study. Clumped isotope temperatures are also based on being substituted into the and OH sites of hydroxylapatite, while REE concentrations occur on the site of hydroxylapatite (Kohn et al. 1999), suggesting that geochemical alteration may not be the same at either bonding site. Based on these factors, and the REE index not accounting for the carbonate contamination of UF 312730 (Fig. 1c; Online Resource 2), we argue that use of the REE index alone does not always provide a robust indicator of diagenesis or lack thereof.
Suárez and Passey (2014) showed that clumped isotope temperatures found in fossil bone can be intermediate between ambient temperatures during the time of fossilization and body temperature. However, predominantly bone was used in their analysis, which is more prone to diagenetic alteration compared to dental tissues such as enamel and dentine (MacFadden et al. 2010). Their dataset also did not include internal orthodentine from fossil xenarthrans, which is less prone to diagenesis than external dentine but not as resistant as enamel. Cathodoluminescence performed by Larmon et al. (2019) demonstrated that internal orthodentine in an Eremotherium tooth exhibits very few cracks where mineral impregnation can occur and they concluded that little diagenesis had occurred. Moreover, Suárez and Passey (2014) did not compare phosphate concentrations of their samples compared to extant mammalian bone, whereas our previous studies (Eagle et al. 2010, 2011) have shown that the temperatures obtained from the teeth of extinct vertebrates using clumped isotope mass spectrometry reflect the core body temperature of the organism, and not an intermediate between the body temperature and the temperature of the environment.
The novel data generated herein offer several intriguing implications for ground sloth integument, thermal energetics, ecology, and behavior. However, it is acknowledged that the results may not be definitive. The presence of low body temperatures determined by clumped isotope paleothermometry is suggestive that ground sloths were heterothermic, akin to their closest living relatives. However, further analyses with a larger sample size consisting of specimens from various depositional environments are needed to verify this hypothesis. Given that some of the cortical bone samples did not come from the same individuals as the respective tooth samples, this could impact the values of REE ratios due to potential differences in REE uptake
during fossilization. While samples were taken from other sections of the fossil teeth used for the geochemical analyses, these were not used due to restrictions of time and funding for this study. Further analyses of these samples could better indicate evidence of a past heterothermic condition as the inner orthodentine could have recorded different isotopic signals, and thus different body temperature estimates.
The modeling simulations presented here assumed consistent fur depth and density across the entire body, whereas these parameters can vary about the body surface area and this factor should be taken into account when considering the full implications of the results. Other factors such as skin thickness, hair shaft thickness, underfur, and the presence/ absence of hollow hair shafts should additionally be taken into account with future analyses. Admittedly, we tested only average adult body masses, and our results could vary if tested while accounting for intraspecific variations in body masses. Pliocene climates were excluded from this analysis due to Pleistocene climate variables being easier to use for extrapolation of reasonable monthly values from modern climate data. The inclusion of Pliocene Epoch climate data would give more insight into how integument and associated fur coverage could have varied in earlier populations of both Eremotherium and Megatherium. Lastly, geologic ages of localities sampled for the microclimate simulations cover a broad range of geologic time encompassing the entire span of the Pleistocene from years ago, and more precise dates were not available in the Paleobiology Database (Online Resource 1). More precise estimates of paleoclimate and dating of ground sloth localities could have greatly improved the microclimate simulations.

Conclusions

Clumped isotope palaeothermometry of inner orthodentine provides convincing evidence for extinct ground sloths having body temperatures akin to extant xenarthrans and potentially being heterothermic, regardless of absolute differences in body sizes. These findings combined with novel biophysical modeling results in Niche Mapper demonstrate that lower body temperatures and basal metabolisms, in addition to environmental variables, impacted ground sloth geographical distribution, thermoregulation ability, and integument form. The “hairless model of integument” previously proposed for Megatherium is untenable when more accurate estimates of xenarthran basal metabolism, while the closely related Eremotherium has been suggested to be sparsely furred in neotropical climates. Giant ground sloths likely had a variety of integument appearances from dense, thinner fur coverage in the neotropics to dense, thicker fur coverage in cold environments. Smaller ground sloths would
have needed dense, thick fur supplemented with behavioral thermoregulation to endure cold climates throughout the year. The modeling data also point to surprising adaptability of xenarthran physiology in both the largest and smaller ground sloths. Moreover, the predicted integument models for the largest taxa are believed to be reliable unless refuted by a future specimen with extensive preserved integument. Prospective future models must examine other taxa (e.g., Megalonyx, Thalassocnus, Lestodon, etc.) and include histological studies to help determine how low the basal metabolic rates of ground sloths could have been beneficial, therefore providing further insight into the physiology of these and other fossil xenarthrans and how their metabolic rates impacted their integument appearance, growth rates, ecology, and behavior.
Supplementary information The online version contains supplementary material available at https://doi.org/10.1007/s10914-0 24-09743-2.
Acknowledgements We immensely thank V. Sinkonen, Y. Kutucu, and Universidad de la República for supplying us with 3D ground sloth models for this research and adding curves to them. We also thank W. Simpson (Field Museum of Chicago) and B. MacFadden and R. Hulbert (University of Florida) for providing fossil samples for geochemical analyses. R.J. Flores was supported by a Basic Energy Sciences grant DE-FG02-8361ER16402 and A.K. Tripati was supported by a grant from the Heising-Simons Foundation 2022-3314. Funding for the analytical infrastructure of the Boise State Isotope Geology Laboratory used for measuring trace elements by LA-ICPMS in this study was provided by National Science Foundation grants EAR-1735889 and EAR-1920336. Thanks to M. Kohn (Boise State University) for feedback on interpreting the rare earth element data. In addition, H.G. McDonald, J. Tejada, R. McAfee, M. Curiel, and S. Lee are thanked for their insight, feedback and critical comments.
Author contributions M.D.D. developed the concepts and approach, collected and analyzed data, prepared figures, and drafted and revised the manuscript; W.P.P. developed code and facilitated data analysis; P.D.M. developed code and facilitated data analysis; D.M.L. performed data mapping; R.J.F. analyzed data; A.K.T. analyzed data; R.A.E. developed analysis approach and interpreted data; D.M.S. developed analysis approach and interpreted data; M.T.B. developed the concepts and approach, and drafted and revised the manuscript.
Funding Basic Energy Sciences grant DE-FG02-8361ER16402; Heis-ing-Simons Foundation grant 2022-3314; National Science Foundation grants EAR-1735889 and EAR-1920336.
Data availability All data collected and analyzed for this study are included in this published article in adherence with disclosure policy of the journal. The authors are sharing new data alongside the results of this report; Niche Mapper input has been reposited to Zenodo (htt ps://doi.org/10.5281/zenodo.10626431) and clumped isotope data has been reposited to EarthChem (https://doi.org/10.60520/IEDA/11312 7).

Declarations

Competing interests The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.o rg/licenses/by/4.0/.

References

Amiot R, Lécuyer C, Buffetaut E, Escarguel G, Fluteau F, Martineau F (2006) Oxygen isotopes from biogenic apatites suggest widespread endothermy in Cretaceous dinosaurs. Earth Planet Sci Lett 246(1-2):41-54. https://doi.org/10.1016/j.epsl.2006.04.018
Anderson NT, Kelson JR, Kele S, Daëron M, Bonifacie M, Horita J, et al (2021) A unified clumped isotope thermometer calibration ( ) using carbonate-based standardization. Geophy Res Lett 48(7):e2020GL092069. https://doi.org/10.1029/2020G L092069
Barbosa FHDS, Alves-Silva L, Liparini A, Porpino KDO (2023) Reviewing the body size of some extinct Brazilian Quaternary xenarthrans. J Quat Sci https://doi.org/10.1002/jqs. 3560
Bargo MS, Vizcaíno SF, Archuby FM, Blanco RE (2000) Limb bone proportions, strength and digging in some Lujanian (Late Pleis-tocene-Early Holocene) mylodontid ground sloths (Mammalia, Xenarthra). J Vertebr Paleontol 20(3):601-610. https://doi.org/10 .1671/0272-4634(2000)020[0601:LBPSAD]2.0.CO;2
Bargo MS, De Iuliis G, Vizcaíno SF (2006) Hypsodonty in Pleistocene ground sloths. Acta Paleontol Pol 51(1):53-61.
Bernasconi SM, Daëron M, Bergmann KD, Bonifacie M, Meckler AN, Affek HP et al (2021) InterCarb: A community effort to improve interlaboratory standardization of the carbonate clumped isotope thermometer using carbonate standards. Geochem Geophys Geosyst 22(5):e2020GC009588 https://doi.org/10.1029/2020GC009 588.
Billet G, Germain D, Ruf I, de Muizon C, Hautier L (2013) The inner ear of Megatherium and the evolution of the vestibular system in sloths. J Anat 223(6):557-567. https://doi.org/10.1111/joa. 12114
Blanco RE, Czerwonogora A (2003) The gait of Megatherium Cuvier 1796. Senckenb Biol 83: 61-68.
Bocherens H, Cotte M, Bonini RA, Straccia P, Scian D, Soibelzon L, Prevosti FJ (2017). Isotopic insight on paleodiet of extinct Pleistocene megafaunal xenarthrans from Argentina. Gondwana Res 48:7-14. https://doi.org/10.1016/j.gr.2017.04.003
Boeskorov GG, Mashchenko EN, Plotnikov VV, Shchelchkova MV, Protopopov AV, Solomonov NG (2016). Adaptation of the woolly mammoth Mammuthus primigenius (Blumenbach, 1799) to habitat conditions in the glacial period. Contemp Probl Ecol 9:54453. https://doi.org/10.1134/S1995425516050024
Boyles JG, Levesque DL, Nowack J, Wojciechowski MS, Stawski C, Fuller A, et al (2019) An oversimplification of physiological principles leads to flawed macroecological analyses. Ecol Evol 9(21):12020-12025. https://doi.org/10.1002/ece3.5721
Brassey CA, Gardiner JD (2015) An advanced shape-fitting algorithm applied to quadrupedal mammals: improving volumetric mass estimates. R Soc Open Sci 2(8):150302. https://doi.org/10.1098 /rsos. 150302
Campbell KL, Roberts JE, Watson LN, Stetefeld J, Sloan AM, Signore AV, Howatt JW, Tame JR, Rohland N, Shen TJ, Austin JJ (2010). Substitutions in woolly mammoth hemoglobin confer biochemical properties adaptive for cold tolerance. Nat Genet 42(6):53640. https://doi.org/10.1038/ng. 574
Cartelle C, Bohórquez GA (1986). Presença de ossículos dérmicos em Eremotherium laurillardi (Lund) Cartelle & Bohórquez 1982 (Edentanta, Megatheriidae). Iheringía Ser Geol 11:9-14.
Chavagnac V, Milton JA, Green DRH, Breuer J, Bruguier O, Jacob DE, et al. (2007). Towards the development of a fossil bone geochemical standard: An inter-laboratory study. Anal Chim Acta 599(2):177-190
Christiansen P, Fariña RA (2003) Mass estimation of two fossil ground sloths (Mammalia, Xenarthra, Mylodontidae). Senckenberg Biol 83(1):95-101
Clarke A, Rothery P (2008) Scaling of body temperature in mammals and birds. Funct Ecol 22(1):58-67. https://doi.org/10.1111/j. 136 5-2435.2007.01341.x
Cliffe RN, Haupt RJ, Avey-Arroyo JA, Wilson RP (2015) Sloths like it hot: ambient temperature modulates food intake in the brownthroated sloth (Bradypus variegatus). PeerJ 3:e875. https://doi.o rg/10.7717/peerj. 875
Cliffe RN, Scantlebury DM, Kennedy SJ, Avey-Arroyo J, Mindich D, Wilson RP (2018) The metabolic response of the Bradypus sloth to temperature. PeerJ 6:e5600. https://doi.org/10.7717/peerj. 5600
Cliffe RN, Haupt RJ, Kennedy S, Felton C, Williams HJ, Avey-Arroyo J, Wilson R (2023) The behaviour and activity budgets of two sympatric sloths; Bradypus variegatus and Choloepus hoffmanni. PeerJ 11:e15430. https://doi.org/10.7717/peerj. 15430
Collins RL (1933) Mylodont (ground sloth) dermal ossicles from Colombia, South America. J Wash Acad Sci 23(9):426-429.
Conte GL, Lopes Le, Mine AH, Trayler RB, Kim SL (2024) SPORA, a new silver phosphate precipitation protocol for oxygen isotope analysis of small, organic-rich bioappatite samples. Chem Geol 65:122000. https://doi.org/10.1016/.chemgeo.2024. 122000
Copploe JV II, Blob RW, Parrish JHA, Butcher MT (2015) In vivo strains in the femur of the nine-banded armadillo (Dasypus novemcinctus). J Morphol 276:889-899. https://doi.org/10.1002 /jmor. 20387
Croft DA, Gelfo JN, López GM (2020) Splendid innovation: the extinct South American native ungulates. Ann Rev Earth Planet Sci 48:259-290.
Dantas, MAT, Omena ÉC, da Silva JLL, Sial A (2021) Could Eremotherium laurillardi (Lund, 1842) (Megatheriidae, Xenarthra) be an omnivore species? Anu Inst Geocienc 44:1-5. https://doi.org/10. 11137/1982-3908_2021_44_36492
Dantas MAT, Cherkinsky A, Bocherens H, Drefahl M, Bernardes C, de Melo França L (2017) Isotopic paleoecology of the Pleistocene megamammals from the Brazilian Intertropical Region: Feeding ecology ( ), niche breadth and overlap. Quat Sci Rev 170:152-163. https://doi.org/10.1016/j.quascirev.2017.06.030
Deka RJ, Sarma NK (2015) Studies on feeding behaviour and daily activities of Rhinoceros unicornis in natural and captive condition of Assam. Indiana J Anim Res 49(4):542-545.
Delsuc, F, Kuch M, Gibb GC, Karpinski, E, Hackenberger D, Szpak P, Billet G (2019) Ancient mitogenomes reveal the evolutionary history and biogeography of sloths. Curr Biol 29(12):2031-2042. https://doi.org/10.1016/j.cub.2019.05.043
Eagle RA, Schauble EA, Tripati AK, Tütken T, Hulbert RC, Eiler JM (2010) Body temperatures of modern and extinct vertebrates from bond abundances in bioapatite. Proc Natl Acad Sci USA 107(23):10377-10382. https://doi.org/10.1073/pnas. 0911115107
Eagle RA, Tütken T, Martin TS, Tripati AK, Fricke HC, Connely M, et al (2011) Dinosaur body temperatures determined from isotopic ( ) ordering in fossil biominerals. Science 333(6041):443445. https://doi.org/10.1126/science. 1206196
Eagle RA, Enriquez M, Grellet-Tinner G, Pérez-Huerta A, Hu D, Tütken T, et al (2015) Isotopic ordering in eggshells reflects body temperatures and suggests differing thermophysiology in two Cretaceous dinosaurs. Nat Commun 6(1):8296. https://doi.org/1 0.1038/ncomms9296
Ehlers J, Gibbard PL, Hughes PD (2018) Quaternary glaciations and chronology. In: Menzies J, van der Meer JJM (eds) Past Glacial Environments. Elsevier, Amsterdam, pp 77-101.
Fariña RA (2002) Megatherium, the hairless: appearance of the great Quaternary sloths (Mammalia; Xenarthra). Ameghiniana 39:241-244.
Fariña RA, Blanco RE (1996) Megatherium, the stabber. Proc R Soc B263(1377):1725-1729.
Fariña RA, Vizcaíno SF, Bargo MS (1998) Body mass estimations in Lujanian (late Pleistocene-early Holocene of South America) mammal megafauna. Mastozool Neotrop 5(2):87-108.
Feldhamer GA, Drickamer LC, Vessey SH, Merritt JF, Krajewski C (2007) Mammalogy: Adaptation, Diversity, and Ecology. John Hopkins University Press, Baltimore.
Fricke HC, Rogers R (2000) Multiple taxa and multiple locality approach to providing oxygen isotope evidence for endothermic homeothermy in theropod dinosaurs. Geology 28(9):799-802.
Ghosh P, Adkins J, Affek H, Balta B, Guo W, Schauble EA, et al (2006) bonds in carbonate minerals: a new kind of paleothermometer. Geochim Cosmochim Acta 70(6):1439-1456. https://d oi.org/10.1016/j.gca.2005.11.014
González-Bernardo E, Russo LF, Valderrábano E, Fernández Á, Penteriani V (2020) Denning in brown bears. Ecol Evol 10(13):68446862. https://doi.org/10.1002/ece3.6372
González-Guarda E, Domingo L, Tornero C, Pino M, Fernández MH, Sevilla P, et al (2017). Late Pleistocene ecological, environmental and climatic reconstruction based on megafauna stable isotopes from northwestern Chilean Patagonia. Quat Sci Rev 170:188202. https://doi.org/10.1016/j.quascirev.2017.06.035
Gorvet MA, Wakeling JM, Morgan DM, Segura DH, Avey-Arroyo J, Butcher MT (2020) Keep calm and hang on: EMG activation in the forelimb musculature of three-toed sloths (Bradypus variegatus). J Exp Biol 223(14):jeb218370. https://doi.org/10.1242/j eb. 218370
Gravett N, Bhagwandin A, Sutcliffe R, Landen K, Chase MJ, Lyamin OI, et al (2017) Inactivity/sleep in two wild free-roaming African elephant matriarchs-Does large body size make elephants the shortest mammalian sleepers? PL oS ONE 12(3):e0171903. https ://doi.org/10.1371/journal.pone. 0171903
Green, JL, Kalthoff, DC (2015) Xenarthran dental microstructure and dental microwear analyses, with new data for Megatherium americanum (Megatheriidae). J Mammal 96(4):645-657. https://doi.o rg/10.1093/mammal/gyv045
Griffiths ML, Eagle RA, Kim SL, Flores R, Becker MA, Maisch IV HM, et al (2023) Thermal physiology of extinct megatooth sharks. Proc Natl Acad Sci USA 120(27):e2218153120. https://d oi.org/10.1073/pnas. 2218153120
Hausman LA (1929) The ovate bodies of the hair of Nothrotherium shastense. Am J Sci 106:331-333.
Hayssen V (2011) Choloepus hoffmanni (Pilosa: Megalonychidae). Mammal Species 43(873):37-55. https://doi.org/10.1644/873.1
Ho TY (1967) Relationship between amino acid contents of mammalian bone collagen and body temperature as a basis for estimation of body temperature of prehistoric mammals. Comp Biochem Physiol 22:113-119.
Holden PB, Edwards NR, Rangel TF, Pereira EB, Tran GT, Wilkinson RD (2019) PALEO-PGEM v1. 0: a statistical emulator of Plio-cene-Pleistocene climate. Geosci Model Dev 12(12):5137-5155. https://doi.org/10.5194/gmd-12-5137-2019
Hutchinson JR, Schwerda D, Famini DJ, Dale RH, Fischer MS, Kram R (2006) The locomotor kinematics of Asian and African elephants:
changes with speed and size. J Exp Biol 209(19):3812-3827. http s://doi.org/10.1242/jeb. 02443
Idso SB, Jackson RD (1969) Thermal radiation from the atmosphere. J Geophys Res 74(23):5397-5403.
Jochum KP, Weis U, Schwager B, Stoll B, Wilson SA, Haug GH, et al. (2016). Reference values following ISO guidelines for frequently requested rock reference materials. Geostands Geoanal Res 40(3):333-350.
John CM, Bowen D (2016) Community software for challenging isotope analysis: First applications of ‘Easotope’ to clumped isotopes. Rapid Commun Mass Spectrom 30:2285-2300. https://doi .org/10.1002/rcm. 7720
Kleiber M (1932) Body size and metabolism. Hilgardia 6(11):315-353.
Knight F, Connor C, Venkataramanan R, Asher R (2020) Body temperatures, life history, and skeletal morphology in the nine-banded armadillo (Dasypus novemcinctus). PCI Ecol. https://doi.org/10 .17863/CAM. 50971
Kohn MJ, Schoeninger MJ, Barker WW (1999) Altered states: effects of diagenesis on fossil tooth chemistry. Geochim Cosmochim Acta 63(18):2737-47. https://doi.org/10.1016/S0016-7037(99)0 0208-2
Krmpotic CM, Carlini AA, Galliari FC, Favaron P, Miglino MA, Scarano AC, Barbeito CG (2014). Ontogenetic variation in the stratum granulosum of the epidermis of Chaetophractus vellerosus (Xenarthra, Dasypodidae) in relation to the development of cornified scales. Zoology 117(6):392-397 https://doi.org/10.101 6/j.zool.2014.06.003.
Larmon JT, McDonald HG, Ambrose S, DeSantis LR, Lucero LJ (2019) A year in the life of a giant ground sloth during the Last Glacial Maximum in Belize. Sci Adv 5(2): 1-9. https://doi.org/1 0.1126/sciadv.aau1200
Lécuyer C, Balter V, Martineau F, Fourel F, Bernard A, Amiot R, Gardien V, Otero O, Legendre S, Panczer G, Simon L, Martini, R (2010) Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at . Geochim Cosmochim Acta 74(7):2072-2081. https://doi.org/10.1016/j.gca.2009.12.02 4
Levesque DL, Marshall KE (2021) Do endotherms have thermal performance curves? J Exp Biol 224(3):jeb141309. https://doi.org/1 0.1242/jeb. 141309
Lindsey EL, Reyes EXL, Matzke GE, Rice KA, McDonald HG (2020) A monodominant late-Pleistocene megafauna locality from Santa Elena, Ecuador: Insight on the biology and behavior of giant ground sloths. Paleogeogr Paleoclimatol Paleoecol 544:109599. https://doi.org/10.1016/j.palaeo.2020.109599
Lovelace DM, Hartman SA, Mathewson PD, Linzmeier BJ, Porter WP (2020) Modeling Dragons: Using linked mechanistic physiological and microclimate models to explore environmental, physiological, and morphological constraints on the early evolution of dinosaurs. PLS ONE 15(5):e0223872. https://doi.org/10.17026/d ans-238-9pjs
MacFadden BJ, DeSantis LR, Hochstein JL Kamenov GD (2010) Physical properties, geochemistry, and diagenesis of xenarthran teeth: prospects for interpreting the paleoecology of extinct species. Paleogeogr Paleoclimatol Paleoecol 291(3-4):180-189. http s://doi.org/10.1016/j.palaeo.2010.02.021
Martin FM (2018) Cueva del Milodón. The hunting grounds of the Patagonian panther. Quat Int 466:212-222. https://doi.org/10.101 6/j.quaint.2016.05.005
McBee K, Baker RJ (1982) Dasypus novemcinctus. Mammal Species 162:1-9.
McCullough EC, Porter WP (1971) Computing clear day solar radiation spectra for the terrestrial ecological environment. Ecology 52:1008-1015.
McDonald HG (2005) Paleoecology of extinct xenarthrans and the Great American Biotic Interchange. Bull Fla State Mus Nat Hist 45(4):313-333.
McDonald HG (2022) Paleoecology of the Eextinct Shasta Ground Sloth, Nothrotheriops shastensis (Xenarthra, Nothrotheridae): The Physical environment. N M Mus Nat Sci Bull 88:33-43
McDonald HG, Lundelius Jr EL (2009) The giant ground sloth Eremotherium laurillardi (Xenarthra. Megatheriidae) in Texas. Bull Mus North Ariz 65:407-421.
McKamy AJ, Young MW, Mossor AM, Young, JW, Avey-Arroyo JA, Granatosky MC, Butcher MT (2023) Pump the Brakes! The hindlimbs of three-toed sloths control suspensory locomotion. J Exp Biol 226(8):jeb245622. https://doi.org/10.1242/jeb. 245622
McMillan BR, Cottam MR, Kaufman DW (2000) Wallowing behavior of American bison (Bison bison) in tallgrass prairie: an examination of alternate explanations. Am Mid Nat 144(1):159-167.
McNab, BK (1985) Energetics, population biology, and distribution of xenarthrans, living and extinct. In: Montgomery GG (ed) The Evolution and Ecology of Armadillos, Sloths, and Vermilinguas. Smithsonian Institution Press, Washington, D.C. pp 219-232.
Mine AH, Waldeck A, Olack G, Hoerner ME, Alex S, Colman AS (2017) Microprecipitation and analysis of phosphate for paleoclimate and biogeochemistry research. Chem Geol 460:1-4 https://doi.org/10.1016/j.chemgeo.2017.03.032.
Mitchell D, Snelling EP, Hetem RS, Maloney SK, Strauss WM, Fuller A (2018) Revisiting concepts of thermal physiology: predicting responses of mammals to climate change. J Anim Ecol 87(4):956-973. https://doi.org/10.1111/1365-2656.12818
Myhrvold CL, Stone HA, Bou-Zeid E (2012) What is the use of elephant hair? PL oS ONE 7:e47018. https://doi.org/10.1371/journa 1.pone. 0047018
Nance WB, Taylor SR (1976) Rare earth element patterns and crustal evolution-I. Australian post-Archean sedimentary rocks. Geochim Cosmochim Acta 40:1539-1551.
Naples VL (1990) Morphological changes in the facial region and a model of dental growth and wear pattern development in Nothrotheriops shastensis. J Vert Paleontol 10(3):372-389.
Nehme C, Todisco D, Breitenbach SF, Couchoud I, Marchegiano M, Peral M, et al (2022) Holocene hydroclimate variability along the Southern Patagonian margin (Chile) reconstructed from Cueva Chica speleothems. Glob Planet Change 222:104050. https://do i.org/10.1016/j.gloplacha.2023.104050
Nelson A, Engelman RK, Croft DA (2023) How to weigh a fossil mammal? South American notoungulates as a case study for estimating body mass in extinct clades. J Mammal Evol 30(3):773-809 h ttps://doi.org/10.1007/s10914-023-09669-1
Pant SRA, Goswami A, Finarelli JA (2014) Complex body size trends in the evolution of sloths (Xenarthra: Pilosa). BMC Evol Biol 14:184. https://doi.org/10.1186/s12862-014-0184-1
Pauli JN, Peery MZ, Fountain ED, Karasov WH (2016) Arboreal folivores limit their energetic output, all the way to slothfulness. Am Nat 188(2):196-204.
Poinar HN, Hofreiter M, Spaulding WG, Martin PS, Stankiewicz BA, Bland H, Evershed RP, Possnert G, Paabo S (1998) Molecular coproscopy: dung and diet of the extinct ground sloth Nothrotheriops shastensis. Science 128(5375):402-406.
Presslee S, Slater GJ, Pujos F, Forasiepi AM, Fischer R, Molloy K, et al. (2019). Palaeoproteomics resolves sloth relationships. Nat Ecol Evol 3(7):1121-1130. https://doi.org/10.1038/s41559-019-0 909-z
Prevosti FJ, Martin FM (2013) Paleoecology of the mammalian predator guild of Southern Patagonia during the latest Pleistocene: ecomorphology, stable isotopes, and taphonomy. Quat Int 305:74-84. https://doi.org/10.1016/j.quaint.2012.12.039
Pujos F, Salas R (2004) A systematic reassessment and paleogeographic review of fossil Xenarthra from Peru. Bull Inst Fr Et Andines 33:331-377.
Ridewood WG (1901). On the structure of the hairs of Mylodon listai and other South American Edentata. J Cell Sci 2(175):393-411. h ttps://doi.org/10.1242/jcs.s2-44.175.393
Santos LE, Ajala-Batista L, Carlini AA, de Araujo Monteiro-Filho EL (2024) Mylodon darwinii (Owen, 1840): hair morphology of an extinct sloth. Zoomorphology 2024:1-9.
Seitz VP, Puig S (2018) Aboveground activity, reproduction, body temperature and weight of armadillos (Xenarthra, Chlamyphoridae) according to atmospheric conditions in the central Monte (Argentina). Mammal Biol 88:43-51.
Seymour R S, Hu Q, Snelling E P, White C R (2019). Interspecific scaling of blood flow rates and arterial sizes in mammals. J Exp Biol 222(7):jeb199554. https://doi.org/10.1242/Jen. 199554
Stears K, McCauley DJ, Finlay JC, Mpemba J, Warrington IT, Mutayoba BM, et al (2018) Effects of the hippopotamus on the chemistry and ecology of a changing watershed. Proc Nat Acad Sci USA 115(22):E5028-E5037. https://doi.org/10.1073/pnas. 1800407115
Stelmock J, Dean FC (1986) Brown bear activity and habitat use, Denali National Park-1980. Int Conf Bear Res Manag 6:155-167.
Suárez MB, Passey BH (2014) Assessment of the clumped isotope composition of fossil bone carbonate as a recorder of subsurface temperatures. Geochim Cosmochim Acta 140:142-159. https:/ dol.org/10.1016|)- gca.2014.05.026
Sunquist ME, Montgomery GG (1973) Activity patterns and rates of movement of two-toed and three-toed sloths (Choloepus hoffmanni and Bradypus infuscatus). J Mammal 54(4):946-954.
Superina M, Abba AM (2014) Zaedyus pichiy (Cingulata: Dasypodidae). Mammal Species 46(905):1-10. https://doi.org/10.1644/90 5.1
Superina M, Boily P (2007) Hibernation and daily torpor in an armadillo, the pichi (Zaedyus pichiy). Comp Biochem Physiol A Mol Int Physiol 148(4):893-898. https://doi.org/10.1016/j.cbpa. 2007 . 09.005
Tagliavento M, Davies AJ, Bernecker M, Staudigel PT, Dawson RR, Dietzel M, et al (2023) Evidence for heterothermic endothermy and reptile-like eggshell mineralization in Troodon, a non-avian maniraptoran theropod. Proc Natl Acad Sci 120(15):e2213987120. https://doi.org/10.1073/pnas. 2213987120
Taube E, Vié JC, Fournier P, Genty C, Duplantier JM (1999) Distribution of two sympatric species of sloths (Choloepus didactylus and Bradypus tridactylus) along the Sinnamary River, French Guiana 1. Biotropica 31(4):686-691.
Tejada JV, Flynn JJ, MacPhee R, O’Connell TC, Cerling TE, Bermudez L, et al (2021) Isotope data from amino acids indicate Darwin’s ground sloth was not an herbivore. Sci Reports 11(1):1-10. https://doi.org/10.1038/s41598-021-97996-9
Tejada JV, Antoine PO, Münch P, Billet G, Hautier L, Delsuc F, Condamine FL (2024) Bayesian total-evidence dating revisits sloth phylogeny and biogeography: a cautionary tale on morphological clock analyses. Syst Biol 73(1):125-139. https://doi.org/10.1093 /sysbio/syad069
Toledo N, Boscaini A, Pérez LM (2021) The dermal armor of mylodontid sloths (Mammalia, Xenarthra) from Cueva del Milodon (Ultima Esperanza, Chile). J Morphol 282(4):612-627. https://do i.org/10.1002/jmor. 21333
Tomassini RL, Montalvo C I, Garrone MC, Domingo L, Ferigolo J, Cruz LE, et al (2020) Gregariousness in the giant sloth Lestodon (Xenarthra): multi-proxy approach of a bonebed from the Last Maximum Glacial of Argentine Pampas. Sci Reports 10(1):19055. https://doi.org/10.1038/s41598-020-67863-0
Tregear RT (1965) Hair density, wind speed, and heat loss in mammals. J Appl Physiol 20(4):796-801.
Varela L, Tambusso PS, Patiño SJ, Di Giacomo M, & Fariña RA (2018) Potential distribution of fossil xenarthrans in South America during the late Pleistocene: co-occurrence and provincialism. J Mamm Evol 25:539-550. https://doi.org/10.1007/s10914-017-9 406-9
Varela L, Clavijo L, Tambusso PS, Fariña RA (2023). A window into a late Pleistocene megafauna community: Stable isotopes show niche partitioning among herbivorous taxa at the Arroyo del Vizcaíno site (Uruguay). Quat Sci Rev 317:108286. https://doi.org/1 0.1016/j.quascirev.2023.108286
Varela L, Tambusso S, Fariña R. (2024). Femora nutrient foramina and aerobic capacity in giant extinct xenarthrans. PeerJ 12:e17815. ht tps://doi.org/10.7717/peerj. 17815
Wang Y, Porter WP, Mathewson PD, Miller PA, Graham RW, Williams JW (2018) Mechanistic modeling of environmental drivers of
woolly mammoth carrying capacity declines on St. Paul Island. Ecology 99(12):2721-2730. https://doi.org/10.1002/ecy. 2524
Wang Y, Widga C, Graham RW, McGuire JL, Porter WP, D. Wårlind, D, Williams JW (2020) Caught in a bottleneck: Habitat loss for woolly mammoths in central North America and the ice-free corridor during the last deglaciation. Glob Ecol Biogeogr 30(2):527542. https://doi.org/10.1111/geb. 13238
White CR, Seymour RS (2004) Does basal metabolic rate contain a useful signal? Mammalian BMR allometry and correlations with a selection of physiological, ecological, and life-history variables. Physiol Biochem Zool 77(6):929-941.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.