We have developed a simple, steady-state, one-dimensional second-order closure model to obtain continuous profiles of turbulent fluxes and mean concentrations of non-conserved scalars in a convective boundary layer without shear. As a basic tool we first set up a model for conserved species with standard parameterizations. This leads to formulations for profiles of the turbulent diffusivity and the ratio of temperature-scalar covariance to the flux of the passive scalar. The model is then extended to solving, in terms of profiles of mean concentrations and fluxes, the NO (x) -O(3) triad problem. The chemical reactions involve one first-order reaction, the destruction of NO(2) with decay time tau, and one second-order reaction, the destruction of NO and O(3) with the reaction constant k. Since the fluxes of the sum concentrations of NO (x) = NO + NO(2) and O(3) + NO(2) turn out to be constant throughout the boundary layer, the problem reduces to solving two differential equations for the concentration and the flux of NO(2). The boundary conditions are the three surface fluxes and the fluxes at the top of the boundary layer, the last obtained from the entrainment velocity, and the concentration differences between the free troposphere and the top of the boundary layer. The equations are solved in a dimensionless form by using 1/(k tau) as the concentration unit, the depth h of the boundary layer as the length unit, the convective velocity scale w (*) as the velocity unit, and the surface temperature flux divided by w (*) as the temperature unit. Special care has been devoted to the inclusion of the scalar-scalar covariance between the concentrations of O(3) and NO. Sample calculations show that the fluxes of the reactive species deviate significantly from those of non-reactive species. Further, the diffusivities, defined by minus the flux divided by the concentration gradient may become negative for reactive species in contrast to those of non-reactive species, which in the present model are never negative.