Udongo wa kilimo unapopewa fosforasi kwa miaka mingi kwa kiasi kikubwa kuliko mimea inavyoweza kuchukua, sehemu ya fosforasi isiyotumika inaweza kujikusanya kwenye udongo na kuunda hifadhi ya muda mrefu. Katika fasihi ya kisayansi, mzigo huu wa fosforasi uliokusanyika kutokana na matumizi ya zamani mara nyingi huitwa legacy phosphorus au kwa kifupi legacy P. Hifadhi hii mwanzoni inaweza kuonekana kama chanzo cha kirutubisho cha kilimo; lakini mkusanyiko unapozidi sana mahitaji ya mimea, kutokuwa na ufanisi wa kiuchumi na hatari ya upotevu wa fosforasi kwenda kwenye maji ya uso yanayozunguka hujitokeza kwa pamoja.

Utafiti huu ulijaribu kutatua tatizo muhimu la mbinu ili kuweza kulinganisha mkusanyiko wa fosforasi katika udongo wa kilimo wa Korea Kusini kwa kiwango cha kimataifa. Korea Kusini hutumia mbinu maalumu ya nchi hiyo ya modified Lancaster P katika uchanganuzi wa kawaida wa udongo. Katika hifadhidata nyingi duniani, Olsen P au Mehlich-3 P hutumiwa. Kwa kuwa mbinu hizi hutoa matokeo tofauti ya nambari zinapotumika kwenye udongo uleule, haiwezekani kulinganisha moja kwa moja data ya kitaifa ya Korea na ya nchi nyingine.

Kwa hiyo, watafiti kwanza walifanya kwa pamoja uchanganuzi wa modified Lancaster, Olsen na Mehlich-3 kwenye udongo 118 wa kilimo wa Korea Kusini na kuunda mahusiano ya ubadilishaji kati ya mbinu hizo. Baadaye, modeli hizo zilitumika kwenye udongo 6727 wa kilimo katika hifadhidata ya kitaifa ya ufuatiliaji wa udongo ya Korea Kusini.

Katika ulinganisho uliosanifishwa kwa kutumia hifadhidata za kimataifa zilizotumika katika utafiti, makadirio ya wastani wa Olsen P yalihesabiwa kuwa 233 ± 106 mg P/kg katika udongo wa kitalu nyumba, 158 ± 77 mg P/kg katika bustani za matunda, 147 ± 86 mg P/kg katika udongo mwingine wa mashamba, na 45 ± 37 mg P/kg katika udongo wa mashamba ya mpunga. Katika maeneo ya kimataifa yaliyolinganishwa, wastani wa Ulaya ulikuwa 38 ± 33 mg P/kg, Marekani 18 ± 22 mg P/kg, China 12 ± 27 mg P/kg na Australia 8 ± 9 mg P/kg.

Matokeo haya hayamaanishi kwamba “kila shamba duniani lilipimwa na Korea ndiyo ya kwanza bila shaka”. Matokeo yanaonyesha, ndani ya mipaka ya hifadhidata za kimataifa, modeli za ubadilishaji na hali za uchukuaji sampuli zilizotumika katika utafiti, kwamba Korea Kusini ina mkusanyiko wa juu sana wa fosforasi inayopatikana kwa mimea, hasa katika mifumo ya uzalishaji wa kina ya kitalu nyumba, bustani na mashamba.

Legacy fosforasi ni nini?

Fosforasi ni kirutubisho muhimu kwa mimea. Hata hivyo, ikiwa kiasi cha fosforasi kinachochukuliwa na mmea katika msimu mmoja wa uzalishaji ni kidogo kuliko kiasi kinachoongezwa kwenye udongo kwa mbolea au nyenzo za kikaboni, mizania chanya ya fosforasi hutokea katika mfumo. Hali hii ikiendelea kwa miaka, fosforasi hujikusanya kwenye udongo.

Si fosforasi yote inayoongezwa kwenye udongo hubaki daima katika umbo lililoyeyushwa ambalo mimea inaweza kulichukua kwa urahisi. Ayoni za fosfati zinaweza kujifunga kwenye oksidi za chuma na alumini, kutuama pamoja na kalsiamu au kushikiliwa katika awamu nyingine ngumu. Kwa hiyo, kiasi cha jumla cha fosforasi na kiasi cha fosforasi inayopatikana kwa mimea kinachopimwa na kipimo fulani cha udongo si dhana moja.

