Articolo metodologico

Rinnovo di superficie: un metodo avanzato micrometeorologici per la misurazione e l'elaborazione del campo scala dell'energia densità di flusso di dati

DOI:

10.3791/50666

⸱

12 dicembre 2013

 ,  ,  ,  ,  , 
1Crops Pathology and Genetics Research Unit, United States Department of Agriculture-Agricultural Research Service, 2Department of Viticulture & Enology, University of California, Davis, 3Faculty of Agricultural Sciences, University of Chile, 4Atmospheric Science, University of California, Davis, 5URS Corporation Australia Pty. Ltd.

In questo articolo

Sommario

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Il rinnovo della superficie è un metodo micrometeorologico che viene utilizzato sempre più spesso per determinare i flussi di energia, ma la sua complessità tecnica lo rende inaccessibile a un vasto pubblico. Descriviamo i passaggi necessari per impostare e calibrare una stazione di campo di rinnovo della superficie, per acquisire ed elaborare i dati e per interpretare correttamente i risultati.

Abstract

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

I metodi micrometeorologici avanzati sono diventati sempre più importanti nelle scienze del suolo, delle colture e dell'ambiente. Per molti scienziati senza una formazione formale in scienze atmosferiche, queste tecniche sono relativamente inaccessibili. Il rinnovo della superficie e altri metodi di misurazione del flusso richiedono una comprensione della meteorologia dello strato limite e una formazione approfondita nella strumentazione e nei molteplici programmi di gestione dei dati. Per migliorare l'accessibilità di queste tecniche, descriviamo la teoria alla base delle misure di rinnovo della superficie, dimostriamo come configurare una stazione di campo per il rinnovo della superficie con calibrazione della covarianza parassita e utilizziamo il nostro programma di data logger chiavi in mano open source per eseguire l'acquisizione e l'elaborazione dei dati di flusso. Il nuovo programma chiavi in mano restituisce all'utente una semplice tabella di dati con i flussi corretti e i parametri di controllo della qualità ed elimina la necessità per i ricercatori di spostarsi tra più programmi di elaborazione per ottenere i dati finali del flusso. Un esempio di dati generati da queste misurazioni dimostra come viene misurato l'uso dell'acqua delle colture con questa tecnica. Le informazioni di output sono utili ai coltivatori per prendere decisioni di irrigazione in una varietà di ecosistemi agricoli. Queste stazioni sono attualmente impiegate in numerosi esperimenti sul campo da parte dei ricercatori del nostro gruppo e del California Department of Water Resources nelle seguenti colture: riso, vino e uva passa, erba medica, mandorlo, noce, pesca, limone, avocado e mais.

Introduzione

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

La scarsità d'acqua minaccia la sostenibilità agricola nelle regioni aride. Questa situazione probabilmente peggiorerà con il cambiamento delle condizioni climatiche e la crescente competizione tra gli enti agricoli, municipali, industriali e di conservazione per le limitate riserve idriche. Per affrontare questo dilemma in corso, i coltivatori sono alla continua ricerca di modi per migliorare l'efficienza dell'irrigazione e la tecnologia che stima meglio l'uso dell'acqua delle colture (ad esempio l'evapotraspirazione, ET) in tempo reale e migliori metodi agrometeorologici aiuteranno sicuramente questi sforzi. Le attuali tecnologie utilizzate dai coltivatori per stimare l'ET dipendono dall'uso dell'evapotraspirazione di riferimento (ET0) e di un coefficiente empirico della coltura (Kc), entrambi altamente suscettibili di errori di stima. Gli agrometeorologi misurano anche l'ET per valutare i trattamenti sperimentali per migliorare l'efficienza dell'uso dell'acqua delle colture1 e per parametrizzare le strategie regionali di allocazione dell'acqua2, ma questi metodi sono altamente tecnici e costosi. Sono attualmente in corso sforzi per tradurre gli attuali metodi di misurazione ET basati sulla ricerca in tecnologie convenienti e di facile utilizzo per i coltivatori.

