Earth
Monsone stagionale
Segui una costa dell’Asia meridionale nel corso dell’anno. Lo stesso sole scalda la terra poco profonda e il mare profondo, ma la terra risponde in pochi giorni e il mare con due mesi di ritardo, come dicono il guadagno in formula chiusa 1/sqrt(1+(ωτ)²) e il ritardo atan(ωτ). Il loro contrasto muove un vento di rovesciamento; verso terra porta l’umidità dell’aria marina, e il calore latente della pioggia rinforza proprio la circolazione che l’ha portata. Così il monsone non cresce dolcemente: appena la forzante della terra diventa positiva, la pioggia salta a un diluvio in una settimana, nubi scure arrivano dal mare e si ammassano contro le colline, le risaie si allagano e diventano verdi; resiste fino all’autunno, finché la forzante scende sotto la piega R_c = −M²/(4Dκ). Blocca la stagione e la costa è bistabile: un’irruzione umida la fa passare all’umido, abbassare il riscaldamento o l’umidità oltre la soglia in formula chiusa la secca per sempre. Il bilancio dell’umidità chiude a meno degli arrotondamenti.

Cosa succede
Due scatole sotto lo stesso riscaldamento stagionale H = S1·cos(ω(t − t_p)), S1 = 85 W/m², massimo al giorno 172, che si rilassano con B = 15 W/m²K verso T_e = 27 °C: la terra, 6 m di suolo (C_L = 1.5×10^7 J/m²K, τ = 11.6 giorni) con un riscaldamento aggiuntivo ΔQ = −76 W/m², e uno strato rimescolato oceanico di 50 m (τ = 158 giorni) di area quattro volte maggiore. Il contrasto x = T_L − T_O muove il vento bilanciato dall’attrito U = κx, κ = gH_a τ_d/(T_0 L_x) = 1.41 m/s per K (H_a = 1 km, τ_d = 1 giorno, L_x = 2000 km). Verso terra importa l’umidità dell’aria marina U·W_s/L_x (W_s = 55 kg/m²), verso il mare esporta quella della terra; l’umidità W sopra la terra precipita in un giorno, P = W/τ_P, e libera L_v·P sopra la terra; il rovesciamento porta D|U|x (D = 12) dalla scatola calda a quella fredda; il mare paga il calore latente di ciò che fornisce. La pioggia riempie il suolo e le risaie, che si svuotano in 30 giorni. A mare bloccato la terra si equilibra su Dκx² − Mx − R = 0, M = L_vκW_s/L_x − B = 82 W/m²K: il monsone si attiva quando R diventa positiva ed esiste fino a R_c = −M²/(4Dκ) = −99 W/m². A stagione bloccata entrambe le scatole si equilibrano su (1+r)Dκx² − ((1+r)L_vκW_s/L_x − B)x − ΔQ = 0: radici secca, instabile e umida tra ΔQ_c = −134 W/m² e zero. Le colline ripartiscono la stessa pioggia, (1 + h/h0)/(2 + h/h0) al lato sopravvento, h0 = 1 km. RK4 a 3 ore con i bilanci integrati con gli stessi pesi; l’anno comincia il 1° marzo dopo due anni di avviamento.
Da provare
Guarda il vento di terra della stagione secca invertirsi e la pioggia saltare in una settimana; confronta l’inizio con il passaggio per lo zero della forzante della terra e il ritiro con R_c; spegni il calore latente o la circolazione; rendi più profondi il mare o il suolo; blocca la stagione e falla passare all’umido con un’irruzione umida, poi abbassa il riscaldamento aggiuntivo della terra o l’umidità oltre la soglia e torna indietro.
Il modello e i suoi limiti
Due scatole, una di terra e una di mare, con un riscaldamento stagionale fisso e un’ora del giorno fissa: niente rotazione, nessuna fascia piovosa in movimento, nessuna evaporazione o riciclo sulla terra, e l’umidità dell’aria marina è un’impostazione invece di seguire la temperatura del mare. Il vento è l’equilibrio stazionario del contrasto con l’attrito. La piega R_c è quella dell’equilibrio della terra a mare bloccato, perciò il ritiro stagionale arriva alcune settimane dopo che la forzante è scesa sotto di essa; le radici a stagione bloccata sono esatte. Le colline si limitano a ripartire la stessa pioggia tra la pianura sopravvento e il lato sottovento. La scena non è in scala; nubi, polvere e fulmini sono decorazione guidata dall’intensità di pioggia.
Equazioni del modulo
Equazioni C_L·dT_L/dt = H(t) + ΔQ − B·(T_L − T_e) − D·|U|·x + L_v·P, H = S1·cos(ω·(t − t_p)), x = T_L − T_O C_O·dT_O/dt = H(t) − B·(T_O − T_e) + r·(D·|U|·x − L_v·F), U = κ·x, κ = g·H_a·τ_d/(T_0·L_x) dW/dt = F − W/τ_P, F = U·W_s/L_x (onshore) or U·W/L_x (offshore); imported = rain + exported + ΔW gain 1/√(1 + (ωτ)²), lag atan(ωτ); R_c = −M²/(4·D·κ), M = L_v·κ·W_s/L_x − B; held: (1+r)·D·κ·x² − ((1+r)·L_v·κ·W_s/L_x − B)·x − ΔQ = 0