Legacy P hueleza hifadhi hii ya fosforasi iliyokusanyika kutokana na fosforasi iliyoingizwa hapo awali na ambayo imebaki katika mfumo wa udongo. Sehemu yake inaweza kupatikana kwa mimea siku zijazo; sehemu nyingine inaweza kubaki imefungwa kwa nguvu.

Kwa nini kuangalia fosforasi ya jumla pekee hakutoshi?

Kuwa na kiasi kikubwa cha fosforasi ya jumla kwenye udongo hakumaanishi kuwa fosforasi hiyo yote inaweza kuchukuliwa na mimea. Kwa hiyo, katika maamuzi ya uwekaji mbolea ya kilimo, mara nyingi hupimwa hazina ya fosforasi ya udongo inayoweza kuondolewa kwa mmumunyo maalumu wa kikemikali.

Vipimo vya Olsen, Mehlich-3 na Modified Lancaster viliundwa kwa kusudi hili. Lakini kwa kuwa vimiminika vya uchimbaji na hali za kazi ni tofauti, havitoi kiasi sawa cha fosforasi kutoka kwenye udongo uleule. Kwa hiyo, udongo unaopimwa kuwa 100 mg P/kg kwa mbinu moja si lazima utoe 100 mg P/kg kwa mbinu nyingine.

Kwa Nini Data ya Fosforasi ya Korea Kusini Haikuweza Kulinganishwa Moja kwa Moja na Dunia?

Katika ufuatiliaji wa kawaida wa udongo wa kilimo nchini Korea Kusini, mbinu ya Modified Lancaster hutumiwa. Mbinu hii ni mfumo maalumu wa uchimbaji uliorekebishwa kwa hali za Korea Kusini na hautumiki katika hifadhidata nyingi za kimataifa za fosforasi ya udongo. Kwa upande mwingine, Olsen P hutumika sana katika tafiti za kimataifa za kilimo na katika ulinganisho wa kiwango kikubwa.

Kwa hiyo, kuweka moja kwa moja matokeo ya Lancaster ya hifadhidata ya kitaifa ya Korea sambamba na matokeo ya Olsen P ya Ulaya, China, Australia au Marekani si sahihi kimethodolojia.

Watafiti kwanza walilinganisha mbinu zote tatu kwenye udongo uleule 118 ili kuunda daraja la pamoja la upimaji.

Katika udongo 118 wa utafiti, kiasi cha wastani cha fosforasi kilichochimbuliwa kilikuwa 360 ± 364 mg P/kg kwa mbinu ya Mehlich-3, 186 ± 164 mg P/kg kwa Modified Lancaster, na 112 ± 81 mg P/kg kwa mbinu ya Olsen.

Tofauti hii haionyeshi kuwa mbinu moja ni “sahihi” na nyingine ni “isiyo sahihi”. Kila mfumo wa uchimbaji hufikia hazina tofauti za kikemikali za fosforasi kwa viwango tofauti.

Matokeo ya Modified Lancaster Yalibadilishwaje Kuwa Olsen P?

Mahusiano makubwa ya kitakwimu yalipatikana kati ya mbinu tatu za uchimbaji. Mgawo wa uhusiano kati ya Modified Lancaster na Olsen uliripotiwa katika utafiti kama:

Olsen P = 0,79 × Lancaster P + 0,28; R² = 0,83

uliripotiwa.

Hata hivyo, mlinganyo huu si ubadilishaji rahisi wa mstari unaoweza kutumika moja kwa moja kwenye thamani ghafi za mg/kg. Kwa kuwa mgawanyo wa data ya fosforasi ulikuwa umepinda kuelekea kulia, thamani zilifanyiwa log-transformation kabla ya uchanganuzi. Katika matumizi kwenye hifadhidata ya kitaifa, Modified Lancaster P ilibadilishwa kwanza kwenda kwenye kipimo cha log, ubadilishaji wa regressioni ukafanywa, kisha thamani ya Olsen P iliyopatikana ikarudishwa kwenye kipimo chake cha awali.

Maelezo haya ni muhimu; kwa sababu kutumia mlinganyo kwenye thamani ghafi bila kuzingatia log-transformation hakuwezi kurudia makadirio ya kitaifa ya Olsen P yaliyopatikana katika utafiti.

