  This paper complements the empirical justification of the revised scheme in
Part I of this work with a mathematical justification leveraging a
semi-discrete analysis framework for assessing the splitting error of process
coupling methods. The novelty of the framework is that splitting error is
distinguished from the process time integration errors, i.e., the errors caused
by discrete time integration of individual processes, leading to expressions
that are more easily interpreted utilizing existing physical understanding of
the processes that the terms represent. This application of this framework to
dust life cycle in EAMv1 showcases such an interpretation, using the
leading-order splitting error that results from the framework to confirm (i)
that the original EAMv1 scheme artificially strengthens the effect of dry
removal processes, and (ii) that the revised splitting reduces that artificial
strengthening. While the error analysis framework is presented in the context
of the dust life cycle in EAMv1, the framework can be broadly leveraged to
evaluate process coupling schemes, both in other physical problems and for any
number of processes. This framework will be particularly powerful when the
various process implementations support a variety of time integration
approaches. Whereas traditional local truncation error approaches require
separate consideration of each combination of time integration methods, this
framework enables evaluation of coupling schemes independent of particular time
integration approaches for each process while still allowing for the
incorporation of these specific time integration errors if so desired. The
framework also explains how the splitting error terms result from (i) the
integration of individual processes in isolation from other processes, and (ii)
the choices of input state and timestep size for the isolated integration of
processes.
