Sea salt aerosols: Difference between revisions

From Geos-chem
Jump to navigation Jump to search
Line 445: Line 445:
'''''Katerina Lapina wrote:'''''
'''''Katerina Lapina wrote:'''''


:I'd like to report a bug in v9-02 version of the code - I think I was using the latest version of the code, but could you please double-check? Thank you!
:I'd like to report a bug in [[GEOS-Chem v9-02]]—I think I was using the latest version of the code, but could you please double-check? Thank you!


:There is a bug in seasalt_mod.F which affects the dry deposition diagnostics (flux of sea salt). Line 871 of seasalt_mod.F reads  
:There is a bug in seasalt_mod.F which affects the dry deposition diagnostics (flux of sea salt). Line 871 of seasalt_mod.F reads  

Revision as of 14:09, 2 June 2014

Overview

The GEOS-Chem sea salt simulation is in the process of being updated following Jaeglé et al. [2011]. These changes will be implemented in v9-01-03. The reference document for the GEOS-Chem sea salt simulation prior to version v9-01-03 is Alexander et al [2005].

Recent Updates to sea salt simulation

Update molecular weight of sea salt tracers

This update was included in GEOS-Chem v9-02 (public release 03 Mar 2014).

Colette Heald wrote:

We'd like to propose a minor change to input.geos:
The molecular weight of the sea salt tracers should be listed as 31.4 g/mol (not as 36.0 as current). This is consistent with actual average composition of sea salt, and international guidelines from the IAPWS (http://www.iapws.org/relguide/seawater.pdf). It may be helpful to add a note somewhere to the effect that users should use the molar mass of any individual constituent of interest. While everything in the GEOS-Chem sea salt code is consistently calculated in terms of total mass, this clarification may be helpful for users who might use the MW from input.geos to convert the sea salt concentrations to global burdens, for example.

--Melissa Sulprizio 12:36, 9 December 2013 (EST)

SST dependent sea salt emissions

This update was tested in the 1-month benchmark simulation v9-01-03d and approved on 12 Jan 2012.

Sea salt emissions now include both a wind speed and sea surface temperature (SST) dependence. The sea salt source function is based on Gong (2003), which is based on Monahan et al. (1986). The Gong (2003) formulation expresses the density function dF/dr80 (in units of particles m-2 s-1 micrometer-1) as follows:

Gong.png

A and B are parameters depending on r80, the particle radius at RH = 80% (with r80 being close to twice the dry radius of sea salt particles).

Based on a comparison of GEOS-Chem sea salt simulation with coarse mode sea salt mass concentration observations obtained on 6 PMEL cruises, a new SST dependent source function was derived (Jaegle et al., 2011):

SST gong.png

where T is the SST expressed in degrees Celsius (valid temperature range: 0-30C).

This new empirical source function leads to improved agreement of GEOS-Chem relative to sea salt mass concentration observations from cruises and ground-based stations, as well as AOD observations from MODIS and AERONET.

Recommended size range for sea salt: Accumulation mode: 0.01-0.5 microns Coarse mode: 0.5 - 8 microns

Note that in Jaeglé et al. (2011) we used 1 accumulation bin (0.01-0.5) and 2 coarse mode bins (0.5-4; 4-10). Due to the non-linearity of dry deposition, using a single coarse bin 0.5-10 microns leads to an overestimate of the sea salt burden, hence we recommend using 0.5-8 microns.

Updates to sea salt dry deposition

This update was tested in the 1-month benchmark simulation v9-01-03d and approved on 12 Jan 2012.

Over land, sea salt dry deposition velocities are calculated using the Zhang et al. (2001) scheme, which is based on the Slinn (1982) model for vegetated canopies. Over the oceans, we have implemented the Slinn and Slinn (1980) deposition model for natural waters. Following the recommendation of Lewis and Schwartz (2004) we assume RH = 98% in the viscous sublayer (0.1-1mm thick layer above the ocean surface). We integrate the dry deposition velocity over each size bin using a bimodal size distribution for sea salt (see below), which includes growth as a function of local RH (see below).

Overall these changes lead to a factor of 3 increase in dry deposition velocity for coarse mode sea salt and a factor of 2 decrease for accumulation mode sea salt.

Updates to hygroscopic growth

The hygroscopic growth of sea salt aerosols is based on Equation (5) in Lewis and Schwartz (2006), which yields more accurate results at RH>98% than the Gerber (1985) formulation previously used in GEOS-Chem.

Updates to optical properties

The size distribution of accumulation mode sea salt aerosols assumes a dry geometric radius rdg=0.085 micrometers with a geometric standard deviation 2.03 micrometers. This is based on cruises in the remote Pacific Ocean (Quinn et al., 1996). For coarse mode sea salt aerosols we use rdg=0.4 micrometers with a geometric standard deviation of 1.8 micrometers based on Reid et al. (2006).

These size distributions are used in the Mie theory calculation of extinction efficiency. They are also used in calculating the size integrated dry deposition velocity of sea salt aerosols.

Overall impact on distribution of sea salt

Implementing these changes leads to small changes in the mean global burdens of accumulation mode (20% decrease) and coarse model (25% increase) sea salt aerosols. However, the spatial changes are much larger, with a 30-50% decrease at high latitudes and a factor of ~2 increase over tropical regions. See media: GC_seasalt_update.pdf for more info.

Prior Updates to sea salt emissions algorithm

Updated hygroscopic growth factors

Becky Alexander recommends new hygroscopic growth factors for sea salt aerosols. For more information, please see her post on the discussion page.

Modification of size bins for coarse mode aerosols

Lyatt Jaeglé wrote:

I think that we should change the dry size bins for the coarse mode aerosols in input.geos. Instead of using:
   Online SEASALT AEROSOLS : T
    => SALA radius bin [um]: 0.1  0.5
    => SALC radius bin [um]: 0.5 10.0
we should use a smaller upper cut for the dry radius of the coarse mode aerosols (up to 4 microns dry size ==> 8-10 microns radius for wet sea salt).
   Online SEASALT AEROSOLS : T
    => SALA radius bin [um]: 0.1  0.5
    => SALC radius bin [um]: 0.5  4.0
The first two changes lead to a factor of 2 decrease in total sea-salt emissions (from ~8000 Tg/yr to ~4300 Tg/yr using GEOS-4 winds for 2003 2x2.5). The last change leads to another reduction by 40% in emissions. However the total burden of sea-salt aerosols (~12.5 Tg) remains nearly unchanged compared to the old formulation (~13.7 Tg) because of the strong non-linearity of the dry deposition velocity at sizes > 2 microns. Indeed the lifetime of coarse mode sea-salt aerosols (0.5-4um vs 0.5-10um) increases by a factor of almost 3.
Here is a summary of the changes:


Old formulation in GEOS-Chem (2x2.5 GEOS-4 winds 2003)
  0.1-0.5 um 0.5-10 um Total: 0.1-10 um
Emissions (Tg/yr) 106 7865 7970
Dry deposition (Tg/yr) 4.7 5012 5016
Wet deposition (Tg/yr) 102 2859 2955
Burden (Tg) 0.73 13.01 13.74
Lifetime (hours) 60 14 15


New formulation, includes the following changes:
  0.1-0.5 um 0.5-4 um Total: 0.1-4 um
Emissions (Tg/yr) 92 2633 2725
Dry deposition (Tg/yr) 4 689 693
Wet deposition (Tg/yr) 87 1944 2031
Burden (Tg) 0.63 11.92 12.55
Lifetime (hours) 60 40 40


I am also including calculations for a 3rd bin size 4-10 microns
  4-10 um
Emissions (Tg/yr) 1544
Dry deposition (Tg/yr) 1156
Wet deposition (Tg/yr) 387
Burden (Tg) 1.6
Lifetime (hours) 9


While these larger aerosols (4-10 microns dry size) add another 50% to the emissions, they only contribute to 12% of the burden because of their short lifetime. So if we want to stick to 2 size bins, I think that it's fine to neglect these larger aerosols and limit the upper cut of the coarse mode sea-salt aerosols to 4 microns.
The overall sea-salt emissions ~3000 Tg/year is now similar to what other studies found when applying the Monahan formula: Monahan (1986), Spillane et al. (1986), Gong et al. (1998), Penner et al. (2001), etc... This is also within the range recommended by Lewis & Schwartz.
I also tried the Gong (2003) formulation which leads to a factor of ~2 decrease in emissions of accumulation mode aerosols but little change to the coarse mode aerosols. I am in the process of evaluating the sea-salt formulation against comparisons to cruise sea-salt observations from PMEL and find that both Gong (2003) and Monahan tend to overestimate sea-salt emissions at the high wind speeds in mid-latitudes and underestimate emissions in subtropical warmer waters. I am working on updating the Gong formulation based on SST.
Note that the optical properties currently used in the GEOS-Chem (jv_spec.dat) assume log-normal size distributions that lead to effective radii that are too large for coarse mode aerosols: ~9 microns at 50% RH! Based on observed sea-salt size distributions, this should be much smaller ~ 1-2 microns. This correction leads to larger AODs due to sea-salt. I think that Colette and Randall are working on updating jv_spec.dat and I will send them my recommendations.

--Bob Y. 10:56, 23 November 2009 (EST)

Comparing modeled sea salt to observations

Please see this discussion about comparing GEOS-Chem modeled sea salt concentrations to observations on the discussion page.

--Bob Y. 12:17, 25 February 2013 (EST)

References

  1. Alexander, B., R.J. Park, D.J. Jacob, Q.B. Li, R.M. Yantosca, J. Savarino, C.C.W. Lee, and M.H. Thiemens, Sulfate formation in sea-salt aerosols: Constraints from oxygen isotopes, J. Geophys. Res., 110, D10307, 2005. PDF
  2. Gong, S. L.: A parameterization of sea-salt aerosol source func- tion for sub- and super-micron particles, Global Biogeochem. Cy., 17(4), 1097, doi:10.1029/2003GB002079, 2003.
  3. Jaeglé, L., P.K. Quinn, T. Bates, B. Alexander, and J.-T. Lin (2011), Global distribution of sea salt aerosols: New constraints from in situ and remote sensing observations, Atmos. Chem. Phys., 11, 3137-3157, doi:10.5194/acp-11-3137-2011.PDF
  4. Lewis E. R. and Schwartz S. E., Comment on "Size distribution of sea-salt emissions as a function of relative humidity" Atmos. Environ. 40, 588-590 (2006); doi:10.1016/j.atmosenv.2005.08.043
  5. Quinn, P. K., et al.: Chem- ical and optical properties of marine boundary layer aerosol particles of the mid-Pacific in relation to sources and meteorological transport, optical properties of sea salt aerosols, J. Geophys. Res., 102, 23269–23275, 1996.
  6. Reid, J. S., et al.: Reconciliation of coarse mode sea-salt aerosol particle size measurements and parameterizations at a sub- tropical ocean receptor site, J. Geophys. Res., 111, D02202, doi:10.1029/2005JD006200, 2006.
  7. Slinn, W. G. N.: Predictions for particle deposition to vegetative canopies, Atmos. Environ., 16, 1785–1794, 1982.
  8. Slinn, S. A. and Slinn, W. G. N.: Predictions for particle deposition on natural-waters, Atmos. Environ., 14, 1013–1016, 1980.
  9. Zhang, L., Gong, S., Padro, J., and Barrie, L.: A size-segregated particle dry deposition scheme for an atmospheric aerosol mod- ule, Atmos. Environ., 35, 549–560, 2001.

Previous issues that have now been resolved

Double-substitution bug in routine GET_ALK

NOTE: This fix was standardized into GEOS-Chem v8-01-02 and higher versions.

Becky Alexander wrote:

The code in GET_ALK (in routine seasalt_mod.f) as it is now is wrong. I did a substitution twice by mistake, that should have been applied only once. This is calculated for both accumulation and coarse mode seasalt, for both SO2 and HNO3, so there are 4 places in the code that must be fixed.
The correct code should be as follows:
            !----------------------------------
            ! SO2 uptake onto fine particles 
            !----------------------------------
 
            ! calculate gas-to-particle rate constant for uptake of 
            ! SO2 onto fine sea-salt aerosols [Jacob, 2000] analytical solution
            CONST1 = 4.D0/(V*GAMMA_SO2)
            A1     = (RAD1/DG)+CONST1
            B1     = (RAD2/DG)+CONST1
   !-----------------------------------------------------------------------------
   ! Prior to 7/18/08:
   ! Becky Alexander's fix to remove double-substitution (bec, bmy, 7/18/08)
   ! Remove these lines:
   !         TERM1A = (((B1/DG)**2)+(2.0D0*CONST1*B1/DG)+(CONST1**2)) -
   !     &            (((A1/DG)**2)+(2.0D0*CONST1*A1/DG)+(CONST1**2))
   !         TERM2A = 2.D0*CONST1*(((B1/DG)+CONST1)-((A1/DG)+CONST1))
   !         TERM3A = (CONST1**2)*(LOG((B1/DG)+CONST1) -
   !     &            LOG((A1/DG)+CONST1))
   !         KT1    = 4.D0*PI*N1*(DG**2)*(TERM1A - TERM2A + TERM3A)
   !-----------------------------------------------------------------------------
            TERM1A = ((B1**2)/2.0d0) - (((A1**2)/2.0d0)
            TERM2A = 2.D0*CONST1*(B1-A1)
            TERM3A = (CONST1**2)*LOG(B1/A1)
            KT1    = 4.D0*PI*N1*(DG**3)*(TERM1A - TERM2A + TERM3A)

            !----------------------------------
            ! SO2 uptake onto coarse particles 
            !----------------------------------
        
            ! calculate gas-to-particle rate constant for uptake of 
            ! SO2 onto coarse sea-salt aerosols [Jacob, 2000] analytical solution
            CONST2 = 4.D0/(V*GAMMA_SO2)
            A2     = (RAD2/DG)+CONST2
            B2     = (RAD3/DG)+CONST2
   !------------------------------------------------------------------------------
   ! Prior to 7/18/08:
   ! Becky Alexander's fix to remove double-substitution (bec, bmy, 7/18/08)
   ! Remove these lines:
   !         TERM1B = (((B2/DG)**2)+(2.0D0*CONST2*B2/DG)+(CONST2**2)) -
   !     &            (((A2/DG)**2)+(2.0D0*CONST2*A2/DG)+(CONST2**2))
   !         TERM2B = 2.D0*CONST2*(((B2/DG)+CONST2)-((A2/DG)+CONST2))
   !         TERM3B = (CONST2**2)*(LOG((B2/DG)+CONST2) -
   !     &             LOG((A2/DG)+CONST2))
   !         KT2    = 4.D0*PI*N2*(DG**2)*(TERM1B - TERM2B + TERM3B)
   !------------------------------------------------------------------------------
            TERM1B = ((B2**2)/2.0d0) - (((A2**2)/2.0d0)
            TERM2B = 2.D0*CONST2*(B2-A2)
            TERM3B = (CONST2**2)*LOG(B2/A2)
            KT2    = 4.D0*PI*N2*(DG**3)*(TERM1B - TERM2B + TERM3B)
            KT     = KT1 + KT2

            !----------------------------------
            ! HNO3 uptake onto fine particles 
            !----------------------------------

            ! calculate gas-to-particle rate constant for uptake of 
            ! HNO3 onto fine sea-salt aerosols [Jacob, 2000] analytical solution
            CONST1N = 4.D0/(V*GAMMA_HNO3)
            A1N     = (RAD1/DG)+CONST1N
            B1N     = (RAD2/DG)+CONST1N
   !-----------------------------------------------------------------------------
   ! Prior to 7/18/08:
   ! Becky Alexander's fix to remove double-substitution (bec, bmy, 7/18/08)
   ! Remove these lines:
   !         TERM1AN = (((B1N/DG)**2)+(2.0D0*CONST1N*B1N/DG)+(CONST1N**2)) -
   !     &             (((A1N/DG)**2)+(2.0D0*CONST1N*A1N/DG)+(CONST1N**2))
   !         TERM2AN = 2.D0*CONST1N*(((B1N/DG)+CONST1N)-((A1N/DG)+CONST1N))
   !         TERM3AN = (CONST1N**2)*(LOG((B1N/DG)+CONST1N) -
   !     &             LOG((A1N/DG)+CONST1N))
   !         KT1N    = 4.D0*PI*N1*(DG**2)*(TERM1AN - TERM2AN + TERM3AN)
   !-----------------------------------------------------------------------------
            TERM1AN = ((B1N**2)/2.0d0) - (((A1N**2)/2.0d0)
            TERM2AN = 2.D0*CONST1N*(B1N-A1N)
            TERM3AN = (CONST1N**2)*LOG(B1N/A1N)
            KT1N    = 4.D0*PI*N1*(DG**3)*(TERM1AN - TERM2AN + TERM3AN)

            !----------------------------------
            ! HNO3 uptake onto coarse particles 
            !----------------------------------

            ! calculate gas-to-particle rate constant for uptake of 
            ! HNO3 onto coarse sea-salt aerosols [Jacob, 2000] analytical solution
            CONST2N = 4.D0/(V*GAMMA_HNO3)
            A2N     = (RAD2/DG)+CONST2N
            B2N     = (RAD3/DG)+CONST2N
   !-----------------------------------------------------------------------------
   ! Prior to 7/18/08:
   ! Becky Alexander's fix to remove double-substitution (bec, bmy, 7/18/08)
   ! Remove these lines:
   !         TERM1BN = (((B2N/DG)**2)+(2.0D0*CONST2N*B2N/DG)+(CONST2N**2)) -
   !     &             (((A2N/DG)**2)+(2.0D0*CONST2N*A2N/DG)+(CONST2N**2))
   !         TERM2BN = 2.D0*CONST2N*(((B2N/DG)+CONST2N)-((A2N/DG)+CONST2N))
   !         TERM3BN = (CONST2N**2)*(LOG((B2N/DG)+CONST2N) -
   !     &             LOG((A2N/DG)+CONST2N))
   !         KT2N    = 4.D0*PI*N2*(DG**2)*(TERM1BN - TERM2BN + TERM3BN)
   !-----------------------------------------------------------------------------
            TERM1BN = ((B2N**2)/2.0d0) - (((A2N**2)/2.0d0)
            TERM2BN = 2.D0*CONST2N*(B2N-A2N)
            TERM3BN = (CONST2N**2)*LOG(B2N/A2N)
            KT2N    = 4.D0*PI*N2*(DG**3)*(TERM1BN - TERM2BN + TERM3BN)

Please make the fix in your version, or you may download it from ftp://ftp.as.harvard.edu/pub/geos-chem/patches/v8-01-01/seasalt_mod.f_w_getalk_fix.

Also see this document by Becky Alexander which describes the analytical solution in more detail.

Duncan Fairlie replied:

Thanks for taking a look at this. I will look at your code corrections and integrate them into my dust code.
Since we're looking back at the analytical solution, I think the last line should read
   Kt = 4.pi.N.D(cubed)[  ]  ,
the extra factor of D coming from
   dr = D.dx,
and the limits of the integral
  ( r=[a,b] )
become
   X = [a/D+c, b/D+c]
I'll recheck my math, and look back at the code.....

--Bob Y. 16:25, 22 February 2010 (EST)

Minor bug in AD44 diagnostic in routine DRY_DEPOSITION

NOTE: This bug fix was included in GEOS-Chem v9-02 (public release 03 Mar 2014).

Bob Yantosca wrote:

I just found and fixed a minor issue in the DRY_DEPOSITION routine in seasalt_mod.F. This issue only arises when the TURBDAY PBL mixing option is selected, so it probably has not affected most people's simulations.
I was getting numerical differences (probably due to REAL*4 roundoff) when we saved directly into the AD44 array within the loop over (I,J,L). It is more robust method to save into a temporary REAL*8 array (ND44_TMP) within the parallel loop and then just sum in the vertical in a separate loop.
   !
   ! !LOCAL VARIABLES:
   !
         ... 

         REAL*8             :: ND44_TMP(IIPAR,JJPAR,LLPAR)

         ...
  
         ! Zero temporary array for drydep diagnostic
         IF ( ND44 > 0 ) ND44_TMP = 0d0

   !$OMP PARALLEL DO
   !$OMP+DEFAULT( SHARED )
   !$OMP+PRIVATE( I, J, AREA_CM2, OLD, NEW, FLUX )
   !$OMP+PRIVATE( L, F_UNDER_TOP   , FREQ ) ! (jaegle 5/11/11)
   !$OMP+SCHEDULE( DYNAMIC )

         ! Note, for grid-independent GEOS-Chem we have to loop over
         ! all levels and then just shunt if we are in the stratosphere
         ! (bmy, 11/20/12)
         DO L = 1, LLPAR
         DO J = 1, JJPAR
         DO I = 1, IIPAR

            ...

            ! Only apply drydep to boxes w/in the PBL
            IF ( F_UNDER_TOP > 0d0 ) THEN

               ...

                  !===========================================================
                  ! ND44 diagnostic: sea salt drydep loss [molec/cm2/s]
                  !===========================================================
                  IF ( ND44 > 0 ) THEN
               
                     ! Convert drydep loss from [kg/s] to [molec/cm2/s]
                     FLUX = ( OLD - NEW ) / DTCHEM 
                     FLUX = FLUX * XNUMOL(IDTSALA) / AREA_CM2 
            
                     ! Save in ND44 temporary array
                     ND44_TMP(I,J,L) = ND44_TMP(I,J,L) + FLUX
                  ENDIF
   
                  ! Update tracer array
                  TC(I,J,L) = NEW 
               ENDIF
            ENDIF
         ENDDO
         ENDDO
         ENDDO
   !$OMP END PARALLEL DO

         !====================================================================
         ! ND44 diagnostic: save into AD44 array, summing in the vertical
         !====================================================================
         IF ( ND44 > 0 ) THEN
   !$OMP PARALLEL DO
   !$OMP+DEFAULT( SHARED )
   !$OMP+PRIVATE( I, J )
            DO J = 1, JJPAR
            DO I = 1, IIPAR
               AD44(I,J,IDDEP(N),1) = SUM( ND44_TMP(I,J,:) ) 
            ENDDO
            ENDDO
   !$OMP END PARALLEL DO
         ENDIF

--Bob Y. 11:58, 28 November 2012 (EST)

Bug in ND44 drydep diagnostic for sea salt aerosols

NOTE: This bug is being validated with 1-month benchmark simulation GEOS-Chem v10-01d.

Katerina Lapina wrote:

I'd like to report a bug in GEOS-Chem v9-02—I think I was using the latest version of the code, but could you please double-check? Thank you!
There is a bug in seasalt_mod.F which affects the dry deposition diagnostics (flux of sea salt). Line 871 of seasalt_mod.F reads
   AD44 (I, J, IDDEP(N), 1)  =  SUM (  ND44_TMP(I, J, :)  )  )
while it should be
   AD44 (I, J, IDDEP(N), 1)  =  AD44 (I, J, IDDEP(N), 1)  + SUM (  ND44_TMP(I, J, :)  )  )
That is, AD44 currently saves only the "dry_deposition" component of the flux and ignores the "wet_settling" component. This part of the code executes when non local PBL mixing is off.

--Bob Y. 10:03, 2 June 2014 (EDT)

Retrieved from "https://wiki.seas.harvard.edu/geos-chem/index.php?title=Sea_salt_aerosols&oldid=17427"