Jednym z powodów przeprowadzenia niniejszego eksperymentu była chęć zbadania czynników, które mogły przyczynić się do słabych wyników w poprzednim doświadczeniu, w którym porównywano straty mocznika w spływie powierzchniowym dla kilku form nawozów mineralnych i organicznych zawierających mocznik. Wszystkie warianty zastosowano do gleb, które zostały nasycone wodą, a następnie pozostawione do osięgnięcia pojemności polowej. Wyniki dla pięciu powtórzeń w wariancie z granulatem mocznika wykazały stężenia urea-N w spływie w zakresie od 1-12 mg/L. Tak duża zmienność między powtórzeniami w warunkach kontrolowanych była niedopuszczalna i zaburzyła wyniki eksperymentu. Silna dodatnia korelacja między całkowitą objętością spływu a stężeniem urea-N w spływie sugerowała, że czynnikami sprawczymi były warunki fizyczne, takie jak zagęszczenie gleby lub zmienna wilgotność początkowa wynikająca z różnych warunków drenażu i wysychania.
Aby zbadać przyczynę tak ekstremalnych różnic w stężeniach mocznika w spływach powierzchniowych, wszystkie skrzynki w niniejszym doświadczeniu zostały starannie wypełnione równymi masami jednorodnie wymieszanego gleby pylastej, jak przedstawiono na rysunkach 1 i 2 , w celu zminimalizowania zmienności warunków fizycznych. Aby uzyskać odpowiednio 50, 60, 70, 80, 90 i 100% przybliżonej pojemności polowej, określonej poprzez nawilżenie, a następnie wysuszenie w suszarce niewielkiej ilości przesianej gleby, obliczono masę wody niezbędną do nawilżenia gleby do odpowiadających im wstępnych wilgotności gleby wynoszących 14, 17, 19, 22, 25 i 27%. Wodę tę dodano do skrzynek i pozostawiono do ustabilizowania O/N. Symulacja opadów przebiegała zgodnie z dokładnym protokołem opisanym powyżej i przedstawionym na rysunku 3-5. Do zapewnienia natężenia opadów wynoszącego 3,2 cm/hr przez okres 40 min wykorzystano dyszę 17 wsq Full Jet 3/8 HH (tabela 1), co odpowiada naturalnemu zjawisku opadowemu, które powszechnie występuje corocznie na wschodnim brzegu zatoki Chesapeake w stanie Maryland.
Otrzymane całkowite objętości spływu, ładunki oraz stężenia ważone przepływem zestawiono w Tabeli 2. Stwierdzono istotną dodatnią zależność między całkowitą objętością spływu a wcześniejszym stanem wilgotności gleby (Rycina 6). Gleby bardziej wilgotne miały mniejszą zdolność do magazynowania wody i niższe wskaźniki infiltracji, co skutkowało większymi objętościami spływu. Odnotowano istotną ujemną zależność między czasem do wystąpienia spływu a wcześniejszym stanem wilgotności gleby (Rycina 7). W glebach bardziej suchych woda infiltrowała przez dłuższy czas, zanim stały się one wilgotne w pobliżu powierzchni, co doprowadziło do wystąpienia spływu. Oczekiwaniem była dodatnia zależność między całkowitym ładunkiem urea-N w spływie a całkowitą objętością spływu (Rycina 8). W badaniach hydrologicznych powszechnie wiadomo, że objętość przepływu jest zazwyczaj silnym predyktorem całkowitego ładunku. Sposób, w jaki stężenie zachowuje się w odpowiedzi na zdarzenie spływowe, jest mniej przewidywalny. Stężenie ważone przepływem obliczono poprzez zsumowanie ładunków dla każdego 2 min poboru spływu i podzielenie przez całkowitą objętość spływu. Jest to równoważne stężeniu w pojedynczej próbce spływu pobranej na koniec 40 min okresu opadów. W niniejszym badaniu stwierdzono istotną dodatnią zależność między stężeniem ważonym przepływem w spływie a wcześniejszym stanem wilgotności gleby (Rycina 9). Biorąc pod uwagę dodatnie zależności liniowe między objętością spływu a wcześniejszą wilgotnością gleby oraz między stężeniem ważonym przepływem a wcześniejszym stanem wilgotności, spodziewano się istotnej dodatniej zależności między całkowitym ładunkiem urea-N a wcześniejszym stanem wilgotności gleby. Jednakże ta istotna zależność była najlepiej opisana równaniem wykładniczym (Rycina 10).
W celu zwizualizowania utraty urea-N w spływie powierzchniowym w czasie, wartości poszczególnych stężeń w interwałach 2 min oraz ładunki skumulowane w jednej powtórce skrzynki glebowej reprezentującej każdy z wcześniejszych stanów wilgotności naniesiono na wykres w 40 min odstępie czasu opadów (Rysunek 11). Chociaż stężenia w spływie mogą zmieniać się w czasie dość nieregularnie (np. w przypadku wilgotności 90%), stężenia zazwyczaj są początkowo wysokie i maleją z upływem czasu. Ładunki skumulowane w czasie są funkcjami znacznie bardziej wygładzonymi i ilustrują omówione wcześniej istotne zależności. Czas do rozpoczęcia spływu jest dłuższy, stężenia urea-N w spływie są niższe, a ładunki skumulowane są mniejsze w przypadku gleb bardziej suchych. Chociaż mocznik ulega szybkiej hydrolizie w glebach, gdy opady występują w ciągu kilku godzin od aplikacji powierzchniowej, znaczna część N nadal występuje w formie mocznika i ulega utracie ze spływem powierzchniowym. Mocznik jest cząsteczką obojętną i nie jest silnie sorbowany na powierzchniach cząstek gleby. Gdy woda infiltruje bardziej suche gleby w początkowej fazie opadu, przenosi rozpuszczony mocznik w głąb gleby, z dala od strefy spływu powierzchniowego. Kiedy spływ faktycznie się rozpoczyna, obecność mocznika jest mniejsza, a stężenia w spływie są niższe. Z praktycznego punktu widzenia mocznik prawie zawsze byłby stosowany w suchszych warunkach, ponieważ sprzęt rolniczy nie mógłby poruszać się po glebach osiągających pełną pojemność wodną.