Il rinnovamento della superficie (SR) è un metodo agrometeorologico utilizzato per misurare l'ET delle colture. La SR si basa sull'analisi del bilancio energetico delle particelle d'aria che risiedono in modo effimero all'interno della chioma della coltura durante il turbolento processo di scambio3, misurato con la stazione mostrata in Figura 1 e illustrato teoricamente in Figura 2 e nel Filmato 1. Le particelle d'aria si manifestano come forme simili a rampe nei dati delle serie temporali di temperatura turbolenta e l'ampiezza e il periodo delle rampe vengono utilizzati per calcolare la densità del flusso (Figure 2 e 3). Con il metodo SR, l'ET per una data superficie di coltura è determinato calcolando la densità del flusso di calore latente (LE) come residuo della seguente equazione del bilancio energetico
LE=R n-G-H,

Qui, LE è la densità del flusso di energia associata al cambiamento di fase dell'acqua dalla superficie di una coltura, Rn è la radiazione netta, G è la densità del flusso di calore del suolo (cioè l'energia condotta dentro o fuori dal terreno) e H è la densità del flusso di calore sensibile (cioè la densità del flusso di energia dalla superficie all'aria o viceversa che provoca un cambiamento di temperatura). Rn è un numero positivo quando il flusso netto è verso il basso (energia aggiunta alla superficie), LE e H sono numeri positivi quando il flusso è verso l'alto (energia aggiunta all'aria) e G è un numero positivo quando il flusso è verso il basso (energia aggiunta al suolo). LE (MJ/m2sec) viene quindi diviso per il calore latente di evaporazione (L=2,45 MJ/kg) per ottenere la densità di flusso di massa del vapore acqueo dalla superficie (cioè ET).

Le misurazioni di R, n e G sono relativamente semplici e poco costose. Le misure dirette di H sono più complesse e richiedono l'acquisizione di dati ad alta frequenza. Il metodo più comune per determinare H è con la covarianza parassita; tuttavia, l'anemometro sonico richiesto per questo metodo è costoso, complesso e, di conseguenza, non ampiamente utilizzato da agronomi, orticoltori o ingegneri per determinare H e LE. D'altra parte, l'H derivato dalla tecnica SR è ottenuto con un metodo più semplice e meno costoso, che utilizza termocoppie a filo sottile per misurare le temperature dell'aria ad alta frequenza all'interfaccia superficie-atmosfera. Nonostante la semplicità dell'SR, le misure di corrente richiedono ancora la calibrazione rispetto alla stima della covarianza parassita di H di un anemometro sonico per ottenere una densità di flusso di calore ragionevole.

Il successo dell'implementazione di una torre di flusso SR per le misure ET può essere una sfida scoraggiante, specialmente per molti ricercatori agricoli senza una formazione formale in scienze atmosferiche. Le misure di SR e eddy covariance richiedono competenze tecniche sofisticate sia nella programmazione di data logger per eseguire attività con strumentazione complessa sia nella scrittura di programmi informatici per post-elaborare i dati grezzi di turbolenza in flussi significativi. Qui descriviamo come configurare una stazione di campo, installare i sensori e utilizzare il nostro nuovo programma di registrazione dati chiavi in mano per le procedure di raccolta e post-elaborazione dei dati.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Protocollo

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

  1. Eseguire un'analisi della direzione del vento e dell'impronta prima dell'installazione sul campo; questo passaggio è importante per assicurarsi che la torre sia posizionata correttamente per catturare il flusso rappresentativo nel campo di studio.
    1. Prima di stabilire la stazione di campo, è necessario valutare la direzione del vento prevalente nel grafico di studio e determinare la quantità di recupero (distanza di bolina) necessaria per la propria particolare applicazione. La distanza del recupero dipende dall'altezza della superficie da misurare e dall'altezza degli strumenti4. Generalmente, il recupero è di 100 metri sottovento per ogni metro sopra il punto medio del baldacchino in cui si trovano i sensori.
  2. In base alle esigenze della configurazione specifica, collegare i sensori al data logger e preparare il più possibile l'attrezzatura prima di recarsi sul campo. Installare il programma nel logger prima dell'implementazione sul campo e verificare che il logger funzioni correttamente utilizzando i passaggi descritti di seguito.
    1. Collegare un laptop al data logger utilizzando la porta del cavo RS-232 (non isolata) sulla parte anteriore del datalogger, aprire il programma di interfaccia per comunicare con il data logger e connettersi al data logger. Quindi caricare il programma di registrazione dati chiavi in mano sul registratore di dati. I cavi standard USB e i telefoni cellulari possono essere utilizzati anche per comunicare con il datalogger.
      1. Il programma è scritto per una specifica suite di strumenti, vale a dire un anemometro sonico 3D, due piastre di flusso di calore, una termocoppia per terreno, una termocoppia per aria a filo sottile e un radiometro a rete. È possibile adattare il programma ad altre configurazioni di sensori, ma l'utente deve testare attentamente sul campo il programma del logger dopo qualsiasi modifica. Se al programma vengono aggiunte molte nuove righe di codice complesso, potrebbe non funzionare correttamente; Soprattutto quando l'anemometro sonico viene campionato in modalità seriale piuttosto che analogica. Il programma è disponibile per il download all'indirizzo: <sites.google.com/site/tmshapland>.
    2. Controlla i valori dei dati per tutti i sensori per vedere se tutto funziona correttamente.
      1. NOTA: Se una di queste variabili produce 'NAN' nella tabella di stato del programma di interfaccia del logger, è necessario controllare il sensore corrispondente e la sua connessione al logger. Dopo aver apportato le correzioni necessarie, tornare al computer e assicurarsi che il segnale non produca più 'NAN'.
  3. La configurazione descritta di seguito è per un sistema SR con una termocoppia a filo sottile e il cablaggio è descritto in modo che corrisponda ai comandi del registratore di dati incorporati nel programma chiavi in mano descritto di seguito.
    1. Una volta identificata una posizione appropriata per la tua stazione nel campo di studio, installa la torre (CM6- Figura 1). Stabilizzare i piedi della torre utilizzando i picchetti forniti, piantare la lunga asta di messa a terra in rame nel terreno con una mini-mazza in una posizione che non interferisca con le piastre di flusso di calore per il suolo - vedi sotto) e collegare l'asta di messa a terra alla torre utilizzando il cavo di calibro spesso e i connettori forniti.
    2. Fissare l'involucro del data logger all'asse principale della torre utilizzando i bulloni a U forniti dal produttore (Figura 1).
    3. Installare le due piastre di flusso di calore del terreno (cioè G1 e G2) con punto bianco a una profondità di 5 cm nello spazio interfilare lontano dagli emettitori di irrigazione a goccia (Figura 4). L'installazione deve garantire un buon contatto termico con il terreno e i cavi di prolunga che riportano al datalogger devono essere interrati per evitare danni dovuti allo strappo dal terreno.
      1. Collegare i fili della piastra di flusso termico ai canali 3 e 4 del data logger come sensori di ingresso differenziali (G1: filo nero 3H, filo rosso 3L - e G2: filo nero 4H, filo rosso 4L). I codici elencati in questo documento corrispondono alle etichette per i canali del data logger e verranno utilizzati in tutto il resto del manoscritto. I dettagli per tutti i sensori sono elencati nella Tabella 1.
      2. Nota: La termopila all'interno della piastra di flusso di calore del terreno misura il gradiente di temperatura attraverso un materiale con proprietà termiche note per arrivare alla densità del flusso di energia. Ogni piastra è calibrata dal produttore. Poiché le piastre di flusso di calore del suolo non possono essere posizionate direttamente sopra il terreno senza intercettare la radiazione solare, devono essere sepolte sotto la superficie del suolo. La densità del flusso di energia alla profondità delle piastre di flusso di calore del terreno non è la stessa alla densità del flusso di energia in superficie perché una parte dell'energia è immagazzinata nello strato di terreno sopra la piastra, quindi una serie di sensori di temperatura del suolo cablati in parallelo sono installati nello strato sopra le piastre.
    4. Installare la termocoppia del terreno (1 cavo con 4 sonde cablate in parallelo, denominate Ts) per coprire il volume del terreno sopra la piastra di flusso di calore del suolo per tenere conto della variazione dell'accumulo di calore sopra le piastre (Figura 4).
      1. Collegare la termocoppia del terreno al datalogger sul canale 5 come sensore differenziale (Ts: filo viola 5H, filo rosso 5L).
      2. Nota: La temperatura dello strato del suolo viene misurata inserendo sonde di temperatura del suolo con un angolo che va dal basso verso la parte superiore dello strato del suolo (circa 0,04-0,01 m sotto la superficie del suolo). È importante che le sonde di temperatura non siano inserite troppo vicino alla superficie. Se il terreno si rompe, le radiazioni e/o il flusso di calore sensibile da o verso i sensori possono causare dati falsi. In generale, è meglio non posizionare i sensori di temperatura a una distanza inferiore a 0,01 m dalla superficie del suolo; E un posizionamento più profondo è necessario in terreni che presentano grandi crepe a causa dell'essiccazione del suolo5.
    5. Installare il radiometro a rete (denominato Rn-Figura 5) sul braccio trasversale sopra il tettuccio e rivolto a sud.
      1. Collegare il radiometro di rete al data logger sul canale 2 come sensore differenziale (Rn: filo bianco 2H, filo verde 2L, schermo (nudo) segnale di massa)
    6. Installare la termocoppia a filo sottile (denominata TC1- Figura 6) a un'estremità del braccio trasversale e sopra la chioma (l'altezza dipenderà dalle piante che crescono in un particolare campo e dal recupero richiesto); la temperatura dell'aria sarà campionata a 10 Hz. Le termocoppie utilizzate per queste misurazioni hanno tipicamente un diametro di 76 μm.
      1. Collegare la termocoppia a filo sottile al registratore di dati sul canale 1 come sensore differenziale in cui il filo viola va a 1H, il filo rosso a 1L, il filo libero a massa del segnale.
      2. NOTA: Vedere i dettagli di seguito sulla manutenzione della termocoppia una volta che il sistema è attivo e funzionante.
    7. Installare l'anemometro sonico (Figura 7) all'altra estremità del braccio trasversale (cioè di fronte alla termocoppia a filo sottile) utilizzando il raccordo di giunzione Nu-rail per misurare le velocità del vento tridimensionali e la temperatura sonica a 10 Hz per calcolare il flusso di calore sensibile alla covarianza parassita. Installare l'anemometro sonico in modo che il centro dell'area di misurazione si trovi alla stessa altezza della termocoppia a filo sottile.
      1. Nota: l'anemometro sonico richiede la riconfigurazione dalle impostazioni di fabbrica prima dell'implementazione sul campo per emettere i segnali desiderati (per ulteriori informazioni, vedere il Instructions.pdf di riconfigurazione RM Young Model 81000RE disponibile all'indirizzo sites.google.com/site/tmshapland).
      2. La densità del flusso di calore sensibile non calibrata ottenuta dall'analisi del rinnovamento superficiale viene quindi calibrata determinando la pendenza della regressione lineare dei minimi quadrati della covarianza parassita rispetto alla densità del flusso di calore sensibile SR.
    8. Collegare l'alimentatore al data logger.
    9. Collegare un laptop al data logger utilizzando la porta del cavo RS-232 (non isolata) sulla parte anteriore del data logger. Caricare il programma di data logger chiavi in mano se questo passaggio non è già stato eseguito in laboratorio.
      1. Il programma è scritto per una specifica suite di strumenti, vale a dire un anemometro sonico, due piastre di flusso di calore, una termocoppia per terreno, una termocoppia per aria a filo sottile e un radiometro a rete. È possibile adattare il programma ad altre configurazioni di sensori, ma l'utente deve testare attentamente sul campo il programma del logger dopo qualsiasi modifica. Se al programma vengono aggiunte molte nuove righe di codice complesso, potrebbe non funzionare correttamente; Soprattutto quando l'anemometro sonico viene campionato in modalità seriale piuttosto che analogica. Il programma è disponibile per il download all'indirizzo: sites.google.com/site/tmshapland. Ulteriori dettagli sul programma e sul suo utilizzo sono disponibili nella nostra recente pubblicazione6.
    10. Controlla i valori dei dati per tutti i sensori per vedere se tutto funziona correttamente.
      1. NOTA: Se una di queste variabili produce 'NAN' nella tabella di stato del programma di interfaccia del logger, è necessario controllare o sostituire il sensore corrispondente. Dopo aver effettuato le regolazioni del sensore o del suo cablaggio, tornare al computer e assicurarsi che il segnale non produca più un "NAN".
  4. Mantenere la stazione utilizzando i seguenti passaggi una volta che la stazione è funzionante.
    1. Svitare il giunto Nu-Rail e abbassare il palo verticale principale che contiene il radiometro di rete.
    2. Piegare il filo isolato della termocoppia alla sua giunzione con il cilindro metallico. Immergi con cura i fili della termocoppia esposti (cioè la punta della termocoppia) nel succo di limone, facendo attenzione a non urtare i fili scoperti contro il barattolo. Pulisci il succo di limone dal sensore immergendo i fili nel barattolo di acqua deionizzata (DI). Raddrizzare nuovamente il filo della termocoppia alla sua giunzione con il cilindro metallico.
    3. Spruzzare il radiometro a rete con acqua deionizzata. Asciugalo con salviette assorbenti. Sollevare il montante verticale fino alla posizione originale, serrare nuovamente le viti di fermo Nu-Rail e livellare il radiometro netto.
    4. Scatta fotografie della torre per registrare eventuali cambiamenti nella chioma delle piante nel tempo.
    5. Elaborare i dati di flusso del bilancio energetico raccolti; la tabella di output del programma di data logger chiavi in mano contiene i già trasformati Rn, G, H e LE. Consulta il nostro recente lavoro per maggiori dettagli sull'uscita6,7 del programma di registrazione dati chiavi in mano.

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Risultati

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Acquisire le componenti R, n e G dell'equazione del bilancio energetico è relativamente semplice e poco costoso. Le misure di H sono più complesse e richiedono tracce di temperatura ad alta frequenza misurate con una termocoppia a filo sottile per il rinnovo della superficie o un anemometro sonico per la covarianza parassita. Un esempio di traccia di temperatura misurata con una termocoppia a filo sottile è rappresentata nella Figura 9 e mostra la necessità di un'ana...

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Discussione

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Il metodo di rinnovamento della superficie rappresenta un'opportunità tangibile per sviluppare una tecnica autonoma ed economica per quantificare l'uso dell'acqua delle colture in tempo reale. I recenti progressi nell'elaborazione dei segnali, nella programmazione dei data logger, nella calibrazione e nella gestione dei dati hanno messo a fuoco questo obiettivo in modo più chiaro. Questo manoscritto e il nostro programma di registrazione dati chiavi in mano6 sviluppato di recente rendono i metodi micrometeorol...

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Dichiarazioni

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Gli autori dichiarano di non avere attualmente interessi finanziari concorrenti.

Ringraziamenti

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Il supporto parziale per questa ricerca è stato fornito da J. Lohr Vineyards & Wines, dal National Grape and Wine Institute, da una sovvenzione NIFA Specialty Crops Research Initiative ad AJM e dal finanziamento USDA-ARS CRIS (Research Project #5306-21220-004-00).

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Materiali

Elenco dei materiali utilizzati in questo articolo
NomeAziendaNumero di catalogoCommenti
Registratore di datiCampbell Scientific, Inc.CR1000
Campbell Scientific, Inc.ENC12/14-SC
Torre (6 piedi) con messa a terra e parafulminiCampbell Scientific, Inc.CM6
Braccio trasversale (4 piedi)Campbell Scientific, Inc.CM204
Raccordo incrociato Nu-railCampbell Scientific, Inc.17953
Alimentatore (12 V) con regolatore e batteriaCampbell Scientific, Inc.PS100
(12 V)Campbell Scientific, Inc.CH100
Campbell Scientific, Inc.6186
Cavo di prolunga per termocoppiaCampbell Scientific, Inc.Termocoppie FWC-20
Campbell Scientific, Inc.FW3
3D con collo lungoRM Young Company81000 RE
8 cavo conduttore per anemometroRM Young Company18660
Piastre di flusso di calore del suoloREBS, Inc.HFT3.1
Termocoppie per terrenoCampbell Scientific, Inc.Radiometro TCAV
Net con cavoKipp and Zonen, Inc.NR Lite 2
Supporto da palo per radiometroKipp and Zonen, Inc.L-CMB1
Custodia per datalogger Regolatore del caricabatterie Cavo di prolunga della batteria Anemometro sonico

Riferimenti

$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Moratiel, R., Martínez-Cob, A. Evapotranspiration of table grape trained to a gable trellis system under netting and black plastic mulching. Irrigation Science. 30 (3), 167-178 (2011).
  2. DWR (Department of Water Resources). California Water Plan Update, Bulletin 160-05. , California Department of Water Resources. Sacramento, CA, U.S.A. (2005).
  3. Snyder, R. L., Spano, D., Duce, P., Paw, U. K. T., Qiu, J., Su, H. B., Watanabe, T., Brunet, Y. Surface renewal analysis: a new method to obtain scalar fluxes without velocity data. Agric. For. Meteorol. 74, 119-137 (1995).
  4. Rosenberg, N. J., Blad, B. L., Verma, S. B. Microclimate: the biological environment. , 2nd, John Wiley and Sons. New York, NY. USA. (1983).
  5. Snyder, R. L., Spano, D., Duce, P., Paw, U. K. T., Anderson, F. E., Falk, M. Surface Renewal Manual. , (2007).
  6. Shapland, T. M., McElrone, A. J., Paw, U. T. K., Snyder, R. L. A turnkey data logger program for field-scale energy flux density measurements using eddy covariance and surface. Ital. J. Agrometeorol. , Forthcoming Forthcoming.
  7. Shapland, T. M., Snyder, R. L., Paw, U. K. T., Lang, D., McElrone, A. J. New tools to improve access to micrometeorological methods for agricultural research scientists. Annual Meeting of the American Society of Agronomy, , (2011).
  8. Van Atta, C. W. Effect of coherent structures on structure functions of temperature in the atmospheric boundary layer. Arch. of Mech. 29, 161-171 (1977).
  9. Snyder, R. L., Spano, D., Paw, U. K. T. Surface Renewal analysis for sensible and latent heat flux density. Boundary-Layer Meteorol. 77, 249-266 (1996).
  10. Spano, D., Duce, P., Snyder, R. L., Paw, U. K. T. Surface renewal estimates of evapotranspiration. Tall canopies. Acta Hort. 449, 63-68 (1997).
  11. Spano, D., Snyder, R. L., Duce, P., Paw, U. K. T. Surface renewal analysis for sensible heat flux density using structure functions. Agric. For. Meteorol. 86, 259-271 (1997).
  12. Shapland, T. M., McElrone, A. J., Snyder, R. L., Paw, U. K. T. Structure function analysis of two-scale scalar ramps. Part I: Theory and modeling. Boundary Layer Meterology. , (2012).
  13. Shapland, T. M., McElrone, A. J., Snyder, R. L., Paw, U. K. T. Structure function analysis of two-scale scalar ramps. Part II: Applications. Boundary Layer Meterology. , (2012).

Accesso limitato. Accedi o avvia una prova gratuita per visualizzare questo contenuto.

Ristampe e permessi

Richiedi il permesso di riutilizzare il testo o le figure di questo articolo JoVE

Richiedi permesso

Tag

Covarianza dei vorticiacquisizione dei dati di flussoprogramma del data loggerconfigurazione della stazione di campoinstallazione dei sensorimisurazione dell evapotraspirazionegestione dell irrigazioneecosistemi agricoli

Articoli correlati