Je, udongo 118 uliwakilisha hifadhidata ya kitaifa?

Ili milinganyo ya ubadilishaji itumike kwenye hifadhidata kubwa ya kitaifa, udongo 118 haupaswi kuwa kundi la sampuli tofauti kabisa na udongo wa nchi nzima. Watafiti walitathmini hili kwa kutumia madaraja ya udongo, pH, mboji, kationi zinazobadilishana, tekstua na sifa nyingine za kifizikia na kikemikali.

Taksonomia ya udongo, sifa kuu za kikemikali, uchanganuzi wa PCA na ulinganisho wa Random Forest ulidokeza kwamba sampuli 118 kwa ujumla zilikuwa na mwingiliano mpana na utofauti wa udongo katika hifadhidata ya kitaifa. Hata hivyo, hii haimaanishi kwamba kila eneo na kila mfumo wa uzalishaji uliwakilishwa kwa msongamano sawa ndani ya sampuli 118.

Ni Kiasi Gani cha Olsen P Kilihesabiwa katika Udongo wa Kilimo wa Korea Kusini?

Hifadhidata ya kitaifa ilijumuisha makundi manne makuu ya matumizi ya ardhi:

• udongo 2110 wa mashamba ya mpunga,
• udongo 1760 wa mashamba mengine,
• udongo 1473 wa bustani za matunda,
• udongo 1384 wa kitalu nyumba.

Hivyo, jumla ya udongo 6727 wa kilimo ulitathminiwa.

Wastani wa Modified Lancaster katika udongo wa kitalu nyumba pia ulikuwa 499 ± 565 mg P/kg, thamani ya juu zaidi kati ya madaraja ya matumizi. Watafiti wanahusisha matumizi makubwa ya mbolea na mboji katika mifumo ya kitalu nyumba na mkusanyiko mkubwa wa fosforasi.

Inamaanisha Nini Kwamba Takribani %90 ya Udongo wa Korea Una “Fosforasi Nyingi Kupita Kiasi”?

Kauli hii hutegemea kizingiti cha agronomia kinachotumika. Utafiti uliainisha udongo uleule kwa mifumo miwili tofauti ya marejeo na matokeo haya mawili hayapaswi kuchanganywa.

Wakati vizingiti vya agronomia vya Korea Kusini vinapotumika

Vizingiti vya Olsen P vilivyobadilishwa kutoka mapendekezo ya Korea yaliyotumika katika utafiti vilikuwa takribani 32–43 mg P/kg kwa mpunga na 52–145 mg P/kg kwa kundi lililounganishwa la upland.

Kulingana na marejeo haya ya kitaifa:

• takribani %38 ya udongo wa mashamba ya mpunga uliainishwa kuwa juu ya kiwango bora au juu kupita kiasi,
• takribani %59 ya udongo katika kundi la upland uliainishwa kuwa juu ya kiwango bora au juu kupita kiasi

.

Wakati vizingiti vya agronomia vya kimataifa vinapotumika

Watafiti pia walilinganisha udongo uleule na viwango vya chini zaidi vya optimum ya Olsen P vinavyotumika katika fasihi ya dunia. Walitumia 15 mg P/kg kwa mpunga na takribani 25–50 mg P/kg kwa mazao mengine makuu.

Mfumo huu wa marejeo ulipotumika, takribani %89,8 ya udongo wa mashamba ya mpunga na takribani %92,2 ya kundi la upland vilibaki juu ya optimum ya kimataifa.

Kwa hiyo, hitimisho kwamba “takribani %90 ya udongo wa kilimo wa Korea una fosforasi nyingi kupita kiasi” lilitolewa kwa mujibu wa vizingiti vya kimataifa vya ulinganisho vilivyotumika katika utafiti. Asilimia hii haitokani moja kwa moja na mfumo wa sasa wa uainishaji wa udongo wa Korea yenyewe.

Kiwango cha Olsen P cha Korea Kusini Kililinganishwaje na Maeneo ya Kimataifa?

Kwa ulinganisho mkuu wa kimataifa, watafiti walitumia hifadhidata ya kimataifa yenye rekodi za Olsen P zilizosanifishwa kutoka vipimo tofauti vya udongo.