Rysunek 1. Schemat skrzynki do pomiaru spływu z zagęszczonej gleby. Metalowa skrzynka (100 cm x 20 cm x 7,5 cm) z 5 cm krawędzią na przednim końcu jest wypełniona zagęszczoną glebą do głębokości 5 cm. Spływ przelewający się przez 5 cm krawędź jest zbierany do przymocowanej rynny, która jest osłonięta przed bezpośrednim opadem deszczu. Dziewięć otworów o średnicy 5 mm umożliwia odpływ wody infiltrującej do gleby ze skrzynek i zapobiega tworzeniu się zastoisk. Przytwierdzony w pobliżu przedniej krawędzi dna rynny króciec umożliwia odpływ wody do lejków i butelek zbiorczych umieszczonych poniżej króćca.

Rycina 2. Materiały do pakowania skrzynki. Około 4 warstwy gazy na dnie skrzynki zapobiegają utracie gleby, jednocześnie umożliwiając swobodny odpływ wody. Przeciwwaga poziomująca, składająca się z pleksiglasu umieszczonego między dwiema drewnianymi deskami, ma szerokość równą szerokości skrzynki (20 cm) i głębokość (2.5 cm) odpowiadającą różnicy między bokami skrzynki (7.5 cm) a górną krawędzią rynienki (5 cm). Opierając deskę o krawędź skrzynki, pleksiglas służy do wyrównania gleby do poziomu rynienki.

Rycina 3. Pozycjonowanie platformy. Ustaw platformę w taki sposób, aby po umieszczeniu skrzynek z zagęszczoną glebą wszystkie miały ten sam kąt nachylenia. W niniejszym badaniu pożądany spadek wynosił 3%. Trzymając deskę w poziomie, ustaw platformę tak, aby dolny koniec skrzynki (przy rynnie) znajdował się 3 cm poniżej końca górnego. Platforma powinna być wypoziomowana w kierunku poprzecznym.

Rysunek 4. Elementy sterujące symulatora opadów, począwszy od źródła wody i przechodząc przez system rur aż do dyszy. (1) Kulowy zawór dźwigienkowy: Jest to zawór szybkiego odcięcia. Dźwignia w linii z rurą oznacza otwarcie; dźwignia pod kątem 90 stopni względem rury oznacza zamknięcie. Należy używać tego zaworu do włączania i wyłączania przepływu bez ingerencji w zawory kontrolujące ciśnienie i natężenie przepływu. Otwierać i zamykać całkowicie. Nie należy próbować używać tego zaworu do regulacji natężenia przepływu. (2) Filtr osadów: Należy okresowo sprawdzać filtr i wymieniać wkład w razie potrzeby, aby zapobiec zapychaniu się osadami. (3) Zawór regulatora ciśnienia: Zawór ten kontroluje ciśnienie w linii od tego punktu w przód. Zbyt wysokie ciśnienie może doprowadzić do uszkodzenia rur, węży lub połączeń. (4) Przepływowy zawór regulacyjny (zawór zasuwowy): Zawór ten służy do precyzyjnego dostrajania przepływu do dyszy w celu uzyskania pożądanego natężenia przepływu i ciśnienia w dyszy. (5) Przepływomierz: Pomiar przybliżonego natężenia przepływu. (6) Manometr: Pomiar przybliżonego ciśnienia w dyszy.

