W pracy zaprezentowano geostatystyczną interpolację plonowania pszenicy ozimej w sąsiedztwie zadrzewień śródpolnych. Interpolację wykonano dla dwóch obszarów w obrębie jednego pola płodozmianowego o nachyleniu ok. 8%. Badania prowadzono przez dwa lata. Punktem wyjścia była analiza funkcji semiwariogramu określającej stopień zdeterminowania wartości plonu w danym punkcie przez wartości w punktach otaczających. Do szacowania plonów w miejscach, w których nie wykonywano pomiarów wykorzystano metodę krigingu. Metoda ta daje najlepsze oszacowanie wartości punktowych lub blokowych oraz pozwala ocenić niepewność przeprowadzonej estymacji w formie mapy błędu standardowego. Ocena zmienności przestrzennej plonowania pszenicy ozimej w sąsiedztwie zadrzewień śródpolnych przeprowadzona przy wykorzystaniu metod geostatystycznych wykazała zróżnicowanie plonowania w obrębie analizowanych powierzchni oraz przestrzenną jego zależność.