Katika ulinganisho huu, wastani wa makundi yote ya matumizi nchini Korea ulikuwa juu kuliko wastani wa maeneo ya kimataifa yaliyolinganishwa. Wastani wa 233 mg P/kg katika udongo wa kitalu nyumba ni takribani mara 6,1 ya wastani wa Ulaya wa 38 mg P/kg.

Hata hivyo, jedwali hili halimaanishi kwamba “hakuna mahali Ulaya penye 233 mg/kg” au kwamba “kila shamba Korea lina kiwango cha juu kuliko kila shamba Ulaya”. Wastani wa kimataifa unaweza kuficha utofauti ndani ya maeneo makubwa ya kijiografia. Ulaya na Amerika Kaskazini pia zina maeneo madogo yenye fosforasi ya juu sana kutokana na ufugaji wa kina au matumizi makubwa ya mbolea kwa muda mrefu.

Mbinu na Matokeo ya Utafiti

Uchukuaji sampuli za udongo

Udongo 118 wa kilimo ulitumika kuunda mahusiano ya ubadilishaji. Kati yake, 34 ulikuwa udongo wa mashamba ya mpunga na 84 ulikuwa udongo wa kilimo wenye hali ya aerobiki. Ndani ya kundi la aerobiki kulikuwa na bustani 22 za matunda, kitalu nyumba 27 na udongo 35 wa mashamba mengine.

Udongo ulichukuliwa kutoka tabaka la juu la 0–15 cm, ukakaushwa hewani na kuchujwa hadi chembe zilizo chini ya 2 mm.

Uchanganuzi wa Modified Lancaster

5 g ya udongo ilitikiswa na 20 mL ya mmumunyo wa uchimbaji kwa dakika 10 kwa 120 rpm. Mmumunyo wa uchimbaji ulikuwa na 0,35 M asidi asetiki, 0,16 M asidi laktiki, 0,03 M floridi ya amoniamu, 0,05 M salfati ya amoniamu na 0,21 M hidroksidi ya sodiamu; pH iliwekwa kuwa 4,25.

Fosfati iliyochimbuliwa ilipimwa kwa kolorimetri ya molibdeti.

Uchanganuzi wa Olsen P

1 g ya udongo mkavu ilitikiswa kwa dakika 30 kwa 120 rpm na 20 mL 0,5 M NaHCO3 iliyowekwa katika pH 8,5. Sampuli zilisentrifugwa kwa 3000 × g kwa dakika 5 na fosfati katika supernatanti ilibainishwa kwa njia ya kolorimetri.

Uchanganuzi wa Mehlich-3

2 g ya udongo mkavu ilichimbuliwa kwa dakika 5 katika joto la kawaida kwa 20 mL ya mmumunyo wa Mehlich-3. Mmumunyo huo una asidi asetiki, nitrati ya amoniamu, floridi ya amoniamu, asidi nitriki na EDTA. Sampuli zilisentrifugwa kwa 3000 × g kwa dakika 10.

Kiwango cha kujaa kwa fosforasi

Ili kutathmini mzigo wa fosforasi wa udongo kwa kuzingatia uwezo wa kufunga wa chuma na alumini, kiwango cha kujaa kwa fosforasi — DPS — kilihesabiwa:

DPS (%) = P / [0,5 × (Al + Fe)] × 100

Katika hesabu hii, thamani za P, Al na Fe ziliwasilishwa kwa mmol/kg.

Uhusiano kati ya mbinu

Matokeo ya Modified Lancaster, Olsen na Mehlich-3 yalionyesha uhusiano mkubwa. Uhusiano wa Modified Lancaster–Mehlich-3 ulikuwa wenye nguvu zaidi, huku R² = 0,83 ikihesabiwa kwa uhusiano wa Modified Lancaster–Olsen.

Watafiti pia walijaribu mahusiano hayo kando katika udongo 91 wenye asidi na udongo 26 wenye alkali na waliripoti kuwa katika mkusanyo wa sampuli uliochunguzwa, mteremko au mkato wa regressioni haukutofautiana kwa kiwango kikubwa kulingana na daraja la pH. Hata hivyo, idadi ndogo ya udongo wenye alkali inaweza kupunguza nguvu ya kitakwimu.

Kwa nini udongo wa mpunga ni tofauti na udongo mwingine wa kilimo?

