Each result cited is universally quantified over the data in its own statement. Integrals against push-forwards are transported by the change-of-variables formula, called "change of variables" below. L2(Ο;Rd) is the space of square-integrable vector fields against ΟβP(R3d), a real Hilbert space by that clause, with norm β₯β
β₯Οβ; the class of a Borel map square-integrable against Ο is written with the same symbol as the map.
Step 1 (Marginals of the gluing). By the definition of a gluing (Gluing Two Couplings over a Common Middle Marginal, and the Composite Coupling Β§glued), (q1β,q2β)#βΟ=Ο12β and (q2β,q3β)#βΟ=Ο23β. For a Borel set BβRd one has q1β1β(B)=(q1β,q2β)β1(pr1β1β(B)), q2β1β(B)=(q1β,q2β)β1(pr2β1β(B)) and q3β1β(B)=(q2β,q3β)β1(pr2β1β(B)); as Ο12ββΞ (Ξ½,Ο) and Ο23ββΞ (Ο,ΞΌ), this gives (q1β)#βΟ=Ξ½, (q2β)#βΟ=Ο and (q3β)#βΟ=ΞΌ.
Step 2 (The three fields over the gluing). Fix representatives of q, Ξ· and ΞΈ, that is, Borel maps RdβRd with finite square integrals against Ξ½, Ο and ΞΌ respectively (Wasserstein Spaces, Random Vectors, Vector Fields and Symmetric Matrices in Every Dimension: Standing Notation Β§fields). The maps qβq1β, Ξ·βq2β and ΞΈβq3β from R3d to Rd are Borel as compositions of Borel maps, and by change of variables and Step 1,
β«R3dββ₯qβq1ββ₯2dΟ=β«Rdββ₯qβ₯2dΞ½<β,
and likewise for Ξ·βq2β against Ο and ΞΈβq3β against ΞΌ. So their classes q^β, Ξ·^β, ΞΈ^ lie in L2(Ο;Rd).
Step 3 (Discrepancies as distances). The function zβ¦β₯q(x)βΞ·(y)β₯2 on Rd+d is Borel and nonnegative by Pairs of Euclidean Points: Coordinate Projections, Pairings, the Product Measure on a Euclidean Space, Borel Norm Functions and Finite Sets Β§functions, applied to the Borel maps qβpr1β and Ξ·βpr2β, and its composition with (q1β,q2β) is β₯qβq1ββΞ·βq2ββ₯2. By change of variables through (q1β,q2β) and the definition of β₯β
β₯Οβ,
β₯q^ββΞ·^ββ₯Ο2β=β«Rd+dββ₯q(x)βΞ·(y)β₯2Ο12β(dz).
In the same way, through (q2β,q3β) and through (q1β,q3β), with (q1β,q3β)#βΟ=Ο13β,
β₯Ξ·^ββΞΈ^β₯Ο2β=β«Rd+dββ₯Ξ·(x)βΞΈ(y)β₯2Ο23β(dz),β₯q^ββΞΈ^β₯Ο2β=β«Rd+dββ₯q(x)βΞΈ(y)β₯2Ο13β(dz).
These three integrals are the discrepancies of the statement, which do not depend on the representatives chosen by The Discrepancy of Two Square-Integrable Vector Fields Along a Coupling of Their Base Measures Β§well-defined.
Step 4 (Conclusion). In L2(Ο;Rd), q^ββΞΈ^=(q^ββΞ·^β)+(Ξ·^ββΞΈ^), so the triangle inequality The Norm Metric of a Real Inner Product Space: Triangle Inequalities, Limits and Continuity Β§triangle gives β₯q^ββΞΈ^β₯Οββ€β₯q^ββΞ·^ββ₯Οβ+β₯Ξ·^ββΞΈ^β₯Οβ. Each norm is the nonnegative square root of the corresponding integral of Step 3, which is the claimed inequality.