Rysunek 5. Pojemniki rozmieszczone na platformie do symulacji opadów deszczu. Umieść od 5 do 6 pojemników w wyznaczonych pozycjach dla każdego zdarzenia symulacji opadów. Unikaj ustawiania pojemnika bezpośrednio pod dyszą, aby zapobiec ściekaniu wody bezpośrednio na powierzchnię pojemnika.

Rycina 6. Całkowita objętość spływu jest dodatnio skorelowana z wcześniejszą zawartością wilgoci w glebie (R2 = 0,64).

Rycina 7. Czas do spływu jest ujemnie skorelowany z początkową zawartością wilgoci w glebie (R2 = 0,48). Powierzchnia wilgotnej gleby szybko ulega nasyceniu. Opady przekraczające przewodność hydrauliczną nasyconej gleby generują spływ.

Rycina 8. Całkowity ładunek urea-N jest dodatnio skorelowany z objętością spływu (R2 = 0.81). Różnice w objętości spływu dominują nad różnicami w stężeniu urea-N w spływie.

Rycina 9. Stężenie mocznika-N ważone przepływem jest dodatnio skorelowane z uprzednią zawartością wilgoci w glebie (R2 = 0.66). Suchsze gleby umożliwiają infiltrację, która wymywa mocznik-N w głąb gleby i z dala od jej powierzchni. W przypadku wystąpienia spływu powierzchniowego, na powierzchni dostępna jest mniejsza ilość mocznika-N do przemieszczenia wraz ze spływem.

Rysunek 10. Całkowity ładunek urea-N jest dodatnio skorelowany z uprzednią zawartością wilgoci w glebie (R2 = 0,74). Pozytywne zależności między całkowitą objętością spływu a uprzednią zawartością wilgoci w glebie oraz między stężeniem urea-N ważonym przepływem a uprzednią zawartością wilgoci łączą się, tworząc zależność wykładniczą (y = 0,2043 e0,0405x).

Rysunek 11. Zależności stężenia Urea-N i ładunku skumulowanego w czasie dla jednej powtórki dla każdej początkowej zawartości wilgoci w glebie. Chociaż stężenie Urea-N nie zawsze jest funkcją gładką w czasie, można zwizualizować istotne zależności omówione wcześniej.
| Rozmiar dyszy | Intensywność | Optymalne ciśnienie | Przepływ | Przepływ 10 s |
| cm/h | psi | gpm | ml |
| 17 wsq pełny strumień 3/8 HH | 3.2 | 6.0 | 1.5 | 940 |
| 24 wsq Full Jet 3/8 HH | 3.3 | 6.0 | 1.8 | 1,140 |
| 30 W Full Jet 1/2 HH | 6.0 | 5.0 | 2.2 | 1,250 |
| 50 w Full Jet 1/2 HH | 7.0 | 4.1 | 3.7 | 2,300 |
Tabela 1. Tabela rozmiarów dysz. Przedstawiono rozmiary dysz zidentyfikowane do użytku z tym symulatorem opadów oraz powiązane z nimi parametry intensywności opadów, ciśnienia i przepływu. Wybór rozmiaru dyszy zależy od pożądanej intensywności opadu. Intensywność i czas trwania opadów odpowiadają zdarzeniu opadowemu o określonym okresie powrotu dla konkretnej lokalizacji badania. W niniejszym badaniu zastosowano dyszę o rozmiarze 17 wsq. Opady o czasie trwania 40 min przy intensywności 3.2 cm/hr są równoważne naturalnemu zdarzeniu opadowemu, które powszechnie występuje w skali roku na wschodnim brzegu zatoki Chesapeake w stanie Maryland.
| Wilgotność gleby | Całkowity spływ | Ważone przepływem | Obciążenie całkowite |
| % | objętość (L) | stężenie | (mg N-mocznikowe) |
| | (mg L⁻¹)-1 azot mocznika) | |
| 27† | 2.96 | 4.99 | 13.66 |
| 27 | 2.87 | 4.37 | 12.55 |
| 25 | 2.52 | 3.57 | 8.62 |
| 25 | 1.81 | 2.33 | 4.21 |
| 22 | 2.52 | 2.18 | 5.50 |
| 22 | 2.47 | 1.54 | 3.81 |
| 19 | 1.99 | 1.72 | 3.41 |
| 19 | 2.35 | 3.70 | 8.68 |
| 17 | 1.91 | 1.69 | 3.22 |
| 17 | 1.66 | 0.90 | 1.50 |
| 14 | 1.51 | 0.78 | 1.18 |
| | | † Powtórzone liczby reprezentują dwie replikacje dla każdego poziomu wilgotności |
Tabela 2. Wstępna zawartość wody w glebie, całkowita objętość spływu, ważone przepływem stężenie urea-N oraz całkowite obciążenie urea-N po symulacji opadów. Powtarzające się liczby reprezentują dwie powtórki dla każdego poziomu wilgotności