Ingawa wastani wa Olsen P ulikuwa chini katika udongo wa mashamba ya mpunga, katika kiwango sawa cha Olsen P kiasi cha jumla cha fosforasi kilikuwa juu kuliko katika udongo mwingine wa kilimo.

Utafiti unaeleza hili kwa sehemu kupitia mafuriko ya mara kwa mara na hali za anaerobiki katika kilimo cha mpunga zinazobadilisha muundo wa oksidi za chuma. Kiasi kikubwa cha oksidi za chuma zenye fuwele zisizokamilika katika udongo wa mpunga kinaweza kufunga fosfati kwa nguvu.

Kwa hiyo, ukiangalia Olsen P pekee, sehemu ya mkusanyiko wa jumla wa fosforasi ya zamani katika udongo wa mpunga inaweza kutoonekana.

Hitimisho zinazoungwa mkono na utafiti

• Mahusiano makubwa ya kitakwimu yanaweza kuanzishwa kati ya Modified Lancaster, Olsen na Mehlich-3.
• Milinganyo ya ubadilishaji iliwezesha hifadhidata ya kitaifa ya Lancaster ya Korea kutathminiwa kwa kutumia thamani sawa za Olsen P.
• Kuna mkusanyiko wa juu sana wa fosforasi inayopatikana kwa mimea katika udongo wa kitalu nyumba, bustani za matunda na mashamba mengine.
• Ndani ya hifadhidata za kimataifa zilizotumika katika utafiti, mifumo ya kilimo ya Korea ilionyesha viwango vya juu sana vya Olsen P.
• Vizingiti vya kimataifa vya agronomia vinapotumika, takribani %90 ya udongo wa kilimo wa Korea iko juu ya kiwango bora cha fosforasi.
• Katika udongo wa mpunga, Olsen P inaweza kuonyesha mkusanyiko wa jumla wa legacy P kuwa wa chini zaidi kwa kiwango kikubwa kuliko katika mifumo mingine ya matumizi.

Hitimisho ambazo utafiti hauungi mkono

• Haiwezi kusemwa kwamba maeneo yote ya kilimo duniani yamepimwa.
• Si kila shamba nchini Korea lina P ya juu kuliko kila shamba katika nchi nyingine.
• Olsen P haiwakilishi fosforasi yote ya jumla katika udongo.
• Olsen P ya juu pekee haitoi ushahidi wa moja kwa moja wa sababu na matokeo kwamba chanzo fulani cha maji kitachafuliwa.
• Milinganyo ya ubadilishaji haiwezi kutumika kiotomatiki kwa data ya Lancaster katika nchi zote na aina zote za udongo.
• Uainishaji wa udongo kuwa “umewekewa mbolea kupita kiasi” hautenganishwi na kizingiti cha agronomia kinachotumika.

Vyanzo vikuu vya kutokuwa na uhakika katika ulinganisho

Wakati sampuli katika hifadhidata ya kitaifa ya Korea zilikusanywa kutoka tabaka la juu la 0–15 cm, hifadhidata za ulinganisho wa kimataifa zinaweza kujumuisha kina tofauti lakini zinazokaribiana za tabaka la juu la udongo. Tofauti hii inaweza kuleta kutokuwa na uhakika, hasa katika mifumo ambako viwango hubadilika kulingana na kina cha udongo.

Ili kulinganisha akiba za fosforasi kwa kg P/ha na hifadhidata ya pili ya kimataifa, matokeo ya Korea ya 0–15 cm yalipanuliwa kwa ekstrapoleshini hadi kina cha 30 cm. Waandishi wanaeleza wazi kwamba kwa kuwa kiwango cha fosforasi kwa kawaida hupungua kadiri kina kinavyoongezeka, mbinu hii inaweza kukadiria akiba halisi kuwa kubwa kuliko ilivyo. Kwa hiyo, tafsiri kuu ya kimataifa ilitegemezwa kwenye viwango vya Olsen P.

Zaidi ya hayo, data ya fosforasi ya udongo inaweza kuwa imepinda sana kuelekea kulia. Sampuli chache zenye thamani za juu kupita kiasi zinaweza kuathiri kwa kiasi kikubwa mahusiano ya regressioni na wastani wa kikanda. Kwa hiyo, milinganyo inayobadilisha mbinu tofauti za uchimbaji haipaswi kuzingatiwa kuwa milinganyo ya ulimwengu mzima iliyo huru kabisa na eneo na sifa za udongo.

Maana ya kisayansi kwa usimamizi

Watafiti wanajadili njia ya kutumia kwa ufanisi zaidi hifadhi iliyopo ya fosforasi kwenye udongo katika mifumo yenye legacy P nyingi, na kupunguza fosforasi mpya inayoongezwa kulingana na uchanganuzi wa udongo au kuiwekea kikomo kwa muda katika hali fulani. Wazo kuu hapa ni kwamba ikiwa tayari kuna hifadhi ya fosforasi kwenye udongo iliyo juu ya mahitaji ya mimea, huenda isiwe lazima kuongeza kiwango kilekile cha fosforasi mpya kila mwaka.

Utafiti huu pekee hautoi pendekezo la dozi ya mbolea kwa mkulima, zao au shamba fulani. Uamuzi wa uwekaji mbolea unapaswa kufanywa kwa kuzingatia aina ya zao, kipimo cha udongo, lengo la mavuno, mapendekezo ya agronomia ya eneo husika na hali maalumu za shamba.

Maelezo ya Chanzo na Mbinu

Utafiti asili:Extreme legacy phosphorus accumulation in agricultural soils of South Korea: A global comparison of plant-available phosphorus

Waandishi: Maria Isabel Bacus; Aleksandra W. Bielnicka; So Yeong Park; Seulbi Lee; Seo Young Ho; Pil Joo Kim; Benjamin L. Turner.

Taasisi: Division of Applied Life Science, Gyeongsang National University; Rural Development Administration, Republic of Korea; Institute of Agriculture and Life Sciences, Gyeongsang National University; College of Resources and Environment, Shandong Agricultural University.

Jarida: Soil & Environmental Health.

Juzuu / nambari ya makala: 4 (2026), 100223.

DOI: 10.1016/j.seh.2026.100223

Mchapishaji: Elsevier B.V., kwa niaba ya Zhejiang University na Zhejiang University Press.

Tarehe ya kupokelewa: 12 Mei 2026.

Tarehe ya marekebisho: 18 Julai 2026.

Tarehe ya kukubaliwa: 18 Julai 2026.

Tarehe ya kuchapishwa mtandaoni: 21 Julai 2026.

Leseni: Creative Commons Attribution 4.0 International — CC BY 4.0.

Ufadhili: National Research Foundation of Korea, RS-2025-00516294. Maria Isabel Bacus, Aleksandra W. Bielnicka, So Yeong Park na Seo Young Ho waliungwa mkono na BK21+ Program; Benjamin L. Turner aliungwa mkono na fellowship ya NRF RS-2023-00217656.

Mgongano wa maslahi: Waandishi hawakuripoti maslahi ya kifedha yanayoshindana.

Upatikanaji wa data: Data zinazounga mkono uundaji wa milinganyo ya regressioni na thamani za Olsen P zilizobadilishwa kwa hifadhidata ya kitaifa zilichapishwa kupitia Figshare.

Vikwazo vikuu vya kimethodolojia: Ubadilishaji wa mbinu ulitengenezwa kwa kutumia udongo 118 wa kilimo wa Korea Kusini. Uhamishaji kwenda kwenye hifadhidata ya kitaifa unategemea dhana kwamba sampuli hizi zinawakilisha udongo wa nchi kwa kiwango cha kutosha. Ulinganisho wa kimataifa hutumia matokeo yaliyosanifishwa ya programu tofauti za ukusanyaji data, na kuna kutokuwa na uhakika kunakotokana na kina cha uchukuaji sampuli, uwakilishi wa kikanda na ubadilishaji wa mbinu. Ekstrapoleshini ya akiba za fosforasi hadi 30 cm inaweza kusababisha makadirio ya juu kuliko halisi.

Kanuni ya uwasilishaji ya Verianla: Maudhui haya si tafsiri kamili ya maandishi ya makala asili ya kitaaluma. Muundo wa majaribio, mantiki ya ubadilishaji, idadi ya sampuli za kitaifa, thamani za fosforasi, ulinganisho wa kimataifa na vikwazo vya utafiti vimehifadhiwa; maudhui ya kisayansi yamepangwa upya ndani ya maelezo huru ya kielimu. Kwa grafu asili za regressioni, ramani na mjadala wa kina wa kitaaluma, utafiti wa chanzo unapaswa kuchunguzwa kupitia